Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,171,5570CA46.3% 26.7 / 66.6 41intergenic (+13/‑80)rrlB/rrfB23S ribosomal RNA of rrnB operon/5S ribosomal RNA of rrnB operon
Reads supporting (aligned to +/- strand):  ref base C (11/11);  new base A (9/10);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.37e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGACGTTGATAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAG  >  NC_000913/4171467‑4171687
                                                                                          |                                                                                                                                        
cgccggtGATAGGCCGGGTGGGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCa                                                                                    <  1:172265/133‑1 (MQ=0)
                                            tAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAg                                        <  2:39561/139‑1 (MQ=14)
                                                   tACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGcgcg                                 <  2:1398705/139‑1 (MQ=18)
                                                        aTGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCggtggt                            >  2:961188/1‑139 (MQ=33)
                                                         tGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTc                           <  1:33957/139‑1 (MQ=33)
                                                         tGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTc                           <  2:1675113/139‑1 (MQ=33)
                                                         tGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTc                           <  2:159303/139‑1 (MQ=32)
                                                         tGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGCTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTc                           <  2:1129961/139‑1 (MQ=18)
                                                            aCCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCa                        >  1:402826/1‑139 (MQ=17)
                                                             ccGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCggtgg                             <  1:1429412/133‑1 (MQ=32)
                                                             ccGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCggtgg                             >  2:1429412/1‑133 (MQ=32)
                                                             ccGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCAc                       >  2:1703017/1‑139 (MQ=33)
                                                             ccGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCggtgg                             >  1:1731017/1‑139 (MQ=11)
                                                             ccGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCAc                       >  2:1567737/1‑139 (MQ=17)
                                                              cGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCAcc                      <  1:652865/139‑1 (MQ=33)
                                                                tGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGGGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTg                    >  2:1190670/1‑139 (MQ=18)
                                                                tGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGcgcg                                 >  1:76134/1‑126 (MQ=33)
                                                                tGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGcgcg                                 >  1:344860/1‑126 (MQ=33)
                                                                tGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGcgcg                                 <  2:344860/126‑1 (MQ=33)
                                                                tGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGcgcg                                 <  2:76134/126‑1 (MQ=33)
                                                                tGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTg                    >  2:1157088/1‑139 (MQ=33)
                                                                tGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTg                    >  2:1562258/1‑139 (MQ=18)
                                                                     cTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGAcccc               <  1:679503/139‑1 (MQ=17)
                                                                     cTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTCTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGcgcg                                 <  1:810094/121‑1 (MQ=16)
                                                                     cTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTCTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGcgcg                                 >  2:810094/1‑121 (MQ=16)
                                                                             tACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGaa       <  1:1653344/139‑1 (MQ=17)
                                                                              aCAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATg            <  2:1612405/133‑1 (MQ=33)
                                                                              aCAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATg            >  1:671204/1‑133 (MQ=33)
                                                                              aCAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATg            >  1:1612405/1‑133 (MQ=33)
                                                                              aCAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATg            <  2:671204/133‑1 (MQ=33)
                                                                              aCAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAAc      >  2:448675/1‑139 (MQ=32)
                                                                              aCAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCAcc                      <  2:1026219/123‑1 (MQ=18)
                                                                              aCAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCAcc                      >  1:1026219/1‑123 (MQ=18)
                                                                              aCAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAAc      >  1:1399709/1‑139 (MQ=17)
                                                                                  cGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGcc          >  1:1614860/1‑131 (MQ=32)
                                                                                  cGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGcc          <  2:1614860/131‑1 (MQ=32)
                                                                                  cGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAg  >  2:225603/1‑139 (MQ=18)
                                                                                   gCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGcgcg                                 <  2:309851/107‑1 (MQ=17)
                                                                                   gCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGcgcg                                 >  1:309851/1‑107 (MQ=18)
                                                                                       aaGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAg                                        <  1:892823/96‑1 (MQ=17)
                                                                                       aaGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAgcgc                                  <  2:255936/102‑1 (MQ=18)
                                                                                          |                                                                                                                                        
CGACGTTGATAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCT‑‑‑‑‑‑GATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAG  >  NC_000913/4171467‑4171687

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑