Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	842,503	G→A	A58T (GCG→ACG)	rlmF →	23S rRNA m(6)A1618 methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	842,503	0	G	A	96.2%	69.3 / NA	26	A58T (GCG→ACG)	rlmF	23S rRNA m(6)A1618 methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (15/10);  total (15/11)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTACACCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGTGAAGGCGCTCAATAAGGCGTTGCTGGCCCATTTTTACGCCGTAGCGAACTGGGATATCCCCGACGGTTTTCTCTGCCCACCAGTACCGGGCCGGGCGGATTATATTCATCACCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGC  >  NC_000913/842263‑842747                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     aaaCGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGCGAATGCTATGTCCGCCCAGACACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTACACCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGTGAAGGCGCTCAATAAGACGTTGCTGGc                                                                                                                                                                                                                                            >  1:250171/1‑251 (MQ=255)   aCGCAGCCAGATATTTACCCTGTCGCGTTATGCGGTTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTACACCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGTGAAGGCGCTCAATAAGACGTTGCTGGccc                                                                                                                                                                                                                                          <  1:184340/251‑1 (MQ=255)     gCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTACACCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGTGAAGGCGCTCAATAAGACGTTGCTGGCCCAt                                                                                                                                                                                                                                        <  1:345866/251‑1 (MQ=255)           aGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTACACCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGTGAAGGCGCTCAATAAGACGTTGCTGGCCCATTTTTAc                                                                                                                                                                                                                                  >  1:315145/1‑251 (MQ=255)               atTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTACACCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGTGAAGGCGCTCAATAAGACGTTGCTGGCCCATTTTTACGccc                                                                                                                                                                                                                              >  2:5606/1‑250 (MQ=255)                        gTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTACACCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGTGAAGGCGCTCAATAAGACGTTGCTGGCCCATTTTTACGCCGTAGCGAACt                                                                                                                                                                                                                     >  1:178161/1‑251 (MQ=255)                                                ccgcTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTACACCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGTGAAGGCGCTCAATAAGACGTTGCTGGCCCATTTTTACGCCGTAGCGAACTGGGATATCCCCGACGGTTTTctct                                                                                                                                                                                             >  1:73708/1‑251 (MQ=255)                                                 cgcTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTACACCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGTGAAGGCGCTCAATAAGACGTTGCTGGCCCATTTTTACGCCGTAGCGAACTGGGATATCCCCGACGGTTTTCTCTg                                                                                                                                                                                            >  1:165960/1‑251 (MQ=255)                                                 cgcTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTACACCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGTGAAGGCGCTCAATAAGACGTTGCTGGCCCATTTTTACGCCGTAGCGAACTGGGATATCCCCGACGGTTTTCTCTg                                                                                                                                                                                            <  1:157668/251‑1 (MQ=255)                                                       ttAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTACACCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGGGAAGGCGCTCAATAAAACGTTGCTGTCCCCTTTTTACCCCGTCACGAACTGGGATATCCCCGTCGTGTTTCTCTGCccacc                                                                                                                                                                                      >  2:271345/1‑251 (MQ=255)                                                                              ccAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTACACCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGTGAAGGCGCTCAATAAGACGTTGCTGGCCCATTTTTACGCCGTAGCGAACTGGGATATCCCCGACGGTTTTCTCTGCCCACCAGTACCGGGCCGGGCGGATtatt                                                                                                                                                               >  2:136158/1‑250 (MQ=255)                                                                                     ccGGGGTTGCAACCGCGCAACCGTCATCACAGCAGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTACACCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGTGGAGGCGCTCAATAAGACGTTGCTGGCCCATTTTTACGCCGTAGCGAACTGGGATATCCCCGACGGTTTTCTCTGCCCACCAGTACCGGGCCGGGCGGATTCTATtcctc                                                                                                                                                         >  2:254048/1‑247 (MQ=255)                                                                                                           gTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTACACCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGTGAAGGCGCTCAATAAGACGTTGCTGGCCCATTTTTACGCCGTAGCGAACTGGGATACCCCCGACGGTTTTCGCTGCCCACCAGTACCGGGCCGGGCGGAGTATATTCATCACCTTGCCGATTTACTGGCAGaa                                                                                                                                  <  2:188996/251‑1 (MQ=255)                                                                                                                                  gCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTACACCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGTGAAGGCGCTCAATAAGACGTTGCTGGCCCATTTTTACGCCGTAGCGAACTGGGATATCCCCGACGGTTTTCTCTGCCCACCAGTACCGGGCCGGGCGGATTATATTCATCACCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGaa                                                                                                           <  2:127014/251‑1 (MQ=255)                                                                                                                                                   cAATCCTGAAGTCAGGCAATTCCTCACGCTTCCACCCGCCGGGGATCAGAGCGTAGACTTTGCCAATCCGCTTGCGGTGAAGTCGCGCAAGAATGCGTTGCTGGCCCATTTTTACGCCGTAGCGAACTGGGATATCCCCGACGGTTTTCTCTTCCCACCAGTACCGGGCCGGGCGGTTTATATTCGTCTCGTGGCCGAGTTGCTGGCAGAAGCAGGCTGAATAATTCCAGCGaaagcgagctgtgtgtatt                                                                                          <  2:173376/251‑18 (MQ=255)                                                                                                                                                                  gCAATTCCTCACGCTTACACCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGAGAAGGCGCTCAATAAGACAATGATGGCCCATTTTTACGCCGTAGCGAACTGGGATATCCCCGACGGTTTTCTCTGCCCACCAGTACCGGGCCGGGCGGATTATATTCATCACCTTGCCGATTTACTGGAGGAAGCAAGCGGAACCATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCg                                                                            >  1:244328/1‑250 (MQ=255)                                                                                                                                                                           cACGCTTACACCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGTGAAGGCGCTCAATAAGACGTTGCTGGCCCATTTTTACGCCGTAGCGAACTGGGATATCCCCGACGGTTTTCTCTGCCCACCAGTACCGGGCCGGGCGGATTATATTCATCACCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTAtt                                                                   <  2:65571/250‑1 (MQ=255)                                                                                                                                                                            aCGCTTACACCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGTGAAGGCGCTCAATAAGACGTTGCTGGCCCATTTTTACGCCGTAGCGAACTGGGATATCCCCGACGGTTTTCTCTGCCCACCAGTACCGGGCCGGGCGGATTATATTCATCACCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTAttt                                                                  >  2:207796/1‑250 (MQ=255)                                                                                                                                                                              gCTTACACCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGTGAAGGCGCTCAATGATACGTTGCGGGCCCAt                                                                                                                                                                                                                                        <  2:215324/82‑1 (MQ=255)                                                                                                                                                                                   cacCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGTGAAGGCGCTCAATAAGACGTTGCTGGCCCATTTTTACGCCGTAGCGAACTGGGATATCCCCGACGGTTTTCTCTGCCCACCAGTACCGGGCCGGGCGGATTATATTCATCACCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCTGCGAATGCCAGCATTCTGGATAGCGGCGTTGGTGCGAACTGTATTTATTCGCTg                                                          <  2:53189/251‑1 (MQ=255)                                                                                                                                                                                                  aaaGCGAAGACTTTGCCAATCCACAGGCGGCGGAGGCGCTCAATAAGACGGTGCTAGCACATTTTTACGCAGTAGCCAACCGGGCTATCCCCGACGGTTTTCTCTGTCCACCAGTACCGGGCCGGGCCGATTATATTCATCACCTTGCCGATTTACTTGCAGAAGCCAGCGGAACAACTCCAGCGCATGCCAGCATTCTGGCTATCGGCGTTGGGGCGAACTGTATTTATCCGCTGATTGGCGTACAtta                                            >  2:231134/1‑248 (MQ=255)                                                                                                                                                                                                       gTAGACTTTGCCAATCCGCTGGCGGTGAAGGCGCTCAATAAGACGTTGCTGGCCCATTTTTACGCCGTAGCGAACTAGGATATCCCCGACGGTTTTCTCTGCCCACCAGTACCGGGCCGGGCGGATTATATTCATCACCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATAtgg                                      >  2:194194/1‑251 (MQ=255)                                                                                                                                                                                                                gCCAATCCGCTGGCGGTGAAGGCGCTCAATAAGACGTTGCTGGCCCATTTTTACGCCGTAGCGAACTGGGATATCCCCGACGGTTTTCTCTGCCCACCAGTACCGGGCCGGGCGGATTATATTCATCACCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGGTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGttt                             <  2:140990/251‑1 (MQ=255)                                                                                                                                                                                                                             cGGTGAAGGCGCTCAATAAGACGTTGCTGGCCCATTTTTACGCCGTAGCGAACTGGGATATCCCCGACGGTTTTCTCTGCCCACCAGTACCGGGCCGGGCGGATTATATTCATCACCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGa                 >  1:36101/1‑250 (MQ=255)                                                                                                                                                                                                                                          cAATAAGACGTTGCTGGCCCATTTTTACGCCGTAGCGAACTGGGAGATCCCCGACGGTTTTCTCTGCCCACCAGTACCGGGCCGGGCGGATTATATTCATCACCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAgg   <  2:117119/251‑1 (MQ=255)                                                                                                                                                                                                                                           aaTAAGACGTTGCTGGCCCATTTTTACGCCGTAGCGAACTGGGATATCCCCGACGGTTTTCTCTGCCCACCAGTACCGGGCCGGGCGGATTATATTCATCACCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGc  >  2:303390/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     AAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTACACCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGTGAAGGCGCTCAATAAGGCGTTGCTGGCCCATTTTTACGCCGTAGCGAACTGGGATATCCCCGACGGTTTTCTCTGCCCACCAGTACCGGGCCGGGCGGATTATATTCATCACCTTGCCGATTTACTGGCAGAAGCAAGCGGAACAATTCCAGCGAATGCCAGCATTCTGGATATCGGCGTTGGTGCGAACTGTATTTATCCGCTGATTGGCGTACATGAATATGGCTGGCGTTTTACCGGTAGTGAAACCAGCAGCCAGGC  >  NC_000913/842263‑842747

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