Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 226,884 T→A 8.1% M327K (ATG→AAG)  frsA → hydrolase, binds to enzyme IIA(Glc)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported226,8840TA 8.1% 170.6 / 4.9 60M327K (ATG→AAG) frsAhydrolase, binds to enzyme IIA(Glc)
Reads supporting (aligned to +/- strand):  ref base T (33/22);  new base A (2/3);  total (35/25)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

TGATTTTAAGTGCCAGCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTGCAAGGATTG  >  Exported/226819‑226951
                                                                 |                                                                   
tGATTTTAAGTGCCAGCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCatg                                                                <  1:513269/71‑1 (MQ=255)
tGATTTTAAGTGCCAGCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCat                                                                 >  1:124862/1‑70 (MQ=255)
tGATTTTAAGTGCCAGCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCat                                                                 >  1:3640510/1‑70 (MQ=255)
  aTTTTAAGTGCCAGCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCatgat                                                              >  1:3055332/1‑71 (MQ=255)
   ttttAAGTGCCAGCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCatgat                                                              >  1:2625742/1‑70 (MQ=255)
          tGCCAGCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGa                                                      >  1:3584917/1‑71 (MQ=255)
          tGCCAGCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGa                                                      <  1:2590749/71‑1 (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:3386003/1‑66 (MQ=255)
                  aCAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGa                                                   >  1:3642173/1‑66 (MQ=255)
                  aCAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGt                                              >  1:1931874/1‑71 (MQ=255)
                  aCAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGt                                              >  1:3575954/1‑71 (MQ=255)
                  aCAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGt                                              >  1:2315821/1‑71 (MQ=255)
                    aGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGAAGCATCATGCATcc                                                        <  1:1588636/59‑1 (MQ=38)
                    aGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGAAGCATCATGCATcc                                                        <  1:566931/59‑1 (MQ=38)
                    aGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGAAGCATCATGCATcc                                                        <  1:439593/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:1379236/61‑1 (MQ=255)
                           ggAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCg                                               <  1:2792925/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:1829862/1‑62 (MQ=255)
                                    tCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTgg                                     >  1:1456278/1‑62 (MQ=255)
                                          cGTTCTGGCCAGTCGTTTGGGGATGCATGATGATTCCGATGAAGCGt                                              <  1:2527099/47‑1 (MQ=255)
                                           gTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTc                     <  1:1050072/71‑1 (MQ=255)
                                           gTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTc                     <  1:2502745/71‑1 (MQ=255)
                                            ttCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCg                                               <  1:2814034/44‑1 (MQ=255)
                                            ttCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCg                                               <  1:3092831/44‑1 (MQ=255)
                                              cTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGc                                  >  1:911384/1‑55 (MQ=255)
                                                 tccAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTa                 <  1:2131884/68‑1 (MQ=255)
                                                 gCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTaa                <  1:2772505/70‑1 (MQ=255)
                                                 gCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTaa                <  1:2335458/70‑1 (MQ=255)
                                                 gCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTaa                >  1:349469/1‑70 (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:2212472/1‑71 (MQ=255)
                                                     gTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTg           >  1:880445/1‑71 (MQ=255)
                                                     gTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTg           >  1:2182875/1‑71 (MQ=255)
                                                     gTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTg           >  1:427046/1‑71 (MQ=255)
                                                     gTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCCTGGAg                                   >  1:269663/1‑47 (MQ=255)
                                                      tCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTg           >  1:794893/1‑70 (MQ=255)
                                                      tCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTg           >  1:1706128/1‑70 (MQ=255)
                                                      tCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTg           >  1:513288/1‑70 (MQ=255)
                                                          ttGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTGCAAg       >  1:597724/1‑70 (MQ=255)
                                                           tGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAg             >  1:840908/1‑63 (MQ=255)
                                                           tGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAg             >  1:1657831/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)
                                                                 |                                                                   
TGATTTTAAGTGCCAGCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTGCAAGGATTG  >  Exported/226819‑226951

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: