Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,442,490	+G	coding (727/2007 nt)	mnmC →	fused 5‑methylaminomethyl‑2‑thiouridine‑forming enzyme methyltransferase and FAD‑dependent demodification enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,442,486	1	.	G	100.0%	147.6 / NA	43	C241W (TGC→TGG)	mnmC	fused 5‑methylaminomethyl‑2‑thiouridine‑forming enzyme methyltransferase and FAD‑dependent demodification enzyme
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (23/20);  total (23/20)

GCCATGGCAAGGTTGGCGCGTCCGGGCGGCACGCTGGCGACATTTACGTCTGCCGGTTTTGTCCGCCGCGGTTTGCAGGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGC‑GGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGG  >  NC_000913/2442349‑2442607                                                                                                                                           |                                                                                                                          gCCATGGCAAGGTTGGCGCGTCCGGGCGGCACGCTGGCGACATTTACGTCTGCCGGTTTTGTCCGCCGCGGTTTGCAGGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGC‑g                                                                                                                          <  2:339127/139‑1 (MQ=255)   cATGGCAAGGTTGGCGCGTCCGGGCGGCACGCTGGCGACATTTACGTCTGCCGGTTTTGTCCGCCGCGGTTTGCAGGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGC‑ggg                                                                                                                        <  2:1447379/139‑1 (MQ=255)                   cgTCCGGGCGGCACGCTGGCGACATTTACGTCTGCCGGTTTTGTCCGCCGCGGTTTGCAGGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACa                                                                                                         >  1:1855975/1‑139 (MQ=255)                      ccGGGCGGCACGCTGGCGACATTTACGTCTGCCGGTTTTGTCCGCCGCGGTTTGCAGGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGGTGATGGa                                                                                                                <  1:849740/129‑1 (MQ=255)                      ccGGGCGGCACGCTGGCGACATTTACGTCTGCCGGTTTTGTCCGCCGCGGTTTGCAGGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGGTGATGGa                                                                                                                >  2:849740/1‑129 (MQ=255)                      ccGGGCGGCACGCTGGCGACATTTACGTCTGCCGGTTTTGTCCGCCGCGGTTTGCAGGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTa                                                                                                      >  2:1357904/1‑139 (MQ=255)                         ggcggcACGCTGGCGACATTTACGTCTGCCGGTTTTGTCCGCCGCGGTTTGCATGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTACCg                                                                                                   <  1:1151306/139‑1 (MQ=255)                         ggcggcACGCTGGCGACATTTACGTCTGCCGGTTTTGTCCGCCGCGGTTTGCAGGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTACCg                                                                                                   <  1:559573/139‑1 (MQ=255)                                                   tGCCGGTTTTGTCCGCCGCGGTTTGCAGGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTaa                                                                         >  2:1853123/1‑139 (MQ=255)                                                     ccGGTTTTGTCCGCCGCGGTTTGCAGGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAAcc                                                                       >  2:1781081/1‑139 (MQ=255)                                                                  ccgcGGTTTGCAGGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGgcagcag                                                          >  2:1108103/1‑139 (MQ=255)                                                                   cgcgGTTTGCAGGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGAc                                                                                                          >  1:1230212/1‑90 (MQ=255)                                                                   cgcgGTTTGCAGGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGAc                                                                                                          <  2:1230212/90‑1 (MQ=255)                                                                             aGGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAg                                               <  1:1357904/139‑1 (MQ=255)                                                                               gACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTACCg                                                                                                   >  1:975485/1‑85 (MQ=255)                                                                               gACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTACCg                                                                                                   <  2:975485/85‑1 (MQ=255)                                                                                 cGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAgcggc                                           >  1:747663/1‑139 (MQ=255)                                                                                  gCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAgcggcg                                          >  1:1187931/1‑139 (MQ=255)                                                                                    cGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTACCg                                                                                                   <  2:1380985/80‑1 (MQ=255)                                                                                    cGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTACCg                                                                                                   >  1:1380985/1‑80 (MQ=255)                                                                                    cGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGAt                                        >  1:76807/1‑139 (MQ=255)                                                                                      gATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTa                                      >  2:1056857/1‑139 (MQ=255)                                                                                            cGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATcggcg                                <  2:76807/139‑1 (MQ=255)                                                                                            cGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATcggcg                                <  1:1056857/139‑1 (MQ=255)                                                                                              aTGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATcggcggg                              >  2:1537301/1‑138 (MQ=255)                                                                                                         tAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATcggcg                                >  1:688182/1‑126 (MQ=255)                                                                                                         tAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATcggcg                                <  2:688182/126‑1 (MQ=255)                                                                                                                 ttGGGCGCAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATcgg                                  <  1:1308283/116‑1 (MQ=255)                                                                                                                 ttGGGCGCAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATcgg                                  >  2:1308283/1‑116 (MQ=255)                                                                                                                 ttGGGCGCAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCgttgt           >  2:215177/1‑139 (MQ=255)                                                                                                                  tGGGCGCAAACGGGAAATGCTTTGCGGGGGTGATGGa                                                                                                                >  1:383308/1‑37 (MQ=255)                                                                                                                  tGGGCGCAAACGGGAAATGCTTTGCGGGGGTGATGGa                                                                                                                <  2:383308/37‑1 (MQ=255)                                                                                                                  tGGGCGCAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATcggcg                                >  1:230990/1‑117 (MQ=255)                                                                                                                  tGGGCGCAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATcggcg                                <  2:230990/117‑1 (MQ=255)                                                                                                                   gggCGCAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTAc                                                                                                     >  1:495686/1‑47 (MQ=39)                                                                                                                   gggCGCAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTAc                                                                                                     <  2:495686/47‑1 (MQ=39)                                                                                                                        cAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTACCGCTCCCCTTCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCt    >  1:1325732/1‑139 (MQ=255)                                                                                                                        cAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCGGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCt    <  1:1537301/139‑1 (MQ=255)                                                                                                                        cAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGc     >  1:1310513/1‑138 (MQ=255)                                                                                                                        cAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGc     <  2:1310513/138‑1 (MQ=255)                                                                                                                        cAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCt    <  1:215177/139‑1 (MQ=255)                                                                                                                        cAAACGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCt    <  2:747663/139‑1 (MQ=255)                                                                                                                          aaCGGGAAATGCTTTGCGGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTgg  <  2:277632/139‑1 (MQ=255)                                                                                                                                     cTTTGCGGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCt    <  2:1112613/126‑1 (MQ=255)                                                                                                                                     cTTTGCGGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCt    >  1:1112613/1‑126 (MQ=255)                                                                                                                                           |                                                                                                                          GCCATGGCAAGGTTGGCGCGTCCGGGCGGCACGCTGGCGACATTTACGTCTGCCGGTTTTGTCCGCCGCGGTTTGCAGGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGC‑GGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGG  >  NC_000913/2442349‑2442607

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