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	92.0%	58.7 / ‑2.1	25	intergenic (+28/+560)	SAUSA300_RS03485/SAUSA300_RS03490	inorganic phosphate transporter/peptidase M23
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (0/23);  total (1/24)
Fisher's exact test for biased strand distribution p-value = 8.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.52e-02

ATTACATTACCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAGTATCACAAACATATGGTGGTATTGAA  >  NC_007793/725809‑725928                                                                                              |                           attacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaa                      <  1:1216389/100‑2 (MQ=255) attacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAA‑‑‑AAACTGAACTTCAATACcaccat                   <  1:337327/100‑4 (MQ=38)  ttacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAGCAAAACTGAACTTCAATACcacc                     <  1:296528/100‑2 (MQ=255)  ttacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcacc                     <  1:740591/100‑2 (MQ=255)  ttacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcacc                     <  1:1547212/100‑2 (MQ=255)  ttacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcacc                     <  1:5535980/100‑2 (MQ=255)  ttacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcacc                     <  1:5316127/100‑2 (MQ=255)  ttacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcacc                     <  1:5218701/100‑2 (MQ=255)   tacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcacca                    <  1:5376851/100‑3 (MQ=255)    acattacCTATTTCGGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccat                   <  1:994712/100‑4 (MQ=255)    acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccat                   <  1:4018427/100‑4 (MQ=255)    acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccat                   <  1:4184094/100‑4 (MQ=255)    acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccat                   <  1:4877304/100‑4 (MQ=255)    acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccat                   <  1:2651304/100‑4 (MQ=255)    acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccat                   <  1:5313711/100‑4 (MQ=255)    acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccat                   <  1:2395040/100‑4 (MQ=255)    acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccat                   <  1:106639/100‑4 (MQ=255)    acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccag                   <  1:2824488/100‑4 (MQ=255)    acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccac                   <  1:1229465/100‑4 (MQ=255)     cattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccata                  <  1:869482/100‑5 (MQ=255)     cattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACCACcaag                  <  1:515602/100‑2 (MQ=255)       ttacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccatatg                <  1:5515359/100‑7 (MQ=39)        tacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAGCAAAACTGAACTTCAATACcaccatacag               <  1:1803784/100‑8 (MQ=255)        tacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAGTATCACAAACata               <  1:449469/100‑1 (MQ=255)                     cctGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAGTATCACAAACATATGGTGGTATTGaa  >  1:2060524/3‑100 (MQ=255)                                                                                              |                           ATTACATTACCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAGTATCACAAACATATGGTGGTATTGAA  >  NC_007793/725809‑725928

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