Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	4,528,968	A→G	74.5%	L242P (CTA→CCA)	yjhH ←	predicted lyase/synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,528,968	0	A	G	74.5%	83.4 / 27.3	47	L242P (CTA→CCA)	yjhH	predicted lyase/synthase
Reads supporting (aligned to +/- strand):  ref base A (10/2);  new base G (24/11);  total (34/13)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACG  >  W3110S.gb/4528900‑4529033                                                                     |                                                                  acTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTagtggt                                                                 >  1:2083930/1‑71 (MQ=255)  cTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGtt                                                                >  1:2812173/1‑71 (MQ=255)    gCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGtttt                                                              >  1:1760252/1‑71 (MQ=255)    gCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGtttt                                                              >  1:2376424/1‑71 (MQ=255)    gCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGtttt                                                              >  1:2692160/1‑71 (MQ=255)    gCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGtttt                                                              >  1:2169991/1‑71 (MQ=255)    gCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGtttt                                                              >  1:2885474/1‑71 (MQ=255)    gCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGtttt                                                              >  1:2900586/1‑71 (MQ=255)    gCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGtttt                                                              >  1:1453985/1‑71 (MQ=255)    gCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGtttt                                                              >  1:1217240/1‑71 (MQ=255)      aTGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGtttttt                                                            <  1:861523/71‑1 (MQ=255)         cTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTAtt                                                         >  1:1527782/1‑71 (MQ=255)         cTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTAtt                                                         >  1:2702742/1‑71 (MQ=255)            tATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTAt                                                          <  1:742802/67‑1 (MQ=255)            tATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTAt                                                          <  1:321330/67‑1 (MQ=255)                tGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAg                                                      >  1:3108565/1‑67 (MQ=255)                tGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAg                                                      >  1:128253/1‑67 (MQ=255)                  aTCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGc                                                <  1:1290269/71‑1 (MQ=255)                        gagaCAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCTCAGCGGt                                          >  1:879198/1‑71 (MQ=255)                        gagaCAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGt                                          >  1:1818681/1‑71 (MQ=255)                        gagaCAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGt                                          >  1:932418/1‑71 (MQ=255)                         agaCAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTc                                         >  1:1110297/1‑71 (MQ=255)                         agaCAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTc                                         >  1:2114729/1‑71 (MQ=255)                          gaCAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCg                                        >  1:3566971/1‑71 (MQ=255)                          gaCAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCg                                        >  1:3564108/1‑71 (MQ=255)                          gaCAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCg                                        >  1:2815641/1‑71 (MQ=255)                          gaCAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCg                                        >  1:3083454/1‑71 (MQ=255)                          gaCAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCg                                        >  1:202998/1‑71 (MQ=255)                           aCAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGc                                       <  1:3655424/71‑1 (MQ=255)                           aCAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGc                                       <  1:269861/71‑1 (MQ=255)                             aaaCGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGc                                       >  1:905900/1‑69 (MQ=255)                             aaaCGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGc                                       >  1:3234020/1‑69 (MQ=255)                              aaCGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAg                                    <  1:2097486/71‑1 (MQ=255)                                    gtTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTc                                                  >  1:2356500/1‑51 (MQ=255)                                    gtTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGcc                              >  1:3564922/1‑71 (MQ=255)                                     tttCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGttttt                                                             <  1:2587070/39‑1 (MQ=255)                                     tttCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGttttt                                                             <  1:2362188/39‑1 (MQ=255)                                       tCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGAt                                  >  1:672595/1‑64 (MQ=255)                                       tCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGAt                                  >  1:523802/1‑64 (MQ=255)                                       tCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGAt                                  >  1:534760/1‑64 (MQ=255)                                         gAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCAc                         <  1:2386307/71‑1 (MQ=255)                                            ggCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGcc                      <  1:2599287/71‑1 (MQ=255)                                                  aaTAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGgcgc                >  1:104026/1‑71 (MQ=255)                                                  aaTAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTAGCAGCGGTCGCCAGATCGCCTTCACGCCAGgcgc                >  1:1612882/1‑71 (MQ=255)                                                      gCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCgg              >  1:2894361/1‑69 (MQ=255)                                                             gTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCa                                     <  1:1695650/39‑1 (MQ=255)                                                              ttGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCa                                     <  1:2203696/38‑1 (MQ=255)                                                               tGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGAc   >  1:158231/1‑71 (MQ=255)                                                                 tagtggTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACg  <  1:1482302/70‑1 (MQ=255)                                                                     |                                                                  ACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACG  >  W3110S.gb/4528900‑4529033

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