Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 1,364,882 G→A 5.9% V142V (GTG→GTA puuC → gamma‑Glu‑gamma‑aminobutyraldehyde dehydrogenase, NAD(P)H‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090481,364,8820GA 5.9% 224.7 / 2.9 68V142V (GTG→GTApuuCgamma‑Glu‑gamma‑aminobutyraldehyde dehydrogenase, NAD(P)H‑dependent
Reads supporting (aligned to +/- strand):  ref base G (37/27);  new base A (2/2);  total (39/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 = 6.61e-01

TGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGCGTGATTGCCGCCA  >  AP009048/1364813‑1364946
                                                                     |                                                                
tgatATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTg                                                                  >  1:3703456/1‑70 (MQ=255)
  atatTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTAGc                                                                <  1:4391032/70‑1 (MQ=255)
   tatTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGa                                                              >  1:4165981/1‑71 (MQ=255)
     ttCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGac                                                             <  1:1869418/70‑1 (MQ=255)
      tCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGacca                                                           >  1:69444/1‑71 (MQ=255)
      tCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGacca                                                           >  1:3646034/1‑71 (MQ=255)
         cggcgcggcgCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGaccac                                                          >  1:198570/1‑69 (MQ=255)
          ggcgcggcgcGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGAccaccac                                                       >  1:3707965/1‑70 (MQ=255)
           gcgcggcgcgCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAg                                                       >  1:1070406/1‑70 (MQ=255)
           gcgcggcgcgCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGt                                                      <  1:339453/71‑1 (MQ=255)
           gcgcggcgcgCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGt                                                      <  1:2845488/71‑1 (MQ=255)
                 cgcgcgcCATTCGCTTGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCAt                                                <  1:1457736/71‑1 (MQ=255)
                 cgcgcgcCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCAt                                                >  1:1341225/1‑71 (MQ=255)
                 cgcgcgcCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCAt                                                <  1:2581137/71‑1 (MQ=255)
                   cgcgcCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCAt                                                >  1:2818056/1‑69 (MQ=255)
                   cgcgcCATTCGCTGGTACGCAGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGa                                              <  1:38094/71‑1 (MQ=255)
                        cATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTg                                          >  1:535933/1‑70 (MQ=255)
                          ttCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGac                                                             <  1:1086785/49‑1 (MQ=255)
                          ttCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGc                                        >  1:4384397/1‑70 (MQ=255)
                          ttCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGc                                        >  1:1854152/1‑70 (MQ=255)
                          ttCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCg                                       <  1:937086/71‑1 (MQ=255)
                          ttCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCg                                       <  1:1759536/71‑1 (MQ=255)
                           tCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCga                                      >  1:373525/1‑71 (MQ=255)
                            cGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTAGCGAccacca                                                        >  1:431214/1‑52 (MQ=255)
                             gCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCgat                                     <  1:2054728/70‑1 (MQ=255)
                             gCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCgat                                     >  1:4586267/1‑70 (MQ=255)
                             gCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCgat                                     >  1:2023945/1‑70 (MQ=255)
                               tGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCgatgat                                  >  1:663505/1‑71 (MQ=255)
                                  tACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATcgt                               <  1:983435/71‑1 (MQ=255)
                                  tACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATcgt                               <  1:555095/71‑1 (MQ=255)
                                    cGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATcgtgc                             >  1:4231604/1‑71 (MQ=255)
                                      ccGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGc                                                   <  1:4059824/47‑1 (MQ=255)
                                      ccGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCAt                                                <  1:3980730/50‑1 (MQ=255)
                                      ccGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATcgtgcgt                           >  1:2582483/1‑71 (MQ=255)
                                        gAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATcgtgcgtg                          >  1:3590011/1‑70 (MQ=255)
                                          aGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAAc                       >  1:4289279/1‑71 (MQ=255)
                                          aGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAAc                       >  1:4284030/1‑71 (MQ=255)
                                          aGCGATCGACAAAGTGTATGGCCAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAAc                       <  1:1345343/71‑1 (MQ=255)
                                           gCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCgatg                                    <  1:118976/57‑1 (MQ=255)
                                             gATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCt                                           <  1:3255352/48‑1 (MQ=255)
                                              aTCGACAAAGTGTATGGCGAAGTAGCGACCACCAGTAGCCATGAGCTGGCGATGATcgtg                              <  1:1189815/60‑1 (MQ=255)
                                               tCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGt                   >  1:3879616/1‑70 (MQ=255)
                                                cGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTc                  >  1:3218741/1‑70 (MQ=255)
                                                cGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCg                 >  1:1744948/1‑71 (MQ=255)
                                                cGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCg                 >  1:3643089/1‑71 (MQ=255)
                                                 gACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAg                                             >  1:1487074/1‑42 (MQ=255)
                                                 gACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGt                   <  1:4519072/68‑1 (MQ=255)
                                                 gACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCgg                <  1:1225278/71‑1 (MQ=255)
                                                  aCAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATcgtg                              >  1:890369/1‑56 (MQ=255)
                                                  aCAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGc               <  1:1596907/71‑1 (MQ=255)
                                                   cAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGa                         <  1:1653350/60‑1 (MQ=255)
                                                   cAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGCg              >  1:3584054/1‑71 (MQ=255)
                                                   cAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGCg              >  1:2973669/1‑71 (MQ=255)
                                                    aaaGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCg                     >  1:3017239/1‑63 (MQ=255)
                                                    aaaGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCTGCGt             <  1:817956/71‑1 (MQ=255)
                                                    aaaGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGCGt             >  1:2333458/1‑71 (MQ=255)
                                                       gtgtATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGCGTGAt          >  1:801925/1‑71 (MQ=255)
                                                        tgtATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGCGTGAtt         <  1:3550479/71‑1 (MQ=255)
                                                         gtATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCgatga                                   <  1:87913/44‑1 (MQ=255)
                                                           aTGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGCGTGATTgcc      >  1:1898447/1‑71 (MQ=255)
                                                            tGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATcgtg                              >  1:2606887/1‑46 (MQ=255)
                                                            tGGCGAAGTAGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTTCGTGAACCGGTCGGCGTGATTgccg     >  1:2355535/1‑71 (MQ=255)
                                                              gCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGCGTGAt          >  1:3517486/1‑64 (MQ=255)
                                                              gCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGCGTGATTgccgc    >  1:685345/1‑70 (MQ=255)
                                                              gCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGCGTGATTgccgc    >  1:318675/1‑70 (MQ=255)
                                                               cGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGCGTGATTGCCGCCa  >  1:2063246/1‑71 (MQ=255)
                                                               cGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGCGTGATTGCCGCCa  <  1:4274172/71‑1 (MQ=255)
                                                                    tGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGCGTGATTgcc      <  1:246374/62‑1 (MQ=255)
                                                                    tGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGCGTGATTg        <  1:1807882/60‑1 (MQ=255)
                                                                     |                                                                
TGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGCGTGATTGCCGCCA  >  AP009048/1364813‑1364946

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: