Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,799,3741.T20.6% 81.9 / 28.9 34intergenic (‑104/+19)pheM/rplTphenylalanyl‑tRNA synthetase operon leader peptide/50S ribosomal subunit protein L20
Reads supporting (aligned to +/- strand):  ref base . (14/13);  new base T (3/4);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGC  >  NC_000913/1799236‑1799520
                                                                                                                                           |                                                                                                                                                  
aaGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGCTAGTCTCCCTCttt                                                                                                                                          <  1:184177/149‑1 (MQ=255)
  gTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGCTAGTCTCCCTCttt                                                                                                                                          <  6:9430/147‑1 (MQ=255)
  gTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGCTAGTCTCCCTCttt                                                                                                                                          >  5:9430/1‑147 (MQ=255)
    aaaaaaaGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAAc                                                                                                                                      >  2:37822/1‑149 (MQ=255)
                ggAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGc                                                                                                                            >  4:114045/1‑147 (MQ=255)
                   aaaTAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCagag                                                                                                                       <  5:109551/149‑1 (MQ=255)
                        gcagcaTTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAgctgct                                                                                                                  >  7:150805/1‑149 (MQ=255)
                            cattcattGCTTGCGTTACCTTTTGGTACTCTACAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGCTAGTCTCCCTCtt                                                                                                                                           >  8:45876/1‑120 (MQ=255)
                            cattcattGCTTGCGTTACCTTTTGGTACTCTACAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGCTAGTCTCCCTCtt                                                                                                                                           <  7:45876/120‑1 (MQ=255)
                               tcattGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCtt                                                                                                           >  8:135548/1‑149 (MQ=255)
                                  ttgcttgcGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCttttt                                                                                                        <  1:37822/149‑1 (MQ=255)
                                    gcttgcGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTc                                                                                                       <  1:133773/148‑1 (MQ=255)
                                    gcttgcGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTc                                                                                                       <  7:135548/148‑1 (MQ=255)
                                     cttgcGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCaa                                                                                                     <  8:150805/149‑1 (MQ=255)
                                          gTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCaa                                                                                                     <  8:49440/144‑1 (MQ=255)
                                          gTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCaa                                                                                                     >  7:49440/1‑144 (MQ=255)
                                                ttttGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGt                                                                                           >  3:7666/1‑148 (MQ=255)
                                                        ctctTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTAc                                                                                  >  5:176197/1‑149 (MQ=255)
                                                         tctTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACt                                                                                 >  6:111842/1‑149 (MQ=255)
                                                              aaaaGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTc                                                                            >  4:1984/1‑149 (MQ=255)
                                                                aaGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTa                                                                                   >  1:225/1‑140 (MQ=255)
                                                                aaGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTa                                                                                   <  2:225/140‑1 (MQ=255)
                                                                       ttATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGc                                                                   >  1:155392/1‑149 (MQ=255)
                                                                        tATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTg                                                                    <  4:7666/147‑1 (MQ=255)
                                                                                aGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGc                                                          <  2:155392/149‑1 (MQ=255)
                                                                                           aTAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTc                                                                            >  6:110232/1‑120 (MQ=255)
                                                                                           aTAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTc                                                                            <  5:110232/120‑1 (MQ=255)
                                                                                                aGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTc                                           >  1:78153/1‑148 (MQ=255)
                                                                                                     aTGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTc                                     <  2:78153/149‑1 (MQ=255)
                                                                                                            gATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACaga                               <  5:111842/148‑1 (MQ=255)
                                                                                                             aTGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACagag                              <  3:1984/148‑1 (MQ=255)
                                                                                                                         ttgaaaGAGGGAGCTAGTCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAAc                                                                                       <  3:138429/76‑1 (MQ=37)
                                                                                                                         ttgaaaGAGGGAGCTAGTCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAAc                                                                                       >  4:138429/4‑79 (MQ=37)
                                                                                                                              agagGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCAtt             >  3:98998/1‑148 (MQ=255)
                                                                                                                                          gCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGc                            <  2:126806/121‑1 (MQ=255)
                                                                                                                                          gCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGc                            >  1:126806/1‑121 (MQ=255)
                                                                                                                                          gCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGc  <  3:114045/147‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                  
AAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGC  >  NC_000913/1799236‑1799520

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