Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 282,417 T→C L47P (CTG→CCG)  yagE → 2‑keto‑3‑deoxy gluconate (KDG) aldolase; CP4‑6 prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913282,4170TC98.5% 222.5 / ‑4.7 65L47P (CTG→CCG) yagE2‑keto‑3‑deoxy gluconate (KDG) aldolase; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (35/29);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

CCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGT  >  NC_000913/282321‑282510
                                                                                                |                                                                                             
ccaccaTTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGTCTGTTCTTCCCggg                                                                                            <  2:1728950/100‑1 (MQ=255)
 caccaTTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGc                                                                                           >  2:1526549/1‑100 (MQ=255)
 caccaTTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGc                                                                                           <  2:228492/100‑1 (MQ=255)
     aTTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCg                                                                                       >  1:1227506/1‑100 (MQ=255)
        tttACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGt                                                                                     >  2:1217671/1‑99 (MQ=255)
          tACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGc                                                                                  <  2:230930/100‑1 (MQ=255)
           aCCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCg                                                                                 <  2:413681/100‑1 (MQ=255)
              gccgACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGt                                                                              <  2:443645/100‑1 (MQ=255)
                  aCGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAATCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTtctc                                                                          <  1:1008598/100‑1 (MQ=255)
                  aCGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTtctc                                                                          <  2:1753364/100‑1 (MQ=255)
                    ggCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTccc                                                                        <  1:1577610/100‑1 (MQ=255)
                    ggCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTccc                                                                        <  2:204541/100‑1 (MQ=255)
                     gCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTccc                                                                        <  2:291850/99‑1 (MQ=255)
                     gCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGAAGAGGTGGCAAGTTCTCCCa                                                                       >  2:719245/1‑100 (MQ=255)
                           tCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCgg                                                                 <  1:907120/100‑1 (MQ=255)
                            cGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGTGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGgc                                                                <  2:1489010/100‑1 (MQ=255)
                             gATAAGCGGTGCAACGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGgcg                                                               <  2:288667/100‑1 (MQ=255)
                             gATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGgcg                                                               >  2:327024/1‑100 (MQ=255)
                             gATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGgcg                                                               >  1:1462419/1‑100 (MQ=255)
                             gATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGgcg                                                               >  2:1403922/1‑100 (MQ=255)
                               tAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGcc                                                             <  2:1051656/100‑1 (MQ=255)
                                aaGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCg                                                            <  1:531333/100‑1 (MQ=255)
                                 aGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGa                                                           >  2:961028/1‑100 (MQ=255)
                                  gCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCAGTGGCGAGTTCTCCCAGCTCGGCGCCGaa                                                          >  2:264268/1‑100 (MQ=255)
                                  gCCGGGCACCCCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGaa                                                          >  2:874962/1‑100 (MQ=255)
                                    cGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGAGCAGCGATGGCAAGTTATACCAACACGGCACCGAaga                                                        >  2:1141800/1‑100 (MQ=255)
                                     gggCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAagag                                                       <  1:480147/100‑1 (MQ=255)
                                      ggCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGc                                                      >  1:1106671/1‑100 (MQ=255)
                                      ggCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGc                                                      <  1:1810589/100‑1 (MQ=255)
                                        cACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCga                                                    <  1:1117050/100‑2 (MQ=255)
                                         aCCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTa                                                   <  2:104123/100‑1 (MQ=255)
                                         aCCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTa                                                   >  1:1318523/1‑100 (MQ=255)
                                         aCCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTa                                                   <  2:104466/100‑1 (MQ=255)
                                           cgccgcGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTaaa                                                 <  1:33855/100‑1 (MQ=255)
                                            gccgcGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAg                                                <  2:721390/100‑1 (MQ=255)
                                              cgcgCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGcc                                              >  2:92431/1‑100 (MQ=255)
                                               gcgcTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCa                                             >  2:1213277/1‑100 (MQ=255)
                                                cgcTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCAt                                            >  2:1019601/1‑100 (MQ=255)
                                                  ctaaTCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCACAGCTTGGCGCCGAAGAGCGTAAAGCCATTg                                          <  1:240964/97‑1 (MQ=255)
                                                  cTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTg                                          >  1:686829/1‑100 (MQ=255)
                                                    gATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGcc                                        <  2:1307264/100‑1 (MQ=255)
                                                    gATCGACGACCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGcc                                        >  2:216721/1‑100 (MQ=255)
                                                       cgacgaTCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGc                                     >  2:1334709/1‑100 (MQ=255)
                                                        gacgaTCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCt                                    >  2:980761/1‑100 (MQ=255)
                                                            aTCTGATCAAAGCAGGCGTTGACGGTCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGc                                <  2:575298/100‑1 (MQ=255)
                                                              cttaTCAAAGCAGGCGTGGACCGACTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTa                              <  1:1659019/97‑1 (MQ=255)
                                                                gATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATc                            >  2:1144400/1‑100 (MQ=255)
                                                                gATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATc                            >  2:729648/1‑100 (MQ=255)
                                                                gATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATc                            >  2:220818/1‑100 (MQ=255)
                                                                gATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATc                            >  2:1089601/1‑100 (MQ=255)
                                                                 aTCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCg                           <  2:179494/100‑1 (MQ=255)
                                                                   cAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGAt                         >  1:61205/1‑100 (MQ=255)
                                                                   cAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGAt                         >  2:1065336/1‑100 (MQ=255)
                                                                   cAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCACCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGAt                         >  1:218792/1‑100 (MQ=255)
                                                                    aaaGCAGGCGTTGACGGCCTGTTCTTCCCGGGAAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATc                        >  1:863467/1‑100 (MQ=255)
                                                                      aGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCAt                      >  1:761187/1‑100 (MQ=255)
                                                                         aGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTc                   <  1:1526538/100‑1 (MQ=255)
                                                                                 aCGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGAtcgtcg           >  1:281038/1‑100 (MQ=255)
                                                                                   ggCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTcgcg         <  1:1650031/100‑1 (MQ=255)
                                                                                    gCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGt        >  1:1785707/1‑100 (MQ=255)
                                                                                     ccTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTg       <  2:1651962/100‑1 (MQ=255)
                                                                                        gTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCg    >  1:341606/1‑100 (MQ=255)
                                                                                          tcttcCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGt  >  2:316587/1‑100 (MQ=255)
                                                                                           cttcCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTcgcccgacgg   >  1:625078/1‑91 (MQ=255)
                                                                                             tcCCGGGCAGGGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTGGCTATGGGACATGTCGATCGTCGCGTGCCgg   >  1:1499859/1‑96 (MQ=255)
                                                                                                |                                                                                             
CCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGT  >  NC_000913/282321‑282510

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