Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	1,380,993	C→T	37.5%	H466Y (CAT→TAT)	ycjT →	predicted hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,380,993	0	C	T	37.5%	30.1 / 38.7	40	H466Y (CAT→TAT)	ycjT	predicted hydrolase
Reads supporting (aligned to +/- strand):  ref base C (4/21);  new base T (6/9);  total (10/30)
Fisher's exact test for biased strand distribution p-value = 1.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

GCGCTACTTCTGGAGACGGCAAAGTTCTGGATTAGCCGCGCGGTGAGAGTTAACGATCGTCTGGAAATTCATGATGTTATTGGGCCAGACGAATATACCGAACATGTCAATAATAATGCATACACCAGCTATAT  >  W3110S.gb/1380924‑1381057                                                                      |                                                                 gcgcTACTTCTGGAGACGGCAAAGTTCTGGATTAGCCGCGCGGTGAGAGTTAACGATCGTCTGGAAATTc                                                                  <  1:3215650/70‑1 (MQ=255) gcgcTACTTCTGGAGACGGCAAAGTTCTGGATTAGCCGCGCGGTGAGAGTTAACGATCGTCTGGAAATTc                                                                  <  1:2453315/70‑1 (MQ=255) gcgcTACTTCTGGAGACGGCAAAGTTCTGGATTAGCCGCGCGGTGAGAGTTAACGATCGTCTGGAAATTc                                                                  <  1:2658511/70‑1 (MQ=255) gcgcTACTTCTGGAGACGGCAAAGTTCTGGATTAGCCGCGCGGTGAGAGTTAACGATCGTCTGGAAATTc                                                                  <  1:75281/70‑1 (MQ=255) gcgcTACTTCTGGAGACGGCAAAGTTCTGGATTAGCCGCGCGGTGAGAGTTAACGATCGTCTGGAAATTc                                                                  <  1:1455880/70‑1 (MQ=255) gcgcTACTTCTGGAGACGGCAAAGTTCTGGATTAGCCGCGCGGTGAGAGTTAACGATCGTCTGGAAATTc                                                                  <  1:3300927/70‑1 (MQ=255) gcgcTACTTCTGGAGACGGCAAAGTTCTGGATTAGCCGCGCGGTGAGAGTTAACGATCGTCTGGAAATTc                                                                  <  1:1035623/70‑1 (MQ=255) gcgcTACTTCTGGAGACGGCAAAGTTCTGGATTACCCGCGCGGTGAGAGTTAACGATCGTCTGGAAATTc                                                                  <  1:2998563/70‑1 (MQ=255) gcgcTACTTATGGAGACGGCAAAGTTCTGGATTAGCCGCGCGGTGAGAGTTAACGATCGTCTGGAAATTc                                                                  <  1:456213/70‑1 (MQ=255)       cTTCTGGAGACGGCAAAGTTCTGGATTAGCCGCGCGGTGAGAGTTAACGATCGTCTGGAAATTCATGATGt                                                           <  1:747374/71‑1 (MQ=255)       cTTCTGGAGACGGCAAAGTTCTGGATTAGCCGCGCGGTGAGAGTTAACGATCGTCTGGAAATTCATGATGt                                                           <  1:2857700/71‑1 (MQ=255)       cTTCTGGAGACGGCAAAGTTCTGGATTAGCCGCGCGGTGAGAGTTAACGATCGTCTGGAAATTCATGATGt                                                           <  1:1128424/71‑1 (MQ=255)       cTTCTGGAGACGGCAAAGTTCTGGATTAGCCGCGCGGTGAGAGTTAACGATCGTCTGGAAATTCATGATGt                                                           <  1:1745274/71‑1 (MQ=255)       cTTCTGGAGACGGCAAAGTTCTGGATTAGCCGCGCGGTGAGAGTTAACGATCGTCTGGAAATTCATGATGt                                                           <  1:590778/71‑1 (MQ=255)       cTTCTGGAGACGGCAAAGTTCTGGATTAGCCGCGCGGTGAGAGTTAACGATCGTCTGGAAATTCATGATGt                                                           <  1:2604820/71‑1 (MQ=255)       cTTCTGGAGACGGCAAAGTTCTGGATTAGCCGCGCGGTGAGAGTTAACGATCGTCTGGAAATTCATGATGt                                                           <  1:2120150/71‑1 (MQ=255)                 cGGCAAAGTTCTGGATTAGGCGCGCGGTGAGAGTTAACGATCGTCTGGAAATTCat                                                                <  1:2433401/56‑1 (MQ=255)                 cGGCAAAGTTCTGGATTAGCCGCGCGGTGAGAGTTAACGATCGTCTGGAAATTCat                                                                <  1:2441557/56‑1 (MQ=255)                 cGGCAAAGTTCTGGATTAGCCGCGCGGTGAGAGTTAACGATCGTCTGGAAATTCat                                                                <  1:382667/56‑1 (MQ=255)                 cGGCAAAGTTCTGGATTAGCCGCGCGGTGAGAGTTAACGATCGTCTGGAAATTCat                                                                <  1:721287/56‑1 (MQ=255)                 cGGCAAAGTTCTGGATTAGCCGCGCGGTGAGAGTTAACGATCGTCTGGAAATTCat                                                                <  1:2844527/56‑1 (MQ=255)                 cGGCAAAGTTCTGGATTAGCCGCGCGGTGAGAGTTAACGATCGTCTGGAAATTCat                                                                <  1:2865091/56‑1 (MQ=255)                  ggCAAAGTTCTGGATTAGCCGCGCGGTGAGAGTTAACGATCGTCTGGAAATTCATGATGTTATTGGGc                                                   >  1:1374341/1‑68 (MQ=255)                       aGTTCTGGATTAGCGGCGCGGTGAGAGTTAACGATCGTCTGGAAATTCATGATGTTATTGGGCCAGACGaa                                           <  1:1932050/71‑1 (MQ=255)                       aGTTCTGGATTAGCCGCGCGGTGAGAGTTAACGATCGTCTGGAAATTCATGATGTTATTGGGCCAGACGaa                                           <  1:3246456/71‑1 (MQ=255)                       aGTTCTGGATTAGCCGCGCGGTGAGAGTTAACGATCGTCTGGAAATTCATGATGTTATTGGGCCAGACGaa                                           <  1:1605639/71‑1 (MQ=255)                                          ggTGAGAGTTAACGATCGTCTGGAAATTTATGATGTTATTGGGCCAGACGAATATACCGAACATGTCaata                        <  1:334096/71‑1 (MQ=255)                                          ggTGAGAGTTAACGATCGTCTGGAAATTTATGATGTTATTGGGCCAGACGAATATACCGAACATGTCaata                        <  1:955219/71‑1 (MQ=255)                                          ggTGAGAGTTAACGATCGTCTGGAAATTTATGATGTTATTGGGCCAGACGAATATACCGAACATGTCaata                        <  1:3239611/71‑1 (MQ=255)                                          ggTGAGAGTTAACGATCGTCTGGAAATTTATGATGTTATTGGGCCAGACGAATATACCGAACATGTCaata                        <  1:3227964/71‑1 (MQ=255)                                          ggTGAGAGTTAACGATCGTCTGGAAATTTATGATGTTATTGGGCCAGACGAATATACCGAACATGTCaata                        <  1:1643347/71‑1 (MQ=255)                                          ggTGAGAGTTAACGATCGTCTGGAAATTTATGATGTTATTGGGCCAGACGAATATACCGAACATGTCaata                        <  1:167735/71‑1 (MQ=255)                                          ggTGAGAGTTAACGATCGTCTGGAAATTTAAGATGTTATTGGGCCAGACGAATATACCGAACATGTCaata                        <  1:2603384/71‑1 (MQ=255)                                           gTGAGAGTTAACGATCGTCTGGAAATTTATGATGTTATTGGGCCAGACGAATATACCGAACATGTCaata                        <  1:508613/70‑1 (MQ=255)                                           gTGAGAGTTAACGATCGTCTGGAAATTTATGATGTTATTGGGCCAGACGAATATACCGAACATGTCaata                        <  1:556103/70‑1 (MQ=255)                                             gagaGTTAACGATCGTCTGGAAATTTATGATGTTa                                                         >  1:991796/1‑35 (MQ=255)                                             gagaGTTAACGATCGTCTGGAAATTTATGATGTTa                                                         >  1:513231/1‑35 (MQ=255)                                             gagaGTTAACGATCGTCTGGAAATTTATGATGTTa                                                         >  1:3299568/1‑35 (MQ=255)                                             gagaGTTAACGATCGTCTGGAAATTTATGATGTTa                                                         >  1:3244736/1‑35 (MQ=255)                                             gagaGTTAACGATCGTCTGGAAATTTATGATGTTa                                                         >  1:1504048/1‑35 (MQ=255)                                             gagaGTTAACGATCGTCTGGAAATTTATGATGTTa                                                         >  1:1122095/1‑35 (MQ=255)                                                  ttAACGATCGTCTGGAAATTCATGATGTTATTGGGCCAGACGAATATACCGAACATGTCAATAATAATGCa                <  1:252277/71‑1 (MQ=255)                                                  ttAACGATCGTCTGGAAATTCATGATGTTATTGGGCCAGACGAATATACCGAACATGTCAATAATAATGCa                <  1:1722198/71‑1 (MQ=255)                                                  ttAACGATCGTCTGGAAATTCATGATGTTATTGGGCCAGACGAATATACCGAACATGTCAATAATAATGCa                <  1:1240547/71‑1 (MQ=255)                                                              tGGAAATTCATGATGTTATTGGGCCAGACGAATATACCGAACATGTCAATAATAATGCATACACCAGCta     <  1:2307680/70‑1 (MQ=255)                                                                gAAATTCATGATGTTATTGGGCCAGACGAATATACCGAACATGTCAATAATAATGCATACACCAGCta     <  1:950966/68‑1 (MQ=255)                                                                gAAATTCATGATGTTATTGGGCCAGACGAATATACCGAACATGTCAATAATAATGCATACACCAGCTatat  >  1:2228553/1‑71 (MQ=255)                                                                gAAATTCATGATGTTATTGGGCCAGACGAATATACCGAACATGTCAATAATAATGCATACACCAGCTatat  >  1:183635/1‑71 (MQ=255)                                                                gAAATTCATGATGTTATTGGGCCAGACGAATATACCGAACATGTCAATAATAATGCATACACCAGCTatat  >  1:1481289/1‑71 (MQ=255)                                                                      |                                                                 GCGCTACTTCTGGAGACGGCAAAGTTCTGGATTAGCCGCGCGGTGAGAGTTAACGATCGTCTGGAAATTCATGATGTTATTGGGCCAGACGAATATACCGAACATGTCAATAATAATGCATACACCAGCTATAT  >  W3110S.gb/1380924‑1381057

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