Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,215,724 T→G W657G (TGG→GGG)  metH → methionine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,215,7240TG100.0% 91.1 / NA 26W657G (TGG→GGG) metHmethionine synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (14/12);  total (14/12)

GAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCC  >  NZ_CP009273/4215640‑4215805
                                                                                    |                                                                                 
gAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGaa                                                                              >  1:86785/1‑90 (MQ=255)
        aCTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATaa                                                                      <  1:64812/90‑1 (MQ=255)
           ggAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACg                                                                   <  1:272805/90‑1 (MQ=255)
            gAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGt                                                                  >  1:276715/1‑90 (MQ=255)
            gAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGt                                                                  >  1:67637/1‑90 (MQ=255)
              gCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCt                                                                >  2:221195/1‑90 (MQ=255)
              gCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCt                                                                >  1:296681/1‑90 (MQ=255)
              gCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCt                                                                >  2:115266/1‑90 (MQ=255)
                ttGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGATCGGGGGAAGTGAATAAGCGGCTgg                                                              >  1:485455/1‑90 (MQ=255)
                 tGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGa                                                             >  1:285522/1‑90 (MQ=255)
                         aaTATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGc                                                     <  1:486266/90‑1 (MQ=255)
                                 ggCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCaaa                                             <  2:145083/90‑1 (MQ=255)
                                 ggCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCaaa                                             <  2:296681/90‑1 (MQ=255)
                                 ggCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCaaa                                             <  2:108260/90‑1 (MQ=255)
                                    agcaAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTCGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGc                                          <  2:548479/90‑1 (MQ=255)
                                                     cGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGa                         >  1:415627/1‑90 (MQ=255)
                                                      gCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAg                        >  2:266988/1‑90 (MQ=255)
                                                       ccAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGc                       <  2:359064/90‑1 (MQ=255)
                                                         aaCGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAg                     <  1:406443/90‑1 (MQ=255)
                                                           cGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGa                   <  1:266988/90‑1 (MQ=255)
                                                                 gcagGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTg                                                               <  2:124057/40‑1 (MQ=39)
                                                                 gcagGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTg                                                               >  1:124057/1‑40 (MQ=255)
                                                                         aGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGccc     >  2:94603/1‑90 (MQ=255)
                                                                          gTGGAGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAg                     <  1:453231/73‑1 (MQ=255)
                                                                          gTGGAGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAg                     >  2:453231/1‑73 (MQ=255)
                                                                            ggCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGcc  >  2:249282/1‑90 (MQ=255)
                                                                                    |                                                                                 
GAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCC  >  NZ_CP009273/4215640‑4215805

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: