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

AAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACG‑CC‑GATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCA‑CCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGG  >  NC_000913/3953419‑3953676
                                                                                                                                   |                                                                                                                                 
aaaaaTGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCGTGCGCCGTGGCGCGCCCCCCGGAAGACCCACG‑CC‑g                                                                                                               >  2:45480/1‑150 (MQ=255)
       ggACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGTCGTGCGCCGTGGGGCGCCACCCGGAAGAACCACG‑CC‑gcattttgg                                                                                                       >  1:54452/1‑143 (MQ=255)
           ggCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGGTCGTTCGCTGTGGCGGGCCCCCCGGAAGAACGCCG‑CC‑GATTTTGGTAAgg                                                                                                   >  1:250354/1‑150 (MQ=255)
           ggCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCGTGCGCTGTGGGGCGCCCCCGGGATGACCGACG‑CC‑GCTTTCGGTGAg                                                                                                    >  2:240588/1‑150 (MQ=255)
                     gTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGGTCGGGCGCTGTGGGGCGCCCCCCGCATGACCGGCC‑CC‑GGTTTTGGGGAGCCggatttcg                                                                                          >  1:264959/1‑143 (MQ=255)
                                       acacTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGCCGTTCCCTGTGGCGCGCCCCCGGCATTACCAACG‑CC‑GGCTTTGGTGAGCCCGTTATTGCGGTtgggggctccgtca                                                                        >  2:217852/1‑138 (MQ=255)
                                                 taaataaaGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACG‑CC‑GATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCA‑CCCAATTTg                                                              <  2:138293/149‑1 (MQ=255)
                                                  aaataaaGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACG‑CC‑GATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCA‑CCCAATTTGTAc                                                           <  1:73488/151‑1 (MQ=255)
                                                    ataaaGTATGCCTAAGTATCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGTCGTGCGCCGGGGCGCGCCCCCCGGAAGACCGACC‑CC‑CACTTTGGTTAGCCGAGTATCGCCGGGGTGGACTCGTTCACCCCAAtttttccg                                                          >  2:260709/1‑146 (MQ=255)
                                                              ccTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGTGTCCTGCGCTGTTGCGCGCCCCCGGGAATACCCACC‑CC‑CGTTTTGGTAAGCCCATTATACCCGCTGTGGACTCGTTCA‑CCCAAttttttacccgtcccgcg                                                >  1:255386/1‑135 (MQ=255)
                                                              ccTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCGCGTTCGCCGTGGCGCGCCCCCCGGAAGAACCACG‑CC‑GCTTTTGGTAAGCCGGTTATCGCGGGTGTGGACCCGTTCT‑CCCCAATTTTACCGGGGCACgcc                                                >  2:94109/1‑148 (MQ=255)
                                                              ccTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCCCGTTCGCCGTGGGGCGCCCCCCGGAGGAACGACG‑CCGGGTTTTGGT‑AGCCGGGTATTGCGGGTGTGGACTCCTTCC‑CCCCATTTTTACCggggggcgtc                                                >  2:94511/1‑143 (MQ=255)
                                                                 aaGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACG‑CC‑GATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCA‑CCCAATTTGTACCGGGTCACGTCCATc                                            <  1:59510/151‑1 (MQ=255)
                                                                 aaGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACG‑CC‑GATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCA‑CCCAATTTGTACCGGGTCACGTCCATc                                            <  1:307924/151‑1 (MQ=255)
                                                                 aaGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACG‑CC‑GATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCA‑CCCAATTTGTACCGGGTCACGTCCATc                                            <  1:288633/151‑1 (MQ=255)
                                                                         ttCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGTCGTGCGCCGGGGGGCGCCCCCGGGAAGACCCACGCCC‑GATTTGGGTAGGCGGATT‑TTGCGGTTGGGGAATTGTTCA‑CCCCATTTTTTCCGGGGCACGGCCCCCTTcgcggg                                    >  2:7158/1‑149 (MQ=255)
                                                                                 ccaccacTCATGGTCGTAATATGGCGGGGGGTCGTGCGCTGTGGCGCGCCCCCCGGATGACCCACG‑CC‑GCTTTTTGTAAGCCGGTTATCGCGGGTGGGGACTCGTTCC‑CCCCCTTTTTTCCCGGGCACGGCCCTCTTCGCGGtctccggga                            >  2:179581/1‑146 (MQ=255)
                                                                                                            gTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACG‑CC‑GATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCA‑CCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCgg  >  2:314517/1‑150 (MQ=255)
                                                                                                                                   |                                                                                                                                 
AAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACG‑CC‑GATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCA‑CCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGG  >  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: ‑