Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,197,853 T→C 100% N129D (AAT→GAT)  ymfE ← e14 prophage; putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,197,8530TC96.4% 265.0 / ‑1.0 83N129D (AAT→GAT) ymfEe14 prophage; putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (36/44);  total (37/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.22e-01

GAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGATA  >  NC_000913/1197756‑1197949
                                                                                                 |                                                                                                
gAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGc                                                                                                <  1:1087164/100‑1 (MQ=255)
gAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGc                                                                                                >  1:296400/1‑100 (MQ=255)
gAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGc                                                                                                >  1:297097/1‑100 (MQ=255)
  aaaTTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCt                                                                                               >  1:1849637/1‑99 (MQ=255)
  aaaTTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTc                                                                                              <  2:11960/100‑1 (MQ=255)
   aaTTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCa                                                                                             <  1:1789935/100‑1 (MQ=255)
   aaTTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCa                                                                                             <  2:544162/100‑1 (MQ=255)
     ttCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCAtt                                                                                           >  1:1323199/1‑100 (MQ=255)
     ttCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCAtt                                                                                           <  1:445622/100‑1 (MQ=255)
     ttCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCAtt                                                                                           <  1:2001/100‑1 (MQ=255)
       cAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTAc                                                                                         >  2:2498064/1‑100 (MQ=255)
         aaaaaTATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACAt                                                                                       >  1:2519220/1‑100 (MQ=255)
            aaTATCCAAATATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGTGTAATATAGATCGCGCATTACATtgt                                                                                    <  1:930386/100‑1 (MQ=255)
             atatCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATtgtg                                                                                   >  2:331184/1‑100 (MQ=255)
              tatCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGa                                                                                  <  1:679665/100‑1 (MQ=255)
                tCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGc                                                                                <  2:1323206/100‑1 (MQ=255)
                 tcAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCt                                                                               <  1:1200360/99‑1 (MQ=255)
                  cAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCtt                                                                              >  2:1556246/1‑100 (MQ=255)
                  cAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCtt                                                                              >  2:464682/1‑100 (MQ=255)
                  cAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCta                                                                              >  2:1200364/1‑99 (MQ=255)
                            ataCTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCtttatt                                                                    <  2:2410844/100‑1 (MQ=255)
                             taCTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCtttattt                                                                   >  2:1939893/1‑100 (MQ=255)
                             taCTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCtttattt                                                                   <  2:479832/100‑1 (MQ=255)
                              aCTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCtttattta                                                                  >  2:774003/1‑100 (MQ=255)
                               cTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCtttattta                                                                  >  2:1290001/1‑99 (MQ=255)
                               cTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTTTTTaa                                                                 >  1:118616/1‑100 (MQ=255)
                                tttATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATATATCGCTCATTACATTGTGAGCTTCATCTTTATTTaa                                                                 <  1:2224599/99‑1 (MQ=255)
                                  ctttCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAAttt                                                              <  2:1958110/98‑1 (MQ=255)
                                    ttCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAAtttag                                                            >  1:1958100/1‑98 (MQ=255)
                                    ttCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTc                                                            <  2:2325086/100‑1 (MQ=255)
                                     tCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCt                                                           >  2:1031092/1‑100 (MQ=255)
                                      ccAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTg                                                          >  2:1154280/1‑100 (MQ=255)
                                      ccAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTg                                                          <  2:792688/100‑1 (MQ=255)
                                       cAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGt                                                         <  2:2299610/100‑1 (MQ=255)
                                        agaagaGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGtt                                                        <  2:1976703/100‑1 (MQ=255)
                                         gaagaGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTg                                                       >  2:542596/1‑100 (MQ=255)
                                         gaagaGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTg                                                       <  2:1610787/100‑1 (MQ=255)
                                         gaagaGTTCAATAGAATGTTTGTCTTCAATTTTTCTTACTTCATGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTg                                                       <  1:943691/100‑1 (MQ=255)
                                            gagTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACt                                                    <  2:1849647/100‑1 (MQ=255)
                                            gagTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACt                                                    <  2:1429434/100‑1 (MQ=255)
                                            gagTTCAATATAATGTTTGTCTTCAATTTTTCTTAATTAAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGagt                                                    >  1:1164502/1‑98 (MQ=255)
                                              gTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTCATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTc                                                   >  1:30414/1‑99 (MQ=255)
                                              gTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTcc                                                  <  1:1275203/100‑1 (MQ=255)
                                              gTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTcc                                                  <  2:557479/100‑1 (MQ=255)
                                              gTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTcc                                                  >  1:134546/1‑100 (MQ=255)
                                               ttCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCa                                                 <  1:2498045/100‑1 (MQ=255)
                                                 cAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGc                                               >  2:2175577/1‑100 (MQ=255)
                                                       aaTGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCg                                         <  1:2003293/100‑1 (MQ=255)
                                                            ttGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATa                                    <  2:2520650/100‑1 (MQ=255)
                                                              gTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAAc                                  <  2:2418049/100‑1 (MQ=255)
                                                              gTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAAc                                  >  2:1023966/1‑100 (MQ=255)
                                                                 ttCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGt                               <  2:1046134/100‑1 (MQ=255)
                                                                 ttCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGt                               >  2:2451480/1‑100 (MQ=255)
                                                                 ttCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGt                               >  2:514467/1‑100 (MQ=255)
                                                                  tCAATTTTTTTTTTTTTTTGGTAATATAAATCGCTCATTACATTGTGAGCTTCATCTGTATTTAATTTTCTGTTGACTCCAGCTCTCCGTTATATCTGtt                              >  1:1011830/1‑100 (MQ=255)
                                                                  tCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCCTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGtt                              >  2:1909667/1‑100 (MQ=255)
                                                                   cAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGttt                             <  2:2219392/100‑1 (MQ=255)
                                                                     aTTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTta                           <  2:175445/100‑1 (MQ=255)
                                                                     aTTTTTCTTACTTCAGGGTAATATAGATCGCTCATTAAATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCAACGTGATAACGGTTTta                           <  1:1509582/100‑1 (MQ=255)
                                                                      tttttCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTtat                          <  1:1664668/100‑1 (MQ=255)
                                                                          tCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAAtt                      <  1:1871294/100‑1 (MQ=255)
                                                                           cTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTa                     <  1:1392265/100‑1 (MQ=255)
                                                                           cTTACCTCAGTGTATCATCGATCGCTAAATACATTGTGACATACATCTTTATTTAATTTACTGCTGACTCGCGATCAAAGTGATAACGGTTTTATGAtt                      <  1:859698/99‑1 (MQ=255)
                                                                             tACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGa                   >  2:2034166/1‑100 (MQ=255)
                                                                             tACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGa                   >  2:1783717/1‑100 (MQ=255)
                                                                              aCTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGAt                  <  2:988096/100‑1 (MQ=255)
                                                                               cTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATg                 >  2:2066896/1‑100 (MQ=255)
                                                                                ttCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGc                <  1:1852142/100‑1 (MQ=255)
                                                                                    gggTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTAt            <  2:595661/100‑1 (MQ=255)
                                                                                      gTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATcc          <  1:1956959/100‑1 (MQ=255)
                                                                                      gTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATcc          <  2:1516558/100‑1 (MQ=255)
                                                                                      gTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATcc          <  2:1616029/100‑1 (MQ=255)
                                                                                       tAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATccc         >  2:1240208/1‑100 (MQ=255)
                                                                                       tAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATccc         <  1:2431103/100‑1 (MQ=255)
                                                                                         atatAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCaa       >  2:2435707/1‑100 (MQ=255)
                                                                                         atatAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCaa       >  2:2463522/1‑100 (MQ=255)
                                                                                           ataGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAg     <  1:1307002/100‑1 (MQ=255)
                                                                                           ataGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAg     <  2:1680176/100‑1 (MQ=255)
                                                                                            taGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAg     >  2:1982672/1‑99 (MQ=255)
                                                                                             aGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTGTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGat   <  1:882038/100‑1 (MQ=255)
                                                                                              gATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGata  >  1:839327/1‑100 (MQ=255)
                                                                                              gATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGaag  >  1:2388674/1‑98 (MQ=255)
                                                                                              gATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGCGATAACGGTTTTATAATTAGATGCTTATCCCAAAGata  >  1:1411708/1‑100 (MQ=255)
                                                                                                 |                                                                                                
GAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGATA  >  NC_000913/1197756‑1197949

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