Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	894,697	T→C	30.9%	C242R (TGC→CGC)	rlmC →	23S rRNA m(5)U747 methyltransferase, SAM‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	894,697	0	T	C	30.9%	8.4 / 7.1	13	C242R (TGC→CGC)	rlmC	23S rRNA m(5)U747 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (6/3);  new base C (2/2);  total (8/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.40e-01

CAGCTGAAAGTTATTACCGTCAATATTCAGCCGGTACATATGGCGATTATGGAAGGGGAGACGGAGATCTACCTGACCGAACAACAGGCACTGGCGGAGCGTTTTAATGACGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACGCCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCTTTGGTTTACACTGCGCGACGCCTGACATGCAGTTAACCGGGATCGAAATTGCATCAGAGGCCATTGCCTGTGCAAAGCAGTCAGCCGCTGAACTGGGCTTAACGCGTTTGCAATTTCAGGCGCTGGACTCCACTCAGTTTGCCACCGCTCAGGGGGATGTGCCGGAGCTGGTGCTGGTTAACCCGCCGCGC  >  CP009273/894463‑894903                                                                                                                                                                                                                                           |                                                                                                                                                                                                               cAGCTGAACGTTATTACCGTCAATATTCAGCCGGTACATATGGCGATTATGGAAGGGGAGACGGAGATCTACCTGACCGAACAACAGGCACTGGCGGAGCGTTTTAATGACGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGAGGTCGCCAGCCAGTTGTACGACACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTATTCTGCGGTGGGGGGAGCATTGGCTTACACTGCGCGACGCCTGa                                                                                                                                                                        >  1:17400/1‑275 (MQ=255)                                                                             ccGAACAACAGGCACTGGCGGAGCGTTTTAATGACGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACGCCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCCGCGGTGTGGGGGGCTTTGGTTTACACTgcg                                                                                                                                                                                  >  2:307930/1‑189 (MQ=255)                                                                                 acaacaGGCACTGGCGGAGCGTTTTAATGACGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACGCCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCCGCGGTGTGGGGGGCTTTGGTTTACACTGCGCGAc                                                                                                                                                                              <  1:307941/189‑1 (MQ=255)                                                                                                        ttAATGACGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACGCCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGAATCTGTTCTGCGGTGTGGGGGGCTTTGGTTTACACTGCGCGACGCCTGACATGCAGTTAACCGGAATCGAACTTGCATCAGAGGCCATTGCCTGTTCAAAGCAGTCAGCCCCTGAACTGGGCTTAACGCGTTTGCAATTTCAGGCGCTGGa                                                                  >  1:70540/1‑274 (MQ=255)                                                                                                                gTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACGCCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTgggg                                                                                                                                                                                                     >  1:54041/1‑135 (MQ=255)                                                                                                                    cGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACGCCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCt                                                                                                                                                                                                 <  2:54039/135‑1 (MQ=255)                                                                                                                                   cgCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACGCCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCTTTGGTTTACACTGCGCGACGCCTGACATGCAGTTAACCGGGATCGAAATTGCATCAGAGGCCATTGCCTGTGCAAAGCAGTCAGCCGCTGAACTGGGCTTAACGCGTTTGCAATTTCAGGCGCTGGACCCCACTCAGTTTGCCACCGCACAg                                         >  1:133808/1‑272 (MQ=255)                                                                                                                                                        aTCCGGCGGTCGCCAGCCAGTTGTACGCCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCTTTGGTTTACACTGCGCGACGCCTGACATGCAGTTAACCGGGATCGAAATTGCATCAGAGGCCATTGCCTGTGCAAAGCAGTCAGCCGCTGAACTGGGCTTAACGCGTTTGCAATTTCAGGCGCTGGAc                                                                 >  2:64505/1‑227 (MQ=255)                                                                                                                                                            ggcggTCGCCAGCCAGTTGTACGCCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCTTTGGTTTACACTGCGCGACGCCTGACATGCAGTTAACCGGGATCGAAATTGCATCAGAGGCCATTGCCTGTGCAAAGCAGTCAGCCGCTGAACTGGGCTTAACGCGTTTGCAATTTCAGGCGCTGGACTCCa                                                             <  1:64507/227‑1 (MQ=255)                                                                                                                                                                                      ccGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCCGCGGTGTGGGGGGCTTTGGTTTACACTg                                                                                                                                                                                    >  2:231542/1‑82 (MQ=255)                                                                                                                                                                                          gcgcgACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCCGCTATGTTAGAGGCTTTGGTTTACACTGcgcg                                                                                                                                                                                <  1:231549/82‑1 (MQ=255)                                                                                                                                                                                                             cTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCTTTGGTTTACACTGCGCGACGCCTGACATGCAGTTAACCGGGATCGAAATTGCATCAGAGGCCATTGCCTGTGCAAAGCAGTCAGCCGCTGAACTGGGCTTAACGCGTTTGCAATTTCAGGCGCTGGACTCCACTCAGTTTGCCACCGCTCAGGGGGATGTGCCGGAGCTGGTGCTGGTTAACccgcc      >  2:89537/1‑233 (MQ=255)                                                                                                                                                                                                                 cGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCTTTGGTTTACACTGCGCGACGCCTGACATGCAGTTAACCGGGATCGAAATTGCATCAGAGGCCATTGCCTGTGCAAAGCAGTCAGCCGCTGAACTGGGCTTAACGCGTTTGCAATTTCAGGCGCTGGACTCCACTCAGTTTGCCACCGCTCAGGGGGATGTGCCGGAGCTGGTGCTGGTTAACCcgccgcgc  <  1:89539/233‑1 (MQ=255)                                                                                                                                                                                                                                           tGCGGTGTGGGGGGCTTTGGTTTACACTGCGCGACGCCTGACATGCAGTTAACCGGGATCGAAATTGCATCAGAGGCCATTGCCTGTGCAAAGCAGTCAGCCGCTGAACTGGGCTTAACGCGTTTGCAATTTCAGGCGCTGGACTCCACTCAGttt                                                     >  1:62273/1‑156 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                               CAGCTGAAAGTTATTACCGTCAATATTCAGCCGGTACATATGGCGATTATGGAAGGGGAGACGGAGATCTACCTGACCGAACAACAGGCACTGGCGGAGCGTTTTAATGACGTACCGCTGTGGATCCGTCCGCAAAGTTTCTTCCAGACTAATCCGGCGGTCGCCAGCCAGTTGTACGCCACCGCGCGCGACTGGGTACGACAGCTGCCGGTTAAACATATGTGGGATCTGTTCTGCGGTGTGGGGGGCTTTGGTTTACACTGCGCGACGCCTGACATGCAGTTAACCGGGATCGAAATTGCATCAGAGGCCATTGCCTGTGCAAAGCAGTCAGCCGCTGAACTGGGCTTAACGCGTTTGCAATTTCAGGCGCTGGACTCCACTCAGTTTGCCACCGCTCAGGGGGATGTGCCGGAGCTGGTGCTGGTTAACCCGCCGCGC  >  CP009273/894463‑894903

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑