Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,953,5500AC56.7% 21.5 / 42.2 30T25P (ACC→CCC) ilvDdihydroxyacid dehydratase
Reads supporting (aligned to +/- strand):  ref base A (1/12);  new base C (17/0);  total (18/12)
Fisher's exact test for biased strand distribution p-value = 1.50e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.53e-06

CAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCC‑AATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTG  >  NC_000913/3953408‑3953688
                                                                                                                                              |                                                                                                                                           
cAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATg                                                                                                                                     <  1:1431/151‑1 (MQ=255)
   taaataCAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGAcc                                                                                                                                  <  1:316547/151‑1 (MQ=255)
         cAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCCTGCGCTGTGGCGCGCCCCCCGGATGACCGACGc                                                                                                                             >  2:553854/1‑150 (MQ=255)
              aaTGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGAttt                                                                                                                       <  2:615222/151‑1 (MQ=255)
                      ggCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGTCGGGCGCTGTGGCGCGCCCCCCGGAGGACCGACCCCGAGTTCGGTAAgg                                                                                                               >  1:389320/1‑150 (MQ=255)
                      ggCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCGGGCGCCGTTGGGCGCCCCCCGGCTGACCGACGCCgcgttctggaagg                                                                                                               >  1:297305/1‑139 (MQ=255)
                      ggCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCCCCGGCGCCCTGGCGCGCCCCCGGCAAGAACGCCCCCGCGTTTGGTAAg                                                                                                                >  1:614294/1‑150 (MQ=255)
                           gcacCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCCCGTGCGCTGTGGCGCGCCCCCGGGATGACCGACGCCGATTTTGGTTAGGCGGtt                                                                                                          >  2:633961/1‑151 (MQ=255)
                                          cGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAAc                                                                                           <  1:552997/151‑1 (MQ=255)
                                                 gacaCTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGTCGTGCGCTGGGGCGCGCCCCCGGGAGGAACGACGCCGGTTTTGGTTAGCCGGTTATCGCGGGTGGGGAATCGTTc                                                                                     >  1:153656/1‑150 (MQ=255)
                                                  acacTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGtt                                                                                      <  1:654612/148‑1 (MQ=255)
                                                   cacTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGGGCGCCCCCCGGCTGACCCACCCCGATTTCGGGGAGCCGGTTATTGCGGTTTTGGAATCGTTCAc                                                                                   >  2:46673/1‑150 (MQ=255)
                                                         gAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGTCGTGCGCTGTGGGGCGCCCCCCGGAAGACCGCCCCCCCTTTCTGTAAGCCGAGTATCGCGGGTGTGGAATCCTTCcccccaatt                                                                            >  1:72306/1‑146 (MQ=255)
                                                             aaataaaGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCC‑AATTTGTAc                                                                       <  2:488108/151‑1 (MQ=255)
                                                                    gTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGTCCTTGGCCGTGGCGCGCCCCCCGAAAGGACGACCCCGCGTTCGGTAAGGCGATTATCGCGGGTGGGGACTCGCTCCCCC‑CATTTTTACCGGGGCa                                                                >  1:452898/1‑151 (MQ=255)
                                                                         ccTAAGTACCGTTCCGCCGCCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCC‑AATTTGTACCGGCTCACGTc                                                            <  1:416237/150‑1 (MQ=255)
                                                                         ccTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGGTCGTGGGCCGTGGGGGGCCCCCGGAATGAACGACGCCGATTTTGGGGAGCCGGTTATTTCGGGTGTGGAATCGGTTACCC‑CAATTTTACCGGGGCGCGTc                                                            >  1:657779/1‑150 (MQ=255)
                                                                         ccTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGGTCGTGCGCTGTGGCGCGCCCCCCGGATGAACCACCCCCCTTTTGGTAAGCCGATTTTTGCGGTTGTGAACTCCTTTCCCC‑CATTTTTACCGGGGCACGTc                                                            >  2:532925/1‑150 (MQ=255)
                                                                         ccTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCCTGCGCTGTGGCGCGCCCCCCGGAAGACCCACGCCCGTTTCGGTTAGCCGGGTATCGCGGGTGTGGAATCGTTCTCCC‑CATTTTTAACGGGGCACGTc                                                            >  1:207465/1‑150 (MQ=255)
                                                                         ccTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCCCGCGCTGCGGCGCGCCCCCCGGACGACCGACGCCGCTTTCGGTAAGCCGATTATCGCGGCTGTGGACACCTTCACCC‑AAATTGTACCggggtacacc                                                            >  2:386450/1‑143 (MQ=255)
                                                                         ccTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGAGGCTCGGGCGCTGGGGCGCGCCCCCCGGATGACCGACGCCGGCTTTGGGGAGCCGGTTATCGCGGGTGTGAAATCGCTCCCCCAAATTTTTCCCgggcccggc                                                             >  1:294062/1‑144 (MQ=255)
                                                                          cTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCC‑AATTTGTACCGGGTCACGTCCa                                                          <  2:14850/151‑1 (MQ=255)
                                                                              gTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGTCGTGCGCCGTGGCGCGCCCCCGGAATGACCGACGCCCGTTTTGGTAAGCCGGTTTTCCCGGGTGGGAAATCGTTCACCC‑CACTTTTTCCgggggcccgccccct                                                       >  1:181914/1‑138 (MQ=255)
                                                                               tACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCC‑AATTTGTACCGGGTCACGTCCATCTgc                                                     <  1:368730/151‑1 (MQ=255)
                                                                                               ccacTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCC‑AATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACt                                     <  2:75130/151‑1 (MQ=255)
                                                                                               ccacTCATGGTCGTAATATGGCGGGGGCTCCTGCGCTGTGGCGCGCCCCCGGAATAACCGACCCCGCTTTTTGGAAGACGGTTATCGCGGGTGTGGAATCCTTCACCC‑AATTTTTATCCGGGCACCTCCCTCTGCGCGATCTCggggaa                                       >  2:485104/1‑145 (MQ=255)
                                                                                                cacTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCC‑AATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTg                                    <  2:27413/151‑1 (MQ=255)
                                                                                                 acTCATGGTCGTAATATGGCGGGGGGTCGGGCGCTGTGGCGCGCCCCCCGGAAGAACGAAGCCGCTTTTGGGGAGGCGGTTTTCGCGGGTGTGGAATCGGTCACCC‑CCCTTTTTTCGGGTCCCGGCCATCCGCGCGGtctccgggaactcg                                   >  1:55695/1‑141 (MQ=255)
                                                                                                                  ggCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCC‑AATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGaa                  >  2:455836/1‑151 (MQ=255)
                                                                                                                                  cTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCC‑AATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTg  <  2:292309/151‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                           
CAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCC‑AATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTG  >  NC_000913/3953408‑3953688

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