Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,123,207	A→G	6.5%	S7P (TCT→CCT)	yraQ ←	predicted permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,123,207	0	A	G	6.5%	175.3 / 2.4	61	S7P (TCT→CCT)	yraQ	predicted permease
Reads supporting (aligned to +/- strand):  ref base A (21/36);  new base G (2/2);  total (23/38)
Fisher's exact test for biased strand distribution p-value = 6.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATAATTGAGCGCGTGAATATTA  >  Exported/2123139‑2123271                                                                     |                                                                 gagaCCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATg                                                                <  1:3053060/71‑1 (MQ=255)      cGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTgaga                                                                  <  1:5503037/64‑1 (MQ=255)      cGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGa                                                            <  1:4480375/70‑1 (MQ=255)         cAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGa                                                        >  1:5151487/1‑71 (MQ=255)          aaTGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGaa                                                              <  1:217848/64‑1 (MQ=255)           aTGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGAc                                                       >  1:2583069/1‑70 (MQ=255)           aTGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGAc                                                       >  1:4178173/1‑70 (MQ=255)              aCAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTgaga                                                                  >  1:5249533/1‑56 (MQ=255)              aCAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGt                                                   <  1:2770750/71‑1 (MQ=255)               cAACGAGAAAGAATAGAACGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGt                                                   <  1:3250064/70‑1 (MQ=255)               cAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGt                                                   >  1:1008084/1‑70 (MQ=255)               cAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGt                                                   >  1:2426632/1‑70 (MQ=255)                aaCGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGaa                                                              <  1:4616937/58‑1 (MQ=255)                 aCGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCa                                                 <  1:767374/70‑1 (MQ=255)                   gagaAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAg                                              <  1:3563324/71‑1 (MQ=255)                      aaagaaaAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCa                                            >  1:1401384/1‑70 (MQ=255)                       aagaaaAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCAtt                                          >  1:1518720/1‑71 (MQ=255)                       aagaaaAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCAtt                                          >  1:3192733/1‑71 (MQ=255)                        agaaaAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTc                                         <  1:3606863/71‑1 (MQ=255)                              gcgcGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGaa                                                              >  1:4076110/1‑44 (MQ=255)                              gcgcGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGaa                                                              >  1:5139188/1‑44 (MQ=255)                              gcgaGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGaa                                                              >  1:3973085/1‑44 (MQ=255)                               cgcgGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGa                                  <  1:921719/71‑1 (MQ=255)                                 cgGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGaa                                 >  1:3954622/1‑70 (MQ=255)                                    gTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATTAAGACTGACCAGTCATAGCATTCCCCAGGAATa                               <  1:5038811/69‑1 (MQ=255)                                    gTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAg                             <  1:808264/71‑1 (MQ=255)                                     tttCCACAACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCAAGGAATAAGt                            <  1:2920774/71‑1 (MQ=255)                                         caccacTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTAt                         <  1:1006002/70‑1 (MQ=255)                                           ccacTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGAt                      <  1:4883417/71‑1 (MQ=255)                                           ccacTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGATATGa                       >  1:3058420/1‑70 (MQ=255)                                            cacTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATa                     <  1:4413309/71‑1 (MQ=255)                                              cTGAATGGGCGTTGCCGCCTGAGATGAAGACTGCCCAGTCATAGCATTCCCCAGGAATAAGTTATGATAAt                   >  1:1255115/1‑71 (MQ=255)                                              cTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATAAt                   <  1:1866234/71‑1 (MQ=255)                                               tGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATa                               <  1:3153554/58‑1 (MQ=255)                                               tGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATa                               <  1:3338260/58‑1 (MQ=255)                                               tGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATa                               <  1:5583223/58‑1 (MQ=255)                                                 aaTGGGCGTTGCCGCCTGAGGTGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATAAtt                  >  1:2285824/1‑69 (MQ=255)                                                 aaTGGGCGTTGCCGCCTGAGGTGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATAAtt                  >  1:1660143/1‑69 (MQ=255)                                                 aaTGGGCGTTGCCGCCTGAGGTGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATAATTg                 <  1:5157565/70‑1 (MQ=255)                                                 aaTGGGCGTTGCCGCCTGAGGTGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATAATTg                 <  1:4809415/70‑1 (MQ=255)                                                 aaTGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCGTAGCATTCCCCAGGAATAAGTTATGATAAtt                  <  1:4439756/69‑1 (MQ=255)                                                 aaTGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATAAtt                  <  1:2296092/69‑1 (MQ=255)                                                 aaTGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATAAtt                  <  1:2741234/69‑1 (MQ=255)                                                 aaTGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATAAtt                  <  1:451593/69‑1 (MQ=255)                                                 aaTGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATAAtt                  >  1:1867392/1‑69 (MQ=255)                                                 aaTGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATAATTg                 <  1:5535808/70‑1 (MQ=255)                                                 aaTGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATAATTGa                <  1:5293539/71‑1 (MQ=255)                                                  aTGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATAATTGAg               >  1:4015141/1‑71 (MQ=255)                                                   tGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATAATTGAgc              >  1:1262209/1‑71 (MQ=255)                                                   tGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATAATTGAg               >  1:1511947/1‑70 (MQ=255)                                                         ttGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATaa                              <  1:36104/49‑1 (MQ=255)                                                         ttGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATa                               <  1:441340/48‑1 (MQ=255)                                                         ttGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATAATTGa                <  1:3075335/63‑1 (MQ=255)                                                         ttGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATAATTGAGCGCGTg         <  1:1422220/70‑1 (MQ=255)                                                         ttGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATAATTGAGCGCGTg         <  1:3207527/70‑1 (MQ=255)                                                           gccgccTGAGATGAAGACTGACGAGTCATAGCATTCCCCAGGAATAAGTTATGATAATTGAGCGCGTGaa       >  1:5585162/1‑70 (MQ=255)                                                           gccgccTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATAATTGAGCGCGTGaa       <  1:5229869/70‑1 (MQ=255)                                                             cgccTGAGATGAAGACTGACCAGTCATAGCATTCCCCGGGAATAAGTTATGATAATTGAGCGCGTGAatat    <  1:583398/71‑1 (MQ=255)                                                             cgccTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATAATTGAGCGCGTGAatat    <  1:4300765/71‑1 (MQ=255)                                                             cgccTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATAATTGAGCGCGTGAatat    <  1:1405141/71‑1 (MQ=255)                                                             cgccTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATAATTGAGCGCGTGAatat    <  1:4435486/71‑1 (MQ=255)                                                               ccTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATAATTGAGCGCGTGAATATTa  >  1:3523974/1‑71 (MQ=255)                                                               ccTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATAATTGAGCGCGTGAATATTa  >  1:2515723/1‑71 (MQ=255)                                                                  gagaTGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATAATTGAGCGCGTGAata     >  1:695028/1‑65 (MQ=255)                                                                     |                                                                 GAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATAATTGAGCGCGTGAATATTA  >  Exported/2123139‑2123271

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