Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 1,197,502 C→T G11R (GGA→AGA)  ymfD ← e14 prophage; putative SAM‑dependent methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA1,197,5020CT100.0% 101.3 / NA 32G11R (GGA→AGA) ymfDe14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/16);  total (16/16)

ACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAG  >  NC_000913_3_pae_tpiA/1197356‑1197647
                                                                                                                                                  |                                                                                                                                                 
aCTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                                                                                                                 <  1:149098/149‑1 (MQ=255)
aCTGATTAATTCATCAGAATATCTAAGTTTTCCGAAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTTAGCGCATATTAACTCTGt                                                                                                                                                 <  2:206031/149‑1 (MQ=255)
    aTTAATTCATCAGAATATCTAAGTTTTCCGCGGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCat                                                                                                                                             >  1:562976/1‑149 (MQ=255)
     ttAATTCATCAGAATATCTAAGTTTTCCGCCGCCAAACTCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCata                                                                                                                                            <  1:386502/149‑1 (MQ=255)
        aTTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAtatta                                                                                                                                         >  1:563508/1‑149 (MQ=255)
         ttCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACCCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAtattat                                                                                                                                        >  1:751409/1‑149 (MQ=255)
          tcatcaGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACCCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATtata                                                                                                                                       <  2:751409/149‑1 (MQ=255)
            atcaGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGt                                                                                                                                     >  1:284116/1‑149 (MQ=255)
             tcaGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGtt                                                                                                                                    <  2:563508/149‑1 (MQ=255)
                      cTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAg                                                                                                                           >  2:21737/1‑149 (MQ=255)
                         aGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTAc                                                                                                                        <  2:284116/149‑1 (MQ=255)
                                cGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTACCAATGTTCATAGTATAGTTCAAAGCAAGTACCATTCaa                                                                                                                 <  2:436763/149‑1 (MQ=255)
                                 gCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAAt                                                                                                                <  2:210350/149‑1 (MQ=255)
                                    gCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCATAGAGATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTaa                                                                                                             >  1:417355/1‑149 (MQ=255)
                                     ccAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAAt                                                                                                            >  1:737440/1‑149 (MQ=255)
                                                         aTTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCGAGAGGGCATCGTATAAGGTTTTGCTGCATTCTCTGCGCGAATATTAACTATGTGCATATTATAGTTCAAAGCAAATACCATTCAATTAATGTTATTATTAGTGAAAAAtt                                                                                        >  2:57457/1‑149 (MQ=255)
                                                                tcAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCTTCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAATACCATACAATTAATGTTATTTGTAGTGAAAAATTCttttat                                                                                 >  2:11294/1‑149 (MQ=255)
                                                                                    gCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTc                                                             >  2:109787/1‑149 (MQ=255)
                                                                                                agagGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACttt                                                 <  1:297312/149‑1 (MQ=255)
                                                                                                    ggCATCGTCTGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAATCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGAGATTTTTCCGACGTGACTTGCAGt                                             >  2:874138/1‑149 (MQ=255)
                                                                                                     gCATCGGATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGtt                                            <  1:674514/149‑1 (MQ=255)
                                                                                                          gTAAGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTggg                                       <  2:402926/149‑1 (MQ=255)
                                                                                                                 gTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAg                                                                                                 >  2:191186/1‑84 (MQ=255)
                                                                                                                 gTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAg                                                                                                 <  1:191186/84‑1 (MQ=255)
                                                                                                                 gTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCttttatttta                                                                             >  1:699001/1‑104 (MQ=255)
                                                                                                                 gTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCttttatttta                                                                             <  2:699001/104‑1 (MQ=255)
                                                                                                                 gTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTAAAAGCAAGTACCATTCAATTAATGTTATTTTTAATGAATAATTCTTGTATTTTATCGTCCGGGATGTTGCGAACTCGACTTTCATGTCTCGGCTTAACt                                >  2:316305/1‑149 (MQ=255)
                                                                                                                      gctgcATTCCCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaa                           >  2:763702/1‑149 (MQ=255)
                                                                                                                                 tGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCttttatt                                                                                >  2:50313/1‑85 (MQ=255)
                                                                                                                                 tGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCttttatt                                                                                <  1:50313/85‑1 (MQ=255)
                                                                                                                                      gAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTCACTGTCAAAGCAAAGATAATACca           <  2:274518/149‑1 (MQ=255)
                                                                                                                                      gAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACca           <  2:43061/149‑1 (MQ=255)
                                                                                                                                               ctctGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGTATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTgat  >  2:105104/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ACTGATTAATTCATCAGAATATCTAAGTTTTCCGCAGCCAAAATCAAGGGCGTGGCCATTTTTCTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAG  >  NC_000913_3_pae_tpiA/1197356‑1197647

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑