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% 150.2 / NA 39A238T (GCG→ACG) yfiFputative methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/18);  total (21/18)

CATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGC  >  NC_000913/2717739‑2717891
                                                                              |                                                                          
catcCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGtcc                                                                          >  2:452573/1‑81 (MQ=255)
    cAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGa                                                                      <  2:452526/81‑1 (MQ=255)
     aGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGAc                                                                     <  1:452527/81‑1 (MQ=255)
        aCGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTcc                                                                  <  1:452528/81‑1 (MQ=255)
        aCGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTc                                                                   >  1:452594/1‑80 (MQ=255)
           ttAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCcagc                                                               <  2:452529/81‑1 (MQ=255)
             aaCAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCcagcag                                                             >  1:452564/1‑81 (MQ=255)
               cAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTg                                                           >  2:452591/1‑81 (MQ=255)
                aaTGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGc                                                          >  1:452580/1‑81 (MQ=255)
                 aTGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGcc                                                         <  1:452530/81‑1 (MQ=255)
                  tgttgtCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCg                                                        >  1:452589/1‑81 (MQ=255)
                     tgtCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCAt                                                     >  2:452561/1‑81 (MQ=255)
                       tCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATcc                                                   <  1:789843/81‑1 (MQ=255)
                         gCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCt                                                  <  1:452531/80‑1 (MQ=255)
                         gCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCt                                                  <  1:452533/80‑1 (MQ=255)
                         gCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTg                                                 <  2:452532/81‑1 (MQ=255)
                          ccAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTg                                                 >  1:452581/1‑80 (MQ=255)
                           cAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGca                                               >  2:452614/1‑81 (MQ=255)
                               aaTCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGcaccac                                           >  2:452578/1‑81 (MQ=255)
                                aTCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGcaccac                                           >  1:452587/1‑80 (MQ=255)
                                aTCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGcaccac                                           >  1:452535/1‑80 (MQ=255)
                                aTCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACg                                          <  2:452534/81‑1 (MQ=255)
                                   ggCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGaca                                       >  1:452579/1‑81 (MQ=255)
                                    gCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGacac                                      >  1:452577/1‑81 (MQ=255)
                                        aaCGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCt                                    <  1:452536/79‑1 (MQ=255)
                                           gTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCAc                               >  1:452583/1‑81 (MQ=255)
                                            tGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACg                              >  1:452595/1‑81 (MQ=255)
                                                       gccTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTgcgc                   <  2:452537/81‑1 (MQ=255)
                                                       gccTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTgcgc                   >  1:452586/1‑81 (MQ=255)
                                                        ccTTCTGCGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTgcgcg                  <  2:452538/81‑1 (MQ=255)
                                                               cGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCt             <  1:452540/79‑1 (MQ=255)
                                                               cGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCt             <  1:452541/79‑1 (MQ=255)
                                                               cGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTgc           >  2:452584/1‑81 (MQ=255)
                                                               cGGTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTg            <  1:452539/80‑1 (MQ=255)
                                                                 gTACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTgcg          >  1:452592/1‑80 (MQ=255)
                                                                  tACGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGca        <  2:452542/81‑1 (MQ=255)
                                                                    cGGATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGcat       <  1:452543/80‑1 (MQ=255)
                                                                     ggATAGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGcatca     >  2:452611/1‑81 (MQ=255)
                                                                         aGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGc  >  1:452582/1‑80 (MQ=255)
                                                                         aGCCGTCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGc  >  2:452596/1‑80 (MQ=255)
                                                                              |                                                                          
CATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGC  >  NC_000913/2717739‑2717891

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