Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,629,792	T→C	41.5%	D467G (GAC→GGC)	guaA ←	GMP synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,629,792	0	T	C	41.5%	18.3 / 47.5	41	D467G (GAC→GGC)	guaA	GMP synthetase
Reads supporting (aligned to +/- strand):  ref base T (13/11);  new base C (3/14);  total (16/25)
Fisher's exact test for biased strand distribution p-value = 2.51e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

CGGCAGATGCGCCCAGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAACACAGTG  >  W3110S.gb/2629725‑2629848                                                                    |                                                         cGGCAGATGCGCCCAGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCAt                                                       <  1:3184255/71‑1 (MQ=255)    cAGATGCGCCGAGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACt                                                    >  1:2348278/1‑71 (MQ=255)         gcgcCCAGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGCCATACTTACGa                                               >  1:4050847/1‑71 (MQ=255)               aGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCg                                         >  1:189758/1‑71 (MQ=255)                 tgtgCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGcc                                       <  1:87456/71‑1 (MQ=255)                 tgtgCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGCCATACTTACGACCATCGcc                                       <  1:354136/71‑1 (MQ=255)                  gtgCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGCCATACTTACGACCAt                                           >  1:136119/1‑66 (MQ=255)                     cGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCa                                     >  1:2822912/1‑69 (MQ=255)                     cGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCa                                     >  1:723322/1‑69 (MQ=255)                     cGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCa                                     >  1:167759/1‑69 (MQ=255)                       gTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATAc                                                     >  1:512188/1‑51 (MQ=255)                       gTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATAc                                                     >  1:3414845/1‑51 (MQ=255)                       gTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAAACCAGTCATAc                                                     >  1:1987824/1‑51 (MQ=255)                           tAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGCCa                                                        <  1:1694232/44‑1 (MQ=255)                           tAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGCCa                                                        <  1:1468616/44‑1 (MQ=255)                            aaaGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCaa                            >  1:742157/1‑71 (MQ=255)                             aaGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTa                                                  >  1:1174362/1‑48 (MQ=255)                             aaGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTa                                                  >  1:2158580/1‑48 (MQ=255)                               gTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACgg                         <  1:3189689/71‑1 (MQ=255)                               gTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACgg                         <  1:742073/71‑1 (MQ=255)                               gTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACgg                         <  1:3561490/71‑1 (MQ=255)                               gTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACgg                         <  1:2852371/71‑1 (MQ=255)                               gTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACg                          >  1:1582216/1‑70 (MQ=255)                               gTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACCCCAACg                          >  1:1098071/1‑70 (MQ=255)                                   aTGGTTTCGACAGCACGCAGAGAGACAACCCAGCCATACTTACGACCATCGCCCATTACGCCAACGGAACg                     <  1:806414/71‑1 (MQ=255)                                   aTGGTTTCGACAGCACGCAGAGAGACAACCCAGCCATACTTACGACCATCGCCCATTACGCCAACGGAACg                     <  1:1178162/71‑1 (MQ=255)                                    tGGTTTCGACAGCACGCAGAGAGACAACCCAGCCATACTTACGACCATCGCCCATTACGCCAACGGAACGt                    <  1:3363032/71‑1 (MQ=255)                                    tGGTTTCGACAGCACGCAGAGAGACAACCCAGCCATACTTACGACCATCGCCCATTACGCCAACGGAACGt                    <  1:2119209/71‑1 (MQ=255)                                       tttCGACAGCACGCAGAGAGACAACCCAGCCATACTTACGACCATCGcc                                       <  1:3191867/49‑1 (MQ=255)                                       tttCGACAGCACGCAGAGAGACAACCCAGCCATACTTACGACCATCGcc                                       <  1:1694409/49‑1 (MQ=255)                                       tttCGACAGCACGCAGAGAGACAACCCAGCCATACTTACGACCATCGcc                                       <  1:1716953/49‑1 (MQ=255)                                       tttCGACAGCACGCAGAGAGACAACCCAGCCATACTTACGACCATCGcc                                       <  1:1886922/49‑1 (MQ=255)                                       tttCGACAGCACGCAGAGAGACAACCCAGCCATACTTACGACCATCGcc                                       <  1:1618252/49‑1 (MQ=255)                                         tCGACAGCACGCAGAGAGACAACCCAGCCATACTTACGACCATCGCCCATTa                                  >  1:3077827/1‑52 (MQ=255)                                                cACGCAGAGAGACAACCCAGCCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAac        <  1:369101/71‑1 (MQ=255)                                                  cGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAacac      <  1:1448420/71‑1 (MQ=255)                                                        gagtCAACCCAGCCATACTTACGACCATCGCCCATTACGCCAACGGAACGt                    <  1:629523/51‑1 (MQ=39)                                                                ccAGTCATACTTACGACCATCGCCCATTACGCCAACgg                         <  1:475052/38‑1 (MQ=255)                                                                   gTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAACACAGTg  <  1:2293045/58‑1 (MQ=255)                                                                   gTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAACACAGTg  <  1:1052956/58‑1 (MQ=255)                                                                   gTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAACACAGTg  <  1:293597/58‑1 (MQ=255)                                                                    |                                                         CGGCAGATGCGCCCAGTGTGCGGTCATAAAGTCGATGGTTTCGACAGCACGCAGAGAGACAACCCAGTCATACTTACGACCATCGCCCATTACGCCAACGGAACGTACCGGCAGGAACACAGTG  >  W3110S.gb/2629725‑2629848

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