Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,441,814	T→C	G17G (GGT→GGC)	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,441,814	0	T	C	100.0%	127.6 / NA	36	G17G (GGT→GGC)	mnmC	fused 5‑methylaminomethyl‑2‑thiouridine‑forming enzyme methyltransferase and FAD‑dependent demodification enzyme
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/19);  total (17/19)

ACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATCCACATCCTCTGTTTGTGGTAG  >  NC_000913/2441666‑2441950                                                                                                                                                     |                                                                                                                                         aCGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGCac                                                                                                                                        >  1:187942/1‑151 (MQ=255)    ccGAACAAGCCCGATCAGCCATTTAATAGAGAAATTTGCACAGCCTAACACACATCGCTAAGATCGACCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGCac                                                                                                                                        <  2:73160/148‑1 (MQ=255)      gAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGCACACCTg                                                                                                                                   >  2:311347/1‑151 (MQ=255)        aCAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGCACACCTGtt                                                                                                                                 >  1:191619/1‑151 (MQ=255)          aaGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGCACACCTGttt                                                                                                                                <  1:124641/150‑1 (MQ=255)                      ccATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGCACACCTGTTTCCCGAGAtttt                                                                                                                     <  1:509926/149‑1 (MQ=255)                           tAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGCACACCTGTTTCCCGAGATTTTGACGAt                                                                                                               <  1:484787/150‑1 (MQ=255)                           tAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGCACACCTGTTTCCCGAGATTTTGACGATg                                                                                                              >  2:121051/1‑151 (MQ=255)                            aaTAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGCACACCTGTTTCCCGAGATTTTGACGATGt                                                                                                             >  1:293152/1‑151 (MQ=255)                              tAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGCACACCTGTTTCCCGAGATTTTGACGATGTCt                                                                                                           >  1:175134/1‑151 (MQ=255)                               agagAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGCACACCTGTTTCCCGAGATTTTGACGATGTCt                                                                                                           >  1:310055/1‑150 (MQ=255)                                gagaAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGCACACCTGTTTCCCGAGATTTTGACGATGTCTAt                                                                                                         >  1:628555/1‑151 (MQ=255)                                    aaTTTGCGCAGCCTTACCCACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGCACACCTGTTTCCCGAGATTTTGACGATGTCTAttttt                                                                                                     <  1:424188/151‑1 (MQ=255)                                      tttGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGCACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTc                                                                                                    <  2:575254/150‑1 (MQ=255)                                        tGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGCACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCa                                                                                                  >  2:569562/1‑150 (MQ=255)                                        tGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGCACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCa                                                                                                  >  2:304186/1‑150 (MQ=255)                                                       acaTCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGCACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTgg                                                                                   >  1:199402/1‑150 (MQ=255)                                                           cGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGCACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAagag                                                                              <  1:159948/151‑1 (MQ=255)                                                                  aTCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGCACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGtt                                                                       >  2:249502/1‑151 (MQ=255)                                                                        ccACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGCACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGttt                                                                 <  1:395665/151‑1 (MQ=255)                                                                          aCCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGCACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGttttt                                                               <  1:195408/151‑1 (MQ=255)                                                                          aCCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGCACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGtt                                                                  <  1:203713/148‑1 (MQ=255)                                                                             gCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGCACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTg                                                            >  2:447785/1‑151 (MQ=255)                                                                                 gTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGCACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTggg                                                          <  1:158722/149‑1 (MQ=255)                                                                                      aCGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGCACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAAc                                                   >  2:436958/1‑151 (MQ=255)                                                                                                 aCGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGCACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGt                                        <  1:369968/151‑1 (MQ=255)                                                                                                                    aCAACCTGCCAACCTCGAATTTAATGCTGAGGGCACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATc                       >  2:189607/1‑149 (MQ=255)                                                                                                                    aCAACCTGCCAACCTCGAATTTAATGCTGAGGGCACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATCca                     <  2:323902/151‑1 (MQ=255)                                                                                                                     cAACCTGCCAACCTCGAATTTAATGCTGAGGGCACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATCcac                    <  2:532417/151‑1 (MQ=255)                                                                                                                        ccTGCCAACCTCGAATTTAATGCTGAGGGCACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATCCACATc                 <  1:404618/151‑1 (MQ=255)                                                                                                                          tGCCAACCTCGAATTTAATGCTGAGGGCACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATCCACATCct               >  2:594298/1‑151 (MQ=255)                                                                                                                           gCCAACCTCGAATTTAATGCTGAGGGCACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATCCACATCct               <  1:522365/150‑1 (MQ=255)                                                                                                                                   cGAATTTAATGCTGAGGGCACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATCCACATCCTCTGTTtgt       <  1:281193/150‑1 (MQ=255)                                                                                                                                    gAATTTAATGCTGAGGGCACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATCCACATCCTCTGTTTGTgg     <  2:630249/151‑1 (MQ=255)                                                                                                                                      aTTTAATGCTGAGGGCACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATCCACATCCTCTGTTTGTGGt    >  1:281214/1‑150 (MQ=255)                                                                                                                                        ttAATGCTGAGGGCACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATCCACATCCTCTGTTTGTGGTAg  <  1:379225/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                         ACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAGACGCGTTATGTTTTTCTGGGAGGCAACCAATTAGAGGTACGCTTTCCTGAGCATCCACATCCTCTGTTTGTGGTAG  >  NC_000913/2441666‑2441950

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