Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,778,494 C→T 8.2% pseudogene (67/639 nt) glvG → ECK3673:JW3658:b3681; predicted 6‑phospho‑beta‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,778,4940CT 8.2% 191.0 / 4.8 61pseudogene (67/639 nt)glvGECK3673:JW3658:b3681; predicted 6‑phospho‑beta‑glucosidase
Reads supporting (aligned to +/- strand):  ref base C (25/31);  new base T (3/2);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGACCAAATTCTCAGTGGTTGTCGCAGGCGGTGGAAGCACCTTTACGCCAGGCATCGTGTTGATGCTCCTGGCGAATCAGGACCGTTTCCCGCTTCGTGCACTGAAATTTTATGATAACGATGGTGCGCGGC  >  AP009048/3778427‑3778560
                                                                   |                                                                  
aaTGACCAAATTCTCAGTGGTTGTCGCAGGCGGTGGAAGCACCTTTACGCCAGGCATCGTGTTGATGCTcc                                                                 >  1:1898692/1‑71 (MQ=255)
 aTGACCAAATTCTCAGTGGTTGTCGCAGGCGGTGGAAGCACCTTTACGCCAGGCATCGTGTTGATGCTCCt                                                                <  1:2736713/71‑1 (MQ=255)
    aCCAAATTCTCAGTGGTTGTCGCAGGCGGTGGAAGCACCTTTACGCCAGGCATCGTGTTGATGCTCCTGGc                                                             >  1:309695/1‑71 (MQ=255)
      cAAATTCTCAGTGGTTGTCGCAGGCGGTGGAAGCACCTTTACGCCAGGCATCGTGTTGATGCTCCTg                                                               <  1:3843516/67‑1 (MQ=255)
      cAAATTCTCAGTGGTTGTCGCAGGCGGTGGAAGCACCTTTACGCCAGGCATCGTGTTGATGCTCCTGGCGa                                                           <  1:1388653/71‑1 (MQ=255)
       aaaTTCTCAGTGGTTGTCGCAGGCGGTGGAAGCACCTTTACGCCAGGCATCGTGTTGATGCTCCTGGCGaa                                                          >  1:3870015/1‑71 (MQ=255)
        aaTTCTCAGTGGTTGTCGCAGGCGGTGGAAGCACCTTTACGCCAGGCATCGTGTTGATGCTCCTGGCGaa                                                          <  1:527602/70‑1 (MQ=255)
        aaTTCTCAGTGGTTGTCGCAGGCGGTGGAAGCACCTTTACGCCAGGCATCGTGTTGATGCTCCTGGCGAAt                                                         >  1:1408686/1‑71 (MQ=255)
        aaTTCTCAGTGGTTGTCGCAGGCGGTGGAAGCACCTTTACGCCAGGCATCGTGTTGATGCTACTGGCGAAt                                                         >  1:455130/1‑71 (MQ=255)
         aTTCTCAGTGGTTGTCGCAGGCGGTGGAAGCACCTTTACGCCAGGCATCGTGTTGATGCTCCTGGCGAATc                                                        <  1:1266428/71‑1 (MQ=255)
         aTTCTCAGTGGTTGTCGCAGGCGGTGGAAGCACCTTTACGCCAGGCATCGTGTTGATGCTCCTGGCGAATc                                                        <  1:3561705/71‑1 (MQ=255)
            ctcAGTGGTTGTCGCAGGCGGTGGAAGCACCTTTACGCCAGGCATCGTGTTGATGCTCCTGGCGAATCAg                                                      <  1:658843/70‑1 (MQ=255)
            ctcAGTGGTTGTCGCAGGCGGTGGAAGCACCTTTACGCCAGGCATCGTGTTGATGCTCCTGGCGAATCAg                                                      <  1:201951/70‑1 (MQ=255)
               aGTGGTTGTCGCAGGCGGTGGAAGCACCTTTACGCCAGGCATCGTGTTGATGCTCCTGGCGAATCAGGa                                                    <  1:2558937/69‑1 (MQ=255)
                gTGGTTGTCGCAGGCGGTGGAAGCACCTTTACGCCAGGCATCGTGTTGATGCTCCTGGCGAATCAGGACCg                                                 <  1:1724931/71‑1 (MQ=255)
                 tGGTTGTCGCAGGCGGTGGAAGCACCTTTACGCCAGGCATCGTGTTGATGCTCCTGGCGAATCAGGACCGt                                                <  1:2653779/71‑1 (MQ=255)
                 tGGTTGTCGCAGGCGGTGGAAGCACCTTTACGCCAGGCATCGTGTTGATGCTCCTGGCGAATCAGGACCGt                                                >  1:2516170/1‑71 (MQ=255)
                   gTTGTCGCAGGCGGTGGAAGCACCTTTACGCCAGGCATCGTGTTGATGCTCCTGGCGaa                                                          <  1:4154274/59‑1 (MQ=255)
                   gTTGTCGCAGGCGGTGGAAGCACCTTTACGCCAGGCATCGTGTTGATGCTCCTGGCGAATCAGGACCGtt                                               <  1:1727721/70‑1 (MQ=255)
                   gTTGTCGCAGGCGGGGGAAGCACCTTTACGCCAGGCATCGTGTTGATGCTCCTGGCGAATCAGGACCGttt                                              <  1:2156077/71‑1 (MQ=255)
                     tGTCGCAGGCGGTGGAAGCACCTTTACGCCAGGCATCGTGTTGATGCTCCTGGCGAATCAGGACCGTTTcc                                            >  1:3544446/1‑71 (MQ=255)
                      gTCGCAGGCGGTGGAAGCACCTTTACGCCAGGCATCGTGTTGATGCTCCTGGCGAATCAGGACCGTTTcc                                            >  1:3129210/1‑70 (MQ=255)
                       tCGCAGGCGGTGGAAGCACCTTTACGCCAGGCATCGTGTTGATGCTcc                                                                 <  1:545278/48‑1 (MQ=255)
                        cGCAGGCGGTGGAAGCACCTTTACGCCAGGCATCGTGTTGATGCTCCTGGCGAATCAgg                                                     >  1:1596985/1‑59 (MQ=255)
                         gCAGGCGGTGGAAGCACCTTTACGCCAGGCATCGTGTTGATGCTCCTGGCGAATCAGGACCGTTTCCCGCt                                        <  1:1310193/71‑1 (MQ=255)
                           aGGCGGTGGAAGCACCTTTACGCCAGGCATCGTGTTGATGCTCCTGGCGAATCAGGACCGTTTCCCGCTTc                                      <  1:2140554/71‑1 (MQ=255)
                             gCGGTGGAAGCACCTTTACGCCAGGCATCGTGTTGATGCTCCTGGCGAATCAGGACCGTTTCCCGCTTCGt                                    >  1:1211236/1‑71 (MQ=255)
                               ggTGGAAGCACCTTTACGCCAGGCATCGTGTTGATGCTCCTGGCGAATCAGGACCGtt                                               >  1:1105490/1‑58 (MQ=255)
                                  ggAAGCACCTTTACGCCAGGCATCGTGTTGATGCTcc                                                                 <  1:2058395/37‑1 (MQ=255)
                                  ggAAGCACCTTTACGCCAGGCATCGTGTTGATGCTCCTGGCGAATCAGGAc                                                   <  1:4479620/51‑1 (MQ=255)
                                    aaGCACCTTTACGCCAGGCATCGTGTTGATGCTCCTGGCGAATCAGGACCGTTTCCCGCTTCGTGCa                                 <  1:2004886/67‑1 (MQ=255)
                                     aGCACCTTTACGCCAGGCATCGTGTTGATGCTCCTGGCGAATCAGGACCGTTTCCCGCTTCGTGCACTGa                             >  1:2914553/1‑70 (MQ=255)
                                        aCCTTTACGCCAGGCATCGTGTTGATGTTCCTGGCGAATCAGGACCGTTTCCCGCTTCGTGCACTg                              <  1:3323331/66‑1 (MQ=255)
                                        aCCTTTACGCCAGGCATCGTGTTGATGTTCCTGGCGAATCAGGACCGTTTCCCGCTTCGTGCACTGAAAtt                         <  1:298745/71‑1 (MQ=255)
                                         ccTTTACGCCAGGCATCGTGTTGATGCTCCTGGCGAATCAGGACCGTTTCCCGCTTCGTGCACTGAAAttt                        <  1:3965243/71‑1 (MQ=255)
                                          cTTTACGCCAGGCATCGTGTTGATGCTCCTGGCGAATCAGGACCGTTTCCCGCTTCGTGCACTGAAAttt                        >  1:2220093/1‑70 (MQ=255)
                                           tttACGCCAGGCATCGTGTTGATGCTCCTGGCGAATCAGGACCGt                                                <  1:3642453/45‑1 (MQ=255)
                                            ttACGCCAGGCATCGTGTTGATGCTCCTGGCGAATCAGGACCGTTTCCCGCTTCGTGCACTGAAATTTTa                      >  1:1684361/1‑70 (MQ=255)
                                             tACGCCAGGCATCGTGTTGATGTTCCTGGCGAATCAGGACCGTTTCCCGCTTCGTGCACTGAAATTTTAt                     >  1:2612496/1‑70 (MQ=255)
                                             tACGCCAGGCATCGTGTTGATGTTCCTGGCGAATCAGGACCGTTTCCCGCTTCGTGCACTGAAATTTTAt                     >  1:1684566/1‑70 (MQ=255)
                                             tACGCCAGGCATCGTGTTGATGCTCCTGGCGAATCAGGACCGTTTCCCGCTTCGTGCACTGAAATTTTa                      <  1:475719/69‑1 (MQ=255)
                                             tACGCCAGGCATCGTGTTGATGCTCCTGGCGAATCAGGACCGTTTCCCGCTTCGTGCACTGAAATTTTAt                     >  1:3709366/1‑70 (MQ=255)
                                              aCGCCAGGCATCGTGTTGATGCTCCTGGCGAATCAGGACCGTTTCCCGCTTCGTGCACTGAAATTTTATg                    <  1:2766333/70‑1 (MQ=255)
                                                  cAGGCATCGTGTTGATGCTCCTGGCGAATCAGGACCGTTTCCCGCtt                                       >  1:2817744/1‑47 (MQ=255)
                                                  cAGGCATCGTGTTGATGCTCCTGGCGAATCAGGACCGTTTCCCGCTTCGTGCACTGAAATTTTATGATAAc               <  1:497414/71‑1 (MQ=255)
                                                   aGGCATCGTGTTGATGCTCCTGGCGAATCAGGACCGTTTCCCGCTTCg                                     >  1:3385392/1‑48 (MQ=255)
                                                    ggCATCGTGTTGATGCTCCTGGCGAATCAGGACCGTTTCCCGCTTCGTGCACTGAAATTTTATGATAACGa             <  1:1939853/71‑1 (MQ=255)
                                                    ggCATCGTGTTGATGCTCCTGGCGAATCAGGACCGTTTCCCGCTTCGTGCACTGAAATTTTATGATAACGa             >  1:1483549/1‑71 (MQ=255)
                                                    ggCATCGTGTTGATGCTCCTGGCGAATCAGGACCGTTTCCCGCTTCGTGCACTGAAATTTTATGATAACGa             >  1:2798170/1‑71 (MQ=255)
                                                     gCATCGTGTTGATGCTCCTGGCGATTCAGGACCGTTTCCCGCTTCGTGCACTGAAATTTTATGATAACGAt            <  1:4110421/71‑1 (MQ=255)
                                                     gCATCGTGTTGATGCTCCTGGCGAATCAGGACCGTTTCCCGCTTCGTGCACTGAAATTTTATGATAACGAt            >  1:1595087/1‑71 (MQ=255)
                                                        tCGTGTTGATGCTCCTGGCGAATCAGGACCGTTTCCCGCTTCGTg                                   <  1:1581328/45‑1 (MQ=255)
                                                        tCGTGTTGATGCTCCTGGCGAATCAGGACCGTTTCCCGCTTCGTGCACTGAAATTTTATGATAACGATGGt         >  1:1656647/1‑71 (MQ=255)
                                                          gtgtTGATGCTCCTGGCGAATCAGGAGCGTTTCCCGCTTCGTGCACTGAAATTTTATGATAACGATGGTg        <  1:2463343/70‑1 (MQ=255)
                                                          gtgtTGATGCTCCTGGCGAATCAGGACCGTTTCCCGCTTCGTGCACTGAAATTTTa                      >  1:3379259/1‑56 (MQ=255)
                                                             ttGATGCTCCTGGCGAATCAGGACCGTTTCCCGCTTCGTGCACTGAAATTTTATGATAACGATGGTgcgc     <  1:4470799/70‑1 (MQ=255)
                                                              tGATGCTCCTGGCGAATCAGGACCGTTTCCCGCTTCGTGCACTGAAATTTTATGATAACGATgg          >  1:2154295/1‑64 (MQ=255)
                                                              tGATGCTCCTGGCGAATCAGGACCGTTTCCCGCTTCGTGCACTGAAATTTTATGATAACGATGGTGcgcg    >  1:1835135/1‑70 (MQ=255)
                                                               gATGTTCCTGGCGAATCAGGACCGTTTCCCGCTTCGTGCACTGAAATTTTATGATAACGATGGTGCGCgg   >  1:3651446/1‑70 (MQ=255)
                                                               gATGCTCCTGGCGAATCAGGACCGTTTCCCGCTTCGTGCACTGAAATTTTATGATAACGATGGTGCGCgg   >  1:382990/1‑70 (MQ=255)
                                                                aTGCTCCTGGCGAATCAGGACCGTTTCCCGCTTCGTGCACTGAAATTTTATGATAACGATGGTGCGCGGc  <  1:3643251/70‑1 (MQ=255)
                                                                   |                                                                  
AATGACCAAATTCTCAGTGGTTGTCGCAGGCGGTGGAAGCACCTTTACGCCAGGCATCGTGTTGATGCTCCTGGCGAATCAGGACCGTTTCCCGCTTCGTGCACTGAAATTTTATGATAACGATGGTGCGCGGC  >  AP009048/3778427‑3778560

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