Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,951,2030AT68.4% 16.7 / 10.4 19I124I (ATA→ATTyecDisochorismatase family protein YecD
Reads supporting (aligned to +/- strand):  ref base A (1/5);  new base T (3/10);  total (4/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCGCCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGG  >  NC_000913/1950969‑1951447
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    
cAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCAGTATGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTGACGGTACGGATCTGGAGTTGCAATTACGCCGCCGAGGTATCGATACAATAGTGTTATGTGGGAtct                                                                                                                                                                                                                                      <  1:88254/251‑1 (MQ=255)
                                                                                                      ngAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCACTGGAGCGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCAGCCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGagcagc                                                                                                                                >  1:18763/2‑251 (MQ=255)
                                                                                                                tggCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCGCCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTTGCTTTAATCTGGTGATTGc                                                                                                                                                                   <  1:50241/206‑1 (MQ=255)
                                                                                                                                 aTTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCGCCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTTGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTAcc                                                                                                     <  1:162152/251‑1 (MQ=255)
                                                                                                                                                  gCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCGCCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCgg                                                                                                                                                                                                                         <  1:16067/118‑1 (MQ=255)
                                                                                                                                                                                                      cgatgtgcggttACGCCGCCGAGGTATCGATACAAGTGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATTATGTACATCGGTCTGCCGc                                 <  1:133693/240‑1 (MQ=255)
                                                                                                                                                                                                        nnnnnnggnnnnnccgncgAGGTATCGATACAATTGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGagc                                                                                                                                   <  2:105299/137‑1 (MQ=255)
                                                                                                                                                                                                             cggttACGCCGCCGAGTTATCGATACAATTGTGTTATGTGGGATCTCGACCAATAGCGGTGGTGAATCCTCCGCCCGCAATGCCTGGGAACTGGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAAGAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCTTAGCGGGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCg                         <  1:65340/248‑1 (MQ=255)
                                                                                                                                                                                                             cggttACGCCGCCGAGGTATCGATACAATTGTGTTATGTgg                                                                                                                                                                                                                                           <  1:64752/38‑1 (MQ=255)
                                                                                                                                                                                                               gttACGCCGCCGAGGTATCGATACAATTGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGc                        >  1:82473/2‑250 (MQ=255)
                                                                                                                                                                                                                   cgccgactaGGTATCGATACAATTGTGGTATGTGGGATCTCGACCAATATCGGTGTTGAATCCCCCGCCCGCAATGCCAGGGATCGTGGCTTTAAGCTGGTGATGGCCGAAGACGCCATTAGTGCCTCCCGTGCCGAGCAACACAATAACAGCATGAATCATAACGAGCCGCGCATCGACCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATc                                           <  1:176469/219‑1 (MQ=255)
                                                                                                                                                                                                                       gcnnAGGTATCGANACAATTGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGNAATGCCTGGNNACNTGGCTNTNNNCNNNTGNTTGCCGAAGACGCCTGGAGNGCNNCCAGCGNCGAGCAACACAATAACAgcattaatc                                                                                                              <  2:38737/156‑9 (MQ=255)
                                                                                                                                                                                                                       gccgAGGTATCGATACAATTGTGTTATGTGGGCTCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAg                                                                                                                                           >  1:550/1‑127 (MQ=255)
                                                                                                                                                                                                                       gccgAGGTATCGATACAATTGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAgcattaatc                                                                                                              >  1:38737/1‑148 (MQ=255)
                                                                                                                                                                                                                          gAGGTATCGATACAAGTGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGGGGAAGAGATCCTTCACGCGTTATGATTTACAGCGTTCTGCCGCAATGGTCGCATCCTAAATggn            <  1:18911/251‑2 (MQ=255)
                                                                                                                                                                                                                            ggTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTg          <  1:138799/251‑1 (MQ=255)
                                                                                                                                                                                                                                     tcAATTGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCGACCCGCTCATCGCCCGTGTGCGTAGCGGGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTggg  <  1:132636/249‑1 (MQ=255)
                                                                                                                                                                                                                                     aCAATTGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTTGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTggg  <  1:73708/250‑1 (MQ=255)
                                                                                                                                                                                                                                     aCAATTGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTggg  <  1:7045/250‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    
CAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCGCCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGG  >  NC_000913/1950969‑1951447

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