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

CAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTG  >  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)
                                                             aaataaaGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTAc                                                                       <  2:488108/151‑1 (MQ=255)
                                                                    gTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGTCCTTGGCCGTGGCGCGCCCCCCGAAAGGACGACCCCGCGTTCGGTAAGGCGATTATCGCGGGTGGGGACTCGCTCCCCCCATTTTTACCGGGGCa                                                                >  1:452898/1‑151 (MQ=255)
                                                                         ccTAAGTACCGTTCCGCCGCCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGCTCACGTc                                                            <  1:416237/150‑1 (MQ=255)
                                                                         ccTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGGTCGTGGGCCGTGGGGGGCCCCCGGAATGAACGACGCCGATTTTGGGGAGCCGGTTATTTCGGGTGTGGAATCGGTTACCCCAATTTTACCGGGGCGCGTc                                                            >  1:657779/1‑150 (MQ=255)
                                                                         ccTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGGTCGTGCGCTGTGGCGCGCCCCCCGGATGAACCACCCCCCTTTTGGTAAGCCGATTTTTGCGGTTGTGAACTCCTTTCCCCCATTTTTACCGGGGCACGTc                                                            >  2:532925/1‑150 (MQ=255)
                                                                         ccTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCCTGCGCTGTGGCGCGCCCCCCGGAAGACCCACGCCCGTTTCGGTTAGCCGGGTATCGCGGGTGTGGAATCGTTCTCCCCATTTTTAACGGGGCACGTc                                                            >  1:207465/1‑150 (MQ=255)
                                                                         ccTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCCCGCGCTGCGGCGCGCCCCCCGGACGACCGACGCCGCTTTCGGTAAGCCGATTATCGCGGCTGTGGACACCTTCACCCAAATTGTACCggggtacacc                                                            >  2:386450/1‑143 (MQ=255)
                                                                         ccTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGAGGCTCGGGCGCTGGGGCGCGCCCCCCGGATGACCGACGCCGGCTTTGGGGAGCCGGTTATCGCGGGTGTGAAATCGCTCCCCCAAATTTTTCCCgggcccggc                                                             >  1:294062/1‑144 (MQ=255)
                                                                          cTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCa                                                          <  2:14850/151‑1 (MQ=255)
                                                                              gTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGTCGTGCGCCGTGGCGCGCCCCCGGAATGACCGACGCCCGTTTTGGTAAGCCGGTTTTCCCGGGTGGGAAATCGTTCACCCCACTTTTTCCgggggcccgccccct                                                       >  1:181914/1‑138 (MQ=255)
                                                                               tACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTgc                                                     <  1:368730/151‑1 (MQ=255)
                                                                                               ccacTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACt                                     <  2:75130/151‑1 (MQ=255)
                                                                                               ccacTCATGGTCGTAATATGGCGGGGGCTCCTGCGCTGTGGCGCGCCCCCGGAATAACCGACCCCGCTTTTTGGAAGACGGTTATCGCGGGTGTGGAATCCTTCACCCAATTTTTATCCGGGCACCTCCCTCTGCGCGATCTCggggaa                                       >  2:485104/1‑145 (MQ=255)
                                                                                                cacTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTg                                    <  2:27413/151‑1 (MQ=255)
                                                                                                 acTCATGGTCGTAATATGGCGGGGGGTCGGGCGCTGTGGCGCGCCCCCCGGAAGAACGAAGCCGCTTTTGGGGAGGCGGTTTTCGCGGGTGTGGAATCGGTCACCCCCCTTTTTTCGGGTCCCGGCCATCCGCGCGGtctccgggaactcg                                   >  1:55695/1‑141 (MQ=255)
                                                                                                                  ggCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGaa                  >  2:455836/1‑151 (MQ=255)
                                                                                                                                  cTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTg  <  2:292309/151‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                           
CAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTG  >  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: