Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,014,008 T→C 54.2% intergenic (‑5/+2) exbD ← / ← exbB membrane spanning protein in TonB‑ExbB‑ExbD complex/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,0080TC54.2% 5.0 / 38.5 51intergenic (‑5/+2)exbD/exbBmembrane spanning protein in TonB‑ExbB‑ExbD complex/membrane spanning protein in TonB‑ExbB‑ExbD complex
Reads supporting (aligned to +/- strand):  ref base T (6/17);  new base C (18/9);  total (25/26)
Fisher's exact test for biased strand distribution p-value = 5.40e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

CACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGATGCGCAGCGGCGCTGGCTTCCAGATCCAG  >  minE/2013944‑2014075
                                                                |                                                                   
cACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta                                                                 <  1:2033038/69‑1 (MQ=255)
cACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta                                                                 <  1:1916164/69‑1 (MQ=255)
cACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta                                                                 <  1:1856769/69‑1 (MQ=255)
 aCGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtattat                                                                >  1:607003/1‑69 (MQ=255)
          tCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGt                                                       <  1:552137/69‑1 (MQ=255)
                   tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta                                                                 >  1:2973099/1‑50 (MQ=255)
                   tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta                                                                 >  1:3005809/1‑50 (MQ=255)
                   tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta                                                                 >  1:2285271/1‑50 (MQ=255)
                   tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta                                                                 >  1:2196410/1‑50 (MQ=255)
                   tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta                                                                 >  1:2185360/1‑50 (MQ=255)
                   tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta                                                                 >  1:3416119/1‑50 (MQ=255)
                   tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta                                                                 >  1:377021/1‑50 (MQ=255)
                   tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta                                                                 >  1:521292/1‑50 (MQ=255)
                   tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta                                                                 >  1:1645949/1‑50 (MQ=255)
                   tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta                                                                 >  1:877743/1‑50 (MQ=255)
                   tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta                                                                 >  1:577900/1‑50 (MQ=255)
                   tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta                                                                 >  1:1307452/1‑50 (MQ=255)
                   tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatt                                                                  >  1:3139784/1‑49 (MQ=255)
                              gtcgAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAAc                                   >  1:322360/1‑69 (MQ=255)
                              gtcgAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAAc                                   >  1:2319648/1‑69 (MQ=255)
                              gtcgAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAAc                                   >  1:1322779/1‑69 (MQ=255)
                                cgAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCg                                 <  1:2439881/69‑1 (MQ=255)
                                cgAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCg                                 <  1:3407145/69‑1 (MQ=255)
                                cgAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCg                                 <  1:2249768/69‑1 (MQ=255)
                                cgAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCg                                 <  1:692533/69‑1 (MQ=255)
                                cgAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCg                                 <  1:3423653/69‑1 (MQ=255)
                                 gAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCg                                 <  1:1876626/68‑1 (MQ=255)
                                  aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACg                                      >  1:1626561/1‑62 (MQ=255)
                                  aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                     >  1:3162929/1‑63 (MQ=255)
                                  aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                     >  1:3131535/1‑63 (MQ=255)
                                  aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                     >  1:2556039/1‑63 (MQ=255)
                                  aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                     >  1:2843410/1‑63 (MQ=255)
                                  aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                     >  1:2850078/1‑63 (MQ=255)
                                  aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                     >  1:2886301/1‑63 (MQ=255)
                                  aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                     >  1:1001794/1‑63 (MQ=255)
                                  aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                     >  1:2551921/1‑63 (MQ=255)
                                  aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                     >  1:1658140/1‑63 (MQ=255)
                                  aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                     >  1:234267/1‑63 (MQ=255)
                                  aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                     >  1:2236872/1‑63 (MQ=255)
                                  aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                     >  1:1341580/1‑63 (MQ=255)
                                  aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                     >  1:1347871/1‑63 (MQ=255)
                                  aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                     >  1:1436332/1‑63 (MQ=255)
                                  aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                     >  1:1750292/1‑63 (MQ=255)
                                  aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                     >  1:51551/1‑63 (MQ=255)
                                  aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa                                     >  1:1694793/1‑63 (MQ=255)
                                  aGGTTTTCGTTAAGATGCATCGCCATCGGAAATTATc                                                               >  1:3025965/1‑37 (MQ=37)
                                   ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa                               <  1:311711/68‑1 (MQ=255)
                                   ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa                               <  1:379593/68‑1 (MQ=255)
                                   ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa                               <  1:1140723/68‑1 (MQ=255)
                                   ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa                               <  1:1189063/68‑1 (MQ=255)
                                   ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa                               <  1:2679952/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:137561/68‑1 (MQ=255)
                                   ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa                               <  1:1662003/68‑1 (MQ=255)
                                       ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa                                                      <  1:2889443/41‑1 (MQ=25)
                                       ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa                                                      <  1:76530/41‑1 (MQ=25)
                                       ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa                                                      <  1:122149/41‑1 (MQ=25)
                                       ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa                                                      <  1:1563156/41‑1 (MQ=25)
                                       ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa                                                      <  1:467569/41‑1 (MQ=25)
                                       ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa                                                      <  1:2139418/41‑1 (MQ=25)
                                       ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa                                                      <  1:2143550/41‑1 (MQ=25)
                                       ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa                                                      <  1:243333/41‑1 (MQ=25)
                                       ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa                                                      <  1:2955377/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)
                                                               atatTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGATGCGCAGCGGCGCTGGCTTCCAGATCCAg  <  1:182134/69‑1 (MQ=255)
                                                               atatTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGATGCGCAGCGGCGCTGGCTTCCAGATCCAg  <  1:1266925/69‑1 (MQ=255)
                                                                tattatCCTGCGCGTAATTTTTGTGCGACACTAACCg                                 <  1:2591879/37‑1 (MQ=255)
                                                                |                                                                   
CACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGATGCGCAGCGGCGCTGGCTTCCAGATCCAG  >  minE/2013944‑2014075

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: