Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,096,599	A→C	8.2%	F114V (TTT→GTT) ‡	garK ←	glycerate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,096,599	0	A	C	8.2%	138.8 / 7.4	61	F114V (TTT→GTT) ‡	garK	glycerate kinase I
Reads supporting (aligned to +/- strand):  ref base A (29/27);  new base C (2/3);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.18e-01

CACGAGTGGATCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCA  >  Exported/2096534‑2096668                                                                  |                                                                      cACGAGTGGATCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATaaa                                                                       >  1:4661834/1‑66 (MQ=255)  aCGAGTGGATCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGt                                                                 <  1:1169047/71‑1 (MQ=255)     aGTGGATCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGttt                                                               >  1:3821196/1‑70 (MQ=255)       tGGATCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGc                                                            <  1:522963/71‑1 (MQ=255)       tGGATCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGc                                                            <  1:4657247/71‑1 (MQ=255)        ggATCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGc                                                            <  1:979523/70‑1 (MQ=255)        ggATCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGc                                                            >  1:2143747/1‑70 (MQ=255)         gATCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGcc                                                           <  1:5144728/70‑1 (MQ=255)           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)                   ttCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTg                                                             >  1:1095900/1‑58 (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)                       gCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGt                                                                 <  1:675312/50‑1 (MQ=255)                       gCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGt                                                                 <  1:1759723/50‑1 (MQ=255)                       gCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATccc                                            >  1:4005276/1‑71 (MQ=255)                            tACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAAc                                       <  1:602647/71‑1 (MQ=255)                              ccAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGcc                                                           <  1:814367/49‑1 (MQ=255)                               cAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTgg                                    <  1:3979976/71‑1 (MQ=255)                                aGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGcc                                                           <  1:4865336/47‑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:3957323/1‑52 (MQ=255)                                              ggCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCaa                                        >  1:4515923/1‑52 (MQ=255)                                              ggCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAAc                                       >  1:2219952/1‑53 (MQ=255)                                              ggCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGc                     >  1:2411976/1‑71 (MQ=255)                                               gctgccgccATCTCAATCACCGCCGTTTTCCCGTCGCCGGTAAGGCCATAAAAAGCATTCACTTTTTCGcc                    >  1:786951/4‑71 (MQ=255)                                               gctgccgccATCTCAATCACCGCCGTTTTCCCGTCGCCGGTAAGGCCATAAAAAGCATTCACTTTTTCGcc                    >  1:697326/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: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)                                                                  |                                                                      CACGAGTGGATCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCA  >  Exported/2096534‑2096668

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