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% 67.7 / NA 24G11R (GGA→AGA) ymfDe14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (10/14);  total (10/14)

CTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTT  >  NC_000913/1197437‑1197566
                                                                 |                                                                
cTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                                <  2:391892/68‑1 (MQ=255)
cTCGCAAAGATATCAAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                                <  2:62208/68‑1 (MQ=255)
  cGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                                              <  1:197352/68‑1 (MQ=255)
  cGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                                              <  2:299008/68‑1 (MQ=255)
    cAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCat                                                            <  2:314005/68‑1 (MQ=255)
       aGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAtatt                                                         <  2:163120/68‑1 (MQ=255)
            tCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGt                                                    <  1:398219/68‑1 (MQ=255)
               agagGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                                >  1:485279/1‑53 (MQ=255)
               agagGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                                <  2:485279/53‑1 (MQ=255)
               agagGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAtatt                                                         >  2:539844/1‑60 (MQ=255)
               agagGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAtatt                                                         <  1:539844/60‑1 (MQ=255)
               agagGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCa                                                 >  1:626842/1‑68 (MQ=255)
                 agGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                                              >  2:32982/1‑53 (MQ=255)
                 agGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                                              <  1:32982/53‑1 (MQ=255)
                  gggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAg                                              >  2:613439/1‑68 (MQ=255)
                  gggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAg                                              >  2:614073/1‑68 (MQ=255)
                  gggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAg                                              >  2:85722/1‑68 (MQ=255)
                       tCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGt                                         <  1:85722/68‑1 (MQ=255)
                       tCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGt                                         <  2:464199/68‑1 (MQ=255)
                            tGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCAt                                    <  1:133325/68‑1 (MQ=255)
                                gTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCaa                                >  2:77167/1‑68 (MQ=255)
                                 ttttGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                                              >  1:479166/1‑37 (MQ=255)
                                 ttttGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                                              <  2:479166/37‑1 (MQ=255)
                                                        tatTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAAt        >  1:489673/1‑68 (MQ=255)
                                                              ctctGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCtttt  >  1:567381/1‑68 (MQ=255)
                                                                 |                                                                
CTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTT  >  NC_000913/1197437‑1197566

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