Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_010473 4,272,972 +T intergenic (‑226/‑135) coaA ← / → ECDH10B_RS21215 type I pantothenate kinase/tRNA‑Thr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0104734,272,9721.T86.4% 71.8 / 0.7 22intergenic (‑226/‑135)coaA/ECDH10B_RS21215type I pantothenate kinase/tRNA‑Thr
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base T (8/11);  total (9/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

AAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGC‑‑GCATTGTTGAGCACAATGATGT‑‑GAAAAAGTGTGCTA‑‑ATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGTGACCAATAAAACAAATTATGCAATTTTTTAGTTGCATGAACTCGCATGTCTCCATAGAATGCGCGCTACTTGATGCCGACTTAGCTCAGTAGGTAGAGCAACTGACTTGTAATCAGTAGGTCACC  >  NC_010473/4272780‑4273157
                                                                                                                                                                                                        |                                                                                                                                                                                               
aaaGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGC‑‑GCATTGTTGAGCACAATGATGTTGaaaaa                                                                                                                                                                                          <  1:9968S42/200‑1 (MQ=255)
      atatTTGCCATGTCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGGATGAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAAGAATTATGGTTTAGCAGGAGC‑‑GCATTGTTGAGCACAATGATGTTGAAAAAgtgt                                                                                                                                                                                     >  1:17671S7/1‑200 (MQ=255)
                  gcgcGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGC‑‑GCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTA‑‑ATCTGCCCTc                                                                                                                                                                     <  1:30243S69/200‑1 (MQ=255)
                             ttCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGC‑‑GCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTA‑‑ATCTGCCCTCCGTTCGGCTg                                                                                                                                                          <  1:14883S16/200‑1 (MQ=255)
                               cTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGC‑‑GCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTA‑‑ATCTGCCCTCCGTTCGGCTGttt                                                                                                                                                       <  1:8226S50/200‑1 (MQ=255)
                                     aaTGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGAGCATTGTTGAGCACAGTGATGTTGAAAAAGCGTGCTA‑‑ATCTGCCCTCCGTTGGGGCTGTttctt                                                                                                                                                    <  1:24799S28/200‑1 (MQ=255)
                                             aaTGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGC‑‑GCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTA‑‑ATCTGCCCTCCGTTCGGCTGTTTCT‑‑ATCGTGTCGCAt                                                                                                                                       <  1:26655S35/200‑1 (MQ=255)
                                              aTGGGCAGGAAGGGCAACACCAGATGGGGAGGAATAAT‑‑‑‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGC‑‑GCATTGTTGAGCACAATGATGTTGAAAAAGTGGTGCGACTTTCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCa                                                                                                                                        <  1:1252S29/200‑1 (MQ=255)
                                                  gCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGC‑‑GCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTA‑‑ATCTGCCCTCCGTTCGGCTGTTCTCATCGTGTCGCATAAAAt                                                                                                                                  >  1:8044S8/1‑200 (MQ=255)
                                                                                           aaaaaTCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGC‑‑GCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTA‑‑ATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGTGACCAATAAAACAAATTATGCAATTTTTTAGtt                                                                                               >  1:11559S6/1‑200 (MQ=255)
                                                                                             aaaTCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGC‑‑GCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTA‑‑ATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGTGACCAATAAAACAAATTATGCAATTTTTTAGTTGc                                                                                             >  1:18082S7/1‑200 (MQ=255)
                                                                                                  tctGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGC‑‑GCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTA‑‑ATCTGCCCTCCGTTCGGCTGT‑‑CTTCATCGTGTCGCATAAAATGTGACCAATAAAACAAATTATGCAATTTTTTAGTTGCATGAACt                                                                                      <  1:2780S107/200‑1 (MQ=255)
                                                                                                                 tCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGC‑‑GCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTA‑‑ATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGTGACCAATAAAACAAATTATGCAATTTTTTAGTTGCATGAACTCGCATGTCTCCAt                                                                        <  1:2520S55/200‑1 (MQ=255)
                                                                                                                                          ttCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGC‑‑GCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTA‑‑ATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGTGACCAATAAAACAAATTATGCAATTTTTTAGTTGTATGAACCGCATGTCTCCATAGAATGCGCGCTACTTGATGCCGACt                                              <  1:24687S21/200‑1 (MQ=255)
                                                                                                                                                  gATAAAAAGAATTATGGTTTAGCAGGAGC‑‑GCATTGTTGAGCACAATGATGT‑‑TGAAAAGTGTGCTA‑‑ATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGTGACCAATAAAACAAATTATGCAATTTTTTAGTTGCATGAACTCGCATGTCTCCATAGAATGCGCGCTACTTGATGCCGACTTAGCTCAg                                      <  1:27321S21/200‑1 (MQ=255)
                                                                                                                                                     aaaaaGAATTATGGTTTAGCAGGAGC‑‑GCATTGTTGAGCACAATGATGTTTGAAAAAGTGTGCTG‑‑ATCTGCCCTCCGTTCGGCTGTTTTGCATCGTGTCGCATAAAATGTGACCAATAAAACAAATTATGCAATTTTTTAGTTGCATGAACTCGCATGTCTCCATAGAATGCGCGCTACTTGATGCCGACTTAGCTCAgta                                    >  1:25383S4/1‑200 (MQ=255)
                                                                                                                                                       aaaGAATTATGGTTTAGCAGGAGC‑‑GCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTA‑‑ATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGTGACCAATAAAACAAATTATGCAATTTTTTAGTTGCATGAACTCGCATGTCTCCATAGAATGCGCGCTACTTGATGCCGACTTAGCTCAgtagg                                  >  1:5607S29/1‑200 (MQ=255)
                                                                                                                                                        aaGAATTATGGTTTAGCAGGAGC‑‑GCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTA‑‑ATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGTGACCAATAAAACAAATTATGCAATTTTTTGACCATGAACTCAGC‑‑GTCTCCATAGAATGCGCGCTACTTGATGCCGACTTAGCTCAGTAGGTaga                              >  1:11569S6/1‑200 (MQ=255)
                                                                                                                                                                aTGGTTTAGCAGGAGC‑‑GCATTGTTGAGCACAATGATGT‑‑TGAAAAGTGTGCTA‑‑ATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGTGACCAATAAAACAAATTATGCAATTTTTTAGTTGCATGAACTCGCATGTCTCCATAGAATGCGCGCTACTTGATGCCGACTTAGCTCAGTAGGTAGAGCAAc                         <  1:27180S13/200‑1 (MQ=255)
                                                                                                                                                                      tAGCAGGAGC‑‑GCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTA‑‑ATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGTGACCAATAAAACAAATTATGCAATTTTTTAGTTGCATGAACTCGCATGTCTCCATAGAATGCGCGCTACTTGATGCCGACTTAGCTCAGTAGGTAGAGCAACTGACt                    <  1:6455S41/200‑1 (MQ=255)
                                                                                                                                                                                   cATTGTTGAGCACGATGATAT‑‑TAAAAAGTGTGTAGATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGTGACCAATAAAACAAATTATGCAATTTTTTAGTTGCATGAACTCGCATGTCTCCATAGAATGCGCGCTACTTGATGCCGACTTAGCTCAGTAGGTAGAGCAACTGACTTGTAATCAGTAg        >  1:27996S22/1‑200 (MQ=255)
                                                                                                                                                                                          gAGCACAATGATGTTGAAAAAGTGTGCTA‑‑ATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGTGACCAATAAAACAAATTATGCAATTTTTTAGTTGCATGAACTCGCATGTCTCCATAGAATGCGCGCTACTTGATGCCGACTTAGCTCAGTAGGTAGAGCAACTGACTTGTAATCAGTAGGTCAcc  >  1:3206S6/1‑200 (MQ=255)
                                                                                                                                                                                                        |                                                                                                                                                                                               
AAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGC‑‑GCATTGTTGAGCACAATGATGT‑‑GAAAAAGTGTGCTA‑‑ATCTGCCCTCCGTTCGGCTGTTTCTTCATCGTGTCGCATAAAATGTGACCAATAAAACAAATTATGCAATTTTTTAGTTGCATGAACTCGCATGTCTCCATAGAATGCGCGCTACTTGATGCCGACTTAGCTCAGTAGGTAGAGCAACTGACTTGTAATCAGTAGGTCACC  >  NC_010473/4272780‑4273157

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 48 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: