Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,013,566	A→G	K172E (AAA→GAA)	metE →	cobalamin‑independent homocysteine transmethylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,013,566	0	A	G	100.0%	127.8 / NA	37	K172E (AAA→GAA)	metE	cobalamin‑independent homocysteine transmethylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/19);  total (18/19)

GTGCCGGAGTTCGTTAAAGGCCAACAGTTCAAACTGACCTGGACGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGG  >  NC_000913/4013431‑4013698                                                                                                                                        |                                                                                                                                     gTGCCGGAGTTCGTTAAAGGCCAACAGTTCAAACTGACCTGGACGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGGAAgg                                                                                                                                  >  2:57843/1‑139 (MQ=255)  tGCCGGAGTTCGTTAAAGGCCAACAGTTCAAACTGACCTGGACGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGGAAGTg                                                                                                                                 >  1:349360/1‑140 (MQ=255)   gCCGGAGTTCGTTAAAGGCCAACAGTTCAAACTGACCTGGACGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGGAAGTGa                                                                                                                                >  2:49151/1‑140 (MQ=255)   gCCGGAGTTCGTTAAAGACCAACAGTTCAAACTGACCTGGACGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGGAAGTGa                                                                                                                                >  2:237113/1‑140 (MQ=255)     cGGAGTTCGTTAAAGGCCAACAGTTCAAACTGACCTGGACGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGGAAGTGaaa                                                                                                                              <  1:49151/140‑1 (MQ=255)       gAGTTCGTTAAAGGCCAACAGTTCAAACTGACCTGGACGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGGAAGTGAAAgg                                                                                                                            >  2:318150/1‑140 (MQ=255)             gTTAAAGGCCAACAGTTCAAACTGACCTGGACGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGAACa                                                                                                                      <  2:349360/140‑1 (MQ=255)                   ggCCAACAGTTCAAACTGACCTGGACGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGAACAATTTGa                                                                                                                <  2:86761/140‑1 (MQ=255)                    gCCAACAGTTCAAACTGACCTGGACGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGAACAATTTGAc                                                                                                               >  2:61787/1‑140 (MQ=255)                    gCCAACAGTTCAAACTGACCTGGACGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGAACAATTTGAc                                                                                                               <  1:476674/140‑1 (MQ=255)                      cAACAGTTCAAACTGACCTGGACGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGAACAATTTGACCg                                                                                                             >  1:424033/1‑140 (MQ=255)                                      ccTGGACGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAAc                                                                                             >  2:511265/1‑140 (MQ=255)                                                    ctgGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGt                                                                               <  1:237113/140‑1 (MQ=255)                                                    ctgGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGt                                                                               >  2:96320/1‑140 (MQ=255)                                                      ggACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGttt                                                                             <  1:477241/140‑1 (MQ=255)                                                       gACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTa                                                                            <  1:375296/140‑1 (MQ=255)                                                        aCGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTAt                                                                           <  1:404962/140‑1 (MQ=255)                                                           aaGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAg                                                                        >  1:150164/1‑140 (MQ=255)                                                            aGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTAt                                                                           <  2:382039/136‑1 (MQ=255)                                                            aGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTAt                                                                           >  1:382039/1‑136 (MQ=255)                                                              tGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCaa                                                                     >  2:154335/1‑140 (MQ=255)                                                                                           gTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCAt                                        <  2:123941/140‑1 (MQ=255)                                                                                           gTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCAt                                        <  1:265757/140‑1 (MQ=255)                                                                                           gTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCAt                                        <  2:424033/140‑1 (MQ=255)                                                                                              aaaCCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGa                                     >  1:1449/1‑140 (MQ=255)                                                                                               aaCCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAg                                    <  2:150164/140‑1 (MQ=255)                                                                                                     tgctgcTGGGGCCGGTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTa                              <  2:354712/140‑1 (MQ=255)                                                                                                     tgctgcTGGGGCCGGTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTa                              >  2:522754/1‑140 (MQ=255)                                                                                                     tgctgcTGGGGCCGGTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTa                              >  1:266922/1‑140 (MQ=255)                                                                                                     tgctgcTGGGGCCGGTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTa                              <  1:61787/140‑1 (MQ=255)                                                                                                            ggggCCGGTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATtg                       <  1:318150/140‑1 (MQ=255)                                                                                                             gggCCGGTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATtga                      >  1:421926/1‑140 (MQ=255)                                                                                                               gCCGGTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATtgatg                    >  2:131014/1‑140 (MQ=255)                                                                                                                   gTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAAcc                <  2:266922/140‑1 (MQ=255)                                                                                                                   gTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAAcc                <  1:211397/140‑1 (MQ=255)                                                                                                                    ttACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAAccc               >  2:194405/1‑140 (MQ=255)                                                                                                                    ttACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAAccc               >  2:148569/1‑140 (MQ=255)                                                                                                                               tggctgGGGGAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACt    <  1:57843/140‑1 (MQ=255)                                                                                                                                 gctgGGGGAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTgg  >  1:376168/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                     GTGCCGGAGTTCGTTAAAGGCCAACAGTTCAAACTGACCTGGACGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGG  >  NC_000913/4013431‑4013698

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