Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 226,894 T→A 7.3% A330A (GCT→GCA frsA → hydrolase, binds to enzyme IIA(Glc)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported226,8940TA 7.3% 195.8 / 6.5 68A330A (GCT→GCAfrsAhydrolase, binds to enzyme IIA(Glc)
Reads supporting (aligned to +/- strand):  ref base T (33/30);  new base A (2/3);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

TGCCAGCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTGCAAGGATTGCTGGGACGTCGCT  >  Exported/226829‑226964
                                                                 |                                                                      
tGCCAGCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGa                                                                   <  1:2590749/71‑1 (MQ=255)
tGCCAGCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGa                                                                   >  1:3584917/1‑71 (MQ=255)
 gCCAGCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGAt                                                                  >  1:995053/1‑71 (MQ=255)
    aGCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGaa                                                               <  1:1438712/71‑1 (MQ=255)
     gCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAg                                                              >  1:2315251/1‑71 (MQ=255)
     gCAACACGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAg                                                              >  1:3582979/1‑71 (MQ=255)
       aaCAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCg                                                            >  1:2155474/1‑71 (MQ=255)
        aCAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATg                                                                 >  1:3080399/1‑65 (MQ=255)
        aCAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGa                                                                >  1:3642173/1‑66 (MQ=255)
        aCAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGa                                                                >  1:3386003/1‑66 (MQ=255)
        aCAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGt                                                           >  1:3575954/1‑71 (MQ=255)
        aCAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGt                                                           >  1:1931874/1‑71 (MQ=255)
        aCAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGt                                                           >  1:2315821/1‑71 (MQ=255)
          aGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGAAGCATCATGCATcc                                                                     <  1:1588636/59‑1 (MQ=38)
          aGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGAAGCATCATGCATcc                                                                     <  1:439593/59‑1 (MQ=38)
          aGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGAAGCATCATGCATcc                                                                     <  1:566931/59‑1 (MQ=38)
             tGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTgcgc                                                      <  1:1669603/71‑1 (MQ=255)
               ccGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGt                                                    <  1:2830726/71‑1 (MQ=255)
                cGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTg                                                   <  1:26798/71‑1 (MQ=255)
                 ggAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCg                                                            <  1:3401596/61‑1 (MQ=255)
                 ggAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCg                                                            <  1:2792925/61‑1 (MQ=255)
                 ggAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCg                                                            <  1:1379236/61‑1 (MQ=255)
                      tGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTg                                             <  1:3258260/71‑1 (MQ=255)
                        tATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGt                                                    <  1:805079/62‑1 (MQ=255)
                        tATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGaa                                           <  1:2627397/71‑1 (MQ=255)
                         aTCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGt                                                    <  1:1546280/61‑1 (MQ=255)
                          tCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTgg                                                  >  1:1456278/1‑62 (MQ=255)
                          tCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTgg                                                  >  1:1829862/1‑62 (MQ=255)
                                cGTTCTGGCCAGTCGTTTGGGGATGCATGATGATTCCGATGAAGCGt                                                           <  1:2527099/47‑1 (MQ=255)
                                 gTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTc                                  <  1:2502745/71‑1 (MQ=255)
                                 gTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTc                                  <  1:1050072/71‑1 (MQ=255)
                                  ttCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCg                                                            <  1:3092831/44‑1 (MQ=255)
                                  ttCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCg                                                            <  1:2814034/44‑1 (MQ=255)
                                    cTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGc                                               >  1:911384/1‑55 (MQ=255)
                                       tccAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTa                              <  1:2131884/68‑1 (MQ=255)
                                       gCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTaa                             <  1:2335458/70‑1 (MQ=255)
                                       gCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTaa                             >  1:349469/1‑70 (MQ=255)
                                       gCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTaa                             <  1:2772505/70‑1 (MQ=255)
                                       gCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTaa                             <  1:827774/70‑1 (MQ=255)
                                         aaGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATAGCTATTCATTAAAAg                          >  1:2134162/2‑71 (MQ=255)
                                           gTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAg                          >  1:1743050/1‑69 (MQ=255)
                                           gTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTg                        >  1:427046/1‑71 (MQ=255)
                                           gTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTg                        >  1:2212472/1‑71 (MQ=255)
                                           gTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTg                        >  1:880445/1‑71 (MQ=255)
                                           gTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTg                        >  1:2182875/1‑71 (MQ=255)
                                           gTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCCTGGAg                                                >  1:269663/1‑47 (MQ=255)
                                            tCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTg                        >  1:1706128/1‑70 (MQ=255)
                                            tCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTg                        >  1:794893/1‑70 (MQ=255)
                                            tCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTg                        >  1:513288/1‑70 (MQ=255)
                                                ttGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTGCAAg                    >  1:597724/1‑70 (MQ=255)
                                                 tGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAg                          >  1:1657831/1‑63 (MQ=255)
                                                 tGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAg                          >  1:840908/1‑63 (MQ=255)
                                                 tGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTGCAAGGa                  >  1:1645794/1‑71 (MQ=255)
                                                    ggAAGCATCATGCATCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTGCAAGGATTg               >  1:2191541/1‑71 (MQ=255)
                                                    ggAAGCATCATGCATCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTGCAAGGATTg               >  1:973165/1‑71 (MQ=255)
                                                          atgatgCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTa                              <  1:1938323/50‑1 (MQ=255)
                                                          atgatgCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTa                              <  1:2088903/50‑1 (MQ=255)
                                                          atgatgCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTa                              <  1:811573/50‑1 (MQ=255)
                                                          atgatgCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTa                              <  1:3116896/50‑1 (MQ=255)
                                                          atgatgCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTa                              <  1:329626/50‑1 (MQ=255)
                                                          atgatgCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTGCAAGGATTGc              <  1:981095/66‑1 (MQ=255)
                                                          atgatgCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTGCAAGGATTGCTGGGa         >  1:854481/1‑71 (MQ=255)
                                                           tgatgCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTGCAAGGATTGCTGGGAc        <  1:929066/71‑1 (MQ=255)
                                                            gatgCTTCCGATGAAGCTTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTGCAAGGAt                 >  1:2448326/1‑61 (MQ=255)
                                                            gatgCTTCCGATGAAGCTTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTGCAAGGAt                 >  1:2428499/1‑61 (MQ=255)
                                                             atgCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTaaaa                           >  1:1420654/1‑50 (MQ=255)
                                                             atgCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTGCAAGGATTGCTGGGACg       <  1:88352/70‑1 (MQ=255)
                                                                cTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTGCAAGGATTGCTGGGACGTCGc   <  1:482781/71‑1 (MQ=255)
                                                                 ttCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTGCAAGGATTGCTg            <  1:423760/61‑1 (MQ=255)
                                                                 ttCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTGCAAGGATTGCTg            <  1:1593839/61‑1 (MQ=255)
                                                                 ttCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTGCAAGGATTGCTGGGACGTCGCt  >  1:3681248/1‑71 (MQ=255)
                                                                 |                                                                      
TGCCAGCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTGCAAGGATTGCTGGGACGTCGCT  >  Exported/226829‑226964

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: