Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,096,612	A→G	7.2%	D109D (GAT→GAC)	garK ←	glycerate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,096,612	0	A	G	7.2%	192.1 / 4.8	69	D109D (GAT→GAC)	garK	glycerate kinase I
Reads supporting (aligned to +/- strand):  ref base A (38/26);  new base G (2/3);  total (40/29)
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.45e-01

TCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTT  >  Exported/2096544‑2096679                                                                     |                                                                    tCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATc                                                                   >  1:146437/1‑71 (MQ=255)       ttttCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGcc                                                                >  1:2997059/1‑68 (MQ=255)        tttCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCg                                                               >  1:3010229/1‑68 (MQ=255)        tttCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGga                                                             <  1:5280595/70‑1 (MQ=255)         ttCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCg                                                               >  1:2439995/1‑67 (MQ=255)          tCCGCAGGTACCTGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGgaga                                                           <  1:1582598/70‑1 (MQ=255)           ccGCAGGTACCAGCTCCAGCCCACTGGCCGCCG‑CATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATcc                                                        >  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:3593843/1‑71 (MQ=255)                                   tGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCg                                 >  1:4799295/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)                                         ccgccATTTCAATAAACGCGGTTTTGCCATCGCCg                                                               <  1:4194555/35‑1 (MQ=255)                                         ccgccATCTCAATCACCGCCGTTTTCCCGTCGCCGGTAAGGCCATAAAAAGCATTCACTTTTTCGCCCa                             <  1:526871/69‑1 (MQ=255)                                          cgccATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATccccc                                                     <  1:5388770/44‑1 (MQ=255)                                          cgccATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATccccc                                                     <  1:3167975/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:3670230/1‑71 (MQ=255)                                                    aTAAACGCGGTTTTGACATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAAc                >  1:3663530/1‑71 (MQ=255)                                                     tAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAg                            >  1:3800223/1‑58 (MQ=255)                                                       aaCGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCa             <  1:4018471/71‑1 (MQ=255)                                                         cgcgGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACt                                       >  1:4705379/1‑43 (MQ=255)                                                         cgcgGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACt                                       >  1:4415560/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:660612/69‑1 (MQ=255)                                                                ttGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGa      <  1:4347880/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:5274956/1‑71 (MQ=255)                                                                 tGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGt   >  1:1211212/1‑71 (MQ=255)                                                                  gCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACg    >  1:592681/1‑69 (MQ=255)                                                                  gCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGtt  >  1:4407448/1‑71 (MQ=255)                                                                     |                                                                    TCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTT  >  Exported/2096544‑2096679

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: ‑