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
                                                                                                                               |                                                                                                                                  
cccTGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCGACTCcag                                                                                                                        <  1:693643/139‑1 (MQ=255)
   tGCTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCGACTCcagcag                                                                                                                     >  1:21442/1‑139 (MQ=255)
     cTCGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCGACTCCAGCAGTg                                                                                                                   <  2:478976/139‑1 (MQ=255)
       cGCTGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCGACTCCAGCAGTGcc                                                                                                                 <  1:930370/139‑1 (MQ=255)
          tGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCGACTCCAGCAGTGCCGCa                                                                                                              <  2:21442/139‑1 (MQ=255)
          tGGAAGTGGTCACTACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCGACTCCAGCAGTGCCGCa                                                                                                              >  2:326502/1‑139 (MQ=255)
                       tACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAgccgcc                                                                                                                                   >  2:655273/1‑106 (MQ=255)
                       tACGGTGTAACCTGCCTGACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAgccgcc                                                                                                                                   <  1:655273/106‑1 (MQ=255)
                                    gcctgACGGAAATCATCCAGCACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGcc                                                                                    >  2:1181029/1‑139 (MQ=255)
                                                        cACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCGACTCCAGCAGTGCCGCATCCTg                                                                                                         >  1:219814/1‑98 (MQ=255)
                                                        cACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCGACTCCAGCAGTGCCGCATCCTg                                                                                                         <  2:219814/98‑1 (MQ=255)
                                                        cACGTTAACAATGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGca                                                                <  2:893258/139‑1 (MQ=255)
                                                                 aaTGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCAc                                                       <  1:318590/139‑1 (MQ=255)
                                                                    gttgcCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCACCCa                                                    >  2:857665/1‑139 (MQ=255)
                                                                     ttgtCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGcatcat                                                            >  2:144697/1‑130 (MQ=255)
                                                                     ttgtCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGcatcat                                                            <  1:144697/130‑1 (MQ=255)
                                                                     ttgtCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATg                                                           >  2:189772/1‑131 (MQ=255)
                                                                     ttgtCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATg                                                           <  1:189772/131‑1 (MQ=255)
                                                                       gtCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCGACTCCAGCAGTGCCGCATCCTg                                                                                                         >  1:446996/1‑83 (MQ=255)
                                                                       gtCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCGACTCCAGCAGTGCCGCATCCTg                                                                                                         <  2:446996/83‑1 (MQ=255)
                                                                           ccAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCACCCAGGTTATg                                             <  1:926689/139‑1 (MQ=255)
                                                                                                tctGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCACCCAGGTTATGCGGGTTAg                                     >  2:899387/1‑126 (MQ=255)
                                                                                                tctGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCACCCAGGTTATGCGGGTTAg                                     <  1:899387/126‑1 (MQ=255)
                                                                                                   gcgcCGCCTTCTGCGGTACGGATAGCCGCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCACCCAGGTTATGCGGGTTAGATTCGTTTTCCAGTGc                     >  1:888960/1‑139 (MQ=255)
                                                                                                   gcgcCGCCTTCTGCGGTACGGATAGCCGCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCACCCAGGTTATGCGGGTTAGATTCGTTTTCCAGTGc                     >  2:64695/1‑139 (MQ=255)
                                                                                                    cgccgcCTTCTGCGGTACGGATAGCCGCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCtt                                                                                           >  1:787025/1‑68 (MQ=255)
                                                                                                    cgccgcCTTCTGCGGTACGGATAGCCGCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCtt                                                                                           <  2:787025/68‑1 (MQ=255)
                                                                                                    cgccgcCTTCTGCGGTACGGATAGCCGCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCACCCAGGTTATGCGGGTTAGATTCGTTTTCCAGTGcc                    >  1:430293/1‑139 (MQ=255)
                                                                                                     gccgccTTCTGCGGTACGGATAGCCGCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCACCCAGGTTATGCGGGTTAGATTCGTTTTCCAGTGCCa                   <  1:1181029/139‑1 (MQ=255)
                                                                                                         ccTTCTGCGGTACGGATAGCCGCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTg                                                                    <  2:139604/86‑1 (MQ=255)
                                                                                                         ccTTCTGCGGTACGGATAGCCGCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTg                                                                    >  1:139604/1‑86 (MQ=255)
                                                                                                         ccTTCTGCGGTACGGATAGCCGCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCACCCAGGTTATGCGGGTTAg                                     >  2:393641/1‑117 (MQ=255)
                                                                                                         ccTTCTGCGGTACGGATAGCCGCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCACCCAGGTTATGCGGGTTAg                                     <  1:393641/117‑1 (MQ=255)
                                                                                                          cTTCTGCGGTACGGATAGCCGCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCACCCAGGTTATGCGGGTTAGATTCGTTTTCCAGTGCCAGAACg              <  1:326502/139‑1 (MQ=255)
                                                                                                                   tACGGATAGCCGCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCACCCAGGTTATGCGGGTTAGATTCGTTTTc                           <  1:775618/117‑1 (MQ=255)
                                                                                                                   tACGGATAGCCGCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCACCCAGGTTATGCGGGTTAGATTCGTTTTc                           >  2:775618/1‑117 (MQ=255)
                                                                                                                    aCGGATAGCCGCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCACCCAGGTTATGCGGGTTAGATTCGTTTTCCAGTGCCAGAACGCAATCCTgcg    <  1:354837/139‑1 (MQ=255)
                                                                                                                      ggATAGCCGCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCATCATGCCACCCAGGTTATGCGGGTTAGATTCGTTTTCCAGTGCCAGAACGCAATCCTgcgcg  >  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: