Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,207,549 Δ1 bp coding (179/606 nt) ttdB → L(+)‑tartrate dehydratase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,207,5490G.96.4% 119.4 / ‑1.5 28coding (179/606 nt)ttdBL(+)‑tartrate dehydratase subunit beta
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base . (12/15);  total (12/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-01

GCCGATCTCACCTTTGAAAATCTGTCTCACACCCGGAGCGCGTTATGAAAAAGATCCTGACAACCCCGATCAAAGCTGAAGATCTGCAAGATATTCGCGTCGGCGATGTGATCTACCTGACCGGTACGCTGGTGACCTGCCGCGACGTTTGTCACCGCCGTTTGATCGAACTGAAACGTCCGATCCCTTACGATCTCAACGGCAAAGCGATTTTCCACGCTGGCCCCATCGTGCGCAAAAACGGCGACAAATGGGAGATGGTCTCCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCAAAGGTGGTATGGGGCCGCTGACCGAAGAAGGCTGCCAGAAATTCAAGGCGCTACATGTGATTTTCCCGGCAGGCTGCGCGGTGCTGGCGGCAACCCAGGT  >  NC_000913/3207327‑3207774
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                 
cccGATCTCCCCTTTGAAAAGCTGTCTCACCCCGGGAGGGCGTTATGAAAAAGATCCTTACAACCCCGATCAAAGCTGAAGATCTGCAAGATATTCGCGTCGGCGATGTGATCTACCTGACCGGTACGCTGGTGACCTGCCGCGACGTTTGTCACCGCCGTTTGATCGAACTGAAACGTCCGATCCCTTACGATCTCAACGGCAAAGCGATTTTCCACGCTGCCCCATCGTGCGCAAAAACGGCGACAAAt                                                                                                                                                                                                      <  2:198223/250‑1 (MQ=255)
    aTCTCACCTTTGAAAATCTGTCTCACACCCGGAGCGCGTTATGAAAAAGATCCTGACAACCCCGATCAAAGCTGAAGATCTGCAAGATATTCGCGTCGGCGATGTGATCTACCTGACCGGTACGCTGGTGACCTGCCGCGACGTTTGTCACCGCCGTTTGATCGAACTGAAACGTCCGATCCCTTACGATCTCAACGGCAAAGCGATTTTCCACGCTGCCCCATCGTGCGCAAAAACGGCGACAAATggg                                                                                                                                                                                                   >  1:51001/1‑250 (MQ=255)
          ccTTTGAAAATCTGTCTCACACCCGGAGCGCGTTATGAAAAAGATCCTGACAACCCCGATCAAAGCTGAAGATCTGCAAGATATTCGCGTCGGCGATGGGATCTACCGGACCGGTACGCGGGTGACCTGCCGCGACGTTTGTCACCGCCGTTTGATCGAACTGAAACGGCCGATCCCTTACGATCTCAACGGCAAAGCGATTTTTCACGCGGCCCCGTCGTGCGCAAAAACGGCGACAAATGGGAGATGGt                                                                                                                                                                                            <  1:23463/251‑1 (MQ=255)
                       gTCTCACACCCGGAGCGCGTTATGAAAAAGATCCTGACAACCCCGATCAAAGCTGAAGATCTGCAAGATATTCGCGTCGGCGATGTGATCTACCTGACCGGTACGCTGGTGACCTGCCGCGACGTTTGTCACCTCCGTTTGATCGAACTGAAACGTCCGATCCCTTACGATCTCAACGGCAAAGCGATTTTCCACGCTGCCCCATCGTGCGCAAAAACGGCGACAAATGGGAGATGGTCTCCGTCGGCCCg                                                                                                                                                                               <  1:82177/251‑1 (MQ=255)
                                        gtttATGAAAAAGATCCTGACACCCCCGACCAAAGCTAAAGATCTGCAAGATATGCCCGCCGGGGATGGGACCTACCTGCCCGGTACGCGGGTGACCTGCCGCGACGTTGGTCACCTCCGTTTGATCGAACTGAAACGTCCGATCCCTTACGATCTCAACGGCAAAGCGATTTTCCACGCTGCCCCATCGTGCGCAAAAACTGCGACAAATGGGAGATGGTCTCCGTCGGCCCGACAACCAGTATGCGTAt                                                                                                                                                              <  2:99741/249‑1 (MQ=255)
                                                          gACAACCCCGATCAAAGCTGAAGATATGCAAGATATTCGCGTCGGCGATGTGATCTACCTGACCGGTACGCTGGTGACCTGCCGCGACGTTTGTCACCGCCGTTTGATCGAACTGAAACGTCCGATCCCTTACGATCTCAACGGCAAAGCGATTTTCCACGCTGCCCCATCGTGCGCAAAAACGGCGACAAATGGGAGATGGTCTCCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTg                                                                                                                                             <  1:76359/250‑1 (MQ=255)
                                                                            tGAAGATCTGCAAGATATTCGCGTCGGCGATGTGATCTACCTGACCGGTACGCTGGTGACCTGCCGCGACGTTTGTCACCGCCGTTTGATCGAACTGAAACGTCCGATCCCTTACGATCTCAACGGCAAAGCGATTTTCCACGCTGCCCCATCGTGCGCAAAAACGGCGACAAATGGGAGATGGTCTCCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCgg                                                                                                                          >  2:72848/1‑251 (MQ=255)
                                                                                        aGATATTCGCGTCGGCGATGTGATCTACCTGACCGGTACGCTGGTGACCTGCCGCGACGTTTGTCACCGCCGTTTGATCGAACTGAAACGTCCGATCCCTTACGATCTCAACGGCAAAGCGATTTTCCACGCTGCCCCATCGTGCGCAAAAACGGCGACAAATGGGAGATGGTCTCCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACtggg                                                                                                              >  2:31834/1‑250 (MQ=255)
                                                                                               ccgtgcggcattgtgATCTCCCGGACCGGTCCGCTGCTGCCCTGCCGCGACGGTTGTCACCGCCGTTTGATCGAACTGAAACGTCCGATCCCTTACGATCTCAACGGCAAAGCGATGTTCCACGCTGCCCCATCGTGCGCAAAAACGGCGACAAATGGGAGATGGTCTCCGTCGGCCCGACAACCATTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGc                                                                                                       <  2:109238/240‑1 (MQ=255)
                                                                                                 ggTCGGCGAGGGGATCTACCGGACGGGTACGCGGGTTACCGTCCGAGACGTTTGTCACCGCCGTATGATCGAACTGAAACGTCCGATCCCTTACGATCTCAACGGCAAAGCGATTTTCCACGCTGCCCCATCGTGCGCAAAAACGGCGACAAATGGGAGATGGTCTCCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAAACCGGCGTGAAACTGGTGGTTGGCaa                                                                                                     <  2:189603/250‑1 (MQ=255)
                                                                                                 cgTCGGCGATGTGATCTACCTGACCGGTACGCTGGTGACCTGCCGCGACGTTTGTCACCGCCGTTTGATCGAACTGAAACGTCCGATCCCTTACGATCTCAACGGCAAAGCGATTTTCCACGCTGCCCCATCGTGCGCAAAAACGGCGACAAATGGGAGATGGTCTCCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGGTGGCa                                                                                                      >  1:192406/1‑250 (MQ=255)
                                                                                                         aTGTGATCTACCTGACCGGTACGCTGGTGACCTGCCGCGACGTTTGTCACCGCCGTTTGATCGAACTGAAACGTCCGATCCCTTACGATCTCAACGGCAAAGCGATTTTCCACGCTGCCCCATCGTGCGCAAAAACGGCGACAAATGGGAGATGGTCTCCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCAAAggtggt                                                                                              >  1:109432/1‑250 (MQ=255)
                                                                                                            tgATCTACCTGACCGGTACGCTGGTGACCTGCCGCGACGTTTGTCACCGCCGTTTGATCGAACTGAAACGTCCGATCCCTTACGATCTCAACGGCAAAGCGATTTTCCACGCTGCCCCATCGTGCGCAAAAACGGCGACAAATGGGAGATGGTCTCCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCAAAGGTGGTATgg                                                                                          >  1:178516/1‑251 (MQ=255)
                                                                                                                 tACCGGACCGGTACCCTGGTGCCCTGCCGCGCCGTTTGTCCCCGCCGTTTGATCGAACTGAAACGTCCGATCCCTTACGATCTCAACGGCAAAGCGATTTTCCACGCTGCCCCATCGTGCGCAAAAACGGCGACAAATGGGAGATGGTCTCCGTCGGCCCGACAACCAGTATGTGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCAAAGGTGGTATGGGGcc                                                                                      <  2:123941/250‑1 (MQ=255)
                                                                                                                 acccTGACCGGTACGCTGGTGCCCTGCCGCTACGTTTGTCACCGCCGTTTGATCGAACTGAAACGTCCGATCCCTTACGATCTCAACGGCAAAGCGATTTTCCACGCTGCCCCATCGTGCGCAAAAACGGCGACAAATGGGAGATGGTCTCCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCAAAGGTGGTATGGGGCCg                                                                                     <  2:123684/249‑1 (MQ=255)
                                                                                                                             aCGCTGGTGACCTGCCGCGACGTTTGTCACCGCCGTTTGATCGAACTGAAACGTCCGATCCCTTACGATCTCAACGGCAAAGCGATTTTCCACGCTGCCCCATCGTGCGCAAAAACGGCGACAAATGGGAGATGGTCTCCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTTAAACTGGTGGTTGGCAAAGGTGGTATGGGGCCGCTGACCgaagaa                                                                         <  1:20464/251‑1 (MQ=255)
                                                                                                                                  ggTGACCTGCCGCGACGTTTGTCACCGCCGTTTGATCGAACTGAAACGTCCGATCCCTTACGATCTCAACGGCAAAGCGATTTTCCACGCTGCCCCATCGTGCGCaaa                                                                                                                                                                                                                   >  1:32496/1‑108 (MQ=255)
                                                                                                                                  ggTGACCTGCCGCGACGTTTGTCACCGCCGTTTGATCGAACTGAAACGTCCGATCCCTTACGATCTCAACGGCAAAGCGATTTTCCACGCTGCCCCATCGTGCGCaaa                                                                                                                                                                                                                   <  2:32496/108‑1 (MQ=255)
                                                                                                                                              cgACGTTTGTCACCGCCGTTTGATCGAACTGAAACGTCCGATCCCTTACGATCTCAACGGCAAAGCGATTTTCCACGCTGCCCCATCGTGCGCAAAAACGGCGACAAATGGGAGATGGTCTCCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCAAAGGTGGTATGGGGCCGCTGACCGAAGAAGGCTGCCAGAAATTCaa                                                        <  2:82530/251‑1 (MQ=255)
                                                                                                                                                  gTTTGTCACCGCCGTTTGATCGAACTGAAACGTCCGATCCCTTACGATCTCAACGGCAAAGCGATTTTCCACGCTGCCCCATCGTGCGCAAAAACGGCGACAAATGGGAGATGGTCTCCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCAAAGGTGGTATGGGGCCGCTGACCGAAGAAGGCTGCCAGAAATTCAAGgcg                                                    <  1:247851/251‑1 (MQ=255)
                                                                                                                                                   tttGTCACCGCCGTTTGATCGACCTGAAACGTCCGATCACTTACGATCTCAACGGCAAAGCGATTTTCCACGCTGCCCCATCGCGCGCACAAACGGCGACAAATGGGAGATGGTCTCCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCAAAGGTGGGATGGGGCCGCTGACCGAAGAAGGCTGCCAGAAATTCAAGgcgc                                                   >  2:226583/1‑251 (MQ=255)
                                                                                                                                                      gTCACCGCCGTTTGATCGAACTGAAACGTCCGATCCCTTACGATCTCAACGGCAAAGCGATTTTCCACGCTGCCCCATCGTGCGCAAAAACGGCGACAAATGGGAGATGGTCTCCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCAAAGGTGGGATGGGGCCGCTGACCGAAGAAGGCTGCCAGAAATTCAAGGCGCTAc                                                >  1:28408/1‑251 (MQ=255)
                                                                                                                                                                  tGATCGAACTGAAACGTCCGATCCCTTACGATCTCAACGGCATAGCGATTTTCCACGCTGCCCCATCGTGCGCAAAAACGGCGACAAATGGGAGATGGTCTCCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCAAAGGTGGTATGGGGCCGCTGACCGAAGAAGGCGGCCAGAAATTCAAGGCGCTACATGTGATTTTc                                     <  2:152477/250‑1 (MQ=255)
                                                                                                                                                                   gATCGAACTGAAACGTCCGATCCCTTACGATCTCAACGGCAAAGCGATTTTCCACGCTGCCCCATCGTGCGCAAAAACGGCGACAAATGGGAGATGGTCTCCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCAAAGGTGGTATGGGGCCGCTGACCGAAGAAGGCTGCCAGAAATTCAAGGCGCTACATGTGATTTTccc                                   >  1:250078/1‑251 (MQ=255)
                                                                                                                                                                    aTCGAACTGACACGTCCGATCCCTTACGATCTCAACGGCAAAGCGATTTTCCACGCTGCCCCATCGTGCGCAAAAACGGCGACAAATGGGAGATGGTCTCCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCAAAGGTGGTATGGGGCCGCTGACCGAAGAAGGCTGCCAGAAATTCAAGGCGCTACATGTGATTTTCCCg                                  >  2:155305/1‑251 (MQ=255)
                                                                                                                                                                                    cGATCCCTTACGATCTCAACGGCAAAGCGATTTTCCACGCTGCCCCATCGTGCGCAAAAACGGCGACAAATGGGAGATGGTCTCCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCAAAGGTGGTATGGGGCCGCTGACCGAAGAAGGCTGCCAGAAATTCAAGGCGCTACATGTGATTTTCCCGGCAGGCTGCGCGGTGc                  <  1:51324/251‑1 (MQ=255)
                                                                                                                                                                                             aCGATCTCAACGGCAAATCGATTTTCCACGCTGCCCCATCGTGCGCAACAACGGCGACAAATGGGAGATGGTCTCCTTCGGCCCGCCACCCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCAAAGGTGGTATGGGGCCGCTGACCGAAGAAGGCTGCCAGAAATTCAAGGCGCTACCTGTGATTTTCCCGGCAGGCTGCGCGGTGCTGGCGGCa          >  2:128349/1‑250 (MQ=255)
                                                                                                                                                                                                    cAACGGCAAAGCGATTTTCCACGCTGCCCCATCGTGCGCAAAAACGGCGACAAATGGGAGATGGTCTCCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCAAAGGTGGTATGGGGCCGCTGACCGAAGAAGGCTGCCAGAAATTCAAGGCGCTACATGTGATTTTCCCGGCAGGCTGCGCGGTGCTGGCGGCAACCCAGGt  <  1:104548/251‑1 (MQ=255)
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                 
GCCGATCTCACCTTTGAAAATCTGTCTCACACCCGGAGCGCGTTATGAAAAAGATCCTGACAACCCCGATCAAAGCTGAAGATCTGCAAGATATTCGCGTCGGCGATGTGATCTACCTGACCGGTACGCTGGTGACCTGCCGCGACGTTTGTCACCGCCGTTTGATCGAACTGAAACGTCCGATCCCTTACGATCTCAACGGCAAAGCGATTTTCCACGCTGGCCCCATCGTGCGCAAAAACGGCGACAAATGGGAGATGGTCTCCGTCGGCCCGACAACCAGTATGCGTATGGAAAGTTTTGAACGTGAATTTATTGAGCAGACCGGCGTGAAACTGGTGGTTGGCAAAGGTGGTATGGGGCCGCTGACCGAAGAAGGCTGCCAGAAATTCAAGGCGCTACATGTGATTTTCCCGGCAGGCTGCGCGGTGCTGGCGGCAACCCAGGT  >  NC_000913/3207327‑3207774

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