Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 481,654 T→A 7.7% L16Q (CTG→CAG)  tolQ → membrane spanning protein in TolA‑TolQ‑TolR complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE481,6540TA 7.7% 144.3 / 4.1 52L16Q (CTG→CAG) tolQmembrane spanning protein in TolA‑TolQ‑TolR complex
Reads supporting (aligned to +/- strand):  ref base T (26/22);  new base A (2/2);  total (28/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

TATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCTTGGGCCATTATTAT  >  minE/481587‑481714
                                                                   |                                                            
tATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTgg                                                            >  1:407996/1‑70 (MQ=255)
tATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTgg                                                            >  1:456836/1‑70 (MQ=255)
  ttGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTa                                                         <  1:1365291/71‑1 (MQ=255)
  ttGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCAgg                                                            <  1:2036339/68‑1 (MQ=255)
  ttGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCAGAAGGCTAGCCTTCAggaac                                                         <  1:1870057/71‑4 (MQ=38)
     tCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAAc                                                      <  1:3599697/71‑1 (MQ=255)
       gcgGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACtt                                                    <  1:1164664/71‑1 (MQ=255)
       gcgGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACt                                                     <  1:1353986/70‑1 (MQ=255)
          gAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTAt                                                  >  1:2243492/1‑70 (MQ=255)
               tAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGtt                                            >  1:2141630/1‑71 (MQ=255)
               tAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGt                                             <  1:2637487/70‑1 (MQ=255)
                 aGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGttga                                          >  1:2963653/1‑71 (MQ=255)
                 aGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGttga                                          >  1:1536707/1‑71 (MQ=255)
                 aGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCAgg                                                            >  1:3318943/1‑53 (MQ=255)
                 aGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCAgg                                                            >  1:3232713/1‑53 (MQ=255)
                   cAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGttgatt                                        >  1:275294/1‑71 (MQ=255)
                    aGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGttgattt                                       <  1:477282/71‑1 (MQ=255)
                     gTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGttgatttt                                      >  1:3906445/1‑71 (MQ=255)
                      tgactgacATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGttgattttg                                     <  1:3584850/71‑1 (MQ=255)
                        actgacATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGttgattttgat                                   >  1:2577582/1‑71 (MQ=255)
                        actgacATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGttgattttga                                    <  1:2379655/70‑1 (MQ=255)
                          tgacATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTaaa                                                       <  1:81429/49‑1 (MQ=255)
                            acATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTaa                                                        <  1:238127/46‑1 (MQ=255)
                            acATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGt                               <  1:1656096/71‑1 (MQ=255)
                              aTGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGttt                             >  1:2418515/1‑71 (MQ=255)
                              aTGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGttt                             >  1:3476690/1‑71 (MQ=255)
                                 aaTATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTTGTTAAACTTATCATGttgattttgat                                   <  1:2699102/62‑1 (MQ=255)
                                 aaTATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGttgatt                                        >  1:1814924/1‑57 (MQ=255)
                                   tatCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTAt                                                  <  1:2747008/45‑1 (MQ=255)
                                   tatCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTAt                                                  <  1:3160771/45‑1 (MQ=255)
                                   tatCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCaa                        <  1:588917/71‑1 (MQ=255)
                                     tCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATc                      >  1:2840926/1‑71 (MQ=255)
                                       cTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCg                     >  1:2500728/1‑70 (MQ=255)
                                          gATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGttttt                           >  1:3036915/1‑61 (MQ=255)
                                          gATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATc                 >  1:932068/1‑71 (MQ=255)
                                          gATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATc                 >  1:3575667/1‑71 (MQ=255)
                                          gATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATc                 >  1:2544461/1‑71 (MQ=255)
                                          gATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATc                 >  1:1175781/1‑71 (MQ=255)
                                             ttGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATc                                                 >  1:1671283/1‑36 (MQ=255)
                                             ttGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGttgatt                                        <  1:3877902/45‑1 (MQ=255)
                                             ttGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGttgatt                                        <  1:1853649/45‑1 (MQ=255)
                                              tGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCTTg              <  1:3296298/70‑1 (MQ=255)
                                                 tCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCTTGGGc           >  1:3490098/1‑70 (MQ=255)
                                                 tCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCTTGGGc           >  1:167930/1‑70 (MQ=255)
                                                  ccTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCTTGGGCCa         >  1:195472/1‑71 (MQ=255)
                                                   cTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGttttt                           <  1:1550926/52‑1 (MQ=255)
                                                        ggCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATc                 <  1:3519691/57‑1 (MQ=255)
                                                          cTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCaa                        <  1:878704/48‑1 (MQ=255)
                                                          cTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCTTGGGCCat        <  1:2609113/64‑1 (MQ=255)
                                                           tAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATc                      >  1:3893511/1‑49 (MQ=255)
                                                           tAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCTTGGGCCATtattat  >  1:662787/1‑69 (MQ=255)
                                                             gCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCa                   >  1:3801503/1‑50 (MQ=255)
                                                                   |                                                            
TATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCTTGGGCCATTATTAT  >  minE/481587‑481714

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