Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,792,318 C→T 11.1% I234I (ATC→ATT gabT → 4‑aminobutyrate aminotransferase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,792,3180CT11.1% 220.4 / 13.5 81I234I (ATC→ATTgabT4‑aminobutyrate aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (32/40);  new base T (5/4);  total (37/44)
Fisher's exact test for biased strand distribution p-value = 7.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

AGGCGGTTTCTACGCCTCGTCGCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCTGTTTGCGATGGAGCAGATG  >  minE/1792249‑1792387
                                                                     |                                                                     
aGGCGGTTTCTACGCCTCGTCGCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCa                                                                      <  1:3694155/71‑1 (MQ=255)
    ggTTTCTACGCCTCGTCGCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCt                                                                  <  1:4571069/71‑1 (MQ=255)
       ttCTACGCCTCGTCGCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGAt                                                               >  1:4443432/1‑71 (MQ=255)
        tCTACGCCTCGTCGCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTTACGAGCACGGGATCATGCTGAtt                                                              >  1:2877950/1‑71 (MQ=255)
            cGCCTCGTCGCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGcc                                                           >  1:4371433/1‑70 (MQ=255)
            cGCCTCGTCGCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGc                                                            <  1:1034862/69‑1 (MQ=255)
               cTCGTCGCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACTAGCACGGGATCATGCTGATTGc                                                            <  1:5473312/66‑1 (MQ=255)
               cTCGTCGCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATg                                                       <  1:5214815/71‑1 (MQ=255)
               cTCGTCGCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATg                                                       <  1:3023233/71‑1 (MQ=255)
                tcgtcgCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCa                                                                      <  1:405858/55‑1 (MQ=255)
                tcgtcgCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCa                                                                      <  1:5311388/55‑1 (MQ=255)
                tcgtcgCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGAtt                                                              <  1:876693/63‑1 (MQ=255)
                tcgtcgCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGAtt                                                              <  1:1491870/63‑1 (MQ=255)
                tcgtcgCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGAtt                                                              <  1:4722991/63‑1 (MQ=255)
                tcgtcgCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGAtt                                                              <  1:349032/63‑1 (MQ=255)
                  gtcgCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATg                                                       <  1:2151589/68‑1 (MQ=255)
                  gtcgCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGaa                                                     <  1:4960442/70‑1 (MQ=255)
                  gtcgCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGa                                                      <  1:2090819/69‑1 (MQ=255)
                    cgCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGAt                                                               >  1:1352033/1‑58 (MQ=255)
                    cgCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGAt                                                               >  1:1566072/1‑58 (MQ=255)
                    cgCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGAt                                                               >  1:2296925/1‑58 (MQ=255)
                    cgCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGAt                                                               >  1:4995256/1‑58 (MQ=255)
                    cgCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGAt                                                               >  1:3393540/1‑58 (MQ=255)
                    cgCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGAt                                                               >  1:3119042/1‑58 (MQ=255)
                    cgCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGAAGTg                                                  >  1:1732872/1‑71 (MQ=255)
                     gCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGcc                                                           <  1:1719065/61‑1 (MQ=255)
                     gCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGcc                                                           <  1:3290203/61‑1 (MQ=255)
                     gCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGcc                                                           <  1:2223772/61‑1 (MQ=255)
                     gCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGAAGTg                                                  >  1:5565677/1‑70 (MQ=255)
                     gCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGAAGTg                                                  <  1:5211219/70‑1 (MQ=255)
                     gCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGAAGTg                                                  >  1:1320184/1‑70 (MQ=255)
                     gCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGAAGTGc                                                 <  1:568096/71‑1 (MQ=255)
                       cAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGc                                                            <  1:2647816/58‑1 (MQ=255)
                       cAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGc                                                            <  1:278195/58‑1 (MQ=255)
                        aGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCa                                                >  1:4922956/1‑69 (MQ=255)
                         gCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCa                                                >  1:5103033/1‑68 (MQ=255)
                           cTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATg                                                       <  1:4135747/59‑1 (MQ=255)
                              tATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCAGAgcgg                                          <  1:1002385/69‑1 (MQ=255)
                               aTGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCAGAgcggcg                                        <  1:3271303/70‑1 (MQ=255)
                                tGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGa                                                                >  1:4482616/1‑45 (MQ=255)
                                tGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGAAg                                                    <  1:2577810/57‑1 (MQ=255)
                                   aGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCAGAGCGGCGCgggg                                   >  1:3037624/1‑71 (MQ=255)
                                    gCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCAGAgc                                            >  1:1679210/1‑61 (MQ=255)
                                     cGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCaga                                              >  1:5051213/1‑58 (MQ=255)
                                     cGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCaga                                              >  1:2879982/1‑58 (MQ=255)
                                     cGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCaga                                              >  1:4511673/1‑58 (MQ=255)
                                       tttACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGcc                                                           <  1:4685251/43‑1 (MQ=255)
                                       tttACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGcc                                                           <  1:4032040/43‑1 (MQ=255)
                                       tttACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGcc                                                           <  1:1842095/43‑1 (MQ=255)
                                         tACGCGCTCTGTGTGACGAGCACGGGATCATGCTGAtt                                                              >  1:4306828/1‑38 (MQ=255)
                                         tACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGa                                                         >  1:759671/1‑43 (MQ=255)
                                                tctGTGTGACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGc                      >  1:1664302/1‑71 (MQ=255)
                                                 ctGTGTGACGAGCACGGGATCATGCTGATTGCCGATg                                                       <  1:5180370/37‑1 (MQ=255)
                                                 ctGTGTGACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCAGAGCGGCGCGGGGCGTAc                              >  1:4993916/1‑62 (MQ=255)
                                                 ctGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCaaag                                             >  1:725134/1‑47 (MQ=12)
                                                 ctGTGCGACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCt                     <  1:4652432/71‑1 (MQ=255)
                                                  tgtgtgACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCTg                    >  1:4684396/1‑71 (MQ=255)
                                                   gtgcgACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAgcggct                                        <  1:2046466/50‑2 (MQ=12)
                                                    tgtgACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCTg                    <  1:1308500/69‑1 (MQ=255)
                                                    tgtgACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCTg                    <  1:849431/69‑1 (MQ=255)
                                                    tgtgACGAGCACGGGATCATGCTGAATGCCGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCt                     <  1:4720948/68‑1 (MQ=255)
                                                     gtgACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCagag                                             <  1:4913888/43‑1 (MQ=255)
                                                     gtgACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCagag                                             <  1:3701008/43‑1 (MQ=255)
                                                     gtgACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCagag                                             <  1:4194455/43‑1 (MQ=255)
                                                     gcgACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAgcggc                                         >  1:4235612/3‑47 (MQ=12)
                                                     gcgACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGttt                 >  1:4831320/3‑71 (MQ=255)
                                                      tgACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCAGAGCGGCGCggg                                    >  1:1961799/1‑51 (MQ=255)
                                                      tgACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCAGAGCGGCGCggg                                    >  1:532730/1‑51 (MQ=255)
                                                      tgACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCAGAGCGGCGCggg                                    >  1:4733063/1‑51 (MQ=255)
                                                      tgACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCAGAGCGGCGCggg                                    >  1:2040951/1‑51 (MQ=255)
                                                       gatgAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGc               <  1:893522/68‑1 (MQ=255)
                                                        aCGAGCACGGGATCATGCTGATTGCCGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCTGTTTGc               <  1:2704639/70‑1 (MQ=255)
                                                        aCGAGCACGGGATCATGCTGATTGCCGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCTGTTTGc               <  1:782166/70‑1 (MQ=255)
                                                          gAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGtt                  >  1:5575742/1‑65 (MQ=255)
                                                            gCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGCCATgg          <  1:5162724/71‑1 (MQ=255)
                                                             cACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCgg                           >  1:2652280/1‑53 (MQ=255)
                                                             cACGGGATCATGCTGATTGCCGATGAAGTGCAGAgcg                                           >  1:2712028/1‑37 (MQ=255)
                                                             cACGGGATCATGCTGATTGCCGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCTGTTTGCGATg           >  1:2924971/1‑69 (MQ=255)
                                                                gggATCATGCTGATTGCCGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCTGTTTGCGATGGAGCa      <  1:2439660/71‑1 (MQ=255)
                                                                    tCATGCTGATTGCCGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCTg                    <  1:3937224/53‑1 (MQ=255)
                                                                    tCATGCTGATTGCCGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCTGTTTGCGATGGAGCAGATg  >  1:606700/1‑71 (MQ=255)
                                                                     |                                                                     
AGGCGGTTTCTACGCCTCGTCGCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCTGTTTGCGATGGAGCAGATG  >  minE/1792249‑1792387

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