Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	2,518,015	A→G	A46A (GCT→GCC)	mntH ←	manganese/divalent cation transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,518,015	0	A	G	98.6%	231.3 / ‑4.7	69	A46A (GCT→GCC)	mntH	manganese/divalent cation transporter
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (0/68);  total (1/68)
Fisher's exact test for biased strand distribution p-value = 1.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

AGGTTGGCCCAAACGACAACCCACAGTAGCTGATAGCCGAAGCTAGCACCCGCCTGAATATT  >  W3110S.gb/2517971‑2518032                                             |                  aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:2578471/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:287668/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:2836852/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:2749577/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:2744344/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:2722705/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:2680716/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:2642002/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:258585/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:31023/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:2558071/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:2475183/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:2466262/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:2451302/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:2414398/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:2386235/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:1069352/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:2357191/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:334587/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:412352/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:429972/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:438397/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:44442/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:444700/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:459067/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:484513/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:5709/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:649743/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:762839/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:814803/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:845421/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:913742/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:950994/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:1451697/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:2316994/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:1647932/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:1611185/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:1600778/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:1574520/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:1522152/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:1521578/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:146659/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:1720319/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:1401813/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:1244804/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:1242510/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:11838/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:1182448/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:1091949/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:1088408/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:20382/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:2293726/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:2220761/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:2209954/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:217451/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:2136093/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:2133168/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:20220/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:1993070/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:1923402/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:1892543/62‑1 (MQ=255) aGGTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCC‑CCTGAATAtt  <  1:2485843/61‑1 (MQ=255) aGGTTAGCCCAAACGACAACCAACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:1611136/62‑1 (MQ=255)  ggTTGGCCCAAACGACAACCCACAGTAGCTGATAGCCGAAGCTAGCa                >  1:2292540/1‑47 (MQ=255)  ggTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:2017609/61‑1 (MQ=255)  ggTTAGCCCAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:1741157/61‑1 (MQ=255)          cAAACGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:2837517/53‑1 (MQ=255)           aaaCGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:1844137/52‑1 (MQ=255)            aaCGACAACCCACAATAACTGATAGCCGAAGCTGGCACCCGCCTGAATAtt  <  1:807877/51‑1 (MQ=255)                                             |                  AGGTTGGCCCAAACGACAACCCACAGTAGCTGATAGCCGAAGCTAGCACCCGCCTGAATATT  >  W3110S.gb/2517971‑2518032

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