Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,029,739	0	A	.	32.3%	68.1 / 33.2	32	coding (231/1191 nt)	rlmI	23S rRNA m(5)C1962 methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (10/11);  new base . (7/3);  total (18/14)
Fisher's exact test for biased strand distribution p-value = 2.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.27e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAA‑CGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACTTTTCTGTC  >  NC_000913/1029493‑1029985                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        gACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAAGCa                                                                                                                                                                                                                                                      >  2:280692/1‑250 (MQ=255)   cGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAAGCAATg                                                                                                                                                                                                                                                   <  1:280692/251‑1 (MQ=255)                     tCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGCTCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCGTTTTTGTGCTTGTTGCAAACGGCGGGAAAACAAAGAAATGTCGATAGATTCAGACggg                                                                                                                                                                                                                                 >  2:77841/1‑251 (MQ=255)                              ccGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCTATATTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCGTGTTGCAAACGGCGGGAAAAAAA‑GCAATGTCGATAGACTCAGACGTGTCAAA‑CGTcc                                                                                                                                                                                                                      <  2:87149/251‑1 (MQ=255)                                             gCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAAGCAATGTCGCTAGAATCAGACGGGGCAAAACG‑CCAGACGCGCGCCCgg                                                                                                                                                                                                        >  2:24719/1‑251 (MQ=255)                                                  gCACTAATTAATGCCGTGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAAGAAATGTCGAATGACTCAGACGGGTCaaaaggtcaaaatcacgcccggatttg                                                                                                                                                                                                    >  2:270856/1‑227 (MQ=255)                                                          taatGCCGCGCGCTGATGTTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACTGAAACTTCCGATAGTAAGACCGGGCAGGCCATCAGACGCCCCGGCGATTAAACGATAGCTGTCGAGGCCTGCTTTTGGCGCCAGCCAGTCACGCGATTCTTTTGCTTGTGGCAAACGGCGGGAAAAAAAAGCAATGTCGATAGACTCAGACGGGGCAAA‑CGTCCAGTCGCGCGCCCGGATTTGCGAAGCt                                                                                                                                                                                            <  2:42253/250‑1 (MQ=255)                                                          taatGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAA‑CGTCCAGACGCGCGCCCTGATTTGCGAAGCt                                                                                                                                                                                            <  1:77753/250‑1 (MQ=255)                                                           aatGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAA‑GAAATGTCGATAGACCCAGACGGGTCAAA‑CGTCCCGACGCGCGCCCGGGATTGCGAAGCTggg                                                                                                                                                                                         >  2:51929/1‑250 (MQ=255)                                                                cgcgcgCTGATATTCTGCCCCAGCACTCAGCCGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAA‑CGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAAt                                                                                                                                                                                     >  1:220402/1‑251 (MQ=255)                                                                         atatTCTGCCCCAGCACTCAGAAGTTGCAGCCCCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAA‑GCAATGTCGATAGACTCAGACGGGTCAAA‑CGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCcgcg                                                                                                                                                                           <  1:218877/251‑1 (MQ=255)                                                                            ttCTGCCCAAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAAGCAATGTCGATAGAATCAGACGGGTCAAA‑CGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTg                                                                                                                                                                         >  1:126689/1‑251 (MQ=255)                                                                                   ccAGCACTCAGCAGTGGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAA‑CGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTTCTAACCa                                                                                                                                                                  <  2:210304/251‑1 (MQ=255)                                                                                                  tGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGGAAAAAA‑GCAATGCCGATCGACTCAACCGGGGCAAA‑CCCCCCGACGCGCGCCCAGATTTGCGAAGCTGGCGCAAAAGCGCCGCGGTCTAACCCTTTTCCCttgtgctc                                                                                                                                                   >  2:209158/1‑243 (MQ=255)                                                                                                  tGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAA‑CGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTgatgat                                                                                                                                                    <  1:64677/250‑1 (MQ=255)                                                                                                     agcaCCAGAAAATTCCCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGGCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAA‑CGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATcaac                                                                                                                                                <  2:134075/251‑1 (MQ=255)                                                                                                             aaaaTTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAA‑CGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCg                                                                                                                                         >  1:266062/1‑250 (MQ=255)                                                                                                               aaTTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAA‑GCAATGTCGATAGACTCAGACGGGTCAAA‑CGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATgg                                                                                                                                     >  1:349442/1‑251 (MQ=255)                                                                                                                 ttACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAA‑CGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGt                                                                                                                                    >  1:127865/1‑251 (MQ=255)                                                                                                                              cgatAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAA‑GCAATGTCGATAGACTCAGACGGGTCAAA‑CGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAggctgg                                                                                                                      >  1:116625/1‑251 (MQ=255)                                                                                                                                        aCCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAA‑CGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCtt                                                                                                             >  1:64408/1‑251 (MQ=255)                                                                                                                                             gCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAA‑CGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATg                                                                                                        <  2:64408/251‑1 (MQ=255)                                                                                                                                               aGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGCGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAA‑CGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCg                                                                                                      <  2:88502/251‑1 (MQ=255)                                                                                                                                                ggCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAA‑CGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCgg                                                                                                     <  1:274576/251‑1 (MQ=255)                                                                                                                                                                           gATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAACGCAATGTCGATAGACTCAGACGGGTCAAA‑CGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTCACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCCACGGCCCCGGAAAAAACCCCCggg                                                                          >  2:191241/1‑250 (MQ=255)                                                                                                                                                                                  gTCGAGGCCATCTTTGTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAA‑CGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCGGAAAAGACCCACGGATGGcga                                                                    <  2:117431/250‑1 (MQ=255)                                                                                                                                                                                        gCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGAGGGAAAAAAA‑GCAATGTCGATAGACTAAGACGGGTCAAA‑CGTCCCGACCCCAGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCGGAAAAGACCCCCGGATGGCGACGAAGTa                                                             >  1:368375/1‑250 (MQ=255)                                                                                                                                                                                               ttttGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAA‑GCAATGTCGATAGACTCAGACGGGTCAAA‑CGTCCAGACGCGCGCCCGGATTTTCGAAGCTGGCGAATAAACGCCGCCGGCTAACCATTTTCCCTGATGATCAACAATATCGATGGGTTCAACCAGGCCGGCTTTACCTTCCCTGCGGGCAAAGGCCCCGGAAAAGACCCCCGGATTGCGAACAAGTAAttgtttt                                                     >  2:328363/1‑245 (MQ=255)                                                                                                                                                                                                           cagTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAA‑CGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTGACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCttt                                          <  2:259956/251‑1 (MQ=255)                                                                                                                                                                                                           cagTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAA‑CGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATTCGGGGAACGGCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCTCGCCCttt                                          >  2:23428/1‑251 (MQ=255)                                                                                                                                                                                                                         ttGTGCTTGTTGCAAACGGCGGGAAAAAAA‑GCAATGTCGATAGACTCAGACGGGTCAAA‑CGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAAc                            <  2:271076/250‑1 (MQ=255)                                                                                                                                                                                                                                    gCAAACGGCGGGAAAAAAA‑GCCATGTCGATAGACTCAGCCGGGTCAAA‑CGGCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATCAGCGCCGCGTGCTAACCCTTTTCCCTGATGATCAACAATATCGATGGTTTCCCCGCGGCTGGCTTTACCTTCCATGCGGGCAACGGACCCGGAAAAGACCCCCGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATa                 >  2:358309/1‑250 (MQ=255)                                                                                                                                                                                                                                                   aaaaaGCAATGTCGATAGACTCAGACGGGTCAAA‑AGTCCAGACGCGCGCACGGCTTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACTTTTCTGTc  >  1:151808/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        GACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAA‑CGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACTTTTCTGTC  >  NC_000913/1029493‑1029985

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