Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,920,963	T→C	C247C (TGT→TGC)	rsmF →	16S rRNA m(5)C1407 methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,920,963	0	T	C	100.0%	90.8 / NA	26	C247C (TGT→TGC)	rsmF	16S rRNA m(5)C1407 methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/14);  total (12/14)

GGTGCGTAAAGATCCCGATGCGCTAAAAAACTGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCTGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCAT  >  NC_000913/1920828‑1921098                                                                                                                                        |                                                                                                                                        ggtgCGTAAAGATCCCGATGCGCTAAAAAACTGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGCACCt                                                                                                                                     >  1:106234/1‑140 (MQ=255)        aaaGATCCCGATGCGCTAAAAAACTGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGCACCTTAAACCa                                                                                                                              <  1:100979/140‑1 (MQ=255)             tCCCGATGCGCTAAAAAACTGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGCACCTTAAACCAGgaag                                                                                                                         >  1:548698/1‑140 (MQ=255)              cccGATGCGCTAAAAAACTGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGCACCTTAAACCAGgaaga                                                                                                                        <  1:519399/140‑1 (MQ=255)                          aaaaaCTGGTCACCATAAAGCACTCCGGATATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGCACCTTAAACCAGGAAGAAAACGAAGCCGt                                                                                                            <  1:196000/140‑1 (MQ=255)                                tGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGCACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCt                                                                                                      >  1:661243/1‑140 (MQ=255)                                        aGAAAGCAATCAGGAAATCGCATCTACACAAGGGGGGCTGATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGCACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGa                                                                                              <  2:261329/140‑1 (MQ=255)                                          aaaGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGCACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGaaa                                                                                            <  2:548698/140‑1 (MQ=255)                                                aTCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGCACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACt                                                                                      <  1:6180/140‑1 (MQ=255)                                                 tCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGCACCATAAACCAGAAAGTAAAACAAGCCGTATGCCTATGGCTGAAAGAGACtt                                                                                     >  2:344833/1‑140 (MQ=255)                                                     gAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGCACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTAccc                                                                                 <  1:69770/140‑1 (MQ=255)                                                     gAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGCACCTTAAACCAGGAAGAAAACGAAGCCGTTAGCATGTGGCAGAAAGAGACTTAccc                                                                                 >  1:340010/1‑140 (MQ=255)                                                          cGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGCACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACg                                                                            <  2:608005/140‑1 (MQ=255)                                                                       cGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGCACCTTAAACCAGGAAGAAAACGAAGTCGTTTGCCTGTGGCAGAAAGAGACTAACCCAAACGCAGTAGAGttttt                                                               >  2:326488/1‑140 (MQ=255)                                                                        gggAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGCACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTa                                                              <  2:661243/140‑1 (MQ=255)                                                                         ggAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGCACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGGGACGTATCCCGACGCAGTAGAGTTTTTAc                                                             >  1:217321/1‑140 (MQ=255)                                                                               tATCGACAGCGCCTTGCATGCAGTATGTCCTGGTGGTGCGCTGGTTTACTCGACCTGCACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTg                                                       <  1:344833/140‑1 (MQ=255)                                                                                 tCGACAGCGCCTTTCATGCATTTCGTCCTGGAGGTTCGCCGGTTTACTCGACCTGCACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGc                                                     <  2:217321/140‑1 (MQ=255)                                                                                    aCAGCGCCTTTCATGCCGTACGTCCTGGTGGTAAGCTGGTTTACTCGACCTGCACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGAt                                                  <  2:340010/140‑1 (MQ=255)                                                                                                             tggtggTACGCTGGTTTACTCGACCTGCACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCTGGTGCAAACAAAgcgc                         >  2:695124/1‑140 (MQ=255)                                                                                                              ggtggtACGCTGGTTTACTCGACCTGCACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATGTCCTCCCTGCTGCAAACAAAGCGCt                        >  1:239120/1‑140 (MQ=255)                                                                                                                  gtACGCTGGTTTACTCGACCTGCACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCTGGTGc                                    >  2:399528/1‑124 (MQ=255)                                                                                                                  gtACGCTGGTTTACTCGACCTGCACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCTGGTGc                                    <  1:399528/124‑1 (MQ=255)                                                                                                                     cGCTGGTTTACTCGACCTGCACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCTGGTGCAAACAAAGCGCTGACCgaa                 >  2:331187/1‑140 (MQ=255)                                                                                                                            ttACTCGACCTGCACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCTGGTGCAAACAAAGCGCTGACCGAAGAAGGCt          <  2:106234/140‑1 (MQ=255)                                                                                                                                    ccTGCACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAGCAAAGCGCGGACCAAAGAAGGCTTTTTGCAt  >  2:234961/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                        GGTGCGTAAAGATCCCGATGCGCTAAAAAACTGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCTGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCAT  >  NC_000913/1920828‑1921098

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