Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,206,4350CT37.0% 31.1 / 16.0 27intergenic (‑179/‑28)ttdR/ttdAtranscriptional activator of ttdABT/L‑tartrate dehydratase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (1/16);  new base T (10/0);  total (11/16)
Fisher's exact test for biased strand distribution p-value = 1.30e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAG‑ATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGA  >  NC_000913/3206236‑3206652
                                                                                                                                                                                                        |                                                                                                                                                                                                                          
tAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGa                                                                                                                                                                                      <  1:665323/238‑1 (MQ=255)
                aGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAgg                                                                                                                                                                     <  2:637074/239‑1 (MQ=255)
                                cGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTTCA‑ACAAAAAAAAAAAAATGAGGTGGCATGATGAGGGGCAGAAAAAAAAAACAAGCAGTGAATAAgaagga                                                                                                                                                     >  1:254797/1‑234 (MQ=255)
                                               cgcAGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AAATTTTTCCGACAAAAACAAAAAATAGAGTTTGCATGATGAGCGGAAGTAATAAACAACAAGCAGCGGAAAAGATTATAGAG‑AATATCGCtga                                                                                                                                     >  2:95479/1‑237 (MQ=255)
                                                            gCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AAATTTTTTCGACAAAAAAAAAAAATTGAGTTGGCATGATGAGGTAAATAAATAAGAAAAAAGCAGGGCATAAGTTTACATAGAATTGT‑GCTCACTTTAACACCact                                                                                                                        >  2:533672/1‑237 (MQ=255)
                                                                     tGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCg                                                                                                                                                                                       <  1:28921/168‑1 (MQ=255)
                                                                             ccTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAG‑Att                                                                                                                                             <  1:391733/201‑1 (MQ=255)
                                                                               tGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGGG‑AATTTTTTTAAAAAAATACAAAAAAAGGAGGTGGCCTGATTAGGGGAAGaaataa                                                                                                                                                                             >  2:28921/1‑168 (MQ=255)
                                                                                       cgcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTTTAGACAAAAACAAAAAAAAGAGGTGTCATGATGAGGGAAAGAAAAAAACAAAAGACAGAGAAAAAGATTACAGAG‑ATAGTCGCTAACt                                                                                                                                   >  2:391733/1‑201 (MQ=255)
                                                                                                                   aTTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAgag‑a                                                                                                                                               <  1:471564/161‑1 (MQ=255)
                                                                                                                             tGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AAATTTTTCAGACAAAAAAAAAAAAATGAGGTGGCCAGATGAGGGGAAGGAAAAAAAAAAAAGCAGGGAATAAAGTGACAGAG‑AGTGTctgtta                                                                                                                                     >  2:471564/1‑156 (MQ=255)
                                                                                                                                   ttACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTCCAGACAAATACAAAAACTGGAGTTGCCatga                                                                                                                                                                                             <  1:224545/100‑1 (MQ=255)
                                                                                                                                             aCGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTTTACACAAAAAAAAAAAATGGAGTTGCCATGATGAGGGGAAg                                                                                                                                                                                   >  2:224545/1‑100 (MQ=255)
                                                                                                                                                       aaaaTAAAAAACCACATTTTGAGTGGTAGTTGGGGGCGGTA‑AATTTTTCCAAAAAAATACAAAAACTGGAGTTGCCATGATGAGCAAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAG‑ATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATa                             <  1:95479/239‑1 (MQ=255)
                                                                                                                                                         aaTATAAACCCACATTTTTAGTGGTAGTTTGTGGCGGTA‑AATTTTCCCAGAAAAATACAAAACCTGGAGTTGCAATGATGAGCAAAAGTAATAAGCAACAGGCAGTAAATAAGTTGACAGAG‑ATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACg                           <  1:533672/239‑1 (MQ=255)
                                                                                                                                                                     caTTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTCCAGACa                                                                                                                                                                                                                      <  1:141990/41‑1 (MQ=255)
                                                                                                                                                                     caTTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTCCAGACa                                                                                                                                                                                                                      >  2:141990/1‑41 (MQ=255)
                                                                                                                                                                      aTTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTCCAGACAAATACAAAACCTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAG‑ATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGaa                                                            <  1:100425/193‑1 (MQ=255)
                                                                                                                                                                      aTTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGta                                                                                                                                                                                 <  2:638755/77‑1 (MQ=255)
                                                                                                                                                                      aTTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAG‑ATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAgg                                                                    <  1:735648/185‑1 (MQ=255)
                                                                                                                                                                            aGTGGTAGTTTGTGGCGGTG‑AATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAG‑ATTGTCGCt                                                                                                                                       <  2:733392/112‑1 (MQ=255)
                                                                                                                                                                              tGGTAGTTTGTGGCGGTA‑AATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCAAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAG‑ATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCtgatg                                                                                                  <  1:715933/147‑1 (MQ=255)
                                                                                                                                                                                gTAGTTTGTGGCGGTG‑AATTTTTTCCGACAAAAACAAAAAATGGGGGTGGCCTGATGAGCGAAAATAAAAAacaaca                                                                                                                                                                       >  1:638755/1‑77 (MQ=255)
                                                                                                                                                                                gTAGTTTGTGGCGGTG‑AATTTTTTCCGACAAAAAAAAAAAATGGGGGTGGCCTGATTAGGGAAAAGAATAAACAAAAAGCAGGGAATAAGATTACAGAG‑ATTGTCGCCTACTTTTACGCCATGAATTTTAACCGAATGCCTGATGACGTGGTGGGGTAACTAAAAAAACTCAAAGAGGCCGaaaac                                                         >  2:735648/1‑184 (MQ=255)
                                                                                                                                                                                   tttggggacgttG‑AATTTTCCAAGACAAATAAAAAAGCTGGAGTGCCAATGAAGGGGGAAATTAATAACCAACAGGAATTAAATAAGTGGACAGAG‑TTTGTCGCTAACTTAACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGa  <  2:254797/227‑1 (MQ=255)
                                                                                                                                                                                   gTTTGTGGCGGTG‑AATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAG‑ATTGTCGCt                                                                                                                                       <  1:53071/105‑1 (MQ=255)
                                                                                                                                                                                      tgtgGCGGTG‑AATTTTTTCAGACAAAAACAAAAAATGGAGTTGTCATGATGAGCGAAAAAAATAAACAACAAGCAGTGAATAAGGTGACAGAG‑ATTGTCGCTAACTTTACCg                                                                                                                             >  1:733392/1‑112 (MQ=255)
                                                                                                                                                                                                        |                                                                                                                                                                                                                          
TAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTG‑AATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAG‑ATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAACGTCGTCGATGGGGAAAATTATCTACCATACGATGTTCGACAACATGCAAAAAGCGA  >  NC_000913/3206236‑3206652

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑