Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 265,418 A→G R42R (CGT→CGC)  ykfF ← CP4‑6 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913265,4180AG98.5% 239.1 / ‑4.4 68R42R (CGT→CGC) ykfFCP4‑6 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (30/37);  total (30/38)
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.06e-01

CTGTACGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGC  >  NC_000913/265321‑265514
                                                                                                 |                                                                                                
cTGTACGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATcgcg                                                                                                >  1:1908343/1‑100 (MQ=255)
 tGTACGGATGCCTGAGGTGCGGATATAGCGGTTAAGACATTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGa                                                                                               >  2:10334/1‑100 (MQ=255)
    aCGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAacc                                                                                            <  2:49075/100‑1 (MQ=255)
    aCGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAacc                                                                                            >  2:1703136/1‑100 (MQ=255)
     cGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAacca                                                                                           >  2:387633/1‑100 (MQ=255)
          gCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGa                                                                                      <  2:1647170/100‑1 (MQ=255)
          gCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGa                                                                                      >  1:178295/1‑100 (MQ=255)
           ccTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGAc                                                                                     >  2:1691134/1‑100 (MQ=255)
            cTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACg                                                                                    <  2:319120/100‑1 (MQ=255)
                ggTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAg                                                                                >  2:264183/1‑100 (MQ=255)
                ggTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAg                                                                                >  2:1781031/1‑100 (MQ=255)
                 gTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAg                                                                                <  2:1161303/99‑1 (MQ=255)
                   gCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACCGAAGTGa                                                                             >  2:152213/1‑100 (MQ=255)
                    cGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGAc                                                                            >  2:1183042/1‑100 (MQ=255)
                     ggATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGAc                                                                            >  2:103965/1‑99 (MQ=255)
                         ataGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGccc                                                                       <  1:664559/100‑1 (MQ=255)
                           aGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTg                                                                     >  1:1594139/1‑100 (MQ=255)
                           aGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTg                                                                     >  1:1531304/1‑100 (MQ=255)
                            gCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTgg                                                                    >  1:1777610/1‑100 (MQ=255)
                              ggTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGgtc                                                                  <  2:867079/100‑1 (MQ=255)
                              ggTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGgtc                                                                  >  1:961091/1‑100 (MQ=255)
                               gtgaaGACCTTCACCGTCATCCGGCTCAAAGTTCCATGCCCGCCGGACCAGCCGTGCTTCATTATCGCGAACCACCAGGCGGAAGTTACTTCCCGGgtgg                                                                 >  1:1622272/4‑98 (MQ=19)
                               gTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCTAACCACCAGACGGAAGTGACTGCCCTGgtcg                                                                 >  1:668408/1‑100 (MQ=255)
                                 tAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGgtcgtc                                                               <  2:463712/100‑1 (MQ=255)
                                  aaGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCt                                                              <  2:974837/100‑1 (MQ=255)
                                   aGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCtt                                                             <  1:1731486/100‑1 (MQ=255)
                                     aCCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCg                                                           <  2:1552707/100‑1 (MQ=255)
                                       cTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAg                                                         <  1:1534797/100‑1 (MQ=255)
                                       cTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAg                                                         <  2:1729623/100‑1 (MQ=255)
                                       cTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAg                                                         >  2:1265618/1‑100 (MQ=255)
                                         tCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAgtg                                                       >  2:522270/1‑100 (MQ=255)
                                          cACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAgtgt                                                      <  2:201724/100‑1 (MQ=255)
                                          cACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCAGCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAgtgt                                                      <  2:178296/100‑1 (MQ=255)
                                            ccGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGa                                                    <  1:11735/100‑1 (MQ=255)
                                            ccGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGa                                                    <  1:1793754/100‑1 (MQ=255)
                                             cGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGat                                                   >  1:815569/1‑100 (MQ=255)
                                              ggCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGata                                                  >  2:9507/1‑100 (MQ=255)
                                               gCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGatat                                                 <  1:780296/100‑1 (MQ=255)
                                                  tCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGc                                              <  2:950431/100‑1 (MQ=255)
                                                     ggCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGt                                           >  1:1935554/1‑100 (MQ=255)
                                                      gCTCAAAGTTCCATGCCCGCCAGACAATCCGGCCTTCAGTATCGCGAACCACCAGACGGAACTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTa                                          <  1:1550115/100‑1 (MQ=255)
                                                       cTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTAc                                         <  1:103965/100‑1 (MQ=255)
                                                        tCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACg                                        >  2:444943/1‑100 (MQ=255)
                                                            aGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGt                                    <  1:1390253/100‑1 (MQ=255)
                                                              ttCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAg                                  <  2:541915/100‑1 (MQ=255)
                                                              ttCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAg                                  <  2:541931/100‑1 (MQ=255)
                                                              ttCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAg                                  <  2:630925/100‑1 (MQ=255)
                                                                ccATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTTGTAGTg                                <  1:1058931/100‑1 (MQ=255)
                                                                ccATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTg                                <  1:1494729/100‑1 (MQ=255)
                                                                ccATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTg                                >  2:1434920/1‑100 (MQ=255)
                                                                ccATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTg                                >  1:721108/1‑100 (MQ=255)
                                                                   tGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGAcc                             >  2:1534664/1‑100 (MQ=255)
                                                                   tGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGAcc                             <  2:1531312/100‑1 (MQ=255)
                                                                       cGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcgctt                         >  2:1275649/1‑100 (MQ=255)
                                                                        gCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcgcttg                        <  2:893024/100‑1 (MQ=255)
                                                                           aGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcgcttgcgc                     >  1:61824/1‑100 (MQ=255)
                                                                            gACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcgcttgcgct                    <  2:1908356/100‑1 (MQ=255)
                                                                              ccATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcgcttgcgcttg                  <  2:1392117/100‑1 (MQ=255)
                                                                                aTCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGtc                >  1:216792/1‑100 (MQ=255)
                                                                                 tCCGGCCTTCAGTATCGTGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGtc                <  1:1183040/99‑1 (MQ=255)
                                                                                     gCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCgg           <  2:961092/100‑1 (MQ=255)
                                                                                      ccTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCggt          <  2:1935567/100‑2 (MQ=255)
                                                                                      ccTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCggg          >  2:1458865/1‑100 (MQ=255)
                                                                                      ccTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCggg          <  2:815570/100‑1 (MQ=255)
                                                                                      ccTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCggg          <  1:183449/100‑1 (MQ=255)
                                                                                       cTTCAGTATCGCGAACCACCCGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGt         <  2:61824/100‑1 (MQ=255)
                                                                                             tATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAgg   <  2:216793/100‑1 (MQ=255)
                                                                                             tATCGCGAACCACCAGACGGAAGTGACTACCCTGGTCGTCTTCGAGTCTGATATTACTGTACGTGGTAGTGACCGCTTGCTCTTGTCTCCGGGTGAAAg    <  1:387631/99‑1 (MQ=255)
                                                                                              aTCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGc  >  2:1197245/1‑100 (MQ=255)
                                                                                                 |                                                                                                
CTGTACGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGC  >  NC_000913/265321‑265514

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