Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,594,183	A→G	20.5%	D28D (GAT→GAC)	yihE ←	predicted kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,594,183	0	A	G	20.5%	77.4 / 17.1	39	D28D (GAT→GAC)	yihE	predicted kinase
Reads supporting (aligned to +/- strand):  ref base A (19/12);  new base G (6/2);  total (25/14)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

GACGATCTTCGTCCTGAAATTGATAGACACGGTTTTCATAGCTGTTAAGCGGGGTAAGACCGGAATCCACCCGGATCCCATGCTCAAACAGAGCGTCCATGATGGTATCCGGGTGTAGTGTCTGGAAAG  >  W3110S.gb/3594119‑3594247                                                                 |                                                                 gacgaTCTTCGTCCTGAAATTGATAGACACGGTTTTCATAGCTGTTAAGCGGGGTAAGACCGGAATCCAcc                                                            <  1:3983611/71‑1 (MQ=255)    gaTCTTCGTCCTGAAATTGATAGACACGGTTTTCATAGCTGTTAAGCGGGGTAAGACCGGAATCCACCCg                                                          >  1:2651485/1‑70 (MQ=255)    gaTCTTCGTCCTGAAATTGATAGACACGGTTTTCATAGCTGTTAAGCGGGGTAAGACCGGAATCCACCCg                                                          >  1:510348/1‑70 (MQ=255)    gaTCTTCGTCCTGAAATTGATAGACACGGTTTTCATAGCTGTTAAGCGGGGTAAGACCGGAATCCACCCg                                                          >  1:387658/1‑70 (MQ=255)      tCTTCGTCCTGAAATTGATAGACACGGTTTTCATAGCTGTTAAGCGGGGTAAGACCGGAATCCACCCGGa                                                        >  1:2501874/1‑70 (MQ=255)      tCTTCGTCCTGAAATTGATAGACACGGTTTTCATAGCTGTTAAGCGGGGTAAGACCGGAATCCACCCGGAt                                                       >  1:107794/1‑71 (MQ=255)          cGTCCTGAAATTGATAGACACGGTTTTCATAGCTGTTAAGCGGGGTAAGACCGGAATCCACCCGGATCCCa                                                   >  1:256936/1‑71 (MQ=255)              cTGAAATTGATAGACACGGTTTTCATAGCTGTTAAGCGGGGTAAGACCGGAATCCACCCGGAta                                                      >  1:2329020/1‑63 (MQ=255)              cTGAAATTGATAGACACGGTTTTCATAGCTGTTAAGCGGGGTAAGACCGGAATCCACCCGGATCCCATGCt                                               >  1:3799695/1‑71 (MQ=255)                gAAATTGATAGACACGGTTTTCATAGCTGTTAAGCGGGGTAAGACCAGAGTCCACCCGGATCCCATGCTc                                              >  1:2000670/1‑70 (MQ=255)                 aaaTTGATAGACACGGTTTTCATAGCTGTTAAGCGGGGTAAGACCGGAATCCACCCGGATCCCATGCTCa                                             <  1:3987757/70‑1 (MQ=255)                 aaaTTGATAGACACGGTTTTCATAGCTGTTAAGCGGGGTAAGACCGGAATCCACCCGGATCCCATGCTCa                                             <  1:2245443/70‑1 (MQ=255)                   aTTGATAGACACGGTTTTCATAGCTGTTAAGCGGGGTAAGACCGGAATc                                                                <  1:1262901/49‑1 (MQ=255)                   aTTGATAGACACGGTTTTCATAGCTGTTAAGCGGGGTAAGACCGGAATCCACCCGGATCCCATGCTCaaa                                           >  1:716998/1‑70 (MQ=255)                    ttGATAGACACGGTTTTCATAGCTGTTAAGCGGGGTAAGACCGGAATc                                                                <  1:2028928/48‑1 (MQ=255)                    ttGATAGACACGGTTTTCATAGCTGTTAAGCGGGGTAAGACCGGAATc                                                                <  1:1274193/48‑1 (MQ=255)                               ggTTTTCATAGCTGTTAAGCGGGGTAAGACCGGAATCCACCCGGATCCCATGCTCAAACAGAGCGTCCat                               <  1:2355548/70‑1 (MQ=255)                                gTTTTCATAGCTGTTAAGCGGGGTAAGACCGGAATCCACCCGGATCCCATGCTCAAACAGAGCGTCCatga                             <  1:3022678/71‑1 (MQ=255)                                  tttCATAGCTGTTAAGCGGGGTAAGACCGGAATCCACCCGGATCCCATGCTCAAACAGAGCGTCCatgatg                           >  1:943824/1‑71 (MQ=255)                                  tttCATAGCTGTTAAGCGGGGTAAGACCGGAATCCACCCGGATCCCATGCTCAAACAGAGCGTCCatgatg                           >  1:3744011/1‑71 (MQ=255)                                  tttCATAGCTGTTAAGCGGGGTAAGACCGGAATCCACCCGGATCCCATCCTCAAACAGAGCGTCCatgatg                           >  1:3541424/1‑71 (MQ=255)                                   ttCATAGCTGTTAAGCGGGGTAAGACCGGAATCCACCCGGATCCCATGCTCAAACAGAGCGTCCATGATgg                          >  1:1585592/1‑71 (MQ=255)                                     cATAGCTGTTAAGCGGGGTAAGACCGGAATCCACCCGGATCCCATGCTCAAACAGAGCGTCCATGATGGTa                        >  1:3474807/1‑71 (MQ=255)                                         gCTGTTAAGCGGGGTAAGACCGGAGTCCACCCGGATCCCa                                                   >  1:836021/1‑40 (MQ=255)                                         gCTGTTAAGCGGGGTAAGACCGGAGTCCACCCGGATCCCa                                                   >  1:663169/1‑40 (MQ=255)                                         gCTGTTAAGCGGGGTAAGACCGGAGTCCACCCGGATCCCAt                                                  >  1:3873125/1‑41 (MQ=255)                                         gCTGTTAAGCGGGGTAAGACCGGAGTCCACCCGGATCCCAt                                                  >  1:637582/1‑41 (MQ=255)                                          cTGTTAAGCGGGGTAAGACCGGAGTCCACCCGGATCCCATGCTCAAACAGAGCGTCCATGATGGTATCCgg                   >  1:1698906/1‑71 (MQ=255)                                             ttAAGCGGGGTAAGCCCGGAATCCACCCGGATCCCATGCTCAAACAGAGCTTCCATGATGGTATCCGGgtg                <  1:2926180/71‑1 (MQ=255)                                             ttAAGCGGGGTAAGACCGGAGTCCACCCGGATCCCAt                                                  <  1:2405121/37‑1 (MQ=255)                                             ttAAGCGGGGTAAGACCGGAATCCACCCGGATCCCATGCTCAAACAGAGCGTCCATGATGGTATCCGGgtg                <  1:1927028/71‑1 (MQ=255)                                              tAAGCGGGGTAAGACCGGAATCCACCCGGATCCCATGCTCAAACAGAGCGTCCATGATGGTATCCGGgtg                >  1:3219394/1‑70 (MQ=255)                                               aaGCGGGGTAAGACCGGAATCCACCCGGATCCCATGCTCAAACAGAGCGTCCATGATGGTATCCGGgtg                >  1:3809797/1‑69 (MQ=255)                                               aaGCGGGGTAAGACCGGAATCCACCCGGATCCCATGCTCAAACAGAGCGTCCATGATGGTATCCGGgtg                >  1:315096/1‑69 (MQ=255)                                                aGCGGGGTAAGACCGGAATCCACCCGGATCCCATGCTCAAACAGAGCGTCCATGATGGTATCCGGGTGTa              <  1:4023419/70‑1 (MQ=255)                                                     ggTAAGACCGGAATCCACCCGGATCCCATGCTCAAACAGAGCGTc                                  >  1:1695465/1‑45 (MQ=255)                                                      gTAAGACCGGAGTCCACCCGGATCCCATGCTCAAACAGAGCGTCCatg                              <  1:418231/48‑1 (MQ=255)                                                       tAAGACCGGAATCCACCCGGATCGCATGCTCAAACAGAGCGTCCATGATGGTATCCGGgtg                <  1:1672358/61‑1 (MQ=255)                                                           aCCGGAATCCACCCGGATCCCATGCTCAAACAGAGCGTCCATGATGGTATCCGGGTGTAGTGTCTGGAAAg  >  1:2664160/1‑71 (MQ=255)                                                                 |                                                                 GACGATCTTCGTCCTGAAATTGATAGACACGGTTTTCATAGCTGTTAAGCGGGGTAAGACCGGAATCCACCCGGATCCCATGCTCAAACAGAGCGTCCATGATGGTATCCGGGTGTAGTGTCTGGAAAG  >  W3110S.gb/3594119‑3594247

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