Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,717,817 C→T A238T (GCG→ACG)  yfiF ← putative methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,717,8170CT100.0% 105.4 / NA 30A238T (GCG→ACG) yfiFputative methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/13);  total (17/13)

CATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGC  >  NC_000913/2717739‑2717891
                                                                              |                                                                          
catcCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGtcc                                                                          >  1:436286/1‑81 (MQ=255)
    cAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGa                                                                      >  1:436293/1‑81 (MQ=255)
           ttAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCcagc                                                               <  2:436255/81‑1 (MQ=255)
              aCAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGt                                                            <  1:436256/81‑1 (MQ=255)
                 aTGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGcc                                                         <  1:436257/81‑1 (MQ=255)
                   gttgtCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCg                                                        >  2:436289/1‑80 (MQ=255)
                        cGCCAGTAATCTGCTGAACGTGCTCGGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCt                                                  <  2:436258/81‑1 (MQ=255)
                         gCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCt                                                  <  1:436259/80‑1 (MQ=255)
                             gTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGcac                                              >  1:436298/1‑80 (MQ=255)
                                   ggCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGac                                        <  1:436261/80‑1 (MQ=255)
                                   ggCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGac                                        <  2:436260/80‑1 (MQ=255)
                                       gAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCtt                                   >  1:436294/1‑81 (MQ=255)
                                       gAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCtt                                   >  1:436301/1‑81 (MQ=255)
                                        aaCGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCttt                                  >  1:436307/1‑81 (MQ=255)
                                           gTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCAc                               >  1:436302/1‑81 (MQ=255)
                                           gTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCAc                               >  2:436297/1‑81 (MQ=255)
                                             gCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGc                             >  1:436292/1‑81 (MQ=255)
                                                ctGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCg                           >  2:436295/1‑80 (MQ=255)
                                                 tGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCg                           >  1:436303/1‑79 (MQ=255)
                                                 tGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCg                           <  2:436262/79‑1 (MQ=255)
                                                   cgccgcCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGa                          <  1:436263/78‑1 (MQ=255)
                                                       gccTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTCCACGCCGAAGTgcgc                   <  2:436264/81‑1 (MQ=255)
                                                         cTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTgcgcg                  >  1:436311/1‑80 (MQ=255)
                                                               cGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTg            <  1:436265/80‑1 (MQ=255)
                                                                ggTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTgcg          >  2:436296/1‑81 (MQ=255)
                                                                 gTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTgcg          >  2:436308/1‑80 (MQ=255)
                                                                 gTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTgcg          >  2:436309/1‑80 (MQ=255)
                                                                   aCGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGcat       >  2:436304/1‑81 (MQ=255)
                                                                      gATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGcatcat    <  2:436267/81‑1 (MQ=255)
                                                                       aTAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATg   <  2:436268/81‑1 (MQ=255)
                                                                        tAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGc  >  1:436310/1‑81 (MQ=255)
                                                                              |                                                                          
CATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGC  >  NC_000913/2717739‑2717891

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