Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,238,421	T→C	E290G (GAG→GGG)	mglA ←	methyl‑galactoside ABC transporter ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,238,421	0	T	C	100.0%	154.1 / NA	41	E290G (GAG→GGG)	mglA	methyl‑galactoside ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (23/18);  total (23/18)

TTCGTTGGCATTATGGTTATTGATCTGTTTGCCGTGCAACGTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCTCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAACC  >  NC_000913/2238288‑2238555                                                                                                                                      |                                                                                                                                       ttCGTTGGCATTATGGTTATTGATCTGTTTGCCGTGCAACGTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCCCCCCt                                                                                                                                   <  2:311735/139‑1 (MQ=255)          aTTATGGTTATTGATCTGTTTGCCGTGCAACGTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCCCCCCTTTATGCAGa                                                                                                                          <  2:1068068/139‑1 (MQ=255)                        tCTGTTTGCCGTGCAACGTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATc                                                                                                            >  2:838185/1‑139 (MQ=255)                          tGTTTGCCGTGCAACGTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATcgc                                                                                                          >  1:977559/1‑139 (MQ=255)                            tttGCCGTGCAACGTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGa                                                                                                        >  2:24288/1‑139 (MQ=255)                              tGCCGTGCAACGTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAAt                                                                                                      >  2:869453/1‑139 (MQ=255)                              tGCCGTGCAACGTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAAt                                                                                                      >  2:213359/1‑139 (MQ=255)                                 cGTGCAACGTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGa                                                                                                   >  1:10054/1‑139 (MQ=255)                                   tGCAACGTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACg                                                                                                 >  1:235562/1‑139 (MQ=255)                                   tGCAACGTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACg                                                                                                 >  1:253594/1‑139 (MQ=255)                                   tGCAACGTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACg                                                                                                 >  2:428115/1‑139 (MQ=255)                                             tGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGCGCAg                                                                                       <  2:253594/139‑1 (MQ=255)                                                  ccAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACg                                                                                  <  2:880670/139‑1 (MQ=255)                                                   cAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGt                                                                                 >  1:329118/1‑139 (MQ=255)                                                    aGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTc                                                                                <  1:738266/139‑1 (MQ=255)                                                                      aCCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAgg                                                              <  2:977559/139‑1 (MQ=255)                                                                                  ctcAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTAcc                                                                    <  2:1061305/121‑1 (MQ=255)                                                                                  ctcAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTAcc                                                                    >  1:1061305/1‑121 (MQ=255)                                                                                    caacaaTATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACTTACCTCGAGGATGACTTCGCCCgg                                                >  1:864594/1‑139 (MQ=255)                                                                                    caacaaTATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCgg                                                <  2:664072/139‑1 (MQ=255)                                                                                          tatCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTc                                                                                <  1:88269/101‑1 (MQ=255)                                                                                          tatCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTc                                                                                >  2:88269/1‑101 (MQ=255)                                                                                          tatCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGgtttgtt                                          >  2:579954/1‑139 (MQ=255)                                                                                               gtgcGTTTCGCCCCCCCCAGACGGGCAATACCGAGGATACCCCCTTTATGCAGGTCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCtt                                     >  1:764290/4‑139 (MQ=255)                                                                                               gTACGTTTCGCCCCCACCAGACCGGCAGTACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGAATGACTTCGCCCGGTTTGTTTTCtt                                     >  2:717496/1‑139 (MQ=255)                                                                                               gTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGTGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCtt                                     >  1:731713/1‑139 (MQ=255)                                                                                               gTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCtt                                     <  1:593882/139‑1 (MQ=255)                                                                                               gTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCtt                                     >  2:501537/1‑139 (MQ=255)                                                                                               gTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCtt                                     >  2:593882/1‑139 (MQ=255)                                                                                                 aCGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACgag                                                                                                                 <  2:176127/61‑1 (MQ=255)                                                                                                 aCGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACgag                                                                                                                 >  1:176127/1‑61 (MQ=255)                                                                                                         cccccACCAGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCTTTGTCAGGGaa                           >  1:804595/1‑139 (MQ=255)                                                                                                            ccaccaGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCTTTGTCAGGGAAACg                        <  2:804595/139‑1 (MQ=255)                                                                                                            ccaccaGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCTTTGTCAGGGAAACg                        <  2:10054/139‑1 (MQ=255)                                                                                                            ccaccaGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCTTTGTCAGGGAAACg                        <  1:808616/139‑1 (MQ=255)                                                                                                            ccaccaGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCTTTGTCAGGGAAACg                        <  1:579954/139‑1 (MQ=255)                                                                                                            ccaccaGACCGGCAATACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCTTTGTCAGGGAAACg                        <  1:213359/139‑1 (MQ=255)                                                                                                                          aTACCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCTTTGTCAGGGAAACGCTGGTTAAGAGAAc          <  1:838185/139‑1 (MQ=255)                                                                                                                            aCCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGgtttgttt                                         <  1:994066/106‑1 (MQ=255)                                                                                                                            aCCGAGGATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGgtttgttt                                         >  2:994066/1‑106 (MQ=255)                                                                                                                                  gATCCCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATCACTTCGCCCGGTTTGTTTTCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAAcc  >  1:657590/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                       TTCGTTGGCATTATGGTTATTGATCTGTTTGCCGTGCAACGTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCTCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAACC  >  NC_000913/2238288‑2238555

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