Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,781,720 T→C 25.6% R57R (CGT→CGC yfjD → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,781,7200TC25.6% 115.7 / 38.1 101R57R (CGT→CGCyfjDpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (49/26);  new base C (18/8);  total (67/34)
Fisher's exact test for biased strand distribution p-value = 8.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.28e-01

Only 100 of 107 total aligned reads displayed.

TTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGA  >  minE/1781658‑1781788
                                                              |                                                                    
ttAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTAc                                                                >  1:2387500/1‑69 (MQ=255)
ttAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTAc                                                                >  1:1378778/1‑69 (MQ=255)
    tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc                                                             <  1:675052/68‑1 (MQ=255)
    tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc                                                             <  1:1895089/68‑1 (MQ=255)
    tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc                                                             <  1:194329/68‑1 (MQ=255)
    tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc                                                             <  1:2835426/68‑1 (MQ=255)
    tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc                                                             <  1:3033260/68‑1 (MQ=255)
    tCGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc                                                             <  1:635200/68‑1 (MQ=255)
     cGGCAATAATCTGGTCAATATCCTTGCTTCCGCACTCGGCACTATTGTCGGGATGCGCTTGTATGGc                                                             <  1:3029205/67‑1 (MQ=255)
      ggCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGAt                                                          >  1:2041821/1‑69 (MQ=255)
            aaCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGc                                                    >  1:1583131/1‑69 (MQ=255)
             aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                                                     <  1:3083514/67‑1 (MQ=255)
             aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                                                     >  1:2963043/1‑67 (MQ=255)
             aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                                                     <  1:2883012/67‑1 (MQ=255)
             aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                                                     <  1:419064/67‑1 (MQ=255)
             aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                                                     <  1:2846702/67‑1 (MQ=255)
             aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                                                     >  1:505960/1‑67 (MQ=255)
             aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                                                     >  1:2692847/1‑67 (MQ=255)
             aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                                                     <  1:2505300/67‑1 (MQ=255)
             aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                                                     <  1:552721/67‑1 (MQ=255)
             aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                                                     <  1:2233189/67‑1 (MQ=255)
             aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                                                     <  1:1066511/67‑1 (MQ=255)
             aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                                                     <  1:187804/67‑1 (MQ=255)
             aCCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCggg                                                     >  1:12834/1‑67 (MQ=255)
                 ggTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGc                                               <  1:2381755/69‑1 (MQ=255)
                 ggTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGc                                               <  1:2593498/69‑1 (MQ=255)
                    cAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAAt                                            <  1:2935050/69‑1 (MQ=255)
                    cAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAAt                                            <  1:763862/69‑1 (MQ=255)
                    cAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAAt                                            <  1:1239142/69‑1 (MQ=255)
                      atatTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAAt                                            <  1:2522312/67‑1 (MQ=255)
                       tatTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAAtt                                           >  1:2963841/1‑67 (MQ=255)
                       tatTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTg                                          >  1:2494517/1‑68 (MQ=255)
                        atTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCg                                                   >  1:2816003/1‑58 (MQ=255)
                        atTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGc                                         >  1:2399167/1‑68 (MQ=255)
                        atTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGc                                         >  1:2820023/1‑68 (MQ=255)
                        atTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGc                                         >  1:2087856/1‑68 (MQ=255)
                        atTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCg                                        >  1:2841937/1‑69 (MQ=255)
                        atTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCg                                        >  1:3341696/1‑69 (MQ=255)
                        atTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCg                                        >  1:477606/1‑69 (MQ=255)
                        atTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCg                                        >  1:187557/1‑69 (MQ=255)
                         ttCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGa                                       >  1:3246407/1‑69 (MQ=255)
                         ttCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGa                                       >  1:1280030/1‑69 (MQ=255)
                                     cgcTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCt                             <  1:1527102/67‑1 (MQ=255)
                                     cgcTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCt                             <  1:2904288/67‑1 (MQ=255)
                                     cgcTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCt                             <  1:556798/67‑1 (MQ=255)
                                     cgcTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCt                             <  1:568684/67‑1 (MQ=255)
                                     cgcTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCt                             <  1:39859/67‑1 (MQ=255)
                                         cGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACttt                       >  1:2964842/1‑69 (MQ=255)
                                                aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCg                                        >  1:1897545/1‑45 (MQ=255)
                                                aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACt                                     >  1:1402707/1‑48 (MQ=255)
                                                aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGTCTTTTGTCGt                 >  1:3410559/1‑68 (MQ=255)
                                                aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTc                   >  1:3086236/1‑66 (MQ=255)
                                                aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCg                  >  1:2747033/1‑67 (MQ=255)
                                                aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGt                 >  1:239761/1‑68 (MQ=255)
                                                aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGt                 >  1:773059/1‑68 (MQ=255)
                                                aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGt                 >  1:1976946/1‑68 (MQ=255)
                                                aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGt                 >  1:3237664/1‑68 (MQ=255)
                                                aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGt                 >  1:730417/1‑68 (MQ=255)
                                                aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGt                 >  1:2055274/1‑68 (MQ=255)
                                                aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGt                 >  1:2786024/1‑68 (MQ=255)
                                                aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGt                 >  1:2098523/1‑68 (MQ=255)
                                                aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGt                 >  1:653992/1‑68 (MQ=255)
                                                aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGt                 >  1:1560504/1‑68 (MQ=255)
                                                aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGt                 >  1:2371700/1‑68 (MQ=255)
                                                aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGt                 >  1:2374221/1‑68 (MQ=255)
                                                aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTa                >  1:3359722/1‑69 (MQ=255)
                                                aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTa                <  1:493486/69‑1 (MQ=255)
                                                aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTa                >  1:42187/1‑69 (MQ=255)
                                                aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTa                <  1:489925/69‑1 (MQ=255)
                                                aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTa                >  1:1347952/1‑69 (MQ=255)
                                                aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTa                >  1:1386945/1‑69 (MQ=255)
                                                aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTa                <  1:1417593/69‑1 (MQ=255)
                                                aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTa                <  1:1051496/69‑1 (MQ=255)
                                                aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTa                >  1:2134023/1‑69 (MQ=255)
                                                aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTa                >  1:1921077/1‑69 (MQ=255)
                                                aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTa                >  1:1883298/1‑69 (MQ=255)
                                                aTTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTa                >  1:887564/1‑69 (MQ=255)
                                                 attttgGGATGCGTTTGTACGGCGATGCGGGCGTGGCAAATGCGACTGGTGTGCTGACTTTTGTCGTAc               >  1:1753910/4‑69 (MQ=255)
                                                         aTGCGCTTGTATGGCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGAcct                        >  1:598451/1‑50 (MQ=37)
                                                         aTGCGCTTGTATGGCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGAcct                        >  1:117257/1‑50 (MQ=37)
                                                         aTGCGCTTGTATGGCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGAcct                        >  1:1163817/1‑50 (MQ=37)
                                                         aTGCGCTTGTATGGCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGAcct                        >  1:708668/1‑50 (MQ=37)
                                                         aTGCGCTTGTATGGCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGAcct                        >  1:787040/1‑50 (MQ=37)
                                                         aTGCGCTTGTATGGCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGAcct                        >  1:952090/1‑50 (MQ=37)
                                                         aTGCGCTTGTATGGCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGAcct                        >  1:1350990/1‑50 (MQ=37)
                                                         aTGCGCTTGTATGGCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGAcct                        >  1:3356089/1‑50 (MQ=37)
                                                         aTGCGCTTGTATGGCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGAcct                        >  1:3139369/1‑50 (MQ=37)
                                                         aTGCGCTTGTATGGCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGAcct                        >  1:3061494/1‑50 (MQ=37)
                                                         aTGCGCTTGTATGGCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGAcct                        >  1:1323499/1‑50 (MQ=37)
                                                         aTGCGCTTGTATGGCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGAcct                        >  1:1326008/1‑50 (MQ=37)
                                                         aTGCGCTTGTATGGCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGAcct                        >  1:1469274/1‑50 (MQ=37)
                                                         aTGCGCTTGTATGGCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGAcct                        >  1:1575575/1‑50 (MQ=37)
                                                         aTGCGCTTGTATGGCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGAcct                        >  1:2010940/1‑50 (MQ=37)
                                                         aTGCGCTTGTATGGCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGAcct                        >  1:1935664/1‑50 (MQ=37)
                                                         aTGCGCTTGTATGGCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGAcct                        >  1:1644877/1‑50 (MQ=37)
                                                         aTGCGCTTGTATGGCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGAcct                        >  1:170944/1‑50 (MQ=37)
                                                              tttGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGa  >  1:2217939/1‑69 (MQ=255)
                                                              tttGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGa  >  1:2201517/1‑69 (MQ=255)
                                                              tttGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGa  >  1:1592228/1‑69 (MQ=255)
                                                              tttGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGa  >  1:185610/1‑69 (MQ=255)
                                                              |                                                                    
TTAATCGGCAATAACCTGGTCAATATTCTTGCCTCCGCGCTCGGCACTATTGTTGGGATGCGTTTGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGA  >  minE/1781658‑1781788

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