Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 226,889 G→C 9.1% D329H (GAT→CAT)  frsA → hydrolase, binds to enzyme IIA(Glc)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported226,8890GC 9.1% 158.6 / 6.0 55D329H (GAT→CAT) frsAhydrolase, binds to enzyme IIA(Glc)
Reads supporting (aligned to +/- strand):  ref base G (29/21);  new base C (2/3);  total (31/24)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

TGATTTTAAGTGCCAGCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTGCAAGGATTGCTGGGAC  >  Exported/226819‑226958
                                                                      |                                                                     
tGATTTTAAGTGCCAGCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCatg                                                                       <  1:513269/71‑1 (MQ=255)
  aTTTTAAGTGCCAGCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCatgat                                                                     >  1:3055332/1‑71 (MQ=255)
   ttttAAGTGCCAGCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCatgat                                                                     >  1:2625742/1‑70 (MQ=255)
          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:1931874/1‑71 (MQ=255)
                  aCAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGt                                                     >  1:3575954/1‑71 (MQ=255)
                  aCAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGt                                                     >  1:2315821/1‑71 (MQ=255)
                    aGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGAAGCATCATGCATcc                                                               <  1:566931/59‑1 (MQ=38)
                    aGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGAAGCATCATGCATcc                                                               <  1:1588636/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: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:1829862/1‑62 (MQ=255)
                                    tCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTgg                                            >  1:1456278/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: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:2335458/70‑1 (MQ=255)
                                                 gCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTaa                       <  1:827774/70‑1 (MQ=255)
                                                 gCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTaa                       >  1:349469/1‑70 (MQ=255)
                                                 gCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTaa                       <  1:2772505/70‑1 (MQ=255)
                                                   aaGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATAGCTATTCATTAAAAg                    >  1:2134162/2‑71 (MQ=255)
                                                     gTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAg                    >  1:1743050/1‑69 (MQ=255)
                                                     gTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTg                  >  1:2182875/1‑71 (MQ=255)
                                                     gTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTg                  >  1:2212472/1‑71 (MQ=255)
                                                     gTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTg                  >  1:880445/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)
                                                                    atgatgCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTa                        <  1:811573/50‑1 (MQ=255)
                                                                    atgatgCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTa                        <  1:1938323/50‑1 (MQ=255)
                                                                    atgatgCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTa                        <  1:2088903/50‑1 (MQ=255)
                                                                    atgatgCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTa                        <  1:329626/50‑1 (MQ=255)
                                                                    atgatgCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTa                        <  1:3116896/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:2428499/1‑61 (MQ=255)
                                                                      gatgCTTCCGATGAAGCTTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTGCAAGGAt           >  1:2448326/1‑61 (MQ=255)
                                                                      |                                                                     
TGATTTTAAGTGCCAGCAACAGGTGCCGGAAATGTATCTTGACGTTCTGGCCAGTCGTTTGGGGATGCATGATGCTTCCGATGAAGCGTTGCGCGTGGAGCTGAATCGCTATTCATTAAAAGTGCAAGGATTGCTGGGAC  >  Exported/226819‑226958

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: