Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,014,886	A→G	S612G (AGC→GGC)	metE →	5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,014,886	0	A	G	100.0%	150.1 / NA	43	S612G (AGC→GGC)	metE	5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/19);  total (24/19)

GTTCCCGCGTGAAGATGTCAGCCGTGAAACCATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGAC  >  NC_000913/4014738‑4015032                                                                                                                                                     |                                                                                                                                                   gTTCCCGCGTGAAGATGTCAGCCGTGAAACCATCGCCAAACAGAT‑‑AGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTGGCg                                                                                                                                                 >  1:82465/1‑150 (MQ=255)      cgcgTGAAGATGTCAGCCGTGAAACCATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTGGCGACTg                                                                                                                                             <  1:20396/151‑1 (MQ=255)       gcgTGAAGATGTCAGCCGTGAAACCATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTGGCGACTgg                                                                                                                                            >  2:304570/1‑151 (MQ=255)                 gTCAGCCGTGAAACCATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTGGCGACTGGGATGCGTAtc                                                                                                                                  <  1:438080/151‑1 (MQ=255)                   cAGCCGTGAAACCATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTGGCGACTGGGATGCGTAtctc                                                                                                                                <  2:182401/151‑1 (MQ=255)                     gCCGTGAAACCATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTGGCGACTGGGATGCGTATCTCCa                                                                                                                              <  1:628473/151‑1 (MQ=255)                              ccATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGt                                                                                                                     >  2:305462/1‑151 (MQ=255)                               cATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGTa                                                                                                                    <  2:190976/151‑1 (MQ=255)                                aTCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTGGCGACTGGGATGCGTAGCTCCAGTGGGGCGTag                                                                                                                   <  1:85482/151‑1 (MQ=255)                                      aaaCAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCt                                                                                                             >  1:130638/1‑151 (MQ=255)                                            aTTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGACGGTTTACCGCTGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGt                                                                                                        >  1:714805/1‑150 (MQ=255)                                                     gcgctgCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAAcgccg                                                                                              >  2:6019/1‑151 (MQ=255)                                                     gcgctgCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAAcgc                                                                                                <  1:424459/149‑1 (MQ=255)                                                       gctgCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACgccgcc                                                                                            >  1:346877/1‑151 (MQ=255)                                                          gCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGccgccg                                                                                           <  2:752305/149‑1 (MQ=255)                                                                tgaAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGaa                                                                                    >  2:502359/1‑150 (MQ=255)                                                                 gaAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGaaa                                                                                   <  1:531259/150‑1 (MQ=255)                                                                         cGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGacac                                                                           >  2:206128/1‑150 (MQ=255)                                                                              tGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAAt                                                                     >  1:523721/1‑151 (MQ=255)                                                                                gAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATcc                                                                   >  2:368306/1‑151 (MQ=255)                                                                                gAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATcc                                                                   <  2:742580/151‑1 (MQ=255)                                                                                 aaGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATcc                                                                   >  1:543716/1‑150 (MQ=255)                                                                                  aGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCcac                                                                 <  2:109391/151‑1 (MQ=255)                                                                                  aGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCca                                                                  >  2:688537/1‑150 (MQ=255)                                                                                   gCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCcaca                                                                <  1:144125/151‑1 (MQ=255)                                                                                        tGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTcac                                                           >  1:768538/1‑151 (MQ=255)                                                                                        tGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTcac                                                           >  1:96821/1‑151 (MQ=255)                                                                                             ttGGCATAATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACAtgtg                                                      >  1:742171/1‑151 (MQ=255)                                                                                                gCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTa                                                   <  2:769783/151‑1 (MQ=255)                                                                                                 catcatCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTAt                                                  >  2:298924/1‑151 (MQ=255)                                                                                                    catcCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGc                                               >  1:289868/1‑151 (MQ=255)                                                                                                                     cGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGAcatcat                              <  2:655320/151‑1 (MQ=255)                                                                                                                           tGCGCGAAGGTTTACCGCTGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTc                        >  1:498407/1‑151 (MQ=255)                                                                                                                              gcgAAGGTTTACCGCTGCGTAGTGGCGCCTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGAt                     >  2:538633/1‑151 (MQ=255)                                                                                                                                gAAGGTTTACCGCTGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGa                      <  2:214980/148‑1 (MQ=255)                                                                                                                                gAAGGTTTACCGCTGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGa                      >  1:214948/1‑148 (MQ=255)                                                                                                                                   ggTTTACCGCTGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTgcgg                <  1:663581/151‑1 (MQ=255)                                                                                                                                       tACCGCTGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCt            >  1:137236/1‑151 (MQ=255)                                                                                                                                          cGCTGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGa         <  2:137269/151‑1 (MQ=255)                                                                                                                                          cGCTGCGTCGTGGCGACGGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTgg          >  1:766845/1‑150 (MQ=255)                                                                                                                                           gCTGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGAc        <  2:389010/151‑1 (MQ=255)                                                                                                                                             tGCGTCGTGGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGc      >  2:415416/1‑151 (MQ=255)                                                                                                                                                  cgtGGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGAc  <  2:703071/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GTTCCCGCGTGAAGATGTCAGCCGTGAAACCATCGCCAAACAGATTGCGCTGGCGCTGCGTGATGAAGTGGCCGATCTGGAAGCCGCTGGAATTGGCATCATCCAGATTGACGAACCGGCGCTGCGCGAAGGTTTACCGCTGCGTCGTAGCGACTGGGATGCGTATCTCCAGTGGGGCGTAGAGGCCTTCCGTATCAACGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGAC  >  NC_000913/4014738‑4015032

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