Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 676,582 G→A 43.0% G57G (GGC→GGT ompA ← outer membrane protein A (3a;II*;G;d)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE676,5820GA43.0% 27.5 / 63.4 86G57G (GGC→GGTompAouter membrane protein A (3a;II*;G;d)
Reads supporting (aligned to +/- strand):  ref base G (14/35);  new base A (15/22);  total (29/57)
Fisher's exact test for biased strand distribution p-value = 2.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

TAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAG  >  minE/676516‑676648
                                                                  |                                                                   
tAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAgcg                                                                     >  1:3340944/1‑67 (MQ=255)
    cAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAg                                                                >  1:2034601/1‑68 (MQ=255)
    cAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGt                                                               >  1:749057/1‑69 (MQ=255)
    cAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGt                                                               >  1:1769811/1‑69 (MQ=255)
    cAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGt                                                               >  1:220052/1‑69 (MQ=255)
    cAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGt                                                               >  1:1790136/1‑69 (MQ=255)
    cAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGt                                                               >  1:809132/1‑69 (MQ=255)
    cAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGt                                                               >  1:3231461/1‑69 (MQ=255)
    cAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGt                                                               >  1:1403015/1‑69 (MQ=255)
    cAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGt                                                               >  1:1310969/1‑69 (MQ=255)
    cAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGt                                                               >  1:1302960/1‑69 (MQ=255)
    cAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGt                                                               >  1:2932948/1‑69 (MQ=255)
    cAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGt                                                               >  1:2841507/1‑69 (MQ=255)
    cAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGt                                                               >  1:2907620/1‑69 (MQ=255)
    cAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGt                                                               >  1:1100410/1‑69 (MQ=255)
    cAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGt                                                               >  1:2896776/1‑69 (MQ=255)
        cGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGt                                                          >  1:2035877/1‑69 (MQ=255)
        cGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGt                                                          >  1:846972/1‑69 (MQ=255)
           aaCCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGtttt                                                       <  1:1921744/69‑1 (MQ=255)
           aaCCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGtttt                                                       <  1:1851496/69‑1 (MQ=255)
           aaCCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGtttt                                                       <  1:1679099/69‑1 (MQ=255)
           aaCCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGtttt                                                       <  1:79619/69‑1 (MQ=255)
            aCCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGtttt                                                       <  1:2094238/68‑1 (MQ=255)
            aCCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGtttt                                                       <  1:1908065/68‑1 (MQ=255)
            aCCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTc                                                      >  1:3236695/1‑69 (MQ=255)
              ccATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGTTGGTTTTcgt                                                    <  1:3013130/69‑3 (MQ=255)
              ccATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGTTGGTTTTcgt                                                    <  1:2117824/69‑3 (MQ=255)
              ccATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGTTGGTTTTcgt                                                    <  1:3005944/69‑3 (MQ=255)
              ccATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGTTGGTTTTcgt                                                    <  1:3048548/69‑3 (MQ=255)
              ccATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGTTGGTTTTcgt                                                    <  1:765378/69‑3 (MQ=255)
              ccATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGTTGGTTTTcgt                                                    <  1:1888692/69‑3 (MQ=255)
              ccATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGTTGGTTTTcgt                                                    <  1:1899187/69‑3 (MQ=255)
              ccATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGTTGGTTTTcgt                                                    <  1:1154126/69‑3 (MQ=255)
              ccATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGTTGGTTTTcgt                                                    <  1:698466/69‑3 (MQ=255)
              ccATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGTTGGTTTTcgt                                                    <  1:1827175/69‑3 (MQ=255)
              ccATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGTTGGTTTTcgt                                                    <  1:69934/69‑3 (MQ=255)
              ccATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGTTGGTTTTcgt                                                    <  1:1719681/69‑3 (MQ=255)
              ccATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGTTGGTTTTcgt                                                    <  1:1387688/69‑3 (MQ=255)
              ccATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGTTGGTTTTcgt                                                    <  1:1395701/69‑3 (MQ=255)
               cATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGTTGGTTTTcgt                                                    <  1:1690918/68‑3 (MQ=255)
               cATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGTTGGTTTTcgt                                                    <  1:1165146/68‑3 (MQ=255)
               cATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGTTGGTTTTcgt                                                    <  1:3360547/68‑3 (MQ=255)
               cATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGTTGGTTTTcgt                                                    <  1:3061051/68‑3 (MQ=255)
               cATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGTTGGTTTTcgt                                                    <  1:1866287/68‑3 (MQ=255)
               cATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCTGCACCCAGTTGGTTTTcgt                                                    <  1:3065040/68‑3 (MQ=255)
                   tCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTc                                               >  1:1487613/1‑69 (MQ=255)
                   tCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTc                                               >  1:2313360/1‑69 (MQ=255)
                   tCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTc                                               >  1:303273/1‑69 (MQ=255)
                   tCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTc                                               >  1:383865/1‑69 (MQ=255)
                   tCAAAGACAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTc                                               >  1:4395/1‑69 (MQ=255)
                                 cGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCAttgttgtt                                 <  1:1032542/69‑1 (MQ=255)
                                 cGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCAttgttgtt                                 >  1:515383/1‑69 (MQ=255)
                                 cGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCAttgttgtt                                 >  1:952881/1‑69 (MQ=255)
                                 cGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCAttgttgtt                                 >  1:1397097/1‑69 (MQ=255)
                                 cGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCAttgttgtt                                 <  1:2996575/69‑1 (MQ=255)
                                 cGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCAttgttgtt                                 <  1:3192789/69‑1 (MQ=255)
                                 cGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCAttgttgtt                                 >  1:3116482/1‑69 (MQ=255)
                                 cGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCAttgttgtt                                 <  1:1898340/69‑1 (MQ=255)
                                 cGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATAGGTCGGGCCAttgttgtt                                 <  1:2507194/69‑1 (MQ=255)
                                      tAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTtgatga                            <  1:2862663/69‑1 (MQ=255)
                                      tAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTtgatga                            <  1:1793834/69‑1 (MQ=255)
                                      tAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTtgatga                            <  1:2003673/69‑1 (MQ=255)
                                        aCCTGGTAACCACCAAAAGCACCTGCACCCAGTTGGTTTTcgt                                                    <  1:2916479/43‑3 (MQ=255)
                                         ccTGGTAACCACCAAAAGCACCTGCACCCAGTTGGTTTTcgt                                                    <  1:2563118/42‑3 (MQ=255)
                                          cTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAAcc                        <  1:3087038/69‑1 (MQ=255)
                                          cTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAAcc                        <  1:2866599/69‑1 (MQ=255)
                                          cTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAAcc                        <  1:103237/69‑1 (MQ=255)
                                          cTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAAcc                        <  1:116053/69‑1 (MQ=255)
                                          cTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGGCGGGCCATTGTTGTTGATGAAAcc                        <  1:1207649/69‑1 (MQ=255)
                                           tGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTTCATGGGTCGGGCCATTGTTGTTGATGAAAcc                        >  1:2213659/1‑69 (MQ=255)
                                           tGGTAACCACCAAAAGCACAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAg                      <  1:1904658/69‑1 (MQ=255)
                                            ggTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAg                      <  1:3268724/69‑1 (MQ=255)
                                            ggTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAg                      <  1:3433079/69‑1 (MQ=255)
                                            ggTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAg                      <  1:1407387/69‑1 (MQ=255)
                                            ggTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAg                      <  1:3222272/69‑1 (MQ=255)
                                            ggTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAg                      <  1:2145526/69‑1 (MQ=255)
                                            ggTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGACATTGTTGTTGATGAAACCAg                      <  1:372351/69‑1 (MQ=255)
                                            ggTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGGCGGGCCATTGTTGTTGATGAAACCAg                      <  1:2513864/69‑1 (MQ=255)
                                            ggTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGGCGGGCCATTGTTGTTGATGAAACCAg                      <  1:1929973/69‑1 (MQ=255)
                                               aaCCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACTAgtgt                   <  1:2483159/69‑1 (MQ=255)
                                               aaCCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAgtgt                   <  1:2102136/69‑1 (MQ=255)
                                                  caccaAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCa                 >  1:2936598/1‑68 (MQ=255)
                                                  caccaAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCAt                <  1:117931/69‑1 (MQ=255)
                                                  caccaAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCAt                <  1:975944/69‑1 (MQ=255)
                                                    ccaAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATgg              <  1:3431460/69‑1 (MQ=255)
                                                    ccaAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATgg              <  1:3389133/69‑1 (MQ=255)
                                                    ccaAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATgg              <  1:3161152/69‑1 (MQ=255)
                                                                  gcCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAg  >  1:2702033/1‑67 (MQ=255)
                                                                  |                                                                   
TAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAG  >  minE/676516‑676648

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: