Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,096,620 A→T 7.7% S107T (TCC→ACC)  garK ← glycerate kinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,096,6200AT 7.7% 181.2 / 6.2 65S107T (TCC→ACC) garKglycerate kinase I
Reads supporting (aligned to +/- strand):  ref base A (35/25);  new base T (2/3);  total (37/28)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

TTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTG  >  Exported/2096551‑2096689
                                                                     |                                                                     
tttCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGga                                                                       <  1:5280595/70‑1 (MQ=255)
  tCCGCAGGTACCTGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGgaga                                                                     <  1:1582598/70‑1 (MQ=255)
   ccGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATcc                                                                  >  1:787884/1‑71 (MQ=255)
    cGCATGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATc                                                                   <  1:3388872/70‑1 (MQ=255)
     gCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATccc                                                                 >  1:4005276/1‑71 (MQ=255)
          tACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAAc                                                            <  1:602647/71‑1 (MQ=255)
             cAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTgg                                                         <  1:3979976/71‑1 (MQ=255)
                cTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATcc                                                                  >  1:1091799/1‑59 (MQ=255)
                cTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATc                                                                   >  1:1253572/1‑58 (MQ=255)
                 tCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCaa                                                             >  1:2957591/1‑63 (MQ=255)
                   cAGCCCACTGGCGGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCa                                                   <  1:2768444/71‑1 (MQ=255)
                       ccACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACttt                                               <  1:506739/71‑1 (MQ=255)
                         aCTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGgaga                                                                     <  1:755891/47‑1 (MQ=255)
                         aCTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCa                                                              <  1:2198443/54‑1 (MQ=255)
                          cTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATccccc                                                               <  1:1745855/52‑1 (MQ=255)
                          cTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATccccc                                                               <  1:4667717/52‑1 (MQ=255)
                          cTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTtctc                                            >  1:3926480/1‑71 (MQ=255)
                          cTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCCGAGATccccc                                                               <  1:179563/52‑1 (MQ=255)
                           tGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCg                                           >  1:4799295/1‑71 (MQ=255)
                           tGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCg                                           >  1:3593843/1‑71 (MQ=255)
                            ggCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCaa                                                             >  1:4515923/1‑52 (MQ=255)
                            ggCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCaa                                                             >  1:3957323/1‑52 (MQ=255)
                            ggCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAAc                                                            >  1:2219952/1‑53 (MQ=255)
                            ggCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGc                                          >  1:2411976/1‑71 (MQ=255)
                             gctgccgccATCTCAATCACCGCCGTTTTCCCGTCGCCGGTAAGGCCATAAAAAGCATTCACTTTTTCGcc                                         >  1:697326/4‑71 (MQ=255)
                             gctgccgccATCTCAATCACCGCCGTTTTCCCGTCGCCGGTAAGGCCATAAAAAGCATTCACTTTTTCGcc                                         >  1:786951/4‑71 (MQ=255)
                              ctgccgccATCTCAATCACCGCCGTTTTCCCGTCGCCGGTAAGGCCATAAAAAGCATTCACTTTTTCGcc                                         <  1:5535532/68‑1 (MQ=255)
                              ccgccgccATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTgg                                                         <  1:395542/54‑1 (MQ=255)
                              ccgccgccATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGccc                                        <  1:3151299/71‑1 (MQ=255)
                                 ccgccATCTCAATCACCGCCGTTTTCCCGTCGCCGGTAAGGCCATAAAAAGCATTCACTTTTTCGCCCa                                       <  1:526871/69‑1 (MQ=255)
                                  cgccATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATccccc                                                               <  1:3167975/44‑1 (MQ=255)
                                  cgccATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATccccc                                                               <  1:5388770/44‑1 (MQ=255)
                                  cgccATCTCAATCACCGCCGTTTTCCCGTCGCCGGTAAGGCCATAAAAAGCATTCACTTTTTCGCCCATCg                                    <  1:5037709/71‑1 (MQ=255)
                                   gccATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATc                                                                   >  1:605773/1‑39 (MQ=255)
                                   gccATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAAc                                                            >  1:5606833/1‑46 (MQ=255)
                                    ccATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGc                                  >  1:433689/1‑71 (MQ=255)
                                       tttCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCt                               >  1:5273455/1‑71 (MQ=255)
                                        ttCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTg                              >  1:5572796/1‑71 (MQ=255)
                                         tCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCAtt                                                     <  1:5260805/47‑1 (MQ=255)
                                         tCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTtct                                             >  1:5192297/1‑55 (MQ=255)
                                         tCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTtc                                              >  1:3790962/1‑54 (MQ=255)
                                            aTAAACGCGGTTTTGACATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAAc                          >  1:3663530/1‑71 (MQ=255)
                                            aTAAACGCGGTTTTGACATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAAc                          >  1:3670230/1‑71 (MQ=255)
                                             tAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAg                                      >  1:3800223/1‑58 (MQ=255)
                                               aaCGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCa                       <  1:4018471/71‑1 (MQ=255)
                                                 cgcgGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACt                                                 >  1:4415560/1‑43 (MQ=255)
                                                 cgcgGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACt                                                 >  1:4705379/1‑43 (MQ=255)
                                                  gcgGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTtct                                             >  1:922753/1‑46 (MQ=255)
                                                  gcgGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTtct                                             >  1:1380603/1‑46 (MQ=255)
                                                   cgGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCg                   <  1:2474521/71‑1 (MQ=255)
                                                       tttGCCATCGCCGGAGATCCCCCAACTGGCATTCACttt                                               <  1:2655758/39‑1 (MQ=255)
                                                        ttGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAg                                      >  1:976787/1‑47 (MQ=255)
                                                        ttGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGcc                                 >  1:4816747/1‑52 (MQ=255)
                                                        ttGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGa                <  1:4347880/69‑1 (MQ=255)
                                                        ttGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGa                <  1:660612/69‑1 (MQ=255)
                                                        ttGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGAc               <  1:3019747/70‑1 (MQ=255)
                                                        ttGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGAc               <  1:1676490/70‑1 (MQ=255)
                                                        ttGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACg              >  1:3011316/1‑71 (MQ=255)
                                                        ttGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTATCGCCCAGCGGCCCTGTAACCCAGGCGTGa                <  1:5476708/69‑1 (MQ=255)
                                                        ttGCCAGCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACg              >  1:737594/1‑71 (MQ=255)
                                                         tGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGt             >  1:1211212/1‑71 (MQ=255)
                                                         tGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGt             >  1:5274956/1‑71 (MQ=255)
                                                          gCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACg              >  1:592681/1‑69 (MQ=255)
                                                          gCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGtt            >  1:4407448/1‑71 (MQ=255)
                                                              tCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCg                   <  1:102826/60‑1 (MQ=255)
                                                                gCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTaa                           >  1:2216563/1‑50 (MQ=255)
                                                                   ggAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCt   <  1:1456238/71‑1 (MQ=255)
                                                                   ggAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCt   <  1:4734987/71‑1 (MQ=255)
                                                                    gagaTCCCCCAACTGGCATTCACTTTCTCGCCCATCGGCCCTGt                             <  1:5020226/44‑1 (MQ=255)
                                                                    gagaTCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTg  >  1:3797718/1‑71 (MQ=255)
                                                                     |                                                                     
TTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTG  >  Exported/2096551‑2096689

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: