Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	658,465	T→A	intergenic (+30/‑266)	SAUSA300_RS03130 → / → SAUSA300_RS03135	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	658,465	0	T	A	90.0%	22.1 / ‑5.0	10	intergenic (+30/‑266)	SAUSA300_RS03130/SAUSA300_RS03135	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (0/9);  minor base C (1/0);  total (1/9)
Fisher's exact test for biased strand distribution p-value = 1.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGAATATGCATTCAATTATAGATAAAAAAGTACATATATTCTTAAAGTCTAATAAATAGTTACAAATTTAGTTAGTTTTCAATTGTAATTTAGGGGTGGTAAACAGTGCTTTGTGAATCTAGACAAATTTATAAAAATCCTAAATATCGAGTTATTA  >  NC_007793/658376‑658533                                                                                          |                                                                     atgAATATGCATTCAATTATAGATAAAAAAGTACATATATTCTTAAAGTCTAATAAATAGTTACAAATTTAGTTAGTTTTCAATTGTTAATTAGGggtgg                                                            <  1:5664653/100‑1 (MQ=255)  tgAATATGCATTCAATTATAGATAAAGAAGTACATATATTCTTAAAGTCTAATAAATAGTTACAAATTTAGTTAGATTTCAATTGTTAATTAGGggtgag                                                           <  1:86652/100‑3 (MQ=255)   gAATATGCATTCAATTATAGATAAAAAAGTACATATATTCTTAAAGTCTAATAAATAGTTACAAATTTAGTTAGTTTTCAATTGTTAATTAGGGGTGGTa                                                          <  1:2441953/100‑1 (MQ=255)   gAATATGCATTCAATTATAGATAAAAAAGTACATATATTCTTAAAGTCTAATAAATAGTTACAAATTTAGTTAGTTTTCAATTGTTAATTAGGGGTGGTa                                                          <  1:4457993/100‑1 (MQ=255)    aaTATGCATTCAATTATAGATAAAAAAGTACATATATTCTTAAAGTCTAATAAATAGTTACAAATTTAGTTAGTTTTCAATTGTTAATTAGGGGTGGTaa                                                         <  1:882144/100‑1 (MQ=255)                      gATAAAAAAGTACATATATTCTTAAAGTCTAATAAATAGTTACAAATTTAGTTAGTTTTCAATTGTTAATTAGGGGTGGTAAACAGTGCTTTGTGAATCt                                       <  1:3581861/100‑1 (MQ=255)                      gATAAAAAAGTACATATATTCTTAAAGTCTAATAAATAGTTACAAATTTAGTTAGTTTTCAATTGTTAATTAGGGGTGGTAAACAGTGCTTTGTGAATCt                                       <  1:655507/100‑1 (MQ=255)                      agataaaaaGTACATATATTCTTAAAGTCTAATAAATAGTTACAAATTTAGTTAGTTTTCAATTGTTAATTAGGGGTGGTAAACAGTGCTTTGTGAATCt                                       <  1:3423708/96‑1 (MQ=255)                                                    aataaataGTTACAAATTTAGTTAGTTTTCAATTGTTACTTAGGGGTGGTAAACAGTGCTTTGTGAATCTAGACAAATTTATAAAAATCCTAAATATCGa         >  1:4409385/1‑100 (MQ=255)                                                           aGTTACAAATTTAGTTAGTTTTCAATTGTTAATTAGGGGTGGTAAACAGTGCTTTGTGAATCTAGACAAATTTATAAAAATCCTAAATATCGAGttacag  <  1:4580670/100‑4 (MQ=255)                                                                                          |                                                                     ATGAATATGCATTCAATTATAGATAAAAAAGTACATATATTCTTAAAGTCTAATAAATAGTTACAAATTTAGTTAGTTTTCAATTGTAATTTAGGGGTGGTAAACAGTGCTTTGTGAATCTAGACAAATTTATAAAAATCCTAAATATCGAGTTATTA  >  NC_007793/658376‑658533

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