Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 481,958 C→T 40.6% N117N (AAC→AAT 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,9580CT40.6% 37.5 / 59.3 86N117N (AAC→AATtolQmembrane spanning protein in TolA‑TolQ‑TolR complex
Reads supporting (aligned to +/- strand):  ref base C (30/21);  new base T (24/11);  total (54/32)
Fisher's exact test for biased strand distribution p-value = 3.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATA  >  minE/481905‑482025
                                                     |                                                                   
gcgcCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCGATGCGTATTTCCATGAATCGTGAACTTGaaaa                                                       >  1:1002907/1‑68 (MQ=255)
gcgcCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCGATGCGTATTTCCATGAATCGTGAACTTGaaaa                                                       >  1:2024900/1‑68 (MQ=255)
gcgcCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCGATGCGTATTTCCATGAATCGTGAACTTGaaaa                                                       >  1:2500003/1‑68 (MQ=255)
gcgcCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCGATGCGTATTTCCATGAATCGTGAACTTGAAAAt                                                      >  1:1568615/1‑69 (MQ=255)
gcgcCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCGATGCGTATTTCCATGAATCGTGAACTTGAAAAt                                                      >  1:1021750/1‑69 (MQ=255)
gcgcCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCGATGCGTATTTCCATGAATCGTGAACTTGAAAAt                                                      >  1:1610555/1‑69 (MQ=255)
gcgcCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCGATGCGTATTTCCATGAATCGTGAACTTGAAAAt                                                      >  1:1727044/1‑69 (MQ=255)
gcgcCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCGATGCGTATTTCCATGAATCGTGAACTTGAAAAt                                                      >  1:1924794/1‑69 (MQ=255)
gcgcCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCGATGCGTATTTCCATGAATCGTGAACTTGAAAAt                                                      >  1:1936277/1‑69 (MQ=255)
gcgcCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCGATGCGTATTTCCATGAATCGTGAACTTGAAAAt                                                      >  1:1947945/1‑69 (MQ=255)
gcgcCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCGATGCGTATTTCCATGAATCGTGAACTTGAAAAt                                                      >  1:2058557/1‑69 (MQ=255)
gcgcCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCGATGCGTATTTCCATGAATCGTGAACTTGAAAAt                                                      >  1:2088650/1‑69 (MQ=255)
gcgcCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCGATGCGTATTTCCATGAATCGTGAACTTGAAAAt                                                      >  1:227788/1‑69 (MQ=255)
gcgcCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCGATGCGTATTTCCATGAATCGTGAACTTGAAAAt                                                      >  1:2890701/1‑69 (MQ=255)
gcgcCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCGATGCGTATTTCCATGAATCGTGAACTTGAAAAt                                                      >  1:3177682/1‑69 (MQ=255)
gcgcCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCGATGCGTATTTCCATGAATCGTGAACTTGAAAAt                                                      >  1:451862/1‑69 (MQ=255)
gcgcCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCGATGCGTATTTCCATGAATCGTGAACTTGAAAAt                                                      >  1:58639/1‑69 (MQ=255)
gcgcCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCGATGCGTATTTCCATGAATCGTGAACTTGAAAAt                                                      >  1:696022/1‑69 (MQ=255)
gcgcCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCGATGCGTATTTCCATGAATCGTGAACTTGAAAAt                                                      >  1:985731/1‑69 (MQ=255)
gcgcCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCGATGCGTATTTCCATGAATCGTGAACTTGAAAAt                                                      >  1:993839/1‑69 (MQ=255)
gcgcCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCGATGCGTATTTCCATGAATCGTGAACTTGAAAAt                                                      >  1:1432045/1‑69 (MQ=255)
gcgcCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCGATGCGTATTTCCATGAATCGTGAACTTGAAAAt                                                      >  1:1382968/1‑69 (MQ=255)
gcgcCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCGATGCGTATTTCCATGAATCGTGAACTTGAAAAt                                                      >  1:133292/1‑69 (MQ=255)
gcgcCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCGATGCGTATTTCCATGAATCGTGAACTTGAAAAt                                                      >  1:1229025/1‑69 (MQ=255)
               gTGGAAGGGGCGTCGCGTGCTATTCGTATCTCCATGAACCGTg                                                                 <  1:1457135/43‑1 (MQ=255)
               gTGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTg                                                                 <  1:2738190/43‑1 (MQ=255)
               gTGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTg                                                                 <  1:495698/43‑1 (MQ=255)
               gTGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTg                                                                 <  1:705888/43‑1 (MQ=255)
               gTGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTg                                                                 <  1:1256874/43‑1 (MQ=255)
               gTGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTg                                                                 <  1:2727339/43‑1 (MQ=255)
                             gcgTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACg                                             >  1:2855434/1‑49 (MQ=255)
                             gcgTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCg                                    >  1:2249708/1‑58 (MQ=255)
                             gcgTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTa                          >  1:1456251/1‑68 (MQ=255)
                             gcgTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTAc                         <  1:2457248/69‑1 (MQ=255)
                             gcgTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTAc                         <  1:2900488/69‑1 (MQ=255)
                             gcgTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTAc                         >  1:2851703/1‑69 (MQ=255)
                             gcgTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTAc                         <  1:858241/69‑1 (MQ=255)
                             gcgTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTAc                         <  1:1608991/69‑1 (MQ=255)
                             gcgTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTAc                         <  1:590747/69‑1 (MQ=255)
                                tGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGt                      <  1:399559/69‑1 (MQ=255)
                                tGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGt                      <  1:1559270/69‑1 (MQ=255)
                                tGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGt                      <  1:2610797/69‑1 (MQ=255)
                                tGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGt                      <  1:1094342/69‑1 (MQ=255)
                                tGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGt                      <  1:2090171/69‑1 (MQ=255)
                                    aTGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAAc                                              >  1:3148689/1‑41 (MQ=255)
                                    aTGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACg                                             >  1:523652/1‑42 (MQ=255)
                                    aTGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACg                                             >  1:3152395/1‑42 (MQ=255)
                                    aTGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGc                  <  1:2524585/69‑1 (MQ=255)
                                    aTGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGc                  <  1:2068815/69‑1 (MQ=255)
                                    aTGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGc                  <  1:3389927/69‑1 (MQ=255)
                                       cGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCTTTCCTCGGTACGGTTGGCta                >  1:845889/1‑67 (MQ=255)
                                       cGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTcc               >  1:212531/1‑69 (MQ=255)
                                       cGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTcc               >  1:1890690/1‑69 (MQ=255)
                                       cGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTcc               >  1:712580/1‑69 (MQ=255)
                                       cGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTc                >  1:795096/1‑68 (MQ=255)
                                          atttcCATGAATCGTGAACTTGAAAATCTGGAAACGCACATTCCTTTCCTCGGTACGGTTGGCTCCATc            <  1:502476/66‑1 (MQ=255)
                                          atttcCATGAATCGTGAACTTGAAAATCTGGAAACGCACATTCCTTTCCTCGGTACGGTTGGCTCCATc            <  1:3192139/66‑1 (MQ=255)
                                          atttcCATGAATCGTGAACTTGAAAATCTGGAAACGCACATTCCTTTCCTCGGTACGGTTGGCTCCATc            <  1:1223947/66‑1 (MQ=255)
                                          atttcCATGAATCGTGAACTTGAAAATCTGGAAACGCACATTCCTTTCCTCGGTACGGTTGGCTCCATc            <  1:1561793/66‑1 (MQ=255)
                                          atttcCATGAATCGTGAACTTGAAAATCTGGAAACGCACATTCCTTTCCTCGGTACGGTTGGCTCCATc            <  1:2629577/66‑1 (MQ=255)
                                          atttcCATGAATCGTGAACTTGAAAATCTGGAAACGCACATTCCTTTCCTCGGTACGGTTGGCTCCATc            <  1:259746/66‑1 (MQ=255)
                                          atttcCATGAATCGTGAACTTGAAAATCTGGAAACGCACATTCCTTTCCTCGGTACGGTTGGCTCCATc            <  1:2573301/66‑1 (MQ=255)
                                          atttcCATGAATCGTGAACTTGAAAATCTGGAAACGCACATTCCTTTCCTCGGTACGGTTGGCTCCATc            <  1:1966781/66‑1 (MQ=255)
                                          atttcCATGAATCGTGAACTTGAAAATCTGGAAACGCACATTCCTTTCCTCGGTACGGTTGGCTCCATc            <  1:1305488/66‑1 (MQ=255)
                                           tttcCATGAATCGTGAACTTGAAAATCTGGAAACGCACATTCCTTTCCTCGGTACGGTTGGCTCCATc            <  1:1703836/66‑1 (MQ=255)
                                           tttcCATGAATCGTGAACTTGAAAATCTGGAAACGCACATTCCTTTCCTCGGTACGGTTGGCTCCATc            <  1:2086056/66‑1 (MQ=255)
                                           tctcCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATc            >  1:619324/1‑68 (MQ=255)
                                           tctcCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCa           >  1:1540489/1‑69 (MQ=255)
                                           tctcCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCa           >  1:107776/1‑69 (MQ=255)
                                             tcCATGAACCGTGAACTTGAAAATCTGGAAACGca                                           <  1:2851436/35‑1 (MQ=255)
                                             tcCATGAACCGTGAACTTGAAAATCTGGAAACGca                                           <  1:2303358/35‑1 (MQ=255)
                                               cATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGtt                                  >  1:3233938/1‑42 (MQ=255)
                                               cATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCgg                            >  1:3144061/1‑48 (MQ=255)
                                               cATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCgg                            >  1:961864/1‑48 (MQ=255)
                                               cATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCt                 >  1:2504075/1‑59 (MQ=255)
                                               cATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTc                >  1:1208275/1‑60 (MQ=255)
                                               cATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTc                >  1:1956924/1‑60 (MQ=255)
                                               cATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTc                >  1:2546931/1‑60 (MQ=255)
                                                 tGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCg      >  1:2944583/1‑68 (MQ=255)
                                                 tGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCg      >  1:287538/1‑68 (MQ=255)
                                                 tGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCg      >  1:2782687/1‑68 (MQ=255)
                                                 tGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCg      >  1:2018626/1‑68 (MQ=255)
                                                 tGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCg      >  1:2045241/1‑68 (MQ=255)
                                                 tGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCg      >  1:3343250/1‑68 (MQ=255)
                                                   aaCCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCg      >  1:3329909/1‑66 (MQ=255)
                                                   aaCCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGtat   >  1:1781463/1‑69 (MQ=255)
                                                     ccGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGtata  <  1:1567651/68‑1 (MQ=255)
                                                     ccGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGtata  <  1:1753870/68‑1 (MQ=255)
                                                     ccGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGtata  <  1:1040857/68‑1 (MQ=255)
                                                     ccGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGtata  <  1:1827201/68‑1 (MQ=255)
                                                     ccGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGtata  <  1:1916029/68‑1 (MQ=255)
                                                     ccGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGtata  <  1:1474188/68‑1 (MQ=255)
                                                     ccGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGtata  <  1:429944/68‑1 (MQ=255)
                                                     ccGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGtata  <  1:362538/68‑1 (MQ=255)
                                                     ccGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGtata  <  1:1397792/68‑1 (MQ=255)
                                                     ccGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGtata  <  1:3190534/68‑1 (MQ=255)
                                                     ccGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGtata  <  1:2877557/68‑1 (MQ=255)
                                                     |                                                                   
GCGCCGGAAGCCGTAGTGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATA  >  minE/481905‑482025

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: