Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	725,902	G→A	intergenic (+28/+560)	SAUSA300_RS03485 → / ← SAUSA300_RS03490	inorganic phosphate transporter/peptidase M23

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	725,902	0	G	A	96.8%	88.6 / ‑5.9	31	intergenic (+28/+560)	SAUSA300_RS03485/SAUSA300_RS03490	inorganic phosphate transporter/peptidase M23
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (0/30);  minor base C (0/1);  total (0/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTACATTACCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAGTATCACAAACAT  >  NC_007793/725809‑725914                                                                                              |             attacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAA‑‑‑AAACTGAACTTCAATACcaccag     <  1:2398705/100‑4 (MQ=38)  ttacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATacaaca       <  1:3711469/100‑6 (MQ=255)  ttacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcacc       <  1:5032088/100‑2 (MQ=255)  ttacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcacc       <  1:4646429/100‑2 (MQ=255)  ttacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcacc       <  1:5258553/100‑2 (MQ=255)  ttacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcacc       <  1:5487610/100‑2 (MQ=255)  ttacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcacc       <  1:5961002/100‑2 (MQ=255)  ctacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcacc       <  1:530880/99‑2 (MQ=255)   tacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAGCTAAACTGAACTTCAATACcacca      <  1:1958248/100‑3 (MQ=255)   tacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCACTAacacca      <  1:1483243/100‑3 (MQ=255)   tacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcacca      <  1:3295999/100‑3 (MQ=255)   tacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTA‑TTTTCTAATTGAAAATGAAACTAAACTGAACTTCAATACcaccat     <  1:2318516/100‑4 (MQ=255)    acattacCTATTTCAGCATTGTTAGCAGGTTTGCTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccat     <  1:1368321/100‑4 (MQ=255)    acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTTAACTTCAATACcaccat     <  1:3008905/100‑4 (MQ=255)    acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccat     <  1:3168320/100‑4 (MQ=255)    acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccat     <  1:4704361/100‑4 (MQ=255)    acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccat     <  1:3999983/100‑4 (MQ=255)    acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccat     <  1:3762940/100‑4 (MQ=255)    acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccat     <  1:990413/100‑4 (MQ=255)    acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccat     <  1:3633771/100‑4 (MQ=255)    acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccat     <  1:1137024/100‑4 (MQ=255)    acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccat     <  1:2990621/100‑4 (MQ=255)    acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccat     <  1:2180272/100‑4 (MQ=255)    acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccat     <  1:1922207/100‑4 (MQ=255)    acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccat     <  1:1639188/100‑4 (MQ=255)    acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccat     <  1:1508591/100‑4 (MQ=255)    acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccat     <  1:1281483/100‑4 (MQ=255)    acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccac     <  1:1672141/100‑4 (MQ=255)    acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATCGAAAATAAAACTAAACTGAACTTCAATACcaccat     <  1:4805198/100‑4 (MQ=255)    acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTAATCTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccat     <  1:5566853/100‑4 (MQ=255)       ttacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccatatg  <  1:2818589/100‑7 (MQ=39)                                                                                              |             ATTACATTACCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAGTATCACAAACAT  >  NC_007793/725809‑725914

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑