Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 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_0009131,197,5020CT100.0% 63.5 / NA 22G11R (GGA→AGA) ymfDe14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (6/16);  total (6/16)

CTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTT  >  NC_000913/1197437‑1197565
                                                                 |                                                               
ctctCAAAGATATCTAGTGGTCATCGTATGAGTTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                               <  2:382790/68‑1 (MQ=255)
cTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                               <  2:628596/68‑1 (MQ=255)
cTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                               <  2:60830/68‑1 (MQ=255)
cTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                               <  1:598626/68‑1 (MQ=255)
cTCGCAAAGATATCTAGAGGGCATCGTACGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                               <  2:450301/68‑1 (MQ=255)
  cGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                                             <  2:280406/68‑1 (MQ=255)
  cGCAAAGATATCTAGAGGGCATCGTACGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                                             <  2:489443/68‑1 (MQ=255)
    cAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCat                                                           <  2:174608/68‑1 (MQ=255)
       aGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAtatt                                                        <  1:64982/68‑1 (MQ=255)
               agagGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCa                                                >  2:140018/1‑68 (MQ=255)
               agagGGCATAGGATGAGGTTTTGCAGCATTCTCTGAGCGAATTTTAACTCTGTGCATATTATAGTTCa                                                >  2:507410/1‑68 (MQ=255)
                  gggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAg                                             >  2:633679/1‑68 (MQ=255)
                  gggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAg                                             <  1:134624/68‑1 (MQ=255)
                  gggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAg                                             <  2:411915/68‑1 (MQ=255)
                  gggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAg                                             <  2:208335/68‑1 (MQ=255)
                         gTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTAc                                      >  1:171840/1‑68 (MQ=255)
                          tATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTAcc                                     <  2:138558/68‑1 (MQ=255)
                          tATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTAcc                                     <  2:584400/68‑1 (MQ=255)
                                gTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCat                                                           <  2:383542/40‑1 (MQ=255)
                                gTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCat                                                           >  1:383544/1‑40 (MQ=255)
                                                   gCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGa            <  1:626524/68‑1 (MQ=255)
                                                             aCTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCttt  >  1:181073/1‑68 (MQ=255)
                                                                 |                                                               
CTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTT  >  NC_000913/1197437‑1197565

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