Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,071,466 T→C 32.7% N429N (AAT→AAC 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,4660TC32.7% 48.3 / 32.1 64N429N (AAT→AACygjGputrescine:2‑oxoglutaric acid aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (22/21);  new base C (7/14);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 1.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

AACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAA  >  minE/2071398‑2071534
                                                                    |                                                                    
aaCGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAAt                                                                      <  1:613419/69‑1 (MQ=255)
aaCGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAAt                                                                      <  1:180272/69‑1 (MQ=255)
aaCGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAAt                                                                      <  1:293743/69‑1 (MQ=255)
aaCGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAAt                                                                      <  1:146983/69‑1 (MQ=255)
aaCGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAAt                                                                      <  1:781849/69‑1 (MQ=255)
aaCGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAAt                                                                      <  1:2101472/69‑1 (MQ=255)
 aCGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAAt                                                                      <  1:1997300/68‑1 (MQ=255)
   gAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAAc                                                                   <  1:1353579/69‑1 (MQ=255)
   gAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAAc                                                                   <  1:1311992/69‑1 (MQ=255)
   gAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAAc                                                                   <  1:2194252/69‑1 (MQ=255)
   gAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAAc                                                                   <  1:2052709/69‑1 (MQ=255)
   gAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAAc                                                                   <  1:103686/69‑1 (MQ=255)
            tataACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACg                                                          <  1:146647/69‑1 (MQ=255)
            tataACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACg                                                          <  1:3298587/69‑1 (MQ=255)
            tataACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACg                                                          <  1:319376/69‑1 (MQ=255)
            tataACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACg                                                          <  1:1728540/69‑1 (MQ=255)
            tataACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACg                                                          <  1:1792116/69‑1 (MQ=255)
                aCTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATc                                                       >  1:923113/1‑68 (MQ=255)
                       cAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACACTCAACAACGCTAAAACGATCCGCATTg                                                <  1:2188672/68‑1 (MQ=255)
                       cAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACACTCAACAACGCTAAAACGATCCGCATTg                                                <  1:2011994/68‑1 (MQ=255)
                       cAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACACTCAACAACGCTAAAACGATCCGCATTg                                                <  1:2329145/68‑1 (MQ=255)
                       cAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACACTCAACAACGCTAAAACGATCCGCATTg                                                <  1:1863135/68‑1 (MQ=255)
                       cAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACACTCAACAACGCTAAAACGATCCGCATTg                                                <  1:1612203/68‑1 (MQ=255)
                       cAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACACTCAACAACGCTAAAACGATCCGCATTg                                                <  1:3259106/68‑1 (MQ=255)
                       cAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACACTCAACAACGCTAAAACGATCCGCATTg                                                <  1:537244/68‑1 (MQ=255)
                       cAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACACTCAACAACGCTAAAACGATCCGCATTg                                                <  1:732964/68‑1 (MQ=255)
                       cAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACACTCAACAACGCTAAAACGATCCGCATTg                                                <  1:827110/68‑1 (MQ=255)
                        aGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACACTCAACAACGCTAAAACGATCCGCATTg                                                <  1:3083277/67‑1 (MQ=255)
                        aGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACACTCAACAACGCTAAAACGATCCGCATTg                                                <  1:2291533/67‑1 (MQ=255)
                        aGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACACTCAACAACGCTAAAACGATCCGCATTg                                                <  1:2508389/67‑1 (MQ=255)
                           gagaTGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCa                                                   >  1:579053/1‑61 (MQ=255)
                           gagaTGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAAcc                                            >  1:3352268/1‑68 (MQ=255)
                           gagaTGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAAcc                                            >  1:964077/1‑68 (MQ=255)
                           gagaTGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAAcc                                            >  1:842373/1‑68 (MQ=255)
                           gagaTGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAAcc                                            >  1:634993/1‑68 (MQ=255)
                           gagaTGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAAcc                                            >  1:1505064/1‑68 (MQ=255)
                           gagaTGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAAcc                                            >  1:1512116/1‑68 (MQ=255)
                           gagaTGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAAcc                                            >  1:369900/1‑68 (MQ=255)
                           gagaTGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAAcc                                            >  1:3387116/1‑68 (MQ=255)
                           gagaTGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAAcc                                            >  1:2056497/1‑68 (MQ=255)
                           gagaTGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAAcc                                            >  1:3349300/1‑68 (MQ=255)
                           gagaTGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAAcc                                            >  1:1781713/1‑68 (MQ=255)
                           gagaTGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAAcc                                            >  1:1026080/1‑68 (MQ=255)
                           gagaTGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAAcc                                            >  1:2320659/1‑68 (MQ=255)
                           gagaTGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAAcc                                            >  1:2902372/1‑68 (MQ=255)
                           gagaTGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAAcc                                            >  1:205617/1‑68 (MQ=255)
                              aTGTTCCGCCAGGGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCca                                        >  1:1446023/1‑69 (MQ=255)
                                 ttCCGCCAGCGCGTACTGGTGGCCGGAACACTCAACAACGCTAAAACGATCCGCATTg                                                <  1:942785/58‑1 (MQ=255)
                                   ccGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCcactgac                                   <  1:714369/69‑1 (MQ=255)
                                   ccGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCcactgac                                   <  1:2991447/69‑1 (MQ=255)
                                   ccGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCcactgac                                   <  1:626131/69‑1 (MQ=255)
                                         gcgcGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCactgacactgac                             <  1:1606568/69‑1 (MQ=255)
                                         gcgcGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCactgacactgac                             <  1:2866502/69‑1 (MQ=255)
                                                cTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTg                        >  1:1878232/1‑67 (MQ=255)
                                                cTGGTGGCCGGAACACTCAACAACGCTAAAACGATCCGCATTg                                                <  1:1462862/43‑1 (MQ=37)
                                                  ggtggCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAAc                     <  1:917162/68‑1 (MQ=255)
                                                  ggtggCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAAc                     <  1:1603782/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:2711566/1‑69 (MQ=255)
                                                               cTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTg       >  1:1571769/1‑69 (MQ=255)
                                                                tCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAAc                     <  1:2843872/54‑1 (MQ=255)
                                                                tCAACAACGCTAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTg       >  1:2353232/1‑68 (MQ=255)
                                                                tCAACAACGCTAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTg       >  1:1201321/1‑68 (MQ=255)
                                                                tCAACAACGCTAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTg       >  1:884929/1‑68 (MQ=255)
                                                                tCAACAACGCTAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTg       >  1:1685649/1‑68 (MQ=255)
                                                                tCAACAACGCTAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTg       >  1:1073320/1‑68 (MQ=255)
                                                                tCAACAACGCTAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTg       >  1:3046198/1‑68 (MQ=255)
                                                                tCAACAACGCTAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTg       >  1:2194319/1‑68 (MQ=255)
                                                                    tAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCaa  <  1:618547/69‑1 (MQ=255)
                                                                    |                                                                    
AACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAA  >  minE/2071398‑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: