Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,792,330 C→T 14.0% A238A (GCC→GCT gabT → 4‑aminobutyrate aminotransferase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,792,3300CT14.0% 160.2 / 14.6 64A238A (GCC→GCTgabT4‑aminobutyrate aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (30/25);  new base T (5/4);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

CGCCTCGTCGCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCTGTTTGCGATGGAGCAGATGGGCGTTGCGCCGG  >  Exported/1792261‑1792400
                                                                     |                                                                      
cGCCTCGTCGCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGcc                                                                        >  1:4371433/1‑70 (MQ=255)
   cTCGTCGCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATg                                                                    <  1:3023233/71‑1 (MQ=255)
   cTCGTCGCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATg                                                                    <  1:5214815/71‑1 (MQ=255)
      gtcgCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATg                                                                    <  1:2151589/68‑1 (MQ=255)
      gtcgCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGaa                                                                  <  1:4960442/70‑1 (MQ=255)
      gtcgCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGa                                                                   <  1:2090819/69‑1 (MQ=255)
        cgCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGAAGTg                                                               >  1:1732872/1‑71 (MQ=255)
         gCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGcc                                                                        <  1:2223772/61‑1 (MQ=255)
         gCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGcc                                                                        <  1:3290203/61‑1 (MQ=255)
         gCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGcc                                                                        <  1:1719065/61‑1 (MQ=255)
         gCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGAAGTg                                                               >  1:1320184/1‑70 (MQ=255)
         gCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGAAGTg                                                               <  1:5211219/70‑1 (MQ=255)
         gCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGAAGTg                                                               >  1:5565677/1‑70 (MQ=255)
         gCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGAAGTGc                                                              <  1:568096/71‑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)
                    tGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGAAg                                                                 <  1:2577810/57‑1 (MQ=255)
                       aGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCAGAGCGGCGCgggg                                                >  1:3037624/1‑71 (MQ=255)
                        gCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCAGAgc                                                         >  1:1679210/1‑61 (MQ=255)
                         cGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCaga                                                           >  1:2879982/1‑58 (MQ=255)
                         cGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCaga                                                           >  1:4511673/1‑58 (MQ=255)
                         cGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCaga                                                           >  1:5051213/1‑58 (MQ=255)
                           tttACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGcc                                                                        <  1:4032040/43‑1 (MQ=255)
                           tttACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGcc                                                                        <  1:1842095/43‑1 (MQ=255)
                           tttACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGcc                                                                        <  1:4685251/43‑1 (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:849431/69‑1 (MQ=255)
                                        tgtgACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCTg                                 <  1:1308500/69‑1 (MQ=255)
                                        tgtgACGAGCACGGGATCATGCTGAATGCCGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCt                                  <  1:4720948/68‑1 (MQ=255)
                                         gtgACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCagag                                                          <  1:3701008/43‑1 (MQ=255)
                                         gtgACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCagag                                                          <  1:4194455/43‑1 (MQ=255)
                                         gtgACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCagag                                                          <  1:4913888/43‑1 (MQ=255)
                                         gcgACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAgcggc                                                      >  1:4235612/3‑47 (MQ=12)
                                         gcgACGAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGttt                              >  1:4831320/3‑71 (MQ=255)
                                          tgACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCAGAGCGGCGCggg                                                 >  1:4733063/1‑51 (MQ=255)
                                          tgACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCAGAGCGGCGCggg                                                 >  1:1961799/1‑51 (MQ=255)
                                          tgACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCAGAGCGGCGCggg                                                 >  1:2040951/1‑51 (MQ=255)
                                          tgACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCAGAGCGGCGCggg                                                 >  1:532730/1‑51 (MQ=255)
                                           gatgAGCACGGTATTGTGATGATTGCTGATGAAGTGCAAAGCGGCTTTGCGCGTACCGGTAAGCTGTTTGc                            <  1:893522/68‑1 (MQ=255)
                                            aCGAGCACGGGATCATGCTGATTGCCGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCTGTTTGc                            <  1:782166/70‑1 (MQ=255)
                                            aCGAGCACGGGATCATGCTGATTGCCGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCTGTTTGc                            <  1:2704639/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)
                                                            gCTGATTGCCGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCTGTTTGCGATGGAGCAGa                 >  1:1979243/1‑65 (MQ=255)
                                                            gCTGATTGCCGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCTGTTTGCGATGGAGCAGa                 >  1:757240/1‑65 (MQ=255)
                                                            gCTGATTGCCGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCTGTTTGCGATGGAGCAGa                 >  1:243532/1‑65 (MQ=255)
                                                            gCTGATTGCCGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCTGTTTGCGATGGAGCAGATGGGCg           >  1:1003948/1‑71 (MQ=255)
                                                            gCTGATTGCCGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCTGTTTGCGATGGAGCAGATGGGCg           >  1:23092/1‑71 (MQ=255)
                                                            gCTGATTGCCGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCTGTTTGCGATGGAGCAGATGGGCg           >  1:3120713/1‑71 (MQ=255)
                                                             cTGATTGCCGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCTGTTTGCGATGGAGCAGATGGGCg           <  1:4467188/70‑1 (MQ=255)
                                                             cTGATTGCCGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCTGTTTGCGATGGAGCAGATGGGCGt          <  1:1239879/71‑1 (MQ=255)
                                                              tGATTGCCGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCTGTTTGCGa                          >  1:5492579/1‑54 (MQ=255)
                                                              tGATTGCCGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCTGTTTGCGATGGAGCAGATGGGCGtt         >  1:62374/1‑71 (MQ=255)
                                                                aTTGCCGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCTGTTTGCGATGGAGCAGATGGGCGtt         <  1:4786949/69‑1 (MQ=255)
                                                                 ttGCCGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCTGTTTGCGa                          >  1:5322686/1‑51 (MQ=255)
                                                                     cGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCTGTTTGCGATgg                       <  1:3804862/50‑1 (MQ=255)
                                                                     cGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCTGTTTGCGATgg                       <  1:4517687/50‑1 (MQ=255)
                                                                     cGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCTGTTTGCGATGGAGCAGATGGGCGTTGCGCCgg  <  1:2118988/71‑1 (MQ=255)
                                                                     cGATGAAGTGCAGAGCGGCGCGGGGAGTACCGGCACGCTGTTTGCGATGGAGCAGATGGGCGTTGCGCCgg  <  1:4446443/71‑1 (MQ=255)
                                                                     |                                                                      
CGCCTCGTCGCCAGCCTTTATGCAGCGTTTACGCGCTCTGTGTGACGAGCACGGGATCATGCTGATTGCCGATGAAGTGCAGAGCGGCGCGGGGCGTACCGGCACGCTGTTTGCGATGGAGCAGATGGGCGTTGCGCCGG  >  Exported/1792261‑1792400

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: