Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,096,627	C→A	8.6%	W104C (TGG→TGT) ‡	garK ←	glycerate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,096,627	0	C	A	8.6%	183.1 / 3.5	58	W104C (TGG→TGT) ‡	garK	glycerate kinase I
Reads supporting (aligned to +/- strand):  ref base C (32/21);  new base A (2/3);  total (34/24)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

TACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGT  >  Exported/2096561‑2096692                                                                   |                                                                  tACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAAc                                                               <  1:602647/71‑1 (MQ=255)    cAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTgg                                                            <  1:3979976/71‑1 (MQ=255)        tCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCaa                                                                >  1:2957591/1‑63 (MQ=255)          cAGCCCACTGGCGGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCa                                                      <  1:2768444/71‑1 (MQ=255)              ccACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACttt                                                  <  1:506739/71‑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)                          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:4415560/1‑43 (MQ=255)                                        cgcgGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACt                                                    >  1:4705379/1‑43 (MQ=255)                                         gcgGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTtct                                                >  1:1380603/1‑46 (MQ=255)                                         gcgGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTtct                                                >  1:922753/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:1676490/70‑1 (MQ=255)                                               ttGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGAc                  <  1:3019747/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)                                                             gaTCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGa                   >  1:2952483/1‑55 (MQ=255)                                                              aTCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGt  <  1:1443069/71‑1 (MQ=255)                                                                  cccAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGt                     <  1:851536/48‑1 (MQ=255)                                                                  cccAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGt                     <  1:4788619/48‑1 (MQ=255)                                                                  cccAACTGGCATTCACTTTCTCGCCCAGAGGCCCTGTAACCCAGGCGt                     <  1:3133575/48‑1 (MQ=255)                                                                   |                                                                  TACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGT  >  Exported/2096561‑2096692

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