Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,661,722 T→A 13.5% intergenic (‑98/‑254) yfgF ← / → yfgG predicted inner membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,661,7220TA13.5% 295.8 / 32.8 140intergenic (‑98/‑254)yfgF/yfgGpredicted inner membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (72/49);  new base A (7/12);  total (79/61)
Fisher's exact test for biased strand distribution p-value = 8.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

Only 100 of 140 total aligned reads displayed.

GTATTCAGACAGGCATCCCACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTACCAATTGATATATATCAAACAAAATAACCCTGATTAATGA  >  minE/1661662‑1661787
                                                            |                                                                 
gtATTCAGACAGGCATCCCACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTa                                                           >  1:2739866/1‑69 (MQ=255)
gtATTCAGACAGGCATCCCACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTa                                                           >  1:43544/1‑69 (MQ=255)
   ttCAGACAGGCATCCCACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAg                                                         >  1:2005771/1‑68 (MQ=255)
     cAGACAGGCATCCCACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAgtgt                                                      >  1:2670375/1‑69 (MQ=255)
     cAGACAGGCATCCCACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAgtgt                                                      >  1:631605/1‑69 (MQ=255)
     cAGACAGGCATCCCACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAgtgt                                                      >  1:3167359/1‑69 (MQ=255)
        aCAGGCATCCCACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCat                                                   <  1:3002427/69‑1 (MQ=255)
        aCAGGCATCCCACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCa                                                    <  1:2227266/68‑1 (MQ=255)
        aCAGGCATCCCACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCa                                                    <  1:257307/68‑1 (MQ=255)
        aCAGGCATCCCACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCa                                                    <  1:1929045/68‑1 (MQ=255)
        aCAGGCATCCCACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCa                                                    <  1:3255978/68‑1 (MQ=255)
         cAGGCATCCCACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCat                                                   >  1:1020183/1‑68 (MQ=255)
         cAGGCATCCCACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCat                                                   >  1:1102984/1‑68 (MQ=255)
          aGGCATCCCACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCata                                                  >  1:3323445/1‑68 (MQ=255)
          aGGCATCCCACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCata                                                  >  1:558234/1‑68 (MQ=255)
            gCATCCCACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCatat                                                 >  1:3162356/1‑67 (MQ=255)
            gCATCCCACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCatat                                                 >  1:1741717/1‑67 (MQ=255)
            gCATCCCACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCatat                                                 >  1:2247392/1‑67 (MQ=255)
             cATCCCACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTTTCATATGaa                                              <  1:1219624/69‑1 (MQ=255)
             cATCCCACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGaa                                              <  1:2251853/69‑1 (MQ=255)
             cATCCCACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGaa                                              <  1:2314811/69‑1 (MQ=255)
             cATCCCACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGaa                                              <  1:2174909/69‑1 (MQ=255)
             cATCCCACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGaa                                              <  1:782888/69‑1 (MQ=255)
             cATCCCACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGaa                                              <  1:2778004/69‑1 (MQ=255)
             cATCCCACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGaa                                              <  1:2956875/69‑1 (MQ=255)
             cATCCCACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGaa                                              <  1:546476/69‑1 (MQ=255)
             cATCCCACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGaa                                              <  1:3139566/69‑1 (MQ=255)
             cATCCCACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGaa                                              <  1:177723/69‑1 (MQ=255)
              aTCCCACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGaa                                              <  1:2824811/68‑1 (MQ=255)
                 ccACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGa                                                                >  1:684809/1‑47 (MQ=255)
                 ccACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGa                                                                >  1:2700200/1‑47 (MQ=255)
                 ccACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGa                                                                >  1:2794969/1‑47 (MQ=255)
                 ccACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGa                                                                >  1:2603945/1‑47 (MQ=255)
                 ccACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGa                                                                >  1:2530830/1‑47 (MQ=255)
                 ccACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGa                                                                >  1:955761/1‑47 (MQ=255)
                 ccACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGa                                                                >  1:2205730/1‑47 (MQ=255)
                 ccACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAg                                            >  1:610057/1‑67 (MQ=255)
                 ccACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGt                                           >  1:2279594/1‑68 (MQ=255)
                 ccACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGt                                           >  1:1596541/1‑68 (MQ=255)
                 ccACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGt                                           >  1:3197866/1‑68 (MQ=255)
                 ccACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGt                                           >  1:2189029/1‑68 (MQ=255)
                 ccACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGt                                           >  1:770034/1‑68 (MQ=255)
                 ccACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGt                                           >  1:72060/1‑68 (MQ=255)
                 ccACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGt                                           >  1:1347131/1‑68 (MQ=255)
                   aCCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTAcc                                        <  1:2215753/69‑1 (MQ=255)
                   aCCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTAcc                                        <  1:1165408/69‑1 (MQ=255)
                   aCCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTAcc                                        <  1:1587108/69‑1 (MQ=255)
                   aCCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTAcc                                        <  1:1890400/69‑1 (MQ=255)
                          ttCGAATGATGATTATTCATCACTATAGAGAGTATTGATTCTAAGTGTCATATGAAAGTACCAATTg                                   <  1:3046442/67‑1 (MQ=255)
                          ttCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTACCAATTg                                   <  1:2391458/67‑1 (MQ=255)
                          ttCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTACCAATTg                                   <  1:2387389/67‑1 (MQ=255)
                          ttCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTACCAATTg                                   <  1:2959044/67‑1 (MQ=255)
                           tCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTACCAATTGata                                <  1:606454/69‑1 (MQ=255)
                           tCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTACCAATTGata                                <  1:459407/69‑1 (MQ=255)
                           tCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTACCAATTGata                                >  1:371513/1‑69 (MQ=255)
                           tCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTACCAATTGata                                <  1:356962/69‑1 (MQ=255)
                           tCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTACCAATTGata                                >  1:734586/1‑69 (MQ=255)
                           tCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTACCAATTGata                                >  1:3355631/1‑69 (MQ=255)
                           tCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTACCAATTGata                                >  1:1582118/1‑69 (MQ=255)
                           tCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTACCAATTGata                                >  1:3313588/1‑69 (MQ=255)
                           tCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTACCAATTGata                                <  1:3303142/69‑1 (MQ=255)
                           tCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTACCAATTGata                                <  1:770323/69‑1 (MQ=255)
                           tCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTACCAATTGata                                >  1:1662031/1‑69 (MQ=255)
                           tCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTACCAATTGata                                >  1:1052008/1‑69 (MQ=255)
                           tCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTACCAATTGata                                >  1:201391/1‑69 (MQ=255)
                           tCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTACCAATTGata                                <  1:3001283/69‑1 (MQ=255)
                           tCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTACCAATTGata                                <  1:1881317/69‑1 (MQ=255)
                           tCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTACCAATTGata                                >  1:1855118/1‑69 (MQ=255)
                           tCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTACCAATTGata                                <  1:2955270/69‑1 (MQ=255)
                           tCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTACCAATTGata                                <  1:2704925/69‑1 (MQ=255)
                           tCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTACCAATTGata                                <  1:2914231/69‑1 (MQ=255)
                           tCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTACCAATTGata                                >  1:2827740/1‑69 (MQ=255)
                           tCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTACCAATTGat                                 >  1:2638597/1‑68 (MQ=255)
                           tCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTACCAATTGat                                 >  1:2081315/1‑68 (MQ=255)
                            cGAATGATGATTATTCATCACTATAGAGAGCAATGATTCTAAGTGTCATATGAAAGTTCCAATTGata                                >  1:2303810/1‑68 (MQ=255)
                            cGAATGATGATTATTCATCACTATAGAGAGCAATGATTCTAAGTGTCATATGAAAGTTCCAATTGata                                >  1:2784869/1‑68 (MQ=255)
                            cGAATGATGATTATTCATCACTATAGAGAGCAATGATTCTAAGTGTCATATGAAAGTTCCAATTGata                                >  1:492052/1‑68 (MQ=255)
                            cGAATGATGATTATTCATCACTATAGAGAGCAATGATTCTAAGTGTCATATGAAAGTTCCAATTGata                                >  1:2580227/1‑68 (MQ=255)
                                              caCTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTACCAATTGATATATATCAAACAAAATAAcc             <  1:147193/69‑1 (MQ=255)
                                               aCTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTAcc                                        >  1:1343539/1‑41 (MQ=255)
                                               aCTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTAcc                                        >  1:1953407/1‑41 (MQ=255)
                                               aCTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTAcc                                        >  1:1819856/1‑41 (MQ=255)
                                                 tataGAGAGCATTGATTCTAAGTGTCATATGAAAGTAcc                                        <  1:815661/39‑1 (MQ=255)
                                                   taGAGAGCAATGATTCTAAGTGTCATATGACAGTTCCAATTGATATATATcaaacaa                    <  1:2112424/57‑1 (MQ=255)
                                                   taGAGAGCAATGATTCTAAGTGTCATATGAAAGTTCCAATTGATATATATcaaacaa                    <  1:2812153/57‑1 (MQ=255)
                                                   taGAGAGCAATGATTCTAAGTGTCATATGAAAGTTCCAATTGATATATATcaaacaa                    <  1:2828731/57‑1 (MQ=255)
                                                   taGAGAGCAATGATTCTAAGTGTCATATGAAAGTTCCAATTGATATATATcaaacaa                    <  1:1490972/57‑1 (MQ=255)
                                                   taGAGAGCAATGATTCTAAGTGTCATATGAAAGTTCCAATTGATATATATcaaacaa                    <  1:2376396/57‑1 (MQ=255)
                                                   taGAGAGCAATGATTCTAAGTGTCATATGAAAGTTCCAATTGATATATATcaaacaa                    <  1:2081013/57‑1 (MQ=255)
                                                   taGAGAGCAATGATTCTAAGTGTCATATGAAAGTTCCAATTGATATATATcaaacaa                    <  1:2413524/57‑1 (MQ=255)
                                                   taGAGAGCAATGATTCTAAGTGTCATATGAAAGTTCCAATTGATATATATcaaacaa                    <  1:3043937/57‑1 (MQ=255)
                                                     gagagCATTGATTCTAAGTGTCATATGAAAGTACCAATTGATATATATCAAACAAAATAACCCTGATTa      >  1:3036375/1‑69 (MQ=255)
                                                      agagCAATGATTCTAAGTGTCATATGAAAGTTCCAATTGATATATATcaaacaa                    <  1:312193/54‑1 (MQ=255)
                                                        agCATTGATTCTAAGTGTCATATGAAAGTACCAATTGATATATATCAAACAAAATAACCCTGATTAATg   <  1:664682/69‑1 (MQ=255)
                                                        agCATTGATTCTAAGTGTCATATGAAAGTACCAATTGATATATATCAAACAAAATAACCCTGATTAATg   <  1:365511/69‑1 (MQ=255)
                                                        agCATTGATTCTAAGTGTCATATGAAAGTACCAATTGATATATATCAAACAAAATAACCCTGATTAATg   <  1:1795010/69‑1 (MQ=255)
                                                        agCATTGATTCTAAGTGTCATATGAAAGTACCAATTGATATATATCAAACAAAATAACCCTGATTAATg   <  1:1499720/69‑1 (MQ=255)
                                                          cATTGATTCTAAGTGTCATATGAAAGTACCAATTGATATATATCAAACAAAATAACCCTGATTAATGa  >  1:3071271/1‑68 (MQ=255)
                                                          cATTGATTCTAAGTGTCATATGAAAGTACCAATTGATATATATCAAACAAAATAACCCTGATTAATGa  >  1:319754/1‑68 (MQ=255)
                                                          cATTGATTCTAAGTGTCATATGAAAGTACCAATTGATATATATCAAACAAAATAACCCTGATTAATGa  >  1:2889390/1‑68 (MQ=255)
                                                            |                                                                 
GTATTCAGACAGGCATCCCACCTGACTTCGAATGATGATTATTCATCACTATAGAGAGCATTGATTCTAAGTGTCATATGAAAGTACCAATTGATATATATCAAACAAAATAACCCTGATTAATGA  >  minE/1661662‑1661787

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