Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,717,821 (C)5→4 coding (708/1038 nt) yfiF ← putative methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,717,8170C.100.0% 147.2 / NA 36coding (712/1038 nt)yfiFputative methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (17/19);  total (17/19)

CCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCACCCAGGTTATGCGGGTTAGATTCGTTTTCCAGTGCCAGAACGCAATCCTGCGCG  >  NC_000913/2717690‑2717947
                                                                                                                               |                                                                                                                                  
cccTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCG‑CCCCGACTCcag                                                                                                                        <  1:693643/139‑1 (MQ=255)
   tGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCG‑CCCCGACTCcagcag                                                                                                                     >  1:21442/1‑139 (MQ=255)
     cTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCG‑CCCCGACTCCAGCAGTg                                                                                                                   <  2:478976/139‑1 (MQ=255)
       cGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCG‑CCCCGACTCCAGCAGTGcc                                                                                                                 <  1:930370/139‑1 (MQ=255)
          tGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCG‑CCCCGACTCCAGCAGTGCCGCa                                                                                                              <  2:21442/139‑1 (MQ=255)
          tGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCG‑CCCCGACTCCAGCAGTGCCGCa                                                                                                              >  2:326502/1‑139 (MQ=255)
                       tACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAgccgcc                                                                                                                                   >  2:655273/1‑106 (MQ=255)
                       tACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAgccgcc                                                                                                                                   <  1:655273/106‑1 (MQ=255)
                                    gcctgACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCG‑CCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGcc                                                                                    >  2:1181029/1‑139 (MQ=255)
                                                        cACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCG‑CCCCGACTCCAGCAGTGCCGCATCCTg                                                                                                         >  1:219814/1‑98 (MQ=255)
                                                        cACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCG‑CCCCGACTCCAGCAGTGCCGCATCCTg                                                                                                         <  2:219814/98‑1 (MQ=255)
                                                        cACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCG‑CCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGca                                                                <  2:893258/139‑1 (MQ=255)
                                                                 aaTGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCG‑CCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCAc                                                       <  1:318590/139‑1 (MQ=255)
                                                                    gttgcCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCG‑CCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCACCCa                                                    >  2:857665/1‑139 (MQ=255)
                                                                     ttgtCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCG‑CCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGcatcat                                                            >  2:144697/1‑130 (MQ=255)
                                                                     ttgtCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCG‑CCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGcatcat                                                            <  1:144697/130‑1 (MQ=255)
                                                                     ttgtCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCG‑CCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATg                                                           >  2:189772/1‑131 (MQ=255)
                                                                     ttgtCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCG‑CCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATg                                                           <  1:189772/131‑1 (MQ=255)
                                                                       gtCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCG‑CCCCGACTCCAGCAGTGCCGCATCCTg                                                                                                         >  1:446996/1‑83 (MQ=255)
                                                                       gtCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCG‑CCCCGACTCCAGCAGTGCCGCATCCTg                                                                                                         <  2:446996/83‑1 (MQ=255)
                                                                           ccAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCG‑CCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCACCCAGGTTATg                                             <  1:926689/139‑1 (MQ=255)
                                                                                                tctGCGCCGCCTTCTGCGGTACGGATAGCCG‑CCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCACCCAGGTTATGCGGGTTAg                                     >  2:899387/1‑126 (MQ=255)
                                                                                                tctGCGCCGCCTTCTGCGGTACGGATAGCCG‑CCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCACCCAGGTTATGCGGGTTAg                                     <  1:899387/126‑1 (MQ=255)
                                                                                                   gcgcCGCCTTCTGCGGTACGGATAGCCG‑CCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCACCCAGGTTATGCGGGTTAGATTCGTTTTCCAGTGc                     >  1:888960/1‑139 (MQ=255)
                                                                                                   gcgcCGCCTTCTGCGGTACGGATAGCCG‑CCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCACCCAGGTTATGCGGGTTAGATTCGTTTTCCAGTGc                     >  2:64695/1‑139 (MQ=255)
                                                                                                    cgccgcCTTCTGCGGTACGGATAGCCG‑CCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCtt                                                                                           >  1:787025/1‑68 (MQ=255)
                                                                                                    cgccgcCTTCTGCGGTACGGATAGCCG‑CCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCtt                                                                                           <  2:787025/68‑1 (MQ=255)
                                                                                                    cgccgcCTTCTGCGGTACGGATAGCCG‑CCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCACCCAGGTTATGCGGGTTAGATTCGTTTTCCAGTGcc                    >  1:430293/1‑139 (MQ=255)
                                                                                                     gccgccTTCTGCGGTACGGATAGCCG‑CCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCACCCAGGTTATGCGGGTTAGATTCGTTTTCCAGTGCCa                   <  1:1181029/139‑1 (MQ=255)
                                                                                                         ccTTCTGCGGTACGGATAGCCG‑CCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTg                                                                    <  2:139604/86‑1 (MQ=255)
                                                                                                         ccTTCTGCGGTACGGATAGCCG‑CCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTg                                                                    >  1:139604/1‑86 (MQ=255)
                                                                                                         ccTTCTGCGGTACGGATAGCCG‑CCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCACCCAGGTTATGCGGGTTAg                                     >  2:393641/1‑117 (MQ=255)
                                                                                                         ccTTCTGCGGTACGGATAGCCG‑CCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCACCCAGGTTATGCGGGTTAg                                     <  1:393641/117‑1 (MQ=255)
                                                                                                          cTTCTGCGGTACGGATAGCCG‑CCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCACCCAGGTTATGCGGGTTAGATTCGTTTTCCAGTGCCAGAACg              <  1:326502/139‑1 (MQ=255)
                                                                                                                   tACGGATAGCCG‑CCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCACCCAGGTTATGCGGGTTAGATTCGTTTTc                           <  1:775618/117‑1 (MQ=255)
                                                                                                                   tACGGATAGCCG‑CCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCACCCAGGTTATGCGGGTTAGATTCGTTTTc                           >  2:775618/1‑117 (MQ=255)
                                                                                                                    aCGGATAGCCG‑CCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCACCCAGGTTATGCGGGTTAGATTCGTTTTCCAGTGCCAGAACGCAATCCTgcg    <  1:354837/139‑1 (MQ=255)
                                                                                                                      ggATAGCCG‑CCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCACCCAGGTTATGCGGGTTAGATTCGTTTTCCAGTGCCAGAACGCAATCCTgcgcg  >  2:72076/1‑139 (MQ=255)
                                                                                                                               |                                                                                                                                  
CCCTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCACCCAGGTTATGCGGGTTAGATTCGTTTTCCAGTGCCAGAACGCAATCCTGCGCG  >  NC_000913/2717690‑2717947

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