Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,953,5500AC61.1% 6.3 / 18.7 18T25P (ACC→CCC) ilvDdihydroxyacid dehydratase
Reads supporting (aligned to +/- strand):  ref base A (2/5);  new base C (11/0);  total (13/5)
Fisher's exact test for biased strand distribution p-value = 2.45e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.86e-03

AAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGG  >  NC_000913/3953419‑3953676
                                                                                                                                   |                                                                                                                                 
aaaaaTGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCGTGCGCCGTGGCGCGCCCCCCGGAAGACCCACGCCg                                                                                                               >  2:45480/1‑150 (MQ=255)
       ggACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGTCGTGCGCCGTGGGGCGCCACCCGGAAGAACCACGCCgcattttgg                                                                                                       >  1:54452/1‑143 (MQ=255)
           ggCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGGTCGTTCGCTGTGGCGGGCCCCCCGGAAGAACGCCGCCGATTTTGGTAAgg                                                                                                   >  1:250354/1‑150 (MQ=255)
           ggCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCGTGCGCTGTGGGGCGCCCCCGGGATGACCGACGCCGCTTTCGGTGAg                                                                                                    >  2:240588/1‑150 (MQ=255)
                     gTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGGTCGGGCGCTGTGGGGCGCCCCCCGCATGACCGGCCCCGGTTTTGGGGAGCCggatttcg                                                                                          >  1:264959/1‑143 (MQ=255)
                                       acacTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGCCGTTCCCTGTGGCGCGCCCCCGGCATTACCAACGCCGGCTTTGGTGAGCCCGTTATTGCGGTtgggggctccgtca                                                                        >  2:217852/1‑138 (MQ=255)
                                                 taaataaaGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTg                                                              <  2:138293/149‑1 (MQ=255)
                                                  aaataaaGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTAc                                                           <  1:73488/151‑1 (MQ=255)
                                                    ataaaGTATGCCTAAGTATCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGTCGTGCGCCGGGGCGCGCCCCCCGGAAGACCGACCCCCACTTTGGTTAGCCGAGTATCGCCGGGGTGGACTCGTTCACCCCAAtttttccg                                                          >  2:260709/1‑146 (MQ=255)
                                                              ccTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGTGTCCTGCGCTGTTGCGCGCCCCCGGGAATACCCACCCCCGTTTTGGTAAGCCCATTATACCCGCTGTGGACTCGTTCACCCAAttttttacccgtcccgcg                                                >  1:255386/1‑135 (MQ=255)
                                                              ccTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCGCGTTCGCCGTGGCGCGCCCCCCGGAAGAACCACGCCGCTTTTGGTAAGCCGGTTATCGCGGGTGTGGACCCGTTCTCCCCAATTTTACCGGGGCACgcc                                                >  2:94109/1‑148 (MQ=255)
                                                              ccTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCCCGTTCGCCGTGGGGCGCCCCCCGGAGGAACGACGCCGGGTTTTGGTAGCCGGGTATTGCGGGTGTGGACTCCTTCCCCCCATTTTTACCggggggcgtc                                                >  2:94511/1‑143 (MQ=255)
                                                                 aaGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATc                                            <  1:59510/151‑1 (MQ=255)
                                                                 aaGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATc                                            <  1:307924/151‑1 (MQ=255)
                                                                 aaGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATc                                            <  1:288633/151‑1 (MQ=255)
                                                                         ttCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGTCGTGCGCCGGGGGGCGCCCCCGGGAAGACCCACGCCCGATTTGGGTAGGCGGATTTTGCGGTTGGGGAATTGTTCACCCCATTTTTTCCGGGGCACGGCCCCCTTcgcggg                                    >  2:7158/1‑149 (MQ=255)
                                                                                 ccaccacTCATGGTCGTAATATGGCGGGGGGTCGTGCGCTGTGGCGCGCCCCCCGGATGACCCACGCCGCTTTTTGTAAGCCGGTTATCGCGGGTGGGGACTCGTTCCCCCCCTTTTTTCCCGGGCACGGCCCTCTTCGCGGtctccggga                            >  2:179581/1‑146 (MQ=255)
                                                                                                            gTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCgg  >  2:314517/1‑150 (MQ=255)
                                                                                                                                   |                                                                                                                                 
AAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGG  >  NC_000913/3953419‑3953676

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