Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,837,204 T→C 26.1% E39G (GAA→GGA)  yjeK ← predicted lysine aminomutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,837,2040TC26.1% 78.5 / 31.1 51E39G (GAA→GGA) yjeKpredicted lysine aminomutase
Reads supporting (aligned to +/- strand):  ref base T (13/24);  new base C (8/5);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 1.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

CGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGATCGGTCACAACATCGGCAAGTTGCGTT  >  minE/2837135‑2837272
                                                                     |                                                                    
cgatcgatAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTtc                                                                     >  1:367900/1‑71 (MQ=255)
      atAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTCCCTCcgc                                                               <  1:1861748/71‑1 (MQ=255)
       tAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTCCCTCcgcg                                                              >  1:2488574/1‑71 (MQ=255)
         aaTGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTc                                                            <  1:1087215/71‑1 (MQ=255)
             aGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCtata                                                        >  1:2171334/1‑71 (MQ=255)
                      cacGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTCCCTCCGCGTCTATATTCAAAAGa                                               >  1:539861/1‑71 (MQ=255)
                        cGCAGGGCAAAAAGCTTTTTGGCGCTGCGTAAGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGAc                                              >  1:1971145/1‑70 (MQ=255)
                         gCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACg                                             >  1:3136112/1‑70 (MQ=255)
                           aGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAg                                          <  1:2837239/71‑1 (MQ=255)
                           aGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTCCCTCCGCGTCTATATTCAAAAGACGCAg                                          <  1:1750510/71‑1 (MQ=255)
                             ggCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGAc                                              <  1:3092369/65‑1 (MQ=255)
                             ggCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAg                                          >  1:574489/1‑69 (MQ=255)
                             ggCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGa                                         <  1:594269/70‑1 (MQ=255)
                              gCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGAc                                              >  1:1677525/1‑64 (MQ=255)
                              gCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTCCCTCCGCGTCTATATTCAAAAGACGCAGAAg                                       >  1:1266071/1‑71 (MQ=255)
                               cAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTCCCTCCGCGTCTATATTCAAAAGACGCAGAAGt                                      >  1:4517700/1‑71 (MQ=255)
                                 aaaaGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGtt                                     >  1:1348110/1‑70 (MQ=255)
                                 aaaaGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGtt                                     >  1:905609/1‑70 (MQ=255)
                                  aaaGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTtca                                   <  1:726527/71‑1 (MQ=255)
                                    aGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCaaaa                                                 >  1:1532222/1‑55 (MQ=255)
                                       tttttGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTtca                                   <  1:3394777/66‑1 (MQ=255)
                                       tttttGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTtca                                   <  1:2567828/66‑1 (MQ=255)
                                       tttttGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTtc                                    <  1:2575059/65‑1 (MQ=255)
                                       tttttGGCGCTGCGTCCGGCTAACAGTTTTACCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAg                               >  1:1574120/1‑70 (MQ=255)
                                        ttttGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGa                             >  1:4186769/1‑71 (MQ=255)
                                        ttttGGCGCTGCGTCCGGCTAACAGTTTTCCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGa                             >  1:972868/1‑71 (MQ=255)
                                         tttGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAgg                              <  1:611251/69‑1 (MQ=255)
                                         tttGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAt                            >  1:2103004/1‑71 (MQ=255)
                                          ttGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGAt                            <  1:1748698/70‑1 (MQ=255)
                                          ttGGCGCTGCGTCAGGCTAACAGTTTTCCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGa                             <  1:3422080/69‑1 (MQ=255)
                                            ggCGCTGCGTCCGGCTAACAGTTTTCCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAggt                             >  1:2974499/1‑66 (MQ=255)
                                               gcTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGATCGGTca                      >  1:708063/1‑71 (MQ=255)
                                                cTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGATCGGTcac                     <  1:1309469/71‑1 (MQ=255)
                                                cTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGATCGGTcac                     <  1:2286513/71‑1 (MQ=255)
                                                 tGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCaaaa                                                 <  1:2789719/42‑1 (MQ=255)
                                                 tGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGATCGGTcac                     >  1:2755181/1‑70 (MQ=255)
                                                 tGCGTCCGGCTAACAGTTTTCCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGATCGGTcac                     >  1:3805300/1‑70 (MQ=255)
                                                   cGTCCGGCTAACAGTTTTCCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGATCGGTCacaac                  <  1:3251887/71‑1 (MQ=255)
                                                     tcaggCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGATCGGTCACAACAt                <  1:4224415/68‑1 (MQ=255)
                                                     tCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGATCGGTCacaaca                 >  1:1274591/1‑70 (MQ=255)
                                                     tCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGATCGGTCACAACAt                <  1:4343032/71‑1 (MQ=255)
                                                       cGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGATCGGTCACAACATCg              <  1:3644892/71‑1 (MQ=255)
                                                         gCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGt                                      <  1:3460306/45‑1 (MQ=255)
                                                         gCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGATCGGTcac                     <  1:4547469/62‑1 (MQ=255)
                                                         gCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGATCGGTCACAACATCGGc            <  1:186624/71‑1 (MQ=255)
                                                         gCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGATCGGTCACAACATCGGc            <  1:140642/71‑1 (MQ=255)
                                                          cTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGATCGGTCACAACATCGGc            <  1:143864/70‑1 (MQ=255)
                                                            aaCAGTTTTCCCTCCGCGTCTATATTCAAAAGACGCAg                                          <  1:732662/38‑1 (MQ=255)
                                                             aCAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGATCGGTCACAACATCGGCaa          <  1:1016060/69‑1 (MQ=255)
                                                             aCAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGATCGGTCACAACATCGGCaa          <  1:2539962/69‑1 (MQ=255)
                                                             aCAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGATCGGTCACAACATCGGCaa          <  1:2545455/69‑1 (MQ=255)
                                                               aGTTTTCCCTCCGCGTCTATATTCAAAAGACGCAGAAGt                                      >  1:1916647/1‑39 (MQ=255)
                                                                 tttttCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGATCGGTCACAACATCGGCAAGTTGc     >  1:3748100/1‑70 (MQ=255)
                                                                  tttcCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGATCGGTca                      >  1:4165931/1‑52 (MQ=255)
                                                                  tttcCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGATCGGTCACAACATCGGCAAGTTGCGt   >  1:289099/1‑71 (MQ=255)
                                                                   tttCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGATCGGTCACAACATCGGCAAGTTGCGt   >  1:536022/1‑70 (MQ=255)
                                                                    ttCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGATCGGTCACAACATCGGCAAGTTGCGtt  <  1:791879/70‑1 (MQ=255)
                                                                     |                                                                    
CGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTTCATCAGGATCGGTCACAACATCGGCAAGTTGCGTT  >  minE/2837135‑2837272

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: