Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	842,236	A→C	intergenic (‑180/‑96)	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,236	0	A	C	100.0%	124.6 / NA	32	intergenic (‑180/‑96)	mcbA/rlmF	DUF1471 domain‑containing protein McbA/23S rRNA m(6)A1618 methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (18/14);  total (18/14)

GTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACC  >  NC_000913/842102‑842368                                                                                                                                       |                                                                                                                                     gTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATCGCCAt                                                                                                                                 >  1:276384/1‑140 (MQ=255)    gCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATCGCCATGtt                                                                                                                              <  2:732121/140‑1 (MQ=255)     ccTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATCGCCATGTTa                                                                                                                             <  1:561698/140‑1 (MQ=255)      cTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATGAACCATCGCCATGTTa                                                                                                                             <  1:341577/139‑1 (MQ=255)      cTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATCGCCATGTTa                                                                                                                             >  2:341577/1‑139 (MQ=255)              tGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATCGCCATGTTa                                                                                                                             >  2:550303/1‑131 (MQ=255)              tGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATCGCCATGTTa                                                                                                                             <  1:550303/131‑1 (MQ=255)              tGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATCGCCATGTTAAGGTTTTgg                                                                                                                    >  2:69680/1‑140 (MQ=255)               gCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATCGCCATGTTAAGGTTTTggg                                                                                                                   <  2:426510/140‑1 (MQ=255)                                 tGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATCGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCa                                                                                                 <  2:712976/140‑1 (MQ=255)                                                   tttttGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATCGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGata                                                                                             >  1:440422/1‑126 (MQ=255)                                                   tttttGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATCGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGata                                                                                             <  2:440422/126‑1 (MQ=255)                                                   tttttGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATCGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTg                                                                                    <  1:797496/135‑1 (MQ=255)                                                   tttttGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATCGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTg                                                                                    >  2:797496/1‑135 (MQ=255)                                                                        cAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATCGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTAcgccg                                                          >  2:499270/1‑140 (MQ=255)                                                                        cAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATCGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTAcgccg                                                          <  1:221091/140‑1 (MQ=255)                                                                        cAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATCGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTAcgccg                                                          >  2:410378/1‑140 (MQ=255)                                                                        cAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATCGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTAcgccg                                                          >  2:221091/1‑140 (MQ=255)                                                                        cAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATCGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTAcgccg                                                          >  1:600080/1‑140 (MQ=255)                                                                             tAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATCGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCgctgct                                                     <  2:276384/140‑1 (MQ=255)                                                                                  cAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATCGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAg                                                <  1:69680/140‑1 (MQ=255)                                                                                     aTCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATCGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGgag                                             >  1:576027/1‑140 (MQ=255)                                                                                          tGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATCGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGc                                        >  2:35729/1‑140 (MQ=255)                                                                                               tGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATCGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGt                                   >  1:364329/1‑140 (MQ=255)                                                                                                      tGCTGATTTTTATGATGAATTGCCATAAACCATCGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCa                            >  2:364808/1‑140 (MQ=255)                                                                                                           aTTTTTATGATGAATTGCCATAAACCATCGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAAc                       >  1:426589/1‑140 (MQ=255)                                                                                                           aTTTTTATGATGAATTGCCATAAACCATCGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAAc                       <  2:787031/140‑1 (MQ=255)                                                                                                                tATGATGAATTGCCATAAACCATCGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCgggg                  >  2:219215/1‑140 (MQ=255)                                                                                                                tATGATGAATTGCCATAAACCATCGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCgggg                  >  1:115441/1‑140 (MQ=255)                                                                                                                       aaTTGCCATAAACCATCGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATc           >  2:106632/1‑140 (MQ=255)                                                                                                                             cATAAACCATCGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCa     <  1:35729/140‑1 (MQ=255)                                                                                                                                aaaCCATCGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAAcc  <  1:364808/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                     GTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACC  >  NC_000913/842102‑842368

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