Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,071,472 C→T 29.4% A431A (GCC→GCT ygjG → putrescine:2‑oxoglutaric acid aminotransferase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,071,4720CT29.4% 70.3 / 29.4 71A431A (GCC→GCTygjGputrescine:2‑oxoglutaric acid aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (32/18);  new base T (7/14);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 2.13e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

TATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAA  >  minE/2071410‑2071534
                                                              |                                                              
tataACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACg                                                          <  1:3298587/69‑1 (MQ=255)
tataACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACg                                                          <  1:146647/69‑1 (MQ=255)
tataACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACg                                                          <  1:319376/69‑1 (MQ=255)
tataACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACg                                                          <  1:1792116/69‑1 (MQ=255)
tataACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACg                                                          <  1:1728540/69‑1 (MQ=255)
    aCTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATc                                                       >  1:923113/1‑68 (MQ=255)
           cAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACACTCAACAACGCTAAAACGATCCGCATTg                                                <  1:2188672/68‑1 (MQ=255)
           cAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACACTCAACAACGCTAAAACGATCCGCATTg                                                <  1:827110/68‑1 (MQ=255)
           cAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACACTCAACAACGCTAAAACGATCCGCATTg                                                <  1:2011994/68‑1 (MQ=255)
           cAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACACTCAACAACGCTAAAACGATCCGCATTg                                                <  1:1863135/68‑1 (MQ=255)
           cAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACACTCAACAACGCTAAAACGATCCGCATTg                                                <  1:2329145/68‑1 (MQ=255)
           cAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACACTCAACAACGCTAAAACGATCCGCATTg                                                <  1:732964/68‑1 (MQ=255)
           cAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACACTCAACAACGCTAAAACGATCCGCATTg                                                <  1:1612203/68‑1 (MQ=255)
           cAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACACTCAACAACGCTAAAACGATCCGCATTg                                                <  1:537244/68‑1 (MQ=255)
           cAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACACTCAACAACGCTAAAACGATCCGCATTg                                                <  1:3259106/68‑1 (MQ=255)
            aGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACACTCAACAACGCTAAAACGATCCGCATTg                                                <  1:2508389/67‑1 (MQ=255)
            aGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACACTCAACAACGCTAAAACGATCCGCATTg                                                <  1:3083277/67‑1 (MQ=255)
            aGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACACTCAACAACGCTAAAACGATCCGCATTg                                                <  1:2291533/67‑1 (MQ=255)
               gagaTGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCa                                                   >  1:579053/1‑61 (MQ=255)
               gagaTGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAAcc                                            >  1:634993/1‑68 (MQ=255)
               gagaTGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAAcc                                            >  1:2902372/1‑68 (MQ=255)
               gagaTGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAAcc                                            >  1:3349300/1‑68 (MQ=255)
               gagaTGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAAcc                                            >  1:2320659/1‑68 (MQ=255)
               gagaTGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAAcc                                            >  1:3352268/1‑68 (MQ=255)
               gagaTGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAAcc                                            >  1:3387116/1‑68 (MQ=255)
               gagaTGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAAcc                                            >  1:369900/1‑68 (MQ=255)
               gagaTGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAAcc                                            >  1:205617/1‑68 (MQ=255)
               gagaTGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAAcc                                            >  1:1505064/1‑68 (MQ=255)
               gagaTGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAAcc                                            >  1:842373/1‑68 (MQ=255)
               gagaTGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAAcc                                            >  1:1026080/1‑68 (MQ=255)
               gagaTGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAAcc                                            >  1:1512116/1‑68 (MQ=255)
               gagaTGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAAcc                                            >  1:1781713/1‑68 (MQ=255)
               gagaTGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAAcc                                            >  1:2056497/1‑68 (MQ=255)
               gagaTGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAAcc                                            >  1:964077/1‑68 (MQ=255)
                  aTGTTCCGCCAGGGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCca                                        >  1:1446023/1‑69 (MQ=255)
                     ttCCGCCAGCGCGTACTGGTGGCCGGAACACTCAACAACGCTAAAACGATCCGCATTg                                                <  1:942785/58‑1 (MQ=255)
                       ccGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCcactgac                                   <  1:626131/69‑1 (MQ=255)
                       ccGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCcactgac                                   <  1:2991447/69‑1 (MQ=255)
                       ccGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCcactgac                                   <  1:714369/69‑1 (MQ=255)
                             gcgcGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCactgacactgac                             <  1:2866502/69‑1 (MQ=255)
                             gcgcGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCactgacactgac                             <  1:1606568/69‑1 (MQ=255)
                                    cTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTg                        >  1:1878232/1‑67 (MQ=255)
                                    cTGGTGGCCGGAACACTCAACAACGCTAAAACGATCCGCATTg                                                <  1:1462862/43‑1 (MQ=37)
                                      ggtggCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAAc                     <  1:1603782/68‑1 (MQ=255)
                                      ggtggCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAAc                     <  1:917162/68‑1 (MQ=255)
                                      ggtggCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAAc                     <  1:3070554/68‑1 (MQ=255)
                                          gCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCACTGACACTGACCATTGAACAGtgtg               <  1:2941902/69‑1 (MQ=255)
                                            cGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGtgtg               <  1:340585/68‑1 (MQ=255)
                                                   cTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTg       >  1:420761/1‑69 (MQ=255)
                                                   cTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTg       >  1:1571769/1‑69 (MQ=255)
                                                   cTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTg       >  1:2711566/1‑69 (MQ=255)
                                                    tCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAAc                     <  1:2843872/54‑1 (MQ=255)
                                                    tCAACAACGCTAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTg       >  1:884929/1‑68 (MQ=255)
                                                    tCAACAACGCTAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTg       >  1:1073320/1‑68 (MQ=255)
                                                    tCAACAACGCTAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTg       >  1:1201321/1‑68 (MQ=255)
                                                    tCAACAACGCTAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTg       >  1:3046198/1‑68 (MQ=255)
                                                    tCAACAACGCTAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTg       >  1:1685649/1‑68 (MQ=255)
                                                    tCAACAACGCTAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTg       >  1:2353232/1‑68 (MQ=255)
                                                    tCAACAACGCTAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTg       >  1:2194319/1‑68 (MQ=255)
                                                        tAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCaa  <  1:618547/69‑1 (MQ=255)
                                                         aaCGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCaa  <  1:2793302/68‑1 (MQ=255)
                                                          aCGCCAAAACGTTCCGCATTGAACCGCCACTGACACTGACCATTGAACAGtgtg               >  1:2894181/1‑54 (MQ=255)
                                                          aCGCCAAAACGATCCGCATTGAACCGCGACTGACACTGACCATTGAACAGtgtg               >  1:1469281/1‑54 (MQ=255)
                                                          aCGCCAAAACGATCCGCATTGAACCGCCACTGACACTGAGCATTGAACAGtgtg               >  1:2525192/1‑54 (MQ=255)
                                                          aCGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGtgtg               >  1:2396104/1‑54 (MQ=255)
                                                          aCGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGtgtg               >  1:1705939/1‑54 (MQ=255)
                                                          aCGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGtgtg               >  1:1798099/1‑54 (MQ=255)
                                                          aCGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGtgtg               >  1:181102/1‑54 (MQ=255)
                                                          aCGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGtgtg               >  1:3296773/1‑54 (MQ=255)
                                                          aCGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGtgtg               >  1:58723/1‑54 (MQ=255)
                                                          aCGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGtgtg               >  1:523205/1‑54 (MQ=255)
                                                              |                                                              
TATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAA  >  minE/2071410‑2071534

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: