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

CTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTT  >  NC_000913/1197437‑1197564
                                                                 |                                                              
cTCGCAAAGATATCTAGAGGGCATCGTATGAGTTTTTGCTTCATTCTCTGAGCGAATATTAACTCTGt                                                              <  1:523932/68‑1 (MQ=255)
cTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                              <  1:1099863/68‑1 (MQ=255)
cTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                              <  2:1168187/68‑1 (MQ=255)
cTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                              <  2:64796/68‑1 (MQ=255)
 tCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCtgtg                                                             >  1:1010780/1‑68 (MQ=255)
  cGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                                            <  1:734089/68‑1 (MQ=255)
  cGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                                            <  1:1246767/68‑1 (MQ=255)
  cGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                                            <  1:350219/68‑1 (MQ=255)
  cGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                                            <  2:276713/68‑1 (MQ=255)
  cGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                                            <  2:396552/68‑1 (MQ=255)
    cAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCat                                                          <  1:870945/68‑1 (MQ=255)
    cAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCat                                                          <  1:955877/68‑1 (MQ=255)
       aGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAtatt                                                       <  2:584771/68‑1 (MQ=255)
       aGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAtatt                                                       >  2:136022/1‑68 (MQ=255)
       aGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAtatt                                                       <  1:1092086/68‑1 (MQ=255)
       aGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAtatt                                                       <  2:1263600/68‑1 (MQ=255)
               agagGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                              <  1:166299/53‑1 (MQ=255)
               agagGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                              >  2:1081521/1‑53 (MQ=255)
               agagGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                              >  2:166297/1‑53 (MQ=255)
               agagGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                              <  1:1081525/53‑1 (MQ=255)
               agagGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCa                                               >  2:1012522/1‑68 (MQ=255)
                 agGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                              <  2:843993/51‑1 (MQ=255)
                 agGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                              >  1:843998/1‑51 (MQ=255)
                 agGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAtatt                                                       <  1:1406761/58‑1 (MQ=255)
                 agGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAtatt                                                       >  2:1406756/1‑58 (MQ=255)
                 agGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCaaa                                             >  2:914951/1‑68 (MQ=255)
                  gggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                                            <  1:1336984/52‑1 (MQ=255)
                  gggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                                            >  2:1336980/1‑52 (MQ=255)
                  gggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAg                                            <  1:301234/68‑1 (MQ=255)
                  gggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAg                                            >  2:1118648/1‑68 (MQ=255)
                  gggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAg                                            <  1:1372011/68‑1 (MQ=255)
                  gggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAg                                            >  2:1281524/1‑68 (MQ=255)
                   ggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGc                                           <  1:761111/68‑1 (MQ=255)
                    gCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCa                                          >  2:1395413/1‑68 (MQ=255)
                    gCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCa                                          >  1:785170/1‑68 (MQ=255)
                       tCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGt                                       <  1:309688/68‑1 (MQ=255)
                               ggTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCa                               >  2:1150453/1‑68 (MQ=255)
                               ggTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCa                               >  1:532623/1‑68 (MQ=255)
                               ggTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCa                               <  1:1408849/68‑1 (MQ=255)
                               ggTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCa                               <  1:1395418/68‑1 (MQ=255)
                               ggTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCa                               <  1:1012527/68‑1 (MQ=255)
                                     gctgcATTCTCTGAGCGAATATTTACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAAt                         <  1:178144/68‑1 (MQ=255)
                                     gctgcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAAt                         <  1:265415/68‑1 (MQ=255)
                                     gctgcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAAt                         <  2:620641/68‑1 (MQ=255)
                                       tgcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGt                       <  1:248260/68‑1 (MQ=255)
                                                            aaCTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCtt  >  1:398698/1‑68 (MQ=255)
                                                                 |                                                              
CTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTT  >  NC_000913/1197437‑1197564

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: