Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,968,626 A→G 6.7% S2P (TCC→CCC)  tktA ← transketolase 1, thiamin‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,968,6260AG 6.7% 401.1 / 15.7 134S2P (TCC→CCC) tktAtransketolase 1, thiamin‑binding
Reads supporting (aligned to +/- strand):  ref base A (63/62);  new base G (4/5);  total (67/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

Only 100 of 134 total aligned reads displayed.

GATTTGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTT  >  Exported/1968559‑1968689
                                                                   |                                                               
gATTTGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTAC‑‑GAGGACAttt                                                            >  1:2967230/1‑71 (MQ=37)
   ttGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACAtttt                                                           >  1:3531816/1‑71 (MQ=255)
   ttGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACAtttt                                                           >  1:2489706/1‑71 (MQ=255)
   ttGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACAtttt                                                           >  1:154496/1‑71 (MQ=255)
   ttGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACAtttt                                                           >  1:572487/1‑71 (MQ=255)
   ttGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACAtttt                                                           >  1:4262856/1‑71 (MQ=255)
   ttGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACAtttt                                                           >  1:711201/1‑71 (MQ=255)
   ttGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACAtttt                                                           >  1:4248372/1‑71 (MQ=255)
   ttGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACAtttt                                                           >  1:811146/1‑71 (MQ=255)
   ttGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACAtt                                                             >  1:3684458/1‑69 (MQ=255)
   ttGGCTTTCTGTACTGCGTACATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACAttt                                                            >  1:3451513/1‑70 (MQ=255)
   ttGGCTTTCTGCACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACAtttt                                                           >  1:4691776/1‑71 (MQ=255)
    tGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACAtttt                                                           <  1:2733097/70‑1 (MQ=255)
     ggCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTg                                                          <  1:4385204/70‑1 (MQ=255)
     ggCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTg                                                          <  1:3751829/70‑1 (MQ=255)
     ggCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGa                                                         >  1:3646760/1‑71 (MQ=255)
     ggCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGa                                                         >  1:4287994/1‑71 (MQ=255)
     ggCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGa                                                         >  1:1585285/1‑71 (MQ=255)
      gCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGAc                                                        <  1:2755512/71‑1 (MQ=255)
      gCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGAc                                                        <  1:3080984/71‑1 (MQ=255)
      gCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGAc                                                        <  1:799837/71‑1 (MQ=255)
         ttCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTc                                                      >  1:1614553/1‑70 (MQ=255)
            tGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACt                                                       >  1:1941592/1‑66 (MQ=255)
                 tGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACAtttt                                                           <  1:455628/57‑1 (MQ=255)
                 tGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACAtttt                                                           <  1:4550363/57‑1 (MQ=255)
                  gCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTc                                                      >  1:3208152/1‑61 (MQ=255)
                   cGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTc                                                      >  1:4474991/1‑60 (MQ=255)
                   cGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTc                                                      >  1:3726828/1‑60 (MQ=255)
                   cGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTc                                                      >  1:4479969/1‑60 (MQ=255)
                   cGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTc                                                      >  1:1995757/1‑60 (MQ=255)
                   cGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGgat                                            <  1:3973860/70‑1 (MQ=255)
                   cGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGgat                                            <  1:4128141/70‑1 (MQ=255)
                   cGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGgat                                            <  1:4294856/70‑1 (MQ=255)
                    gTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGgatga                                          >  1:1923928/1‑71 (MQ=255)
                     tCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGgatgat                                         >  1:1081733/1‑71 (MQ=255)
                     tCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGgat                                            >  1:4028300/1‑68 (MQ=255)
                     tCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGgat                                            >  1:4059086/1‑68 (MQ=255)
                        aTGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGGCAtttt                                                           <  1:1706445/50‑1 (MQ=255)
                        aTGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACAtttt                                                           <  1:1278222/50‑1 (MQ=255)
                        aTGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACAtttt                                                           <  1:3854970/50‑1 (MQ=255)
                        aTGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACAtttt                                                           <  1:674255/50‑1 (MQ=255)
                        aTGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACAtttt                                                           <  1:4325689/50‑1 (MQ=255)
                        aTGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACAtttt                                                           <  1:2174618/50‑1 (MQ=255)
                         tGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAg                                                   <  1:2813735/57‑1 (MQ=255)
                           cTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAggg                                   >  1:2865227/1‑71 (MQ=255)
                           cTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAggg                                   >  1:281613/1‑71 (MQ=255)
                                cgcACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGc                              <  1:265137/71‑1 (MQ=255)
                                  cACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGgatgat                                         <  1:149173/58‑1 (MQ=255)
                                  cACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGgatgat                                         <  1:2312085/58‑1 (MQ=255)
                                  cACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGgatgat                                         <  1:4626978/58‑1 (MQ=255)
                                   aCGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCt                           >  1:250295/1‑71 (MQ=255)
                                    cGAATAGCATTGGCAAGCTCTTTACGTGAGGACAtttt                                                           <  1:485042/38‑1 (MQ=255)
                                    cGAATAGCATTGGCAAGCTCTTTACGTGAGGACAtttt                                                           <  1:3368134/38‑1 (MQ=255)
                                    cGAATAGCATTGGCAAGCTCTTTACGTGAGGACAtttt                                                           <  1:3362612/38‑1 (MQ=255)
                                    cGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGAtgatga                                       >  1:1182319/1‑58 (MQ=255)
                                    cGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCtt                          <  1:3585619/71‑1 (MQ=255)
                                    cGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCtt                          <  1:3778256/71‑1 (MQ=255)
                                    cGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCtt                          <  1:3903582/71‑1 (MQ=255)
                                      aaTAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTaa                        >  1:474951/1‑71 (MQ=255)
                                      aaTAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTaa                        >  1:389307/1‑71 (MQ=255)
                                       aTAGCATTGGCAAGCTCTTTACGTGAGGACAtttt                                                           >  1:3835747/1‑35 (MQ=255)
                                       aTAGCATTGGCAAGCTCTTTACGTGAGGACAtttt                                                           >  1:1248371/1‑35 (MQ=255)
                                       aTAGCATTGGCAAGCTCTTTACGTGAGGACAtttt                                                           >  1:1530152/1‑35 (MQ=255)
                                       aTAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAAc                       <  1:4246106/71‑1 (MQ=255)
                                       aTAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAAc                       <  1:96792/71‑1 (MQ=255)
                                       aTAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAAc                       <  1:4756319/71‑1 (MQ=255)
                                       aTAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAAc                       >  1:169834/1‑71 (MQ=255)
                                        tAGCATTGGCAAGCTCTTTACGTGAGGGCATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAAc                       >  1:952934/1‑70 (MQ=255)
                                        tAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAg                                                   <  1:568708/42‑1 (MQ=255)
                                        tAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAg                                                   <  1:2949929/42‑1 (MQ=255)
                                        tAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAggg                                   >  1:545119/1‑58 (MQ=255)
                                            aTTGGCAAGCTCTTTACGTGAGGACGTTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACtt                  >  1:1423155/1‑71 (MQ=255)
                                            aTTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACtt                  <  1:4035588/71‑1 (MQ=255)
                                            aTTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACtt                  <  1:1506708/71‑1 (MQ=255)
                                            aTTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACtt                  <  1:258341/71‑1 (MQ=255)
                                             ttGGCAAGCTCTTTACGTGAGGGCATTTTGACTCCAGATCGgat                                            <  1:4570584/44‑1 (MQ=255)
                                             ttGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGatgatg                                        <  1:3903238/48‑1 (MQ=255)
                                             ttGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGatgatg                                        <  1:3839208/48‑1 (MQ=255)
                                             ttGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGatgatg                                        <  1:4598768/48‑1 (MQ=255)
                                               ggCAAGCTCTTTACGTGAGGGCATTTTGACTCCAGATCGgatga                                          >  1:1006063/1‑44 (MQ=255)
                                               ggCAAGCTCTTTACGTGAGGGCATTTTGACTCCAGATCGgatga                                          >  1:3438764/1‑44 (MQ=255)
                                               ggCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTaa                        >  1:3538208/1‑62 (MQ=255)
                                               ggCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTgac               >  1:1257138/1‑71 (MQ=255)
                                                gCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTgac               <  1:1188549/70‑1 (MQ=255)
                                                 cAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTgacga             <  1:2686952/71‑1 (MQ=255)
                                                 cAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTgacga             <  1:4164535/71‑1 (MQ=255)
                                                 cAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTgacga             <  1:3585436/71‑1 (MQ=255)
                                                    gCTCTTTATGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAg          >  1:2186796/1‑71 (MQ=255)
                                                    gCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAg          <  1:3977476/71‑1 (MQ=255)
                                                      tctTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAgcg        >  1:1010180/1‑71 (MQ=255)
                                                       ctTTACGTGAGGACATTTTGACTCCAGATCGGATGATGaa                                      <  1:798458/40‑1 (MQ=255)
                                                       ctTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACa           <  1:4388587/67‑1 (MQ=255)
                                                        tttACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGc                              >  1:1010065/1‑47 (MQ=255)
                                                        tttACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGc                              >  1:4449358/1‑47 (MQ=255)
                                                          tACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGtt    >  1:1590745/1‑71 (MQ=255)
                                                          tACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGtt    >  1:1372268/1‑71 (MQ=255)
                                                          tACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGtt    >  1:3141758/1‑71 (MQ=255)
                                                          tACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGtt    >  1:1834759/1‑71 (MQ=255)
                                                              tGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTgac               <  1:3677743/56‑1 (MQ=255)
                                                              tGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTgac               <  1:664241/56‑1 (MQ=255)
                                                                   |                                                               
GATTTGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTT  >  Exported/1968559‑1968689

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: