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.

AAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGC  >  NC_000913/1799236‑1799520
                                                                                                                                           |                                                                                                                                                  
aaGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAG‑CTAGTCTCCCTCttt                                                                                                                                          <  1:184177/149‑1 (MQ=255)
  gTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAG‑CTAGTCTCCCTCttt                                                                                                                                          <  6:9430/147‑1 (MQ=255)
  gTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAG‑CTAGTCTCCCTCttt                                                                                                                                          >  5:9430/1‑147 (MQ=255)
    aaaaaaaGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAAc                                                                                                                                      >  2:37822/1‑149 (MQ=255)
                ggAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGc                                                                                                                            >  4:114045/1‑147 (MQ=255)
                   aaaTAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCagag                                                                                                                       <  5:109551/149‑1 (MQ=255)
                        gcagcaTTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAgctgct                                                                                                                  >  7:150805/1‑149 (MQ=255)
                            cattcattGCTTGCGTTACCTTTTGGTACTCTACAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAG‑CTAGTCTCCCTCtt                                                                                                                                           >  8:45876/1‑120 (MQ=255)
                            cattcattGCTTGCGTTACCTTTTGGTACTCTACAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAG‑CTAGTCTCCCTCtt                                                                                                                                           <  7:45876/120‑1 (MQ=255)
                               tcattGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCtt                                                                                                           >  8:135548/1‑149 (MQ=255)
                                  ttgcttgcGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCttttt                                                                                                        <  1:37822/149‑1 (MQ=255)
                                    gcttgcGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTc                                                                                                       <  1:133773/148‑1 (MQ=255)
                                    gcttgcGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTc                                                                                                       <  7:135548/148‑1 (MQ=255)
                                     cttgcGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCaa                                                                                                     <  8:150805/149‑1 (MQ=255)
                                          gTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCaa                                                                                                     <  8:49440/144‑1 (MQ=255)
                                          gTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCaa                                                                                                     >  7:49440/1‑144 (MQ=255)
                                                ttttGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGt                                                                                           >  3:7666/1‑148 (MQ=255)
                                                        ctctTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTAc                                                                                  >  5:176197/1‑149 (MQ=255)
                                                         tctTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACt                                                                                 >  6:111842/1‑149 (MQ=255)
                                                              aaaaGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTc                                                                            >  4:1984/1‑149 (MQ=255)
                                                                aaGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTa                                                                                   >  1:225/1‑140 (MQ=255)
                                                                aaGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTa                                                                                   <  2:225/140‑1 (MQ=255)
                                                                       ttATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGc                                                                   >  1:155392/1‑149 (MQ=255)
                                                                        tATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTg                                                                    <  4:7666/147‑1 (MQ=255)
                                                                                aGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGc                                                          <  2:155392/149‑1 (MQ=255)
                                                                                           aTAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTc                                                                            >  6:110232/1‑120 (MQ=255)
                                                                                           aTAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTc                                                                            <  5:110232/120‑1 (MQ=255)
                                                                                                aGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTc                                           >  1:78153/1‑148 (MQ=255)
                                                                                                     aTGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTc                                     <  2:78153/149‑1 (MQ=255)
                                                                                                            gATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACaga                               <  5:111842/148‑1 (MQ=255)
                                                                                                             aTGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACagag                              <  3:1984/148‑1 (MQ=255)
                                                                                                                         ttgaaaGAGGGAG‑CTAGTCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAAc                                                                                       <  3:138429/76‑1 (MQ=37)
                                                                                                                         ttgaaaGAGGGAG‑CTAGTCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAAc                                                                                       >  4:138429/4‑79 (MQ=37)
                                                                                                                              agagGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCAtt             >  3:98998/1‑148 (MQ=255)
                                                                                                                                          g‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGc                            <  2:126806/121‑1 (MQ=255)
                                                                                                                                          g‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGc                            >  1:126806/1‑121 (MQ=255)
                                                                                                                                          g‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGc  <  3:114045/147‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                  
AAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGC  >  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: ‑