Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,585,874	G→A	D78D (GAC→GAT)	SAUSA300_RS07790 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,585,874	0	G	A	94.4%	44.4 / ‑4.0	18	D78D (GAC→GAT)	SAUSA300_RS07790	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (5/12);  minor base C (0/1);  total (5/13)
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

CTCACCAAATTTGTAATACATTGTTTTAGCCCTCCTTGTTATCATCAATACCGTGAAATTTTTGTGATTTACACATTTGGAGAACATTGACAATGTCTTTATAACTCTTAGTGCTATCCAATAAGGAAGCAAGATCGAAAGTATGACCAATCACAGAACTTG  >  NC_007793/1585780‑1585941                                                                                               |                                                                    gtcacccgttttGTAATACATTGTTTTAGTCCTCCTTGTCGTCATCTATACCGAGAATTTTTTGTGATTTACACATTTGGAGAACATTGACAATATCttt                                                                >  1:3259601/10‑100 (MQ=18)                   cATTGTTTTAGTCCTCCTTGTCGTCATCTATACCGAGAATTTTTTGTGATTTACACATTTGGAGAACATTGACAATATCTTTATAACTCTTAGTGCTATc                                              >  1:2214343/1‑100 (MQ=255)                                gatcttggcgtcatctATACCGAGAATTTTTTGTGATTTACACATTTGGAGAACATTGACAATATCTTTATAACTCTTAGTGCTAGCCAATAAGTAAGCa                                 <  1:4391051/90‑1 (MQ=255)                                gatctTGTCGTCATCTATACCGAGAATTTTTTGTGATTTACACATTTGGAGAACATTGACAATATCTTTATAACTCTTAGTGCTATCCAATAAGTAAGCa                                 <  1:2865473/97‑1 (MQ=255)                                gatctTGTCGTCATCTATACCGAGAATTTTTTGTGATTTACACATTTGGAGAACATTGACAATATCTTTATAACTCTTAGTGCTATCCAATAAGTAAGCa                                 <  1:4950536/97‑1 (MQ=255)                                ctcctTGTCGTCATCTATACCGAGAATTTTTTGTGATTTACACATTTGGAGAACATTGACAATATCTTTATAACTCTTAGTGCTATCCAATAAGTAAGCa                                 <  1:2790828/100‑1 (MQ=255)                                ctcctTGTCGTCATCTATACCGAGAATTTTTTGTGATTTACACATTTGGAGAACATTGACAATATCTTTATAACTCTTAGTGCTATCCAATAAGTAAGCa                                 <  1:3312560/100‑1 (MQ=255)                                ctcctTGTCGTCATCTATACCGAGAATTTTTTGTGATTTACACATTTGGAGAACATTGACAATATCTTTATAACTCTTAGTGCTATCCAATAAGTAAGCa                                 <  1:415331/100‑1 (MQ=255)                                ctcctTGTCGTCATCTATACCGAGAATTTTTTGTGATTTACACATTTGGAGAACATTGACAATATCTTTATAACTCTTAGTGCTATCCAATAAGTAAGCa                                 <  1:5462929/100‑1 (MQ=255)                                ctcctTGTCGTCATCTATACCGAGAATTTTTTGTGATTTAAAAAGGGGGAGAACATTGACAATCTCTTTATAACTCTTAGTGCTATCCAATAAGTAAGCa                                 <  1:228697/100‑1 (MQ=255)                                 atctTGTCGTCATCTATACCGAGAATTTTTTGTGATTTACACATTTGGAGAACATTGACAATATCTTTATAACTCTTAGTGCTATCCAATAAGTAAGCaa                                <  1:641682/98‑1 (MQ=255)                                           catctATACCGAGAATTTTTTGTGATTTACACATTTGGAGAACATTGACAATATCTTTATAACTCTTAGTGCTATCCAATAAGTAAGCAAGATCAAAAGt                      <  1:1905053/100‑1 (MQ=255)                                              ctaTACCGAGAATTTTTTGTGATGTACACATTTGGAGAACATTGACAATATCTTTATAACGCTTAGTGCTATCCAATAAGTAAGCAAGATCAAAAGTATg                   <  1:5793291/98‑1 (MQ=255)                                                aTACCGAGAATTTTTTGTGATTTACACATTTGGAGAACATTGACAATATCTTTATAACTCTTAGTGCTATCCAATAAGTAAGCAAGATCGGAAgagcaca                 >  1:5532923/1‑94 (MQ=255)                                                aTACCGAGAATTTTTTGTGATTTACACATTTGGAGAACATTGACAATATCTTTATAACTCTTAGTGCTATCCAATAAGCAGGCAAGATCGGAAgagcaca                 >  1:4723013/1‑94 (MQ=255)                                                  aCCGAGAATTTTTTGTGATTTACACATTTGGAGAACATTGACAATATCTTTATAACTCTTAGTGCTATCCAATAAGTAAGCAAGATCAAAAGTATGACCa               >  1:4581521/1‑100 (MQ=255)                                                       gAATTTTTTGTGATTTACACATTTGGAGAACATTGACAATATCTTTATAACTCTTAGTGCTATCCAATAAGTAAGCAAGATCAAAAGTATGACCAATcac          <  1:5452040/100‑1 (MQ=255)                                                               tgtgATTTACACATTTGGAGAACATTGACAATATCTTTATAACTCTTAGTGCTATCCAATAAGTAAGCAAGATCAAAAGTATGACCAATCACAGAACTTg  <  1:2408189/100‑1 (MQ=255)                                                                                               |                                                                    CTCACCAAATTTGTAATACATTGTTTTAGCCCTCCTTGTTATCATCAATACCGTGAAATTTTTGTGATTTACACATTTGGAGAACATTGACAATGTCTTTATAACTCTTAGTGCTATCCAATAAGGAAGCAAGATCGAAAGTATGACCAATCACAGAACTTG  >  NC_007793/1585780‑1585941

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