Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,224,682	T→G	D285E (GAT→GAG)	metH →	cobalamin‑dependent methionine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,224,682	0	T	G	100.0%	149.6 / NA	42	D285E (GAT→GAG)	metH	cobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (20/22);  total (20/22)

GGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCT  >  NC_000913/4224604‑4224758                                                                               |                                                                             ggAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAgct                                                                            <  1:839940/81‑1 (MQ=255)  gAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAgctc                                                                           <  1:839941/81‑1 (MQ=255)  gAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAgctc                                                                           <  1:839942/81‑1 (MQ=255)   aGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCg                                                                          <  2:839943/81‑1 (MQ=255)     cTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGAc                                                                        <  1:839944/81‑1 (MQ=255)          aCGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCg                                                                    <  1:839945/80‑1 (MQ=255)           cGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGac                                                                  <  2:839946/81‑1 (MQ=255)           cGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGa                                                                   >  1:840015/1‑80 (MQ=255)             gATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGacac                                                                <  1:839947/81‑1 (MQ=255)                 gCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGAt                                                             <  1:839948/80‑1 (MQ=255)                  cGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATgg                                                           >  1:840004/1‑81 (MQ=255)                   ggAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATgg                                                           >  2:840017/1‑80 (MQ=255)                   ggAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGc                                                          <  1:839949/81‑1 (MQ=255)                     aaTGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCa                                                         >  1:840009/1‑80 (MQ=255)                      aTGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCaaa                                                       <  2:839950/81‑1 (MQ=255)                      aTGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCaa                                                        >  1:839991/1‑80 (MQ=255)                               cACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACg                                              >  2:840028/1‑81 (MQ=255)                                     gcACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATg                                        >  1:840007/1‑81 (MQ=255)                                     gcACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATg                                        >  1:840008/1‑81 (MQ=255)                                     gcACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATg                                        <  1:839951/81‑1 (MQ=255)                                      cACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATgg                                       <  2:839952/81‑1 (MQ=255)                                           gAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGgcgc                                   >  2:840012/1‑80 (MQ=255)                                           gAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCa                                  <  1:839953/81‑1 (MQ=255)                                               gCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCg                              >  2:840002/1‑81 (MQ=255)                                               gCCGGGCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCg                              <  1:839954/81‑1 (MQ=255)                                                  gggCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGt                           >  1:840020/1‑81 (MQ=255)                                                   ggCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGtt                          >  1:840001/1‑81 (MQ=255)                                                   ggCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGtt                          >  1:840024/1‑81 (MQ=255)                                                   ggCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGtt                          <  2:839955/81‑1 (MQ=255)                                                   ggCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGtt                          <  2:839956/81‑1 (MQ=255)                                                    gCTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGttt                         >  2:840016/1‑81 (MQ=255)                                                     cTACCCAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGtttt                        <  2:839957/81‑1 (MQ=255)                                                        cccAACGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTtctc                     <  2:839958/81‑1 (MQ=255)                                                           aaCGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCaa                   >  2:840032/1‑80 (MQ=255)                                                           aaCGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAat                  <  2:839959/81‑1 (MQ=255)                                                           aaCGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAat                  <  2:839960/81‑1 (MQ=255)                                                           aaCGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAat                  <  2:839961/81‑1 (MQ=255)                                                             cGCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAatat                >  2:840005/1‑81 (MQ=255)                                                              gCCTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATAtc               <  1:839962/81‑1 (MQ=255)                                                                cTTTGGTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATAtcgt             >  1:840033/1‑81 (MQ=255)                                                                    ggTGAGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTcggc         >  1:840031/1‑81 (MQ=255)                                                                        aGTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGgct     <  1:839963/81‑1 (MQ=255)                                                                         gTACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGTAAGCGGGTTTTCTCAATATCGTCGGCGgctg    >  2:840025/1‑81 (MQ=255)                                                                          tACGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGgctgc   <  1:839964/81‑1 (MQ=255)                                                                           aCGAGCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGgctgct  >  1:840022/1‑81 (MQ=255)                                                                               |                                                                             GGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCT  >  NC_000913/4224604‑4224758

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