Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	842,301	T→C	intergenic (‑245/‑31)	mcbA ← / → rlmF	DUF1471 domain‑containing protein McbA/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,301	0	T	C	100.0%	83.7 / NA	23	intergenic (‑245/‑31)	mcbA/rlmF	DUF1471 domain‑containing protein McbA/23S rRNA m(6)A1618 methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (11/12);  total (11/12)

TACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGG  >  NC_000913/842171‑842424                                                                                                                                   |                                                                                                                            tACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGCATCTTAcgc                                                                                                                    <  2:63154/140‑1 (MQ=255)   cAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGCATCTTAcgccg                                                                                                                  <  1:208552/140‑1 (MQ=255)                     tGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGCATCTTACGCCGCTGCTTAAAGgag                                                                                                     <  1:405570/135‑1 (MQ=255)                     tGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGCATCTTACGCCGCTGCTTAAAGgag                                                                                                     >  2:405570/1‑135 (MQ=255)                          tGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGCATCTTACGCCGCTGCTTaaa                                                                                                         <  2:360665/126‑1 (MQ=255)                          tGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGCATCTTACGCCGCTGCTTaaa                                                                                                         >  1:360665/1‑126 (MQ=255)                           gCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGCATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTc                                                                                          >  2:272520/1‑140 (MQ=255)                                 tGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGCATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTTCGCCCa                                                                                    >  1:294553/1‑140 (MQ=255)                                 tGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGCATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCa                                                                                    >  1:120988/1‑140 (MQ=255)                                                        cATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGCATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAAcc                                                                              <  1:25825/123‑1 (MQ=255)                                                        cATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGCATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAAcc                                                                              >  2:25825/1‑123 (MQ=255)                                                                 tAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGCATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCAt                                                                    >  1:387708/1‑124 (MQ=255)                                                                 tAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGCATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCAt                                                                    >  1:29404/1‑124 (MQ=255)                                                                 tAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGCATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCAt                                                                    <  2:387708/124‑1 (MQ=255)                                                                 tAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGCATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCAt                                                                    <  2:29404/124‑1 (MQ=255)                                                                         ttAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGCATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGc                                            >  2:352014/1‑140 (MQ=255)                                                                                                   aGCCAGATATTTACCCTGTCGCGTTATGCGGGCATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGt                  <  1:203062/140‑1 (MQ=255)                                                                                                   aGCCAGATATTTACCCTGTCGCGTTATGCGGGCATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGt                  <  2:294553/140‑1 (MQ=255)                                                                                                     ccAGATATTTACCCTGTCGCGTTATGCGGGCATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGCTACGATCTCGc                                 >  2:281012/1‑123 (MQ=255)                                                                                                     ccAGATATTTACCCTGTCGCGTTATGCGGGCATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGCTACGATCTCGc                                 <  1:281012/123‑1 (MQ=255)                                                                                                              tACCCTGTCGCGTTATGCGGGCATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAAc       <  1:44513/140‑1 (MQ=255)                                                                                                               aCCCTGTCGCGTTATGCGGGCATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACt      <  1:352014/140‑1 (MQ=255)                                                                                                                   tGTCGCGTTATGCGGGCATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAgg  >  2:188145/1‑140 (MQ=255)                                                                                                                                   |                                                                                                                            TACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGG  >  NC_000913/842171‑842424

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