Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,014,019 G→A 48.2% R242R (CGC→CGT exbB ← membrane spanning protein in TonB‑ExbB‑ExbD complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,014,0190GA48.2% 1.1 / 47.0 56R242R (CGC→CGTexbBmembrane spanning protein in TonB‑ExbB‑ExbD complex
Reads supporting (aligned to +/- strand):  ref base G (9/20);  new base A (18/9);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 1.53e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

TCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGATGCGCAGCGGCGCTGGCTTCCAGATCCAGGTCAC  >  minE/2013954‑2014080
                                                                 |                                                             
tCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGt                                                            <  1:552137/69‑1 (MQ=255)
                    gtcgAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAAc                                        >  1:322360/1‑69 (MQ=255)
                    gtcgAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAAc                                        >  1:1322779/1‑69 (MQ=255)
                    gtcgAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAAc                                        >  1:2319648/1‑69 (MQ=255)
                      cgAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCg                                      <  1:692533/69‑1 (MQ=255)
                      cgAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCg                                      <  1:3407145/69‑1 (MQ=255)
                      cgAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCg                                      <  1:2439881/69‑1 (MQ=255)
                      cgAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCg                                      <  1:3423653/69‑1 (MQ=255)
                      cgAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCg                                      <  1:2249768/69‑1 (MQ=255)
                       gAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCg                                      <  1:1876626/68‑1 (MQ=255)
                        aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACg                                           >  1:1626561/1‑62 (MQ=255)
                        aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                          >  1:234267/1‑63 (MQ=255)
                        aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                          >  1:2236872/1‑63 (MQ=255)
                        aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                          >  1:2551921/1‑63 (MQ=255)
                        aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                          >  1:2843410/1‑63 (MQ=255)
                        aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                          >  1:2556039/1‑63 (MQ=255)
                        aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                          >  1:2850078/1‑63 (MQ=255)
                        aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                          >  1:2886301/1‑63 (MQ=255)
                        aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                          >  1:3162929/1‑63 (MQ=255)
                        aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                          >  1:1750292/1‑63 (MQ=255)
                        aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                          >  1:3131535/1‑63 (MQ=255)
                        aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                          >  1:1694793/1‑63 (MQ=255)
                        aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                          >  1:1658140/1‑63 (MQ=255)
                        aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                          >  1:1436332/1‑63 (MQ=255)
                        aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                          >  1:1347871/1‑63 (MQ=255)
                        aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                          >  1:1341580/1‑63 (MQ=255)
                        aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                          >  1:51551/1‑63 (MQ=255)
                        aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                          >  1:1001794/1‑63 (MQ=255)
                         ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa                                    <  1:311711/68‑1 (MQ=255)
                         ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa                                    <  1:1140723/68‑1 (MQ=255)
                         ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa                                    <  1:2679952/68‑1 (MQ=255)
                         ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa                                    <  1:379593/68‑1 (MQ=255)
                         ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa                                    <  1:2587099/68‑1 (MQ=255)
                         ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa                                    <  1:2562140/68‑1 (MQ=255)
                         ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa                                    <  1:228451/68‑1 (MQ=255)
                         ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa                                    <  1:1189063/68‑1 (MQ=255)
                         ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa                                    <  1:137561/68‑1 (MQ=255)
                         ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa                                    <  1:1662003/68‑1 (MQ=255)
                             ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa                                                           <  1:243333/41‑1 (MQ=25)
                             ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa                                                           <  1:76530/41‑1 (MQ=25)
                             ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa                                                           <  1:2143550/41‑1 (MQ=25)
                             ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa                                                           <  1:2139418/41‑1 (MQ=25)
                             ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa                                                           <  1:2889443/41‑1 (MQ=25)
                             ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa                                                           <  1:2955377/41‑1 (MQ=25)
                             ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa                                                           <  1:122149/41‑1 (MQ=25)
                             ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa                                                           <  1:467569/41‑1 (MQ=25)
                             ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa                                                           <  1:1563156/41‑1 (MQ=25)
                                       tGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGATGCGCAGCGGCGCTg                     >  1:796282/1‑69 (MQ=255)
                                       tGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGATGCGCAGCGGCGCTg                     >  1:1696644/1‑69 (MQ=255)
                                       tGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGATGCGCAGCGGCGCTg                     >  1:1731907/1‑69 (MQ=255)
                                                     atatTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGATGCGCAGCGGCGCTGGCTTCCAGATCCAg       <  1:1266925/69‑1 (MQ=255)
                                                     atatTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGATGCGCAGCGGCGCTGGCTTCCAGATCCAg       <  1:182134/69‑1 (MQ=255)
                                                      tattatCCTGCGCGTAATTTTTGTGCGACACTAACCg                                      <  1:2591879/37‑1 (MQ=255)
                                                        ttatCCTGCGCGTAATTTTTGTGCGACACGAACCGGATGCGCAGCGGCGCTGGCTTCCAGATCCAGGTc    >  1:1703494/1‑69 (MQ=255)
                                                        ttatCCTGCGCGTAATTTTTGTGCGACACGAACCGGATGCGCAGCGGCGCTGGCTTCCAGATCCAGGTc    >  1:2079072/1‑69 (MQ=255)
                                                          atCCTGCGCGTAATTTTTGTGCGACACGAACCGGATGCGCAGCGGCGCTGGCTTCCAGATCCAGGTCAc  >  1:2591540/1‑69 (MQ=255)
                                                                 |                                                             
TCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGATGCGCAGCGGCGCTGGCTTCCAGATCCAGGTCAC  >  minE/2013954‑2014080

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: