Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,921,492	T→C	W424R (TGG→CGG)	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,921,492	0	T	C	100.0%	176.9 / NA	50	W424R (TGG→CGG)	rsmF	16S rRNA m(5)C1407 methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/25);  total (25/25)

TCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTG  >  NC_000913/1921355‑1921626                                                                                                                                          |                                                                                                                                       tCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCgg                                                                                                                                      >  1:624934/1‑140 (MQ=255) tCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCgg                                                                                                                                      >  1:218726/1‑140 (MQ=255) tCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCgg                                                                                                                                      >  1:225270/1‑140 (MQ=255)       tttCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATcgc                                                                                                                                <  2:1011242/140‑1 (MQ=255)            cGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGgcg                                                                                                                           <  2:225270/140‑1 (MQ=255)            cGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGgcg                                                                                                                           <  2:218726/140‑1 (MQ=255)            cGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGgcg                                                                                                                           <  2:202442/140‑1 (MQ=255)              gTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGcgcg                                                                                                                         <  2:387982/140‑1 (MQ=255)              gTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGcgcg                                                                                                                         <  1:882115/140‑1 (MQ=255)                  gggATTAAACTTTCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGt                                                                                                                     >  1:875324/1‑140 (MQ=255)                              gCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGc                                                                                                         >  1:374842/1‑140 (MQ=255)                              gCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGc                                                                                                         >  1:447589/1‑140 (MQ=255)                                      gCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAg                                                                                                 <  2:548255/140‑1 (MQ=255)                                       cacaACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGt                                                                                                <  2:447589/140‑1 (MQ=255)                                          aacaaAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGc                                                                                             <  2:624934/140‑1 (MQ=255)                                            caaAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCgg                                                                                           <  2:531245/140‑1 (MQ=255)                                               aGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGgatg                                                                                        <  2:875324/140‑1 (MQ=255)                                                ggTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGgatga                                                                                       >  1:748178/1‑140 (MQ=255)                                                 gTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGgatgat                                                                                      <  2:748178/140‑1 (MQ=255)                                                  ttATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGatgatg                                                                                     >  1:553064/1‑140 (MQ=255)                                                      cGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTAtt                                                                                 >  1:197673/1‑140 (MQ=255)                                                        cTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTgg                                                                               <  2:1012462/140‑1 (MQ=255)                                                         tGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGt                                                                              <  2:611844/140‑1 (MQ=255)                                                           gcagcaTGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGc                                                                                             <  2:842087/123‑1 (MQ=255)                                                           gcagcaTGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGc                                                                                             >  1:842087/1‑123 (MQ=255)                                                           gcagcaTGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGgatgat                                                                                      >  1:526577/1‑130 (MQ=255)                                                           gcagcaTGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGgatgat                                                                                      <  2:526577/130‑1 (MQ=255)                                                           gcagcaTGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTa                                                                            <  2:197673/140‑1 (MQ=255)                                                           gcagcaTGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTa                                                                            >  1:952667/1‑140 (MQ=255)                                                                    aGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGc                                                                   <  1:948547/140‑1 (MQ=255)                                                                      cGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCAt                                                                 <  2:374842/140‑1 (MQ=255)                                                                         ttATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCaa                                                              >  1:124955/1‑140 (MQ=255)                                                                             tGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGa                                                          >  2:147775/1‑140 (MQ=255)                                                                                  ttGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGa                                                          >  1:503421/1‑135 (MQ=255)                                                                                  ttGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGa                                                          <  2:503421/135‑1 (MQ=255)                                                                                  ttGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGt                                                     >  1:415960/1‑140 (MQ=255)                                                                                   tGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGtt                                                    >  2:421428/1‑140 (MQ=255)                                                                                   tGCCACACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGtt                                                    >  1:105502/1‑140 (MQ=255)                                                                                          cccGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGa                                                          >  2:503011/1‑127 (MQ=255)                                                                                          cccGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGa                                                          <  1:503011/127‑1 (MQ=255)                                                                                          cccGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCaa                                             <  2:952667/140‑1 (MQ=255)                                                                                                        tGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTcgc                               >  1:450345/1‑140 (MQ=255)                                                                                                         gCATTCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTcgcg                              >  1:1064333/1‑140 (MQ=255)                                                                                                            ttCGAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGAtt                           >  1:629476/1‑140 (MQ=255)                                                                                                               gAGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGaa                        <  2:415960/140‑1 (MQ=255)                                                                                                                aGCTGACACCGCAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGaaa                       >  2:397231/1‑140 (MQ=255)                                                                                                                           cAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGAt                            <  1:68267/124‑1 (MQ=255)                                                                                                                           cAGGAAGCGGAGGAGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGAt                            >  2:68267/1‑124 (MQ=255)                                                                                                                                   ggaggaGCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAAc    <  1:1024570/140‑1 (MQ=255)                                                                                                                                     aggagCGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTg  >  2:17270/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                       TCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTG  >  NC_000913/1921355‑1921626

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