Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,144,974	C→T	A12V (GCT→GTT)	rluC →	23S rRNA pseudouridine(955/2504/2580) synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,144,974	0	C	T	90.9%	23.4 / NA	11	A12V (GCT→GTT)	rluC	23S rRNA pseudouridine(955/2504/2580) synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (7/3);  minor base G (0/1);  total (7/4)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

TTGTTATAGCAAGATGACTTTTACCATTTATCACCCGCTTACTCACAGTTTTTTCACTTCTTGCTGGTGATTGGTTTAATAACCATCAATAAAGTGATCACGCAGTGAAGTAATAAGGTGCAAAAGTAAATATAAGCATAGAAAAATGAGTGGCGCGAATCCTGTTGGCTAATTAGAATCGCCAACCATGAAAACAGAGACTCCATCCGTAAAAATTGTTGCTATCACCGCTGACGAAGCGGGGCAACGTATCGATAACTTTTTGCGTACCCAATTGAAAGGCGTACCAAAAAGTATGATTTACCGTATTTTGCGTAAAGGCGAAGTGCGGGTGAACAAAAAACGTATTAAGCCTGAATATAAACTCGAAGCGGGTGATGAGGTGCGTATTCCACCGGTTCGCGTTGCTGAGCGGGAAGAAGAGGCG  >  NC_000913/1144753‑1145179                                                                                                                                                                                                                              |                                                                                                                                                                                                              ttGTTATAGCAAGATTACTTTTACCATTTATCACCCGACTACTCACAGTTTTTTCACTTCTTGCTGGTGATTGTGTTAATAACCAACAATAAAGTGATCACGCAGTGAAGTAATAAGGTGCAAAAGTAAATATAAGCATAGAAAAATGAGTGGCGCGAATCCTGTTGGCTAATTAGAATCGCCACCCATGAAAACAGAGACTCCATCCGTAAAAATTGTTGTTATCCCCGCTGACTAAGCGGGGCAACGt                                                                                                                                                                                   <  2:6649/250‑1 (MQ=255)                    tttACCATGTATCACCCTCTTACTCACAGTTTTTTCACTTCGTGCTTGTGATTGGTTTAATAACCAGCAAAAAAGTGATCACGCAGTGAAGTAATAAGGTGCAAAAGTAAATATAAGCATGGAAAAATGAGTGGCGCGAATCCTGTTGGCTAATTAGAATCGCCAACCATGAAAACAGAGACGCCATCCGTAAAAATTGTTGTTAGCACCGCTGACGAAGCGGGTCAACGTCTCGATAACTTTTTGCGTa                                                                                                                                                                <  2:13907/250‑1 (MQ=255)                                       ttACTCACAGTTTTTTCACTTCTTGCTGGTGCTTGGTTTAATAACCATCAATAAAGTGATCACGCAGTGAAGTAATAAGGTGCAAAAGTAAATATAAGCATAGAAAAATGAGTGGCGCGAATCCTGTTGGCTAATTAGAATCGCCAACCATGAAAACAGAGACTCCATCCGTAAAAATTGTTGTTATCACCGCTGACGAAGCGGGGCAACGTATCGATAACTTTTTGCGTACCCAATTGAAAGGCGTAcc                                                                                                                                             >  1:66979/1‑250 (MQ=255)                                                      tCACTTTTTGATGGTGATTGGTTTAATAACCATCAATAAAGGGATCCCGCAGTGAGGTAATAAGGTGAAAAAGTAAATATAAGCATAGAAAAATGAGTGGCGCGAATCCTGTTGGCTAATTAGAATCGCCAACCATGAAAACAGAGACTCCATCCGTAAAAATTGTTGTTATCACCGCTGACGAAGCGGGGCCACGTATCGATAACTTTTTGCGTACCCAATTGAAAGGCGTACCAAAAAGTATGATTTAc                                                                                                                             >  1:204891/1‑251 (MQ=255)                                                              tGCTGGTGATTGGTTTAATAACCATCAATAAAATGATCACGCAGTAAAGTAATAAGTTGCAAAAGTAAATAGAAACATAGAAAAATGAGTGGCGCTAATCCTGTTGGCTAAGTAGAATCGCCAACCATGAAAACAGAGACGCCATCCGTAAAAATTGTTGGTAGCACCGCAGACGAAGCGGGGCAACGTATCGATAACTTTTTGCGTTCCCAATTTAAAGGCGTACCAAAAAGTATGATTTACCGTAtt                                                                                                                       <  2:204891/249‑1 (MQ=255)                                                                                taaCCATCAATAAAGGGATCACGCAGTGAAGTTATAAGGTGCAAAAGTAAATATAAGCATAGAAAAATGAGTGGCGCGAATCCTGTTGGCTAATTAGAATCGCCAAACATGAAAACAGAGACTCCATACGTAAAAATTGTTGTTATCACCGCTGACGAAGCGGGGCAACGTATCGATAACTTTTTGCGTACCCAATTGAAAGGCGGACCAAAAAGTATGATTTACCGTATTTTGCGTAAAGGCGAAGTGCg                                                                                                   >  2:203000/1‑251 (MQ=255)                                                                                                                        gCAAAAGTAAATATAAGCATAGAAAAATGAGTGGCGCGAATCCTGTTGGCTAATTAGAATCGCCAACCATGAAAACAGAGACTCCATCCGTAAAAATTGTTGTTATCACCGCTGACGAAGCGGGGCAACGTATCGATAACTTTTTGCGTACCCAATTGAAAGGCGTACCAAAAAGTATGATTTACCGTATTTTGCGTAAAGGCGAAGTGCGGGTGAACAAAAAACGTATTAAGCCTGAATATAAACTCGa                                                            >  2:152207/1‑250 (MQ=255)                                                                                                                                 aaTATAAGCATAGAAAAATGAGTGGCGCGAATCCTGTTGGCTAATTAGAATCGCCAACCATGAAAACAGAGACTCCATCCGTAAAAATTGTTGTTATCACCGCTGACGAAGCGGGGCAACGTATCGATAACTTTTTGCGTACCCAATTGAAAGGCGTACCAAAAAGTATGATTTACCGTATTTTGCGTAAAGGCGAAGTGCGGGTGAACAAAAAACGTATTAAGCCTGAATATAAACTCGAAGCGGGtgat                                                  >  1:201939/1‑251 (MQ=255)                                                                                                                                                 aaTGAGAGGCGCGAATCCTGTTGGCTAATTAGAATCGCCAACCATGAAAACAGAGACTCCATCCGTAAAAATTGTTGTTATCACCGCTGACGAAGCGGGGCAACGTATCGATAACTTTTTGCGTACCCAATTGAAAGGCGTACCAAAAAGTATGATTTACCGTATTTTGCGTAAAGGCGAAGTGCGGGTGAACAAAAAACGTATTAAGCCTGAATATAAACTCGAAGCGGGAGATGAGGAGGGTATTCCAc                                  >  2:210147/1‑251 (MQ=255)                                                                                                                                                                     ttGGCTAATTAGCATCGCCAACCATGAAAACAGAGACTCCATCCGTAAAAATTGTTGTTATCACCGCTGACGAAGCGGGGCAACGTATCGATAACTTTTTGCGTACCCAATTGAAAGGCGTACCAAAAAGTATGATTTACCGTATTTTGCGTAAAGGCGAAGTGCGGGTGAACAAAAAACGTATTAAGCCTGAATATAAACTCGAAGCGGGTGATGAGGTGCGTATTCCACCGGTTCGCGTTGCTGAGCgg              >  1:223225/1‑251 (MQ=255)                                                                                                                                                                                 aaTCGCCAACCATGAAAACAGAGACTCCATCCGTAAAAATTGTTGTTATCACCGCTGACGAAGCGGGGCAACGTATCGATAACTTTTTGCGTACCCAATTGAAAGGCGTACCAAAAAGTATGATTTACCGTATTTTGCGTAAAGGCGAAGTGCGGGTGAACAAAAAACGTATTAAGCCTGAATATAAACTCGAAGCGGGTGATGAGGTGCGTATTCCACCGGGTCGCGTTGCTGAGCGGGAAGAAGAGGCg  <  2:130696/251‑1 (MQ=255)                                                                                                                                                                                                                              |                                                                                                                                                                                                              TTGTTATAGCAAGATGACTTTTACCATTTATCACCCGCTTACTCACAGTTTTTTCACTTCTTGCTGGTGATTGGTTTAATAACCATCAATAAAGTGATCACGCAGTGAAGTAATAAGGTGCAAAAGTAAATATAAGCATAGAAAAATGAGTGGCGCGAATCCTGTTGGCTAATTAGAATCGCCAACCATGAAAACAGAGACTCCATCCGTAAAAATTGTTGCTATCACCGCTGACGAAGCGGGGCAACGTATCGATAACTTTTTGCGTACCCAATTGAAAGGCGTACCAAAAAGTATGATTTACCGTATTTTGCGTAAAGGCGAAGTGCGGGTGAACAAAAAACGTATTAAGCCTGAATATAAACTCGAAGCGGGTGATGAGGTGCGTATTCCACCGGTTCGCGTTGCTGAGCGGGAAGAAGAGGCG  >  NC_000913/1144753‑1145179

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