Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,013,998 T→C 51.8% A2A (GCA→GCG exbD ← 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,013,9980TC51.8% 3.2 / 72.3 83A2A (GCA→GCGexbDmembrane spanning protein in TonB‑ExbB‑ExbD complex
Reads supporting (aligned to +/- strand):  ref base T (20/20);  new base C (34/9);  total (54/29)
Fisher's exact test for biased strand distribution p-value = 6.46e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGATGCGCAGCGGCGCTG  >  minE/2013932‑2014061
                                                                  |                                                               
gATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATCGc                                                               >  1:1137346/1‑69 (MQ=255)
gATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATCGc                                                               >  1:1973233/1‑69 (MQ=255)
gATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATCGc                                                               >  1:1968940/1‑69 (MQ=255)
gATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATCGc                                                               >  1:2429462/1‑69 (MQ=255)
gATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATCGc                                                               >  1:2446569/1‑69 (MQ=255)
gATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATCGc                                                               >  1:2675969/1‑69 (MQ=255)
gATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATCGc                                                               >  1:1664826/1‑69 (MQ=255)
gATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATCGc                                                               >  1:3220841/1‑69 (MQ=255)
gATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATCGc                                                               >  1:1514358/1‑69 (MQ=255)
gATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATCGc                                                               >  1:517858/1‑69 (MQ=255)
gATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATCGc                                                               >  1:722471/1‑69 (MQ=255)
gATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATCGc                                                               >  1:117405/1‑69 (MQ=255)
gATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATCGc                                                               >  1:699122/1‑69 (MQ=255)
gATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATCGc                                                               >  1:594665/1‑69 (MQ=255)
gATAAACGGCGTCACGTTGATATAATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATCGc                                                               >  1:538436/1‑69 (MQ=255)
            cACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta                                                   <  1:1856769/69‑1 (MQ=255)
            cACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta                                                   <  1:1916164/69‑1 (MQ=255)
            cACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta                                                   <  1:2033038/69‑1 (MQ=255)
             aCGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtattat                                                  >  1:607003/1‑69 (MQ=255)
                      tCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGt                                         <  1:552137/69‑1 (MQ=255)
                               tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta                                                   >  1:577900/1‑50 (MQ=255)
                               tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta                                                   >  1:3005809/1‑50 (MQ=255)
                               tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta                                                   >  1:2973099/1‑50 (MQ=255)
                               tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta                                                   >  1:3416119/1‑50 (MQ=255)
                               tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta                                                   >  1:377021/1‑50 (MQ=255)
                               tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta                                                   >  1:2285271/1‑50 (MQ=255)
                               tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta                                                   >  1:2196410/1‑50 (MQ=255)
                               tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta                                                   >  1:521292/1‑50 (MQ=255)
                               tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta                                                   >  1:877743/1‑50 (MQ=255)
                               tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta                                                   >  1:2185360/1‑50 (MQ=255)
                               tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta                                                   >  1:1645949/1‑50 (MQ=255)
                               tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta                                                   >  1:1307452/1‑50 (MQ=255)
                               tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatt                                                    >  1:3139784/1‑49 (MQ=255)
                                          gtcgAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAAc                     >  1:2319648/1‑69 (MQ=255)
                                          gtcgAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAAc                     >  1:1322779/1‑69 (MQ=255)
                                          gtcgAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAAc                     >  1:322360/1‑69 (MQ=255)
                                            cgAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCg                   <  1:2249768/69‑1 (MQ=255)
                                            cgAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCg                   <  1:692533/69‑1 (MQ=255)
                                            cgAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCg                   <  1:2439881/69‑1 (MQ=255)
                                            cgAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCg                   <  1:3423653/69‑1 (MQ=255)
                                            cgAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCg                   <  1:3407145/69‑1 (MQ=255)
                                             gAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCg                   <  1:1876626/68‑1 (MQ=255)
                                              aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACg                        >  1:1626561/1‑62 (MQ=255)
                                              aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                       >  1:2850078/1‑63 (MQ=255)
                                              aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                       >  1:1347871/1‑63 (MQ=255)
                                              aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                       >  1:2886301/1‑63 (MQ=255)
                                              aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                       >  1:1658140/1‑63 (MQ=255)
                                              aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                       >  1:3131535/1‑63 (MQ=255)
                                              aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                       >  1:3162929/1‑63 (MQ=255)
                                              aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                       >  1:2843410/1‑63 (MQ=255)
                                              aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                       >  1:51551/1‑63 (MQ=255)
                                              aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                       >  1:1694793/1‑63 (MQ=255)
                                              aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                       >  1:1436332/1‑63 (MQ=255)
                                              aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                       >  1:2556039/1‑63 (MQ=255)
                                              aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                       >  1:2551921/1‑63 (MQ=255)
                                              aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                       >  1:1341580/1‑63 (MQ=255)
                                              aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                       >  1:234267/1‑63 (MQ=255)
                                              aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                       >  1:1750292/1‑63 (MQ=255)
                                              aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                       >  1:2236872/1‑63 (MQ=255)
                                              aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                       >  1:1001794/1‑63 (MQ=255)
                                              aGGTTTTCGTTAAGATGCATCGCCATCGGAAATTATc                                                 >  1:3025965/1‑37 (MQ=37)
                                               ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa                 <  1:311711/68‑1 (MQ=255)
                                               ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa                 <  1:1662003/68‑1 (MQ=255)
                                               ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa                 <  1:1140723/68‑1 (MQ=255)
                                               ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa                 <  1:1189063/68‑1 (MQ=255)
                                               ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa                 <  1:379593/68‑1 (MQ=255)
                                               ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa                 <  1:228451/68‑1 (MQ=255)
                                               ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa                 <  1:2679952/68‑1 (MQ=255)
                                               ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa                 <  1:2587099/68‑1 (MQ=255)
                                               ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa                 <  1:137561/68‑1 (MQ=255)
                                               ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa                 <  1:2562140/68‑1 (MQ=255)
                                                   ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa                                        <  1:2139418/41‑1 (MQ=25)
                                                   ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa                                        <  1:76530/41‑1 (MQ=25)
                                                   ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa                                        <  1:122149/41‑1 (MQ=25)
                                                   ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa                                        <  1:243333/41‑1 (MQ=25)
                                                   ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa                                        <  1:467569/41‑1 (MQ=25)
                                                   ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa                                        <  1:1563156/41‑1 (MQ=25)
                                                   ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa                                        <  1:2889443/41‑1 (MQ=25)
                                                   ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa                                        <  1:2955377/41‑1 (MQ=25)
                                                   ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa                                        <  1:2143550/41‑1 (MQ=25)
                                                             tGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGATGCGCAGCGGCGCTg  >  1:1731907/1‑69 (MQ=255)
                                                             tGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGATGCGCAGCGGCGCTg  >  1:1696644/1‑69 (MQ=255)
                                                             tGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGATGCGCAGCGGCGCTg  >  1:796282/1‑69 (MQ=255)
                                                                  |                                                               
GATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGATGCGCAGCGGCGCTG  >  minE/2013932‑2014061

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: