Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,971,934 G→A 41.2% S324S (TCC→TCT speA ← biosynthetic arginine decarboxylase, PLP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,971,9340GA41.2% 33.5 / 62.4 89S324S (TCC→TCTspeAbiosynthetic arginine decarboxylase, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base G (17/34);  new base A (19/17);  total (36/53)
Fisher's exact test for biased strand distribution p-value = 8.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATT  >  minE/1971875‑1971999
                                                           |                                                                 
ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                                                          >  1:3010784/1‑69 (MQ=255)
ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                                                          >  1:370606/1‑69 (MQ=255)
ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                                                          >  1:230099/1‑69 (MQ=255)
ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                                                          >  1:380015/1‑69 (MQ=255)
ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                                                          >  1:389537/1‑69 (MQ=255)
ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                                                          >  1:3423026/1‑69 (MQ=255)
ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                                                          >  1:3331415/1‑69 (MQ=255)
ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                                                          >  1:1923788/1‑69 (MQ=255)
ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                                                          >  1:3113155/1‑69 (MQ=255)
ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                                                          >  1:1695219/1‑69 (MQ=255)
ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                                                          >  1:2639000/1‑69 (MQ=255)
ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                                                          >  1:473729/1‑69 (MQ=255)
ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                                                          >  1:529345/1‑69 (MQ=255)
ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                                                          >  1:2784022/1‑69 (MQ=255)
ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                                                          >  1:68579/1‑69 (MQ=255)
ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                                                          >  1:782268/1‑69 (MQ=255)
ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                                                          >  1:1196652/1‑69 (MQ=255)
ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                                                          >  1:902332/1‑69 (MQ=255)
ccAATCGCCCAGATGATGTTGTTGGCGTATTCATTGAGGCCGTAGTTTACGGAGCAGTCAGACTGAGaa                                                          >  1:1229318/1‑69 (MQ=255)
  aaTCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACg                                                        <  1:1523618/69‑1 (MQ=255)
  aaTCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACg                                                        <  1:2124176/69‑1 (MQ=255)
     cGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGt                                                     >  1:2364673/1‑69 (MQ=255)
     cGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGt                                                     >  1:3208619/1‑69 (MQ=255)
              aaTGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATa                                            <  1:2483269/69‑1 (MQ=255)
              aaTGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATa                                            <  1:2487157/69‑1 (MQ=255)
               aTGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATaa                                           <  1:2054038/69‑1 (MQ=255)
               aTGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATaa                                           <  1:1300562/69‑1 (MQ=255)
               aTGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATaa                                           <  1:2939213/69‑1 (MQ=255)
               aTGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATaa                                           <  1:3107842/69‑1 (MQ=255)
               aTGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATaa                                           <  1:853436/69‑1 (MQ=255)
               aTGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATaa                                           <  1:420101/69‑1 (MQ=255)
               aTGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATaa                                           <  1:283719/69‑1 (MQ=255)
               aTGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATa                                            <  1:419807/68‑1 (MQ=255)
               aTGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATa                                            <  1:2373790/68‑1 (MQ=255)
                 gttgttGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATc                                         >  1:2524258/1‑69 (MQ=255)
                 gttgttGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATc                                         >  1:967167/1‑69 (MQ=255)
                 gttgttGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATc                                         >  1:2045850/1‑69 (MQ=255)
                 gttgttGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATc                                         >  1:622461/1‑69 (MQ=255)
                  ttgttgGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCTGACTGCTAACGAGTACCTTCATaa                                           <  1:2907334/66‑1 (MQ=255)
                    gttgGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCa                                       <  1:3430517/68‑1 (MQ=255)
                    gttgGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCa                                       <  1:3241055/68‑1 (MQ=255)
                    gttgGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCa                                       <  1:403994/68‑1 (MQ=255)
                    gttgGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCa                                       <  1:2688682/68‑1 (MQ=255)
                            aTTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGAc                              <  1:1679931/69‑1 (MQ=255)
                                   gAGGCCGTAGTTCACCGAACAGTCGGACTGCTAACGAGTACCTTCATAATCCACGCCCAGAccgccgcc                       <  1:217942/69‑1 (MQ=255)
                                   gAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgccgcc                       <  1:2430607/69‑1 (MQ=255)
                                   gAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgccgcc                       <  1:121414/69‑1 (MQ=255)
                                   gAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgccgcc                       <  1:1514907/69‑1 (MQ=255)
                                   gAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgccgcc                       <  1:2667815/69‑1 (MQ=255)
                                   gAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgccgcc                       <  1:4162/69‑1 (MQ=255)
                                   gAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgccgcc                       <  1:3354413/69‑1 (MQ=255)
                                   gAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgccgcc                       <  1:2111334/69‑1 (MQ=255)
                                            gTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGa               >  1:927662/1‑68 (MQ=255)
                                            gTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGa               >  1:1183401/1‑68 (MQ=255)
                                            gTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGa               >  1:1087814/1‑68 (MQ=255)
                                            gTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGa               >  1:2278498/1‑68 (MQ=255)
                                            gTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGa               >  1:2475707/1‑68 (MQ=255)
                                            gTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGACAGCCGCCGACGTCGa               >  1:1796948/1‑68 (MQ=255)
                                               cACCGAACAGTCGGACTGCGAACGAGTAGCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCa           >  1:1649022/1‑69 (MQ=255)
                                               cACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgccgcc                       >  1:1460357/1‑57 (MQ=255)
                                               cACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgccgcc                       >  1:2519722/1‑57 (MQ=255)
                                               cACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgccgcc                       >  1:515480/1‑57 (MQ=255)
                                               cACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCa           >  1:2626575/1‑69 (MQ=255)
                                               cACCGAACAGGCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgccgcc                       <  1:1597310/57‑1 (MQ=255)
                                                aCCGAACAGTCTGACTGCGAACGAGTACCTTCATa                                            <  1:370074/35‑1 (MQ=255)
                                                   gagcAGTCAGACTGAGAACGGGTGCCTTCATAGTCCACGCCCAGACCGCCGCCGACGTCGAAGCACTg        <  1:854288/65‑1 (MQ=255)
                                                   gagcAGTCAGACTGAGAACGGGTGCCTTCATAGTCCACGCCCAGACCGCCGCCGACGTCGAAGCACTg        <  1:1458448/65‑1 (MQ=255)
                                                   gagcAGTCAGACTGAGAACGGGTGCCTTCATAGTCCACGCCCAGACCGCCGCCGACGTCGAAGCACTg        <  1:1527328/65‑1 (MQ=255)
                                                   gagcAGTCAGACTGAGAACGGGTGCCTTCATAGTCCACGCCCAGACCGCCGCCGACGTCGAAGCACTg        <  1:856614/65‑1 (MQ=255)
                                                   gagcAGTCAGACTGAGAACGGGTGCCTTCATAGTCCACGCCCAGACCGCCGCCGACGTCGAAGCACTg        <  1:531889/65‑1 (MQ=255)
                                                   gagcAGTCAGACTGAGAACGGGTGCCTTCATAGTCCACGCCCAGACCGCCGCCGACGTCGAAGCACTg        <  1:1375786/65‑1 (MQ=255)
                                                   gagcAGTCAGACTGAGAACGGGTGCCTTCATAGTCCACGCCCAGACCGCCGCCGACGTCGAAGCACTg        <  1:3002854/65‑1 (MQ=255)
                                                   gagcAGTCAGACTGAGAACGGGTGCCTTCATAGTCCACGCCCAGACCGCCGCCGACGTCGAAGCACTg        <  1:3027568/65‑1 (MQ=255)
                                                   gagcAGTCAGACTGAGAACGGGTGCCTTCATAGTCCACGCCCAGACCGCCGCCGACGTCGAAGCACTg        <  1:1830043/65‑1 (MQ=255)
                                                   gagcAGTCAGACTGAGAACGGGTGCCTTCATAGTCCACGCCCAGACCGCCGCCGACGTCGAAGCACTg        <  1:2245832/65‑1 (MQ=255)
                                                   gagcAGTCAGACTGAGAACGGGTGCCTTCATAGTCCACGCCCAGACCGCCGCCGACGTCGAAGCACTg        <  1:2406956/65‑1 (MQ=255)
                                                   gagcAGTCAGACTGAGAACGGGTGCCTTCATAGTCCACGCCCAGACCGCCGCCGACGTCGAAGCACTg        <  1:2637931/65‑1 (MQ=255)
                                                   gagcAGTCAGACTGAGAACGGGTGCCTTCATAGTCCACGCCCAGACCGCCGCCGACGTCGAAGCACTg        <  1:335122/65‑1 (MQ=255)
                                                   gagcAGTCAGACTGAGAACGGGTGCCTTCATAGTCCACGCCCAGACCGCCGCCGACGTCGAAGCACTg        <  1:1653784/65‑1 (MQ=255)
                                                   gagcAGTCAGACTGAGAACGGGTGCCTTCATAGTCCACGCCCAGACCGCCGCCGACGTCGAAGCACTg        <  1:2578519/65‑1 (MQ=255)
                                                   gagcAGTCAGACTGAGAACGGGTGCCTTCATAGTCCACGCCCAGACCGCCGCCGACGTCGAAGCACTg        <  1:2618990/65‑1 (MQ=255)
                                                   gagcAGTCAGACTGAGAACGGGTGCCTTCATAGTCCACGCCCAGACCGCCGCCGACGTCGAAGCACTg        <  1:3262574/65‑1 (MQ=255)
                                                    aaCAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGaa      <  1:688541/69‑1 (MQ=255)
                                                    aaCAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGaa      <  1:1393149/69‑1 (MQ=255)
                                                    aaCAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGaa      <  1:836408/69‑1 (MQ=255)
                                                    aaCAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGaa      <  1:2539935/69‑1 (MQ=255)
                                                    aaCAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGaa      <  1:1760149/69‑1 (MQ=255)
                                                        gTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATAtt  <  1:1428845/69‑1 (MQ=255)
                                                        gTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATAtt  <  1:62339/69‑1 (MQ=255)
                                                           |                                                                 
CCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATT  >  minE/1971875‑1971999

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