Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,718,339 A→G W64R (TGG→CGG)  yfiF ← putative methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,718,3390AG100.0% 118.2 / NA 33W64R (TGG→CGG) yfiFputative methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/14);  total (19/14)

TCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTG  >  NC_000913/2718191‑2718475
                                                                                                                                                    |                                                                                                                                        
tCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCGCt                                                                                                                                        >  2:76404/1‑151 (MQ=255)
     gTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCGCTTGCgg                                                                                                                                   >  1:602479/1‑151 (MQ=255)
      tAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCGCTTGCGGt                                                                                                                                  >  2:120691/1‑151 (MQ=255)
             gggTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCGCTTGCGGTCACgggg                                                                                                                           >  2:101356/1‑151 (MQ=255)
                           cGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCGCTTGCGGTCACGGGGCTGACTCTGTTTGt                                                                                                             <  1:211704/151‑1 (MQ=255)
                              tGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCGCTTGCGGTCACGGGGCTGACTCTGTTTGTCAt                                                                                                          >  2:236291/1‑151 (MQ=255)
                                                ggATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCGCTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGc                                                                                        >  2:575097/1‑151 (MQ=255)
                                                        aaaaCTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCGCTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTc                                                                                 >  1:557951/1‑150 (MQ=255)
                                                        aaaaCTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCGCTTGCGGTCACGGGGCTGACTCTGTTTGCCATCGCGGGCGGGGCGACGGCCACCGTCa                                                                                <  1:540639/151‑1 (MQ=255)
                                                          aaCTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCGCTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGc                                                                              >  2:233521/1‑151 (MQ=255)
                                                              tttACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCGCTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACg                                                                           >  2:571556/1‑150 (MQ=255)
                                                                 aCCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCGCTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACg                                                                           >  1:78895/1‑147 (MQ=255)
                                                                        aTGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCGCTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTg                                                                <  2:416852/151‑1 (MQ=255)
                                                                        aTGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCGCTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTg                                                                >  1:33943/1‑151 (MQ=255)
                                                                            cACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCGCTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGAcg                                                            >  1:173092/1‑151 (MQ=255)
                                                                                 tGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCGCTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGAcgcccgc                                                       >  1:397561/1‑151 (MQ=255)
                                                                                   aTCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCGCTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACgcccgccc                                                     <  2:571194/151‑1 (MQ=255)
                                                                                     cGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCGCTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGccccc                                                   <  1:220196/151‑1 (MQ=255)
                                                                                           tttCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCGCTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCAc                                                                            >  1:664336/1‑120 (MQ=255)
                                                                                           tttCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCGCTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCAc                                                                            <  2:664347/120‑1 (MQ=255)
                                                                                            ttCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCGCTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCg                                            >  1:172918/1‑151 (MQ=255)
                                                                                              ccGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCGCTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCgg                                           <  2:119836/150‑1 (MQ=255)
                                                                                              ccGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCGCTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGt                                          >  2:300876/1‑151 (MQ=255)
                                                                                               cGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCGCTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGt                                          >  2:155910/1‑150 (MQ=255)
                                                                                                        atcaCCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCGCTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTtg                                <  2:90910/151‑1 (MQ=255)
                                                                                                           aCCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCGCTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTgg                             <  2:459398/151‑1 (MQ=255)
                                                                                                                          aaCCGTGCGCCACGGCGAATCTTCCCGCTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGaa              <  1:269641/151‑1 (MQ=255)
                                                                                                                              gTGCGCCACGGCGAATCTTCCCGCTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAatcat          <  2:78895/151‑1 (MQ=255)
                                                                                                                               tGCGCCACGGCGAATCTTCCCGCTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAatcatc         >  1:442679/1‑151 (MQ=255)
                                                                                                                                gcgcCACGGCGAATCTTCCCGCTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAatcatc         <  2:619111/150‑1 (MQ=255)
                                                                                                                                  gcCACGGCGAATCTTCCCGCTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAatcatcatc      >  2:36537/1‑151 (MQ=255)
                                                                                                                                     aCGGCGAATCTTCCCGCTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACt   <  2:379188/151‑1 (MQ=255)
                                                                                                                                      cGGCGAATCTTCCCGCTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTg  <  2:347780/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                        
TCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTG  >  NC_000913/2718191‑2718475

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