Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,521,572 G→A 14.0% L157L (CTG→CTA glpX → fructose 1,6‑bisphosphatase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,521,5720GA14.0% 173.5 / 17.2 64L157L (CTG→CTAglpXfructose 1,6‑bisphosphatase II
Reads supporting (aligned to +/- strand):  ref base G (35/20);  new base A (4/5);  total (39/25)
Fisher's exact test for biased strand distribution p-value = 2.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TGAACCTGCCGCTGGCGGATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCAACTCGGCGTA  >  AP009048/3521505‑3521638
                                                                   |                                                                  
tGAACCTGCCGCTGGCGGATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACg                                                                 >  1:4117628/1‑71 (MQ=255)
  aaCCTGCCGCTGGCGGATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACgg                                                                >  1:2935567/1‑70 (MQ=255)
  aaCCTGCCGCTGGCGGATAACCTGCGCAATGTAGCGGCGGCGCTCGGAAAACCGTTGAGCGAACTGACgg                                                                >  1:4122443/1‑70 (MQ=255)
   aCCTGCCGCTGGCGGATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTAACGGTa                                                              <  1:3933328/71‑1 (MQ=255)
     cTGCCGCTGGCGGATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTAACGGTAAc                                                            <  1:2486905/71‑1 (MQ=255)
         cGCTGGCGGATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGAtt                                                        >  1:1495199/1‑71 (MQ=255)
           cTGGCGGATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCt                                                      <  1:3477276/71‑1 (MQ=255)
             ggCGGATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTgg                                                    >  1:794539/1‑71 (MQ=255)
              gCGGATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTgg                                                    >  1:4440279/1‑70 (MQ=255)
                ggATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTa                                                 <  1:2646022/71‑1 (MQ=255)
                ggATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTa                                                 <  1:290807/71‑1 (MQ=255)
                 gATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTaa                                                <  1:2468480/71‑1 (MQ=255)
                 gATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTaa                                                <  1:1821530/71‑1 (MQ=255)
                    aaCCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTAACGGTAACGATTCTg                                                     >  1:96570/1‑63 (MQ=255)
                     aCCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAAcca                                            >  1:1477781/1‑71 (MQ=255)
                        tGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAAccacgc                                         >  1:263074/1‑71 (MQ=255)
                        tGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAAccacgc                                         >  1:111196/1‑71 (MQ=255)
                         gcgcAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGAtt                                                        >  1:214103/1‑55 (MQ=255)
                         gcgcAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTAACGGTAAc                                                            <  1:4402551/51‑1 (MQ=255)
                          cgcAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAAccacgcc                                        <  1:1104704/70‑1 (MQ=255)
                           gcAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAAccacgccac                                      <  1:2973455/71‑1 (MQ=255)
                               tGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCt                                                  >  1:154530/1‑55 (MQ=255)
                               tGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATg                                  >  1:2664155/1‑71 (MQ=255)
                                gTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGAc                                                                  >  1:1387786/1‑38 (MQ=255)
                                gTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAAcca                                            >  1:1176647/1‑60 (MQ=255)
                                   gcggcggcgCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGt                              >  1:1476100/1‑71 (MQ=255)
                                    cggcggcgCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTaa                                                <  1:1492735/52‑1 (MQ=255)
                                    cggcggcgCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGtt                             <  1:4428286/71‑1 (MQ=255)
                                     ggcggcgCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCt                                                      <  1:1420346/45‑1 (MQ=255)
                                      gcggcgCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTaaa                                               >  1:4332671/1‑51 (MQ=255)
                                       cggcgCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTAt                           >  1:684268/1‑70 (MQ=255)
                                        ggcgCTCGGCAAACCGTTGAGCGAACTGCCGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCg                         >  1:2959971/1‑71 (MQ=255)
                                         gcgcTCGGCAAACCGTTGAGCGAACTAACGGTAACGATTCTGGCTAAACCACGCCACGATg                                  <  1:3486864/61‑1 (MQ=255)
                                          cgcTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCt                       >  1:3465992/1‑71 (MQ=255)
                                           gcTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTg                                                     >  1:893384/1‑40 (MQ=255)
                                           gcTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCg                         <  1:3199221/68‑1 (MQ=255)
                                           gcTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCg                         <  1:3979360/68‑1 (MQ=255)
                                           gcTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCt                       >  1:4436518/1‑70 (MQ=255)
                                           gcTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCt                       >  1:1996117/1‑70 (MQ=255)
                                           gcTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTg                      <  1:3249583/71‑1 (MQ=255)
                                            cTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGa                     <  1:524300/71‑1 (MQ=255)
                                             tCGGCAAACCGTTGAGCGAACTAACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGa                     >  1:452432/1‑70 (MQ=255)
                                               ggCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAAt                  >  1:1400573/1‑71 (MQ=255)
                                               ggCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAAt                  >  1:3064762/1‑71 (MQ=255)
                                               ggCAAACCGTTGAGCGAACTAACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGaaa                   >  1:3845019/1‑70 (MQ=255)
                                                gCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATg                 <  1:85417/71‑1 (MQ=255)
                                                 cAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATg                 >  1:1043659/1‑70 (MQ=255)
                                                    aCCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATgcag              >  1:4330867/1‑70 (MQ=255)
                                                     ccGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATgcagca            <  1:4051727/71‑1 (MQ=255)
                                                      cGTTGAGCGAACTGACGGTAACGATTCTGGCTAAAccac                                           >  1:1387208/1‑39 (MQ=255)
                                                        ttGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTa                            <  1:2796129/52‑1 (MQ=255)
                                                         tGAGCGAACTAACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCAACTc        >  1:4178351/1‑71 (MQ=255)
                                                          gAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCg                               >  1:4411042/1‑47 (MQ=255)
                                                          gAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATGGCTGAAATGCAGCAACTCg       >  1:4152280/1‑71 (MQ=255)
                                                          gAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCg                         >  1:3376046/1‑53 (MQ=255)
                                                          gAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCg                         >  1:2747611/1‑53 (MQ=255)
                                                           aGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCAACTCg       <  1:1921430/70‑1 (MQ=255)
                                                            gCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCAACTCGGc     <  1:692506/71‑1 (MQ=255)
                                                            gCGAACTAACGGTAACGATTCTGGCTAAAccacgccac                                      <  1:3795895/38‑1 (MQ=255)
                                                              gAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCAACTCGGCg    >  1:4284084/1‑70 (MQ=255)
                                                              gAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCAACTCGGCg    >  1:4067976/1‑70 (MQ=255)
                                                               aaCTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCAACTCGGCGTa  >  1:134244/1‑71 (MQ=255)
                                                               aaCTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCAACTCGGCGTa  <  1:3662380/71‑1 (MQ=255)
                                                                 cTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCAACTCGGCGTa  >  1:2631815/1‑69 (MQ=255)
                                                                   |                                                                  
TGAACCTGCCGCTGGCGGATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCAACTCGGCGTA  >  AP009048/3521505‑3521638

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: