Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,663,570	Δ1 bp	30.9%	coding (552/882 nt)	yfhR →	predicted peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,663,570	0	T	.	30.9%	104.2 / 61.1	55	coding (552/882 nt)	yfhR	predicted peptidase
Reads supporting (aligned to +/- strand):  ref base T (21/17);  new base . (9/8);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GTATTGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGG  >  W3110S.gb/2663512‑2663632                                                           |                                                               gTATTGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGTGATCGTGAAGGCATACG‑GCGGTGATCCTc                                                    >  1:890849/1‑70 (MQ=255) gTATTGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGTGATCGTGAAGGCATACG‑GCGGTGATCCTc                                                    >  1:3046484/1‑70 (MQ=255)    ttGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGAc                                                 >  1:1050983/1‑71 (MQ=255)    ttGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGAc                                                 >  1:1421807/1‑71 (MQ=255)    ttGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGAc                                                 >  1:3439855/1‑71 (MQ=255)      ggCGGGGCGAATATTCTGGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTc                                               >  1:430744/1‑71 (MQ=255)      ggCGGGGCGAATATTCTGGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTc                                               >  1:3702492/1‑71 (MQ=255)      ggCGGGGCGAATATTCTGGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTc                                               >  1:3511/1‑71 (MQ=255)      ggCGGGGCGAATATTCTGGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTc                                               >  1:2948959/1‑71 (MQ=255)      ggCGGGGCGAATATTCTGGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTc                                               >  1:2287706/1‑71 (MQ=255)      ggCGGGGCGAATATTCTGGATGTTATTGGTCGGGGTGATCGTGAAGGCATACG‑GCGGTGATCCTCGACTc                                               >  1:866491/1‑70 (MQ=255)       gCGGGGCGAATATTCTGGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTc                                               <  1:88203/70‑1 (MQ=255)         ggggCGAATATTCTGGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCca                                             >  1:2922041/1‑70 (MQ=255)            gCGAATATTCTGGATGTTATTGGTCGGGGTGATCGTGAAGGCATACG‑GCGGTGATCCTCGACTCCACAttt                                        <  1:2079231/71‑1 (MQ=255)            gCGAATATTCTGGATGTTATTGGTCGGGGTGATCGTGAAGGCATACG‑GCGGTGATCCTCGACTCCACAttt                                        <  1:2758315/71‑1 (MQ=255)               aaTATTCTGGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTg                                                               <  1:987581/46‑1 (MQ=255)               aaTATTCTGGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGc                                      <  1:1486534/71‑1 (MQ=255)               aaTATTCTGGATGTTATTGGTCGGGGTGATCGTGAAGGCATACG‑GCGGTGATCCTCGACTCCACATTTGcc                                     <  1:180606/71‑1 (MQ=255)                 tatTCTGGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCct                                    >  1:1887559/1‑71 (MQ=255)                 tatTCTGGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGT‑‑GGTGATCCTCGACTCCACATTTGCctct                                  <  1:1806842/71‑1 (MQ=255)                 tatTCTGGATGTTATTGGGCGGGGTGATCGTGAAGGCATACGTG‑GG‑GATCCTCGACTCCACATTTGCctct                                  <  1:1802124/71‑1 (MQ=255)                   ttCTGGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCctct                                  <  1:1577231/71‑1 (MQ=255)                   ttCTGGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCctct                                  <  1:3479575/71‑1 (MQ=255)                           gTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACAtt                                         >  1:2026275/1‑56 (MQ=255)                           gTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACAtt                                         >  1:397096/1‑56 (MQ=255)                           gTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACAtt                                         >  1:2135767/1‑56 (MQ=255)                            ttATTGGTCGGGGTGATCGTGAAGGCATACG‑GCGGTGATCCTCGACTCCACATTTGCctct                                  <  1:3448854/61‑1 (MQ=255)                            ttATTGGTCGGGGTGATCGTGAAGGCATACG‑GCGGTGATCCTCGACTCCACATTTGCctct                                  <  1:3871139/61‑1 (MQ=255)                            ttATTGGTCGGGGTGATCGTGAAGGCATACG‑GCGGTGATCCTCGACTCCACATTTGCctct                                  <  1:3536058/61‑1 (MQ=255)                              aTTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCAt                       >  1:2836561/1‑71 (MQ=255)                                tGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCg                     >  1:894156/1‑71 (MQ=255)                                           aTCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAAt           >  1:1480297/1‑70 (MQ=255)                                           aTCGTGAAGGCATACG‑GCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGa         >  1:966935/1‑71 (MQ=255)                                             cGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGAt        >  1:2921387/1‑71 (MQ=255)                                             cGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGAt        >  1:1012858/1‑71 (MQ=255)                                             cGTGAAGGCATACG‑GCGGTGATCCTCGACTCCACATTTGCctct                                  <  1:2676822/44‑1 (MQ=39)                                               tGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATcc      <  1:1215797/71‑1 (MQ=255)                                               tGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATcc      <  1:3267116/71‑1 (MQ=255)                                                 aaGGCATACG‑GCGGTGATCCTCGACTCCACATTTGCCTCTTa                                >  1:1964927/1‑42 (MQ=39)                                                 aaGGCATACG‑GCGGTGATCCTCGACTCCACATTTGCCTCTTa                                >  1:1224483/1‑42 (MQ=39)                                                 aaGGCATACG‑GCGGTGATCCTCGACTCCACATTTGCCTCTTa                                >  1:2690089/1‑42 (MQ=39)                                                 aaGGCATACG‑GCGGTGATCCTCGACTCCACATTTGCCTCTTa                                >  1:3245891/1‑42 (MQ=39)                                                  aGGCATACG‑GCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCgg  <  1:3274361/71‑1 (MQ=255)                                                   ggCATACG‑GCGGTGATCCTCGACTCCACATTTGCCTCTTa                                >  1:2342042/1‑40 (MQ=39)                                                      aTACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGc                             <  1:2504908/41‑1 (MQ=255)                                                      aTACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGc                             <  1:3679408/41‑1 (MQ=255)                                                      aTACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGc                             <  1:3479510/41‑1 (MQ=255)                                                       tACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGccaaccaa             <  1:3720552/56‑1 (MQ=255)                                                       tACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGccaaccaa             <  1:399035/56‑1 (MQ=255)                                                       tACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGccaaccaa             <  1:63152/56‑1 (MQ=255)                                                       tACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGccaaccaa             <  1:789175/56‑1 (MQ=255)                                                       tACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGccaaccaa             <  1:3190814/56‑1 (MQ=255)                                                         cGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAAt           >  1:2773739/1‑56 (MQ=255)                                                         cGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAAt           >  1:2777950/1‑56 (MQ=255)                                                         cGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAAt           >  1:3602443/1‑56 (MQ=255)                                                           tGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGccaaccaa             <  1:2113827/52‑1 (MQ=255)                                                           tGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGccaaccaa             <  1:2718576/52‑1 (MQ=255)                                                           tGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGccaaccaa             <  1:1716436/52‑1 (MQ=255)                                                           |                                                               GTATTGGCGGGGCGAATATTCTGGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGG  >  W3110S.gb/2663512‑2663632

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