Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,162,506 T→A 48.8% M3K (ATG→AAG)  ydjZ → conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,162,5060TA48.8% ‑0.5 / 102.5 84M3K (ATG→AAG) ydjZconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (7/36);  new base A (21/20);  total (28/56)
Fisher's exact test for biased strand distribution p-value = 1.04e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

Only 100 of 104 total aligned reads displayed.

GGCGGATAACACGCTGGCAACGGTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCAATCGCGAAAAATCTGGTACTACCGTATAACCCTCATCATCCTGTTGTTCGCTAT  >  minE/1162439‑1162566
                                                                   |                                                            
ggcggATAACACGCTGGCAACGGTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGgatgatgatg                                                             >  1:1613017/1‑69 (MQ=255)
ggcggATAACACGCTGGCAACGGTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGgatgatgat                                                              <  1:1595321/68‑1 (MQ=255)
    gATAACACGCTGGCAACGGTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGc                                                         <  1:2773632/69‑1 (MQ=255)
    gATAACACGCTGGCAACGGTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGc                                                         <  1:3233280/69‑1 (MQ=255)
    gATAACACGCTGGCAACGGTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGc                                                         <  1:696847/69‑1 (MQ=255)
      tAACACGCTGGCAACGGTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCaa                                                       >  1:3085594/1‑69 (MQ=255)
        acacGCTGGCAACGGTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCAATc                                                     >  1:1651856/1‑69 (MQ=255)
          acGCTGGCAACGGTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCAATcg                                                    >  1:2785206/1‑68 (MQ=255)
          acGCTGGCAACGGTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCAATcg                                                    >  1:2945717/1‑68 (MQ=255)
          acGCTGGCAACGGTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCAATcg                                                    >  1:1741778/1‑68 (MQ=255)
             cTGGCAACGGTTACCTTTAAAATCACCGAATAAAGCAGGATAACGGATGATGAAGATGCAATCGCGtaa                                               <  1:2112329/69‑1 (MQ=255)
             cTGGCAACGGTTACCTTTAAAATCACCGAATAAAGCAGGATAACGGATGATGAAGATGCAATCGCGtaa                                               <  1:2007976/69‑1 (MQ=255)
             cTGGCAACGGTTACCTTTAAAATCACCGAATAAAGCAGGATAACGGATGATGAAGATGCAATCGCGtaa                                               <  1:1868852/69‑1 (MQ=255)
             cTGGCAACGGTTACCTTTAAAATCACCGAATAAAGCAGGATAACGGATGATGAAGATGCAATCGCGtaa                                               <  1:317715/69‑1 (MQ=255)
             cTGGCAACGGTTACCTTTAAAATCACCGAATAAAGCAGGATAACGGATGATGAAGATGCAATCGCGtaa                                               <  1:403217/69‑1 (MQ=255)
             cTGGCAACGGTTACCTTTAAAATCACCGAATAAAGCAGGATAACGGATGATGAAGATGCAATCGCGtaa                                               <  1:428076/69‑1 (MQ=255)
             cTGGCAACGGTTACCTTTAAAATCACCGAATAAAGCAGGATAACGGATGATGAAGATGCAATCGCGtaa                                               <  1:1441706/69‑1 (MQ=255)
             cTGGCAACGGTTACCTTTAAAATCACCGAATAAAGCAGGATAACGGATGATGAAGATGCAATCGCGtaa                                               <  1:2223879/69‑1 (MQ=255)
             cTGGCAACGGTTACCTTTAAAATCACCGAATAAAGCAGGATAACGGATGATGAAGATGCAATCGCGtaa                                               <  1:1382048/69‑1 (MQ=255)
             cTGGCAACGGTTACCTTTAAAATCACCGAATAAAGCAGGATAACGGATGATGAAGATGCAATCGCGtaa                                               <  1:604920/69‑1 (MQ=255)
             cTGGCAACGGTTACCTTTAAAATCACCGAATAAAGCAGGATAACGGATGATGAAGATGCAATCGCGtaa                                               <  1:694397/69‑1 (MQ=255)
             cTGGCAACGGTTACCTTTAAAATCACCGAATAAAGCAGGATAACGGATGATGAAGATGCAATCGCGtaa                                               <  1:2232104/69‑1 (MQ=255)
             cTGGCAACGGTTACCTTTAAAATCACCGAATAAAGCAGGATAACGGATGATGAAGATGCAATCGCGtaa                                               <  1:83977/69‑1 (MQ=255)
             cTGGCAACGGTTACCTTTAAAATCACCGAATAAAGCAGGATAACGGATGATGAAGATGCAATCGCGtaa                                               <  1:1161609/69‑1 (MQ=255)
             cTGGCAACGGTTACCTTTAAAATCACCGAATAAAGCAGGATAACGGATGATGAAGATGCAATCGCGtaa                                               <  1:1197583/69‑1 (MQ=255)
             cTGGCAACGGTTACCTTTAAAATCACCGAATAAAGCAGGATAACGGATGATGAAGATGCAATCGCGtaa                                               <  1:1187059/69‑1 (MQ=255)
             cTGGCAACGGTTACCTTTAAAATCACCGAATAAAGCAGGATAACGGATGATGAAGATGCAATCGCGtaa                                               <  1:1176286/69‑1 (MQ=255)
             cTGGCAACGGTTACCTTTAAAATCACCGAATAAAGCAGGATAACGGATGATGAAGATGCAATCGCGtaa                                               <  1:774555/69‑1 (MQ=255)
              tGGCAACGGTTACCTTTAAAATCACCGAATAAAGCAGGATAACGGATGATGAAGATGCAATCGCGtaa                                               <  1:3255707/68‑1 (MQ=255)
                   aCGGTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGgatgatgat                                                              <  1:3226587/49‑1 (MQ=255)
                   aCGGTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGgatgatgat                                                              <  1:2086135/49‑1 (MQ=255)
                   aCGGTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGgatgatgat                                                              <  1:2140909/49‑1 (MQ=255)
                   aCGGTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGgatgatgat                                                              <  1:313435/49‑1 (MQ=255)
                   aCGGTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGgatgatgat                                                              <  1:1523825/49‑1 (MQ=255)
                   aCGGTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGgatgatgat                                                              <  1:1246881/49‑1 (MQ=255)
                   aCGGTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGgatgatgat                                                              <  1:1253709/49‑1 (MQ=255)
                   aCGGTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGgatgatgat                                                              <  1:2756421/49‑1 (MQ=255)
                     ggTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCAATCGCGAAAAATCTgg                                        <  1:2420909/69‑1 (MQ=255)
                     ggTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCAATCGCGAAAAATCTgg                                        <  1:39483/69‑1 (MQ=255)
                     ggTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCAATCGCGAAAAATCTgg                                        <  1:2791726/69‑1 (MQ=255)
                     ggTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCAATCGCGAAAAATCTgg                                        <  1:37019/69‑1 (MQ=255)
                     ggTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCAATCGCGAAAAATCTgg                                        <  1:24380/69‑1 (MQ=255)
                     ggTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCAATCGCGAAAAATCTgg                                        <  1:521681/69‑1 (MQ=255)
                     ggTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCAATCGCGAAAAATCTgg                                        <  1:920617/69‑1 (MQ=255)
                     ggTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCAATCGCGAAAAATCTgg                                        <  1:1702353/69‑1 (MQ=255)
                     ggTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCAATCGCGAAAAATCTgg                                        <  1:2196785/69‑1 (MQ=255)
                     ggTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCAATCGCGAAAAATCTgg                                        <  1:1870603/69‑1 (MQ=255)
                     ggTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCAATCGCGAAAAATCTgg                                        <  1:1873194/69‑1 (MQ=255)
                     ggTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCAATCGCGAAAAATCTgg                                        <  1:1987008/69‑1 (MQ=255)
                     ggTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCAATCGCGAAAAATCTgg                                        <  1:1333755/69‑1 (MQ=255)
                     ggTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCAATCGCGAAAAATCTgg                                        <  1:1344828/69‑1 (MQ=255)
                        tACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCAATCGCGAAAAATCTGGtac                                     >  1:2309605/1‑69 (MQ=255)
                        tACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCAATCGCGAAAAATCTGGtac                                     >  1:785219/1‑69 (MQ=255)
                            tttAAAATCGCCGAATAAAGTCAGGATAAAGGATGATGATGATGCAATCGCGAAAAATCTGTtactac                                  <  1:2450103/68‑1 (MQ=255)
                            tttAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCAATCGCGAAAAATCTGGtactac                                  <  1:575851/68‑1 (MQ=255)
                            tttAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCAATCGCGAAAAATCTGGtactac                                  <  1:1847434/68‑1 (MQ=255)
                            tttAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCAATCGCGAAAAATCTGGtactac                                  <  1:3098052/68‑1 (MQ=255)
                            tttAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCAATCGCGAAAAATCTGGtactac                                  <  1:1689982/68‑1 (MQ=255)
                            tttAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCAATCGCGAAAAATCTGGtactac                                  <  1:452173/68‑1 (MQ=255)
                            tttAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCAATCGCGAAAAATCTGGtactac                                  <  1:2451487/68‑1 (MQ=255)
                            tttAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCAATCGCGAAAAATCTGGtactac                                  <  1:79215/68‑1 (MQ=255)
                            tttAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCAATCGCGAAAAATCTGGtactac                                  <  1:1152982/68‑1 (MQ=255)
                            tttAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCAATCGCGAAAAATCTGGtactac                                  <  1:2288312/68‑1 (MQ=255)
                            tttAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCAATCGCGAAAAATCTGGtactac                                  <  1:2375014/68‑1 (MQ=255)
                                           tAAAGCAGGATAACGGATGATGAAGATGCAATCGCGtaa                                               <  1:591162/39‑1 (MQ=16)
                                                        acggATGATGAAGATGCAATCGCGTAAAATCTGGtactac                                  >  1:2004225/3‑40 (MQ=37)
                                                        acggATGATGAAGATGCAATCGCGTAAAATCTGGTACTACCGTATaa                           >  1:1656630/3‑47 (MQ=38)
                                                        acggATGATGAAGATGCAATCGCGTAAAATCTGGTACTACCGTATAAccc                        >  1:1185324/3‑50 (MQ=255)
                                                        acggATGATGAAGATGCAATCGCGTAAAATCTGGTACTACCGTATAACCCt                       >  1:115265/3‑51 (MQ=255)
                                                        acggATGATGAAGATGCAATCGCGTAAAATCTGGTACTACCGTATAACCCt                       >  1:142379/3‑51 (MQ=255)
                                                        acggATGATGAAGATGCAATCGCGTAAAATCTGGTACTACCGTATAACCCt                       >  1:587222/3‑51 (MQ=255)
                                                        acggATGATGAAGATGCAATCGCGTAAAATCTGGTACTACCGTATAACCCt                       >  1:3300012/3‑51 (MQ=255)
                                                        acggATGATGAAGATGCAATCGCGTAAAATCTGGTACTACCGTATAACCCt                       >  1:3299631/3‑51 (MQ=255)
                                                        acggATGATGAAGATGCAATCGCGTAAAATCTGGTACTACCGTATAACCCt                       >  1:1713891/3‑51 (MQ=255)
                                                        acggATGATGAAGATGCAATCGCGTAAAATCTGGTACTACCGTATAACCCt                       >  1:3199882/3‑51 (MQ=255)
                                                        acggATGATGAAGATGCAATCGCGTAAAATCTGGTACTACCGTATAACCCt                       >  1:2419846/3‑51 (MQ=255)
                                                        acggATGATGAAGATGCAATCGCGTAAAATCTGGTACTACCGTATAACCCt                       >  1:2422342/3‑51 (MQ=255)
                                                        acggATGATGAAGATGCAATCGCGTAAAATCTGGTACTACCGTATAACCCt                       >  1:1097678/3‑51 (MQ=255)
                                                        acggATGATGAAGATGCAATCGCGTAAAATCTGGTACTACCGTATAACCCt                       >  1:2472114/3‑51 (MQ=255)
                                                        acggATGATGAAGATGCAATCGCGTAAAATCTGGTACTACCGTATAACCCt                       >  1:2615925/3‑51 (MQ=255)
                                                        acggATGATGAAGATGCAATCGCGTAAAATCTGGTACTACCGTATAACCCt                       >  1:2711591/3‑51 (MQ=255)
                                                        acggATGATGAAGATGCAATCGCGTAAAATCTGGTACTACCGTATAACCCt                       >  1:2189444/3‑51 (MQ=255)
                                                        acggATGATGAAGATGCAATCGCGTAAAATCTGGTACTACCGTATAACCCt                       >  1:2052896/3‑51 (MQ=255)
                                                        acggATGATGAAGATGCAATCGCGTAAAATCTGGTACTACCGTATAACCCt                       >  1:1762789/3‑51 (MQ=255)
                                                        acggATGATGAAGATGCAATCGCGTAAAATCTGGTACTACCGTATAACCCt                       >  1:3157149/3‑51 (MQ=255)
                                                        acggATGATGAAGATGCAATCGCGTAAAATCTGGTACTACCGTATAACCCt                       >  1:306804/3‑51 (MQ=255)
                                                        aaGGATGATGATGATGCAATCGCGAAAAATCTGGTACTACCGTATAACCCTCATCATCCTGTTGTTCGc     <  1:1649112/69‑1 (MQ=255)
                                                          ggATGATGATGATGCAATCGCGAAAAATCTGGTACTACCGTATAACCCTCATCATCCTGTTGTTCGCTa   <  1:2799209/69‑1 (MQ=255)
                                                          ggATGATGATGATGCAATCGCGAAAAATCTGGTACTACCGTATAACCCTCATCATCCTGTTGTTCGCTa   <  1:244892/69‑1 (MQ=255)
                                                          ggATGATGATGATGCAATCGCGAAAAATCTGGTACTACCGTATAACCCTCATCATCCTGTTGTTCGCTa   <  1:1573856/69‑1 (MQ=255)
                                                          ggATGATGATGATGCAATCGCGAAAAATCTGGTACTACCGTATAACCCTCATCATCCTGTTGTTCGCTa   <  1:1302499/69‑1 (MQ=255)
                                                           gatgatgatgatGCAATCGCGAAAAATCTGGTACTACCGTATAACCCTCATCATCCTGTTGTTCGCTAt  <  1:293315/69‑1 (MQ=255)
                                                           gatgatgatgatGCAATCGCGAAAAATCTGGTACTACCGTATAACCCTCATCATCCTGTTGTTCGCTAt  <  1:3299587/69‑1 (MQ=255)
                                                           gatgatgatgatGCAATCGCGAAAAATCTGGTACTACCGTATAACCCTCATCATCCTGTTGTTCGCTAt  <  1:345593/69‑1 (MQ=255)
                                                           gatgatgatgatGCAATCGCGAAAAATCTGGTACTACCGTATAACCCTCATCATCCTGTTGTTCGCTAt  <  1:2881639/69‑1 (MQ=255)
                                                           gatgatgatgatGCAATCGCGAAAAATCTGGTACTACCGTATAACCCTCATCATCCTGTTGTTCGCTAt  <  1:2107526/69‑1 (MQ=255)
                                                           gatgatgatgatGCAATCGCGAAAAATCTGGTACTACCGTATAACCCTCATCATCCTGTTGTTCGCTAt  <  1:531710/69‑1 (MQ=255)
                                                           gatgatgatgatGCAATCGCGAAAAATCTGGTACTACCGTATAACCCTCATCATCCTGTTGTTCGCTAt  <  1:473793/69‑1 (MQ=255)
                                                           gatgatgatgatGCAATCGCGAAAAATCTGGTACTACCGTATAACCCTCATCATCCTGTTGTTCGCTAt  <  1:2804211/69‑1 (MQ=255)
                                                            atgatgatgatgCAATCGCGAAAAATCTGGTACTACCGTATAACCCTCATCATCCTGTTGTTCGc     <  1:789301/65‑1 (MQ=255)
                                                                   |                                                            
GGCGGATAACACGCTGGCAACGGTTACCTTTAAAATCGCCGAATAAAGCCAGGATAAAGGATGATGATGATGCAATCGCGAAAAATCTGGTACTACCGTATAACCCTCATCATCCTGTTGTTCGCTAT  >  minE/1162439‑1162566

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