Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,413,709 T→C E11G (GAA→GGA)  ydaC ← DUF1187 family protein, Rac prophage; putative double‑strand break reduction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,413,7090TC95.7% 222.7 / 0.2 69E11G (GAA→GGA) ydaCDUF1187 family protein, Rac prophage; putative double‑strand break reduction protein
Reads supporting (aligned to +/- strand):  ref base T (3/0);  new base C (26/40);  total (29/40)
Fisher's exact test for biased strand distribution p-value = 6.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

CGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACT  >  NC_000913/1413614‑1413804
                                                                                               |                                                                                               
cGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCAcctccctt                                                                                             >  2:1729711/1‑100 (MQ=255)
    tGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCa                                                                                         >  1:1002275/1‑100 (MQ=255)
     gCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCaa                                                                                        >  1:881134/1‑100 (MQ=255)
      cttcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAAt                                                                                       <  1:1777112/100‑1 (MQ=255)
        tcttTTTTACTTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAg                                                                                     >  2:582762/1‑100 (MQ=255)
         cttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGt                                                                                    >  2:339450/1‑100 (MQ=255)
          ttttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTg                                                                                   <  2:1959926/100‑1 (MQ=255)
                aCCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTa                                                                             >  2:1777570/1‑100 (MQ=255)
                  cTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAAt                                                                           <  2:154576/100‑1 (MQ=255)
                   tGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAAtt                                                                          >  1:1926110/1‑100 (MQ=255)
                       agCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTtat                                                                      <  2:1848551/100‑1 (MQ=255)
                       agCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTtat                                                                      <  1:908066/100‑1 (MQ=255)
                        gCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTtata                                                                     >  2:549212/1‑100 (MQ=255)
                         cATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATAc                                                                    >  1:1688084/1‑100 (MQ=255)
                         cATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATAc                                                                    <  1:1707438/100‑1 (MQ=255)
                         cATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATAc                                                                    <  2:98588/100‑1 (MQ=255)
                          aTTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACa                                                                   <  1:1571620/100‑1 (MQ=255)
                          aTTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACa                                                                   >  1:1857790/1‑100 (MQ=255)
                             ttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGa                                                                <  2:581570/100‑1 (MQ=255)
                             ttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGa                                                                <  2:881131/100‑1 (MQ=255)
                               ttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAac                                                              <  1:101228/100‑1 (MQ=255)
                                tCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAaca                                                             <  2:1217253/100‑1 (MQ=255)
                                 cGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGa                                                                >  2:798879/1‑96 (MQ=255)
                                 cGCATTCTGATTTAGTTAGTTTAGATTTTGAATATCTTATCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAacac                                                            >  2:1259275/1‑100 (MQ=255)
                                  tcATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACACGAACAcc                                                           <  2:1784211/99‑1 (MQ=255)
                                  gCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACAcc                                                           <  1:1489493/100‑1 (MQ=255)
                                   cATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACACGAACAcca                                                          >  1:1784226/1‑99 (MQ=255)
                                   cATTCTGATTTCGTTAGTTTAGATTTTGAACATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCt                                                          <  1:123177/100‑1 (MQ=255)
                                    aTTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTc                                                         >  2:1281890/1‑100 (MQ=255)
                                     ttCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTcc                                                        >  1:456658/1‑100 (MQ=255)
                                      tCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCa                                                       <  2:1408249/100‑1 (MQ=255)
                                             tCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCattattat                                                <  1:1022908/100‑1 (MQ=255)
                                              cGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATATTATACATGAACACCTCCattattatt                                               >  2:1699528/1‑100 (MQ=255)
                                                   gTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAg                                          >  2:602284/1‑100 (MQ=255)
                                                    tttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGt                                         <  2:975858/100‑1 (MQ=255)
                                                        gATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGtt                                     >  2:1001519/1‑100 (MQ=255)
                                                        gATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGtt                                     >  2:794235/1‑100 (MQ=255)
                                                        gATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATAAATCACCACCTCCATTATTATTTcctgcgct                                      >  2:1237336/1‑93 (MQ=255)
                                                          ttttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCg                                   >  1:766698/1‑100 (MQ=255)
                                                           tttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGt                                  >  2:1629663/1‑100 (MQ=255)
                                                            ttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGtt                                 <  1:936052/100‑1 (MQ=255)
                                                            ttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGtt                                 <  1:1281896/100‑1 (MQ=255)
                                                                atatCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCTAGTGGTTCGTTTAtt                             >  1:2015927/1‑100 (MQ=255)
                                                                atatCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTAtt                             <  1:1206407/100‑1 (MQ=255)
                                                                  atCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTcc                           <  1:798882/100‑1 (MQ=255)
                                                                   tCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCa                          <  1:475991/100‑1 (MQ=255)
                                                                   tCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCa                          <  2:73898/100‑1 (MQ=255)
                                                                   tCTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCa                          <  1:1682831/100‑1 (MQ=255)
                                                                    cTTGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCa                          >  2:823308/1‑99 (MQ=255)
                                                                     ttGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATc                        <  1:1404923/100‑1 (MQ=255)
                                                                      tGTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCt                       <  2:496679/100‑1 (MQ=255)
                                                                       gTCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCtt                      >  2:1814483/1‑100 (MQ=255)
                                                                        tCCAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCttt                     <  1:602286/100‑1 (MQ=255)
                                                                         ccAGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTc                    >  1:5256/1‑100 (MQ=255)
                                                                           aGTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGa                  <  2:456657/100‑1 (MQ=255)
                                                                            gTTAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAg                 <  1:985811/100‑1 (MQ=255)
                                                                             ttAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGt                <  2:766695/100‑1 (MQ=255)
                                                                             ttAGTAGGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGt                >  1:979723/1‑100 (MQ=255)
                                                                              tagtagGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTg               <  1:549214/100‑1 (MQ=255)
                                                                                gtagGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGct             >  2:180062/1‑100 (MQ=255)
                                                                                 tagGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGctt            >  2:872371/1‑100 (MQ=255)
                                                                                 tagGAGTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGctt            <  1:1051651/100‑1 (MQ=255)
                                                                                      gTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCtttt       <  2:1952655/100‑1 (MQ=255)
                                                                                      gTGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCtttt       <  1:1777584/100‑1 (MQ=255)
                                                                                       tGCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCttttt      <  1:582764/100‑1 (MQ=255)
                                                                                        gCCACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTc     <  2:1172459/100‑1 (MQ=255)
                                                                                         ccACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTc     <  2:1857773/99‑1 (MQ=255)
                                                                                          cACCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAc   <  1:1814499/100‑1 (MQ=255)
                                                                                           aCCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACt  >  1:1806299/1‑100 (MQ=255)
                                                                                           aCCTCCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACt  <  1:125380/100‑1 (MQ=255)
                                                                                               |                                                                                               
CGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACT  >  NC_000913/1413614‑1413804

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