Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,906,464	T→C	A199A (GCA→GCG)	rlmA ←	23S rRNA m(1)G745 methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,906,464	0	T	C	100.0%	140.4 / NA	39	A199A (GCA→GCG)	rlmA	23S rRNA m(1)G745 methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/25);  total (14/25)

ACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCAATAATGCGACGGCTTCATCACCGCGAAGACGCATCGGATAACACAACTCCGCACTCTGCTGTAATGTAAAACCTTCCAGTTGTTCTGCATGAGGTGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTA  >  NC_000913/1906330‑1906595                                                                                                                                       |                                                                                                                                    aCTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCAATAATGCGACGGCTTCATCACCGCGAAGACGCATCGGATAACACAACTCCGCACTCTGCTGTAATGTAAAACCTTCCAGTTGTTCTGCATGAGGCGCATg                                                                                                                                <  1:981474/140‑1 (MQ=255)            cGCACGCCAGGCAAACGTCGTCATCTGCAATCATGCGACGGCTTCATCGCCGCGAAGACGCATCGGATAACACAACTCCGCACTCTGCTGTAATGTAAAACCTTCCAGTTGTTCTGCATGAGGCGCATGAAGATGTACtt                                                                                                                     <  2:323282/140‑1 (MQ=255)            cGCACGCCAGGCAAACGGCGTCATCTGCAATAATGCGACGGCTTCATCACCGCGAAGACGCATCGGATAACACAACTCCGCACTCTGCTGTAATGTAAAACCTTCCAGTTGTTCTGCATGAGGCGCATGAAGATGTACtt                                                                                                                     <  1:57312/140‑1 (MQ=255)                  ccAGGCAAACGGCGTCATCTGCAATAATGCGACGGCTTCATCACCGCGAAGACGCATCGGATAACACAACTCCGCACTCTGCTGTAATGTAAAACCTTCCAGTTGTTCTGCATGAGGCGCATGAAGATGTACTTCATTGt                                                                                                               <  1:103007/140‑1 (MQ=255)                   cAGGCAAACGGCGTCATCTGCAATAATGCGACGGCTTCATCACCGCGAAGACGCATCGGATAACACAACTCCGCACTCTGCTGTAATGTAAAACCTTCCAGTTGTTCTGCATGAGGCGCATGAAGATGTACTTCATTGTa                                                                                                              >  2:545505/1‑140 (MQ=255)                          aCGGCGTCATCTGCAATAATGCGACGGCTTCATCACCGCGAAGACGCATCGGATAACACAACTCCGCACTCTGCTGTAATGTAAAACCTTCCAGTTGTTCTGCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGc                                                                                                       <  1:228691/140‑1 (MQ=255)                                 cATCTGCAATAATGCGACGGCTTCATCACCGCGAAGACGCATCGGATAACACAACTCCGCACTCTGCTGTAATGTAAAACCTTCCAGTTGTTCTGCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCa                                                                                                <  2:729956/140‑1 (MQ=255)                                          taatGCGACGGCTTCATCACCGCGAAGACGCATCGGATAACACAACTCCGCACTCTGCTGTAATGTAAAACCTTCCAGTTGTTCTGCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCa                                                                                       >  1:561116/1‑140 (MQ=255)                                                 aCGGCTTCATCACCGCGAAGACGCATCGGATAACACAACTCCGCACTCTGCTGTAATGTAAAACCTTCCAGTTGTTCTGCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTcg                                                                                <  2:561116/140‑1 (MQ=255)                                                 aCGGCTTCATCACCGCGAAGACGCATCGGATAACACAACTCCGCACTCTGCTGTAATGTAAAACCTTCCAGTTGTTCTGCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTcg                                                                                >  1:877654/1‑140 (MQ=255)                                                 aCGGCTTCATCACCGCGAAGACGCATCGGATAACACAACTCCGCACTCTGCTGTAATGTAAAACCTTCCAGTTGTTCTGCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTcg                                                                                <  2:842390/140‑1 (MQ=255)                                                 aCGGCTTCATC‑CCGCGAAGACGCATCGGATAACACAACTCCGCACTCTGCTGTAATGTAAAACCTTCCAGTTGTTCTGCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTcgc                                                                               >  2:753939/1‑140 (MQ=255)                                                       tcatcaCCGCGAAGACGCATCGGATAACACAACTCCGCACTCTGCTGTAATGTAAAACCTTCCAGTTGTTCTGCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTcc                                                                          <  2:685392/140‑1 (MQ=255)                                                          atcccGCGAAGACGCATCGGATAACACAACTCCGCACTCTGCTGTAATGTAAAACCTTCCAGTTGTTCTGCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCgg                                                                       <  1:753939/137‑1 (MQ=255)                                                              cgcgAAGACGCATCGGAAAACATAACTCCGCACTCTGCTTTAATGTAAAACCTTCCAGTTGTTCTGCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTg                                                                   <  2:450716/140‑1 (MQ=255)                                                               gcgAAGACGCATCGGATAACACAACTCCGCACTCTGCTGTAATGTAAAACCTTCCAGTTGTTCTGCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTgg                                                                  <  2:684381/140‑1 (MQ=255)                                                               gcgAAGACGCATCGGATAACACAACTCCGCACTCTGCTGTAATGTAAAACCTTCCAGTTGTTCTGCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTgg                                                                  <  2:58446/140‑1 (MQ=255)                                                               gcgAAGACGCATCGGATAACACAACTCCGCACTCTGCTGTAATGTAAAACCTTCCAGTTGTTCTGCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTgg                                                                  <  2:943768/140‑1 (MQ=255)                                                                                      aCTCCGCACTCTGCTGTAATGTAAAACCTTCCAGTTGTTCTGCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGGGc                                           <  2:877654/140‑1 (MQ=255)                                                                                       cTCCGCACTCTGCTGTAATGTAAAACCTTCCAGTTGTTCTGCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGGGCt                                          >  2:975369/1‑140 (MQ=255)                                                                                          cGCACTCTGCTGTAATGTAAAACCTTCCAGTTGTTCTGCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGGGCTTCa                                       >  1:125222/1‑140 (MQ=255)                                                                                           gCACTCTGCTGTAATGTAAAACCTTCCAGTTGTTCTGCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGGGCTTCac                                      <  1:545505/140‑1 (MQ=255)                                                                                               tctGCTGTAATGTAAAACCTTCCAGTTGTTCTGCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGGGCTTCactact                                  <  1:522361/140‑1 (MQ=255)                                                                                                 tgctgTAATGTAAAACCTTCCAGTTGTTCTGCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGGGCTTCACTACTCg                                >  1:737663/1‑140 (MQ=255)                                                                                                 tgctgTAATGTAAAACCTTCCAGTTGTTCTGCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGGGCTTCACTACTCg                                >  2:76919/1‑140 (MQ=255)                                                                                                  gctgTAATGTAAAACCTTCCAGTTGTTCTGCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGGGCTTCACTACTCGt                               >  1:216781/1‑140 (MQ=255)                                                                                                           taaAACCTTCCAGTTGTTCTGCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGGGCTTCACTACTCGTGCTAAttct                      >  2:505001/1‑140 (MQ=255)                                                                                                           taaAACCTTCCAGTTGTTCTGCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGGGCTTCACTACTCGTGCTAAttct                      >  2:849379/1‑140 (MQ=255)                                                                                                           taaAACCTTCCAGTTGTTCTGCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAAATACCCAGGCGCAGGTCTTCACCACTCGTGCTAAttct                      >  1:494561/1‑140 (MQ=255)                                                                                                            aaaaCCTTCCAGTTGTTCTGCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGGGCTTCACTACTCGTGCTAAttctt                     <  2:990880/140‑1 (MQ=255)                                                                                                            aaaaCCTTCCAGTTGTTCTGCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGGGCTTCACTACTCGTGCTAAttctt                     <  2:216781/140‑1 (MQ=255)                                                                                                            aaaaCCTTCCAGTTGTTCTGCATGAAGCGCATGAAGAAGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGGGCTTCACTACTCGTGCTAAttctt                     <  1:416611/140‑1 (MQ=255)                                                                                                             aaaCCTTCCAGTTGTTCTGCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGGGCTTCACTACTCGTGCTAAttcttc                    >  1:147766/1‑140 (MQ=255)                                                                                                                     cAGTTGTTCTGCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCCCCGGGCTGCACAACGCGTGCTGATTCTTCTGCTTTac            >  1:306596/1‑140 (MQ=255)                                                                                                                      aGTTGTTCTGCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGGGCTTCACTACTCGTGCTAATTCTTCTGCTTTaca           <  2:147766/140‑1 (MQ=255)                                                                                                                      aGTTGTTCTGCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGGGCTTCACTACTCGTGCTAATTCTTCTGCTTTaca           <  1:76919/140‑1 (MQ=255)                                                                                                                       gttgttCTGCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGGGCTTCACTACTCGTGCTAATTCTTCTGCTTTacac          <  1:849379/140‑1 (MQ=255)                                                                                                                        ttgttCTGCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACg         <  1:552916/140‑1 (MQ=255)                                                                                                                               gCATGAGGCGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTa  <  1:103121/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                    ACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCAATAATGCGACGGCTTCATCACCGCGAAGACGCATCGGATAACACAACTCCGCACTCTGCTGTAATGTAAAACCTTCCAGTTGTTCTGCATGAGGTGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTA  >  NC_000913/1906330‑1906595

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: ‑