Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,920,205	(T)7→8	intergenic (+62/‑18)	yebT → / → rsmF	MCE domain protein/16S rRNA m(5)C1407 methyltransferase, SAM‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,920,198	1	.	T	100.0%	124.9 / NA	39	intergenic (+55/‑25)	yebT/rsmF	MCE domain protein/16S rRNA m(5)C1407 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (21/18);  total (21/18)

AGGCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑TTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCG  >  NC_000913/1920071‑1920330                                                                                                                                 |                                                                                                                                     aGGCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑ttttttttgcc                                                                                                                           <  1:1273051/139‑5 (MQ=255)  ggCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑ttttttttgc                                                                                                                            <  1:1135230/137‑4 (MQ=255)  ggCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTCCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑ttttttttgc                                                                                                                            >  2:1135230/1‑134 (MQ=255)    cAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGt                                                                                                                         >  1:1361650/1‑139 (MQ=255)     aaaCACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTg                                                                                                                        >  2:475402/1‑139 (MQ=255)             ccTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGc                                                                                                                >  2:56219/1‑139 (MQ=255)              cTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCt                                                                                                               <  2:1422263/139‑1 (MQ=255)                     aGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑ttttttttgc                                                                                                                            <  2:425779/118‑4 (MQ=255)                     aGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑ttttttttgc                                                                                                                            >  1:425779/1‑115 (MQ=255)                       gAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCg                                                                                                             <  2:433151/132‑1 (MQ=255)                       gAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCg                                                                                                             >  1:433151/1‑132 (MQ=255)                       gAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCa                                                                                                      >  1:1518979/1‑139 (MQ=255)                        aaaGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCaa                                                                                                     <  1:475402/139‑1 (MQ=255)                         aaGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAac                                                                                                    <  2:1518979/139‑1 (MQ=255)                              aaCCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACAcc                                                                                               >  2:1408696/1‑139 (MQ=255)                               aCCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCg                                                                                              >  1:1196107/1‑139 (MQ=255)                                         tGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCTGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTccc                                                                                    >  2:1667729/1‑139 (MQ=255)                                         tGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTccc                                                                                    >  1:50212/1‑139 (MQ=255)                                                    gggAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTc                                                                         <  1:56219/139‑1 (MQ=255)                                                     ggAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAacaca                                                                                                 <  1:1483864/114‑1 (MQ=255)                                                     ggAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAacaca                                                                                                 >  2:1483864/1‑114 (MQ=255)                                                     ggAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAaca                                                                                                   >  1:239577/1‑112 (MQ=255)                                                     ggAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAaca                                                                                                   <  2:239577/112‑1 (MQ=255)                                                          tGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGc                                                                              <  2:1077340/128‑1 (MQ=255)                                                          tGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGc                                                                              >  1:1077340/1‑128 (MQ=255)                                                            cgcTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGAcacat                                                                 >  2:102800/1‑138 (MQ=255)                                                             gcTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACaaa                                                                <  2:1361650/139‑1 (MQ=255)                                                                 cccAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAg                                                                                                                          <  1:1712340/77‑1 (MQ=255)                                                                 cccAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAg                                                                                                                          >  2:1712340/1‑77 (MQ=255)                                                                  ccAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATgcgc                                                           <  1:1408696/139‑1 (MQ=255)                                                                      ataatGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑ttttttttgcca                                                                                                                          <  1:1076473/71‑6 (MQ=255)                                                                      ataatGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑ttttttttgcca                                                                                                                          >  2:1076473/1‑66 (MQ=255)                                                                             ccACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGcc                                                >  1:1443358/1‑139 (MQ=255)                                                                                  gCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGa                                           <  2:50212/139‑1 (MQ=255)                                                                                           tGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGa                                                                                 <  2:1402830/92‑1 (MQ=255)                                                                                           tGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGa                                                                                 >  1:1402830/1‑92 (MQ=255)                                                                                             ccTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCAttt                                >  1:941681/1‑139 (MQ=255)                                                                                             ccTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑ttttttttgc                                                                                                                            >  2:1099408/1‑43 (MQ=255)                                                                                             ccTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑ttttttttgc                                                                                                                            <  1:1099408/46‑4 (MQ=255)                                                                                                 cgcCGGAGCGTTTATGCTAAACTGCGCGCCTG‑ttttttttgcc                                                                                                                           <  2:861918/43‑5 (MQ=255)                                                                                                 cgcCGGAGCGTTTATGCTAAACTGCGCGCCTG‑ttttttttgcc                                                                                                                           >  1:861918/1‑39 (MQ=255)                                                                                                       aGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTa                                                                                                                     >  1:1689222/1‑44 (MQ=38)                                                                                                       aGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTa                                                                                                                     <  2:1689222/44‑1 (MQ=255)                                                                                                          gTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTg                   <  2:17685/139‑1 (MQ=255)                                                                                                                         cgcgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAac                                                                                                    >  1:554143/1‑43 (MQ=38)                                                                                                                         cgcgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAac                                                                                                    <  2:554143/43‑1 (MQ=38)                                                                                                                         cgcgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCgg                                                                                  >  1:557561/1‑61 (MQ=255)                                                                                                                         cgcgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCgg                                                                                  <  2:557561/61‑1 (MQ=255)                                                                                                                         cgcgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGcc    >  1:984064/1‑139 (MQ=255)                                                                                                                           cgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCg  <  2:1220753/139‑1 (MQ=255)                                                                                                                                 |                                                                                                                                     AGGCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑TTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCG  >  NC_000913/1920071‑1920330

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