Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,244,760 A→T 41.8% A264A (GCT→GCA yebA ← predicted peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,244,7600AT41.8% 74.6 / 156.2 153A264A (GCT→GCAyebApredicted peptidase
Reads supporting (aligned to +/- strand):  ref base A (49/40);  new base T (15/49);  total (64/89)
Fisher's exact test for biased strand distribution p-value = 1.18e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

Only 100 of 153 total aligned reads displayed.

GCAAGAATCCTTTCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGTTTAC  >  minE/1244696‑1244824
                                                                |                                                                
gCAAGAATCCTTTCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGt                                                              >  1:559941/1‑69 (MQ=255)
gCAAGAATCCTTTCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGt                                                              >  1:730404/1‑69 (MQ=255)
      aTCCTTTCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTaataat                                                        >  1:2860880/1‑69 (MQ=255)
      aTCCTTTCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTaataat                                                        >  1:2093460/1‑69 (MQ=255)
      aTCCTTTCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTaataat                                                        >  1:909918/1‑69 (MQ=255)
      aTCCTTTCGCCAGACAAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTaataat                                                        >  1:1034055/1‑69 (MQ=255)
        ccTTTCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATc                                                       >  1:1867327/1‑68 (MQ=255)
        ccTTTCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATc                                                       >  1:666401/1‑68 (MQ=255)
        ccTTTCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATc                                                       >  1:2432685/1‑68 (MQ=255)
         cTTTCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCtt                                                     <  1:583670/69‑1 (MQ=255)
         cTTTCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCtt                                                     <  1:1042544/69‑1 (MQ=255)
         cTTTCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCtt                                                     <  1:1038660/69‑1 (MQ=255)
         cTTTCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCtt                                                     <  1:330239/69‑1 (MQ=255)
            tCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTa                                                   >  1:503704/1‑68 (MQ=255)
            tCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTAc                                                  >  1:1948166/1‑69 (MQ=255)
            tCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTAc                                                  >  1:2078260/1‑69 (MQ=255)
            tCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTAc                                                  >  1:1673656/1‑69 (MQ=255)
            tCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTAc                                                  >  1:1657130/1‑69 (MQ=255)
            tCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTAc                                                  >  1:1734489/1‑69 (MQ=255)
            tCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTAc                                                  >  1:584212/1‑69 (MQ=255)
            tCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTAc                                                  >  1:689093/1‑69 (MQ=255)
            tCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTAc                                                  >  1:753058/1‑69 (MQ=255)
            tCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGAAATAATCTTTAc                                                  >  1:1448665/1‑69 (MQ=255)
                    ccAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTACCTTCTg                                            <  1:2825418/67‑1 (MQ=255)
                    ccAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTACCTTCTg                                            <  1:2466833/67‑1 (MQ=255)
                    ccAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTACCTTCTg                                            <  1:3239043/67‑1 (MQ=255)
                    ccAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTACCTTCTg                                            <  1:2269874/67‑1 (MQ=255)
                    ccAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTACCTTCTg                                            <  1:3336014/67‑1 (MQ=255)
                    ccAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTACCTTCTg                                            <  1:3402083/67‑1 (MQ=255)
                    ccAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTACCTTCTg                                            <  1:1966162/67‑1 (MQ=255)
                    ccAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTACCTTCTg                                            <  1:940434/67‑1 (MQ=255)
                    ccAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTACCTTCTg                                            <  1:1636168/67‑1 (MQ=255)
                         aCCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCaa                                     >  1:1241028/1‑69 (MQ=255)
                         aCCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCaa                                     >  1:2436333/1‑69 (MQ=255)
                         aCCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCaa                                     >  1:2602487/1‑69 (MQ=255)
                         aCCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCaa                                     >  1:1633723/1‑69 (MQ=255)
                         aCCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCaa                                     >  1:3264868/1‑69 (MQ=255)
                         aCCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCaa                                     >  1:905006/1‑69 (MQ=255)
                         aCCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCaa                                     >  1:1707965/1‑69 (MQ=255)
                         aCCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCaa                                     >  1:2164984/1‑69 (MQ=255)
                         aCCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCaa                                     >  1:2109491/1‑69 (MQ=255)
                         aCCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCaa                                     >  1:605553/1‑69 (MQ=255)
                          ccGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCaaa                                    >  1:321655/1‑69 (MQ=255)
                          ccGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCaaa                                    >  1:1949683/1‑69 (MQ=255)
                            gTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCAAACg                                  <  1:3054522/69‑1 (MQ=255)
                               aCGATCGTAGAACTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTAc                               >  1:1169282/1‑69 (MQ=255)
                               aCGATCGTAGAACTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTAc                               >  1:1697294/1‑69 (MQ=255)
                               aCGATCGTAGAACTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTAc                               >  1:1482159/1‑69 (MQ=255)
                                cGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTACCTTCTg                                            <  1:1484796/55‑1 (MQ=255)
                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:2948949/42‑1 (MQ=255)
                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:3098962/42‑1 (MQ=255)
                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:2925066/42‑1 (MQ=255)
                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:2559859/42‑1 (MQ=255)
                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:1540651/42‑1 (MQ=255)
                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:3233246/42‑1 (MQ=255)
                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:1166556/42‑1 (MQ=255)
                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:328234/42‑1 (MQ=255)
                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:434696/42‑1 (MQ=255)
                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:445238/42‑1 (MQ=255)
                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:127837/42‑1 (MQ=255)
                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:510836/42‑1 (MQ=255)
                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:562129/42‑1 (MQ=255)
                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:222989/42‑1 (MQ=255)
                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:173778/42‑1 (MQ=255)
                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:1853376/42‑1 (MQ=255)
                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:1110124/42‑1 (MQ=255)
                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:1650590/42‑1 (MQ=255)
                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:266073/42‑1 (MQ=255)
                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:2402183/42‑1 (MQ=255)
                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:2441525/42‑1 (MQ=255)
                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:819321/42‑1 (MQ=255)
                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:255404/42‑1 (MQ=255)
                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:2651372/42‑1 (MQ=255)
                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTCCGTaataat                                                        <  1:2405485/42‑1 (MQ=38)
                                             tGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGc                 >  1:1703813/1‑69 (MQ=255)
                                                cATCTTGAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCAAACg                                  <  1:2100809/49‑1 (MQ=255)
                                                    ttCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTc          >  1:1118077/1‑69 (MQ=255)
                                                    ttCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTc          >  1:131986/1‑69 (MQ=255)
                                                    ttCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTc          >  1:1217996/1‑69 (MQ=255)
                                                    ttCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTc          >  1:1878939/1‑69 (MQ=255)
                                                    ttCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTc          >  1:180400/1‑69 (MQ=255)
                                                         cgcgGATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGtt     <  1:2425197/69‑1 (MQ=255)
                                                         cgcgGATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGtt     <  1:1417341/69‑1 (MQ=255)
                                                            ggATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGTTTAc  <  1:113202/69‑1 (MQ=255)
                                                            ggATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGTTTAc  <  1:1740448/69‑1 (MQ=255)
                                                            ggATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGTTTAc  <  1:1135828/69‑1 (MQ=255)
                                                            ggATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGTTTAc  <  1:554350/69‑1 (MQ=255)
                                                            ggATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGTTTAc  <  1:539922/69‑1 (MQ=255)
                                                            ggATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGTTTAc  <  1:344529/69‑1 (MQ=255)
                                                            ggATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGTTTAc  <  1:1458076/69‑1 (MQ=255)
                                                            ggATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGTTTAc  <  1:3286485/69‑1 (MQ=255)
                                                            ggATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGTTTAc  <  1:3265404/69‑1 (MQ=255)
                                                            ggATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGTTTAc  <  1:2512465/69‑1 (MQ=255)
                                                            ggATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGTTTAc  <  1:2501010/69‑1 (MQ=255)
                                                            ggATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGTTTAc  <  1:1544470/69‑1 (MQ=255)
                                                            ggATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCCGGCTTTGCTCACGTTTAc  <  1:1914997/69‑1 (MQ=255)
                                                             gATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGTTTAc  <  1:2590691/68‑1 (MQ=255)
                                                               tAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGtt     >  1:523302/1‑63 (MQ=255)
                                                               tAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGtt     >  1:1733748/1‑63 (MQ=255)
                                                               tAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGtt     >  1:1010475/1‑63 (MQ=255)
                                                                |                                                                
GCAAGAATCCTTTCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGTTTAC  >  minE/1244696‑1244824

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: