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% 128.1 / NA 40R77H (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 (19/21);  total (19/21)

GGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCA  >  NC_000913/1417026‑1417280
                                                                                                                                   |                                                                                                                           
ggAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACt                                                                                                                      <  1:183476/139‑1 (MQ=255)
ggAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACt                                                                                                                      <  2:584057/139‑1 (MQ=255)
                 tcAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGAtct                                                                                                     >  1:929428/1‑139 (MQ=255)
                    aTCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTCGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCAc                                                                                                  <  1:996525/139‑1 (MQ=255)
                        cAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTc                                                                                              >  2:283255/1‑139 (MQ=255)
                               gTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGc                                                                                       <  1:104150/139‑1 (MQ=255)
                               gTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGc                                                                                       >  2:104150/1‑139 (MQ=255)
                                tCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCa                                                                                      <  1:326295/139‑1 (MQ=255)
                                  ttCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAgg                                                                                    <  1:219320/139‑1 (MQ=255)
                                                      tGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCAc                                                                                                  >  1:443872/1‑105 (MQ=255)
                                                      tGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCAc                                                                                                  <  2:443872/105‑1 (MQ=255)
                                                      tGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCAc                                                                                                  >  2:673932/1‑105 (MQ=255)
                                                      tGCCTTGATAGTGAGCGGTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGCTAAAATCGATCTCAc                                                                                                  <  1:673932/105‑1 (MQ=255)
                                                       gCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTc                                                               >  2:988361/1‑139 (MQ=255)
                                                             aTAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGc                                                                                       <  1:813587/109‑1 (MQ=255)
                                                             aTAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGc                                                                                       >  2:813587/1‑109 (MQ=255)
                                                             aTAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCAtt                                                                               >  1:696871/1‑117 (MQ=255)
                                                             aTAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCAtt                                                                               <  2:696871/117‑1 (MQ=255)
                                                               aGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAAc                                                       <  1:988361/139‑1 (MQ=255)
                                                                 tGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAAc                                                       >  1:103315/1‑137 (MQ=255)
                                                                 tGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAAc                                                       <  2:103315/137‑1 (MQ=255)
                                                                  gAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGc                                                    >  1:222621/1‑139 (MQ=255)
                                                                  gAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGc                                                    >  1:771564/1‑139 (MQ=255)
                                                                   aGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCt                                                   >  1:1706/1‑139 (MQ=255)
                                                                       ttGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGc                                                                                       >  1:1002123/1‑99 (MQ=255)
                                                                       ttGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGc                                                                                       <  2:1002123/99‑1 (MQ=255)
                                                                         gTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTaaaa                                             <  2:929428/139‑1 (MQ=255)
                                                                                  tgctCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTc                                    <  2:1706/139‑1 (MQ=255)
                                                                                    ctCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTaaaa                                                                                                            <  2:462535/65‑1 (MQ=255)
                                                                                    ctCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTaaaa                                                                                                            >  1:462535/1‑65 (MQ=255)
                                                                                    ctCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAg                                  >  2:464644/1‑139 (MQ=255)
                                                                                    ctCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAg                                  >  1:985384/1‑139 (MQ=255)
                                                                                                     cATGTCATGGAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTaa                 <  2:222621/139‑1 (MQ=255)
                                                                                                             ggAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGa              >  1:907648/1‑134 (MQ=255)
                                                                                                             ggAGTCTTTGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGa              <  2:907648/134‑1 (MQ=255)
                                                                                                                    ttGCTGAGTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCa  <  1:283255/139‑1 (MQ=255)
                                                                                                                           gTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAgg                                                                                    <  1:793257/50‑1 (MQ=255)
                                                                                                                           gTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAgg                                                                                    >  2:793257/1‑50 (MQ=255)
                                                                                                                           gTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAAc                                                       >  1:941153/1‑79 (MQ=255)
                                                                                                                           gTTGATAGTGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAAc                                                       <  2:941153/79‑1 (MQ=255)
                                                                                                                                   |                                                                                                                           
GGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCA  >  NC_000913/1417026‑1417280

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