Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,417,290 A→T C33S (TGT→AGT)  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,2900AT97.1% 245.0 / ‑1.8 70C33S (TGT→AGT) recERac prophage; exonuclease VIII, 5' to 3' specific dsDNA exonuclease
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (40/28);  total (41/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

CAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCAT  >  NC_000913/1417195‑1417386
                                                                                               |                                                                                                
cAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGtn                                                                                              <  2:887793/100‑2 (MQ=255)
cAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTc                                                                                              <  2:1349586/100‑1 (MQ=255)
  ggTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGa                                                                                            <  1:951424/100‑1 (MQ=255)
  ggTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGa                                                                                            >  1:492412/1‑100 (MQ=255)
  ggTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGa                                                                                            >  1:1505191/1‑100 (MQ=255)
    tCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGAtt                                                                                          <  2:86851/100‑1 (MQ=255)
    tCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGAtt                                                                                          <  2:85633/100‑1 (MQ=255)
     cATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTc                                                                                         >  2:194587/1‑100 (MQ=255)
      aTTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCa                                                                                        >  1:509337/1‑100 (MQ=255)
       ttAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCaa                                                                                       <  2:292413/100‑1 (MQ=255)
           cGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATc                                                                                   >  2:656739/1‑100 (MQ=255)
           cGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATc                                                                                   <  1:846790/100‑1 (MQ=255)
           cGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATc                                                                                   >  1:1155590/1‑100 (MQ=255)
           cGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATc                                                                                   >  1:1233878/1‑100 (MQ=255)
            gACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCg                                                                                  <  2:629276/100‑1 (MQ=255)
              cAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGtt                                                                                >  1:86500/1‑100 (MQ=255)
               aGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGttg                                                                               <  2:1047190/100‑1 (MQ=255)
                ggAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGttgc                                                                              >  2:1319268/1‑100 (MQ=255)
                ggAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGttgc                                                                              >  2:87873/1‑100 (MQ=255)
                     aTTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGcc                                                                         <  2:1419680/100‑1 (MQ=255)
                       tCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCca                                                                       <  1:1120941/100‑1 (MQ=255)
                         gTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCcaca                                                                     <  1:1186570/100‑1 (MQ=255)
                            gCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGga                                                                  >  2:761743/1‑100 (MQ=255)
                               aCAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGgacga                                                               >  1:191050/1‑100 (MQ=255)
                                  gCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGt                                                            >  1:1294796/1‑100 (MQ=255)
                                  gCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGt                                                            >  1:627712/1‑100 (MQ=255)
                                   cTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTc                                                           >  1:664247/1‑100 (MQ=255)
                                     ttAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAg                                                         <  1:595102/100‑1 (MQ=255)
                                     ttAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAg                                                         >  2:154621/1‑100 (MQ=255)
                                     ttAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAg                                                         >  1:1224453/1‑100 (MQ=255)
                                     ttAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAg                                                         >  2:74551/1‑100 (MQ=255)
                                      tAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAgg                                                        <  2:81281/100‑1 (MQ=255)
                                      tAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAgg                                                        <  1:479945/100‑1 (MQ=255)
                                       aaaaTAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGt                                                       >  1:892179/1‑100 (MQ=255)
                                       aaaaTAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGt                                                       >  1:1515437/1‑100 (MQ=255)
                                        aaaTAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGtt                                                      <  2:1350947/100‑1 (MQ=255)
                                            aGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCAcc                                                  <  1:1293666/100‑1 (MQ=255)
                                               tgctCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGga                                               <  1:149187/100‑1 (MQ=255)
                                                gctCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGgat                                              >  1:109962/1‑100 (MQ=255)
                                                gctCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGgat                                              >  1:203251/1‑100 (MQ=255)
                                                gctCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGgat                                              >  1:399309/1‑100 (MQ=255)
                                                 ctCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGgatg                                             <  1:697751/100‑1 (MQ=255)
                                                  tCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGgatga                                            >  1:1016532/1‑100 (MQ=255)
                                                   cAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACCGGACGACGTCAGGTTCACCGgatgat                                           <  1:958483/100‑1 (MQ=255)
                                                   cAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGgatgat                                           >  1:371645/1‑100 (MQ=255)
                                                    aGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGAtt                                          >  1:884976/1‑100 (MQ=255)
                                                        ttttACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGAtttttt                                      <  2:461242/100‑1 (MQ=255)
                                                         tttACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTc                                     <  2:1515427/100‑1 (MQ=255)
                                                          ttACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCg                                    <  2:86500/100‑1 (MQ=255)
                                                          ttACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCg                                    <  1:1288885/100‑1 (MQ=255)
                                                          ttACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCg                                    <  2:492409/100‑1 (MQ=255)
                                                           tACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGc                                   >  1:977949/1‑100 (MQ=255)
                                                             ccTGACTTAACGATCAGGTAGTCCAGCGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCtt                                 >  2:1471837/1‑100 (MQ=255)
                                                                 aCTTAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCg                             <  2:389387/100‑1 (MQ=255)
                                                                    tAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTaa                          >  2:1183441/1‑100 (MQ=255)
                                                                    tAACGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTaa                          >  2:1085183/1‑100 (MQ=255)
                                                                      aCGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACa                        >  2:1216841/1‑100 (MQ=255)
                                                                       cGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAg                       >  2:1536803/1‑100 (MQ=255)
                                                                       cGATCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAg                       >  2:1111392/1‑100 (MQ=255)
                                                                         aTCAGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGa                     <  1:1462/100‑1 (MQ=255)
                                                                            aGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAaga                  >  1:120165/1‑100 (MQ=255)
                                                                            aGGTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAaga                  >  1:1192590/1‑100 (MQ=255)
                                                                             ggTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTGACCGGATGATTTTTTCGCTTTCCGTAACAGGAagag                 >  2:867383/1‑100 (MQ=255)
                                                                             ggTAGTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAagag                 >  2:259879/1‑100 (MQ=255)
                                                                                 gTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTgg              >  2:445390/1‑99 (MQ=255)
                                                                                 gTCCAGAGTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGt             >  1:376912/1‑100 (MQ=255)
                                                                                     agagTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTtg         <  2:1067697/100‑1 (MQ=255)
                                                                                      gagTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTtgt        >  1:973/1‑100 (MQ=255)
                                                                                       agTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTtgtg       <  1:87873/100‑1 (MQ=255)
                                                                                        gTGGCACTGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGc      >  1:189026/1‑100 (MQ=255)
                                                                                            cactGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCAt  <  1:672748/100‑1 (MQ=255)
                                                                                            cactGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTATTGTGCTCAt  >  2:1028675/1‑100 (MQ=255)
                                                                                               |                                                                                                
CAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCAT  >  NC_000913/1417195‑1417386

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