Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,417,157 C→T R77H (CGC→CAC)  recE ← Rac prophage, exonuclease VIII, 5' to 3' specific dsDNA exonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,417,1570CT100.0% 108.2 / NA 35R77H (CGC→CAC) recERac prophage, exonuclease VIII, 5' to 3' specific dsDNA exonuclease
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/19);  total (16/19)

ATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTC  >  NC_000913/1417022‑1417278
                                                                                                                                       |                                                                                                                         
aTTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGtt                                                                                                                        <  1:574750/139‑1 (MQ=255)
   gggAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCAc                                                                                                                     <  1:406871/139‑1 (MQ=255)
       gTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCa                                                                                                                 >  1:162635/1‑139 (MQ=255)
                                  tGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCggg                                                                                      <  2:162635/139‑1 (MQ=255)
                                   gTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGc                                                                                     >  1:496986/1‑139 (MQ=255)
                                    tCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCa                                                                                    <  2:496986/139‑1 (MQ=255)
                                             ttGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAAt                                                                                                         <  1:802206/109‑1 (MQ=255)
                                             ttGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAAt                                                                                                         >  2:802206/1‑109 (MQ=255)
                                              tGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAAc                                                                          >  1:276174/1‑139 (MQ=255)
                                                             cTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGt                                                           >  1:321812/1‑139 (MQ=255)
                                                                gATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGtt                                                                                                                        >  2:721267/1‑75 (MQ=255)
                                                                gATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGtt                                                                                                                        <  1:721267/75‑1 (MQ=255)
                                                                   aGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAAc                                                     <  2:371305/139‑1 (MQ=255)
                                                                   aGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAAc                                                     <  2:470925/139‑1 (MQ=255)
                                                                     tGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGtt                                                                                                                        >  1:348588/1‑70 (MQ=255)
                                                                     tGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGtt                                                                                                                        <  2:348588/70‑1 (MQ=255)
                                                                     tGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGtt                                                                                                                        >  1:31818/1‑70 (MQ=255)
                                                                     tGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGtt                                                                                                                        <  2:31818/70‑1 (MQ=255)
                                                                      gAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGtt                                                                                                                        <  1:257194/69‑1 (MQ=255)
                                                                      gAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGtt                                                                                                                        >  2:257194/1‑69 (MQ=255)
                                                                      gAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCggg                                                                                      >  2:455130/1‑103 (MQ=255)
                                                                      gAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCggg                                                                                      <  1:455130/103‑1 (MQ=255)
                                                                       aGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTGCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCt                                                 >  1:522950/1‑139 (MQ=255)
                                                                                      tgctCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTc                                  <  2:473431/139‑1 (MQ=255)
                                                                                             ggTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGtt                                                                                                                        <  2:581309/46‑1 (MQ=255)
                                                                                             ggTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGtt                                                                                                                        >  1:581309/1‑46 (MQ=255)
                                                                                             ggTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGttttt                           >  2:279999/1‑139 (MQ=255)
                                                                                                 ttAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTAcc                       <  1:279999/139‑1 (MQ=255)
                                                                                                               aTGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGAc                                                                       <  2:386357/77‑1 (MQ=255)
                                                                                                               aTGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGAc                                                                       >  1:386357/1‑77 (MQ=255)
                                                                                                                      cTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTc  <  2:321812/139‑1 (MQ=255)
                                                                                                                         tGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGa            <  2:199983/126‑1 (MQ=255)
                                                                                                                         tGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGa            >  1:199983/1‑126 (MQ=255)
                                                                                                                                 tGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTAc                        >  2:567189/1‑106 (MQ=255)
                                                                                                                                 tGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTAc                        <  1:567189/106‑1 (MQ=255)
                                                                                                                                       |                                                                                                                         
ATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTC  >  NC_000913/1417022‑1417278

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑