Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,632,043	T→C	intergenic (‑81/+28)	SAUSA300_RS08065 ← / ← SAUSA300_RS08070	integrase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,632,043	0	T	C	100.0%	39.6 / NA	17	intergenic (‑81/+28)	SAUSA300_RS08065/SAUSA300_RS08070	integrase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (0/17);  total (0/17)

TTTTTTC‑TTTTTTTCATGTATATCCTCCTTTACTTTTAAAAAGGTTATTTATAAAATAAGTTATAAGAGGATATATGTTCATCAACTCCCGATACCATCGCTAAAGCGTATAGGGAGTTGATTCTTTATTTTTTCATATTATTTTCTTCATAAACTGGTTTA  >  NC_007793/1631946‑1632107                                                                                                   |                                                                 ttttttc‑tttttttcATGTATATCCTCCTTTACTTTTAAAAAGGTTATTTATAAAATAAGTTATAAGAGGATATATGTTCATCAACTCCCGATACCACCg                                                                <  1:5242484/100‑1 (MQ=255) ttttttc‑tttttttcATGTATATCCTCCTTTACTTTTAAAAAGGTTATTTATAAAATAAGTTATAAGAGGATATATGTTCATCAACTCCCGATACCACCg                                                                <  1:462544/100‑1 (MQ=255)   ttttcttttttttCATGTATATCCTCCTTTACTTTTAAAAAGGTTATTTATAAAATAAGTTATAAGAGGATATATGTTCATCAACTCCCGATACCACCGc                                                               <  1:3080082/100‑1 (MQ=255)   ttttc‑tttttttcATGTATATCCTCCTTTACTTTTAAAAAGGTTATTTATAAAATAAGTTATAAGAGGATATATGTTCATCAACTCCCGATACCACCGCt                                                              <  1:2078824/100‑1 (MQ=255)    tttc‑tttttttcATGTATATCCTCCTTTACTTTTAAAAAGGTTATTTATAAAATAAGTTATAAGAGGATATATGTTCATCAACTCCCGATACCACCGCTa                                                             <  1:1109718/100‑1 (MQ=255)    tttc‑tttttttcATGTATATCCTCCTTTACTTTTAAAAAGGTTATTTATAAAATAAGTTATAAGAGGATATATGTTCATCAACTCCCGATACCACCGCTa                                                             <  1:2581044/100‑1 (MQ=255)     ttc‑tttttttcATGTATATCCTCCTTTACTTTTAAAAAGGTTATTTATAAAATAAGTTATAAGAGGATATATGTTCATCAACTCCCGATACCACCGCTaa                                                            <  1:1637170/100‑1 (MQ=255)      cc‑tttttttcATGTATATCCTCCTTTACTTTTAAAAAGGTTATTTATAAAATAAGTTATAAGAGGATATATGTTCATCAACTCCCGATACCACCGCTaaa                                                           <  1:1030305/99‑1 (MQ=255)            ttttcATGTATATCCTCCTTTACTTTTAAAAAGGTTATTTATAAAATAAGTTATAAGAGGATATATGTTCATCAACTCCCGATACCACCGCTAAAGCGta                                                      <  1:3783492/100‑1 (MQ=255)             tttcATGTATATCCTCCTTTACTTTTAAAAAGGTTATTTATAAAATAAGTTATAAGAGGATATATTTTCATCAACTCCCGATACCACCGCTAAAGCGtat                                                     <  1:5855131/100‑1 (MQ=255)             tttcATGTATATCCTCCTTTACTTTTAAAAAGGTTATTTATAAAATAAGTTATAAGAGGATATATGTTCATCAACTCCCGATACCACCGCTAAAGCGtat                                                     <  1:4413367/100‑1 (MQ=255)             tttcATGTATATCCTCCTTTACTTTTAAAAAGGTTATTTATAAAATAAGTTATAAGAGGATATATGTTCATCAACTCCCGATACCACCGCTAAAGCGtat                                                     <  1:5720737/100‑1 (MQ=255)              ttcATGTATATCCTCCTTTACTTTTAAAAAGGTTATTTATAAAATAAGTTATAAGAGGATATATGTCCATCAACTCCCGATACCACCGCTAAAGCGtata                                                    <  1:5091144/100‑1 (MQ=255)                                                            aGTTATAAGAGGATATATGTTCATCAACTCCCGATACCACCGCTAAAGCGTATAGGGAGTTGATTCTTTATTTTTTCATATTATTTTCTTCATAAACTgg      <  1:4088133/100‑1 (MQ=255)                                                              ttATAAGAGGATATATGTTCATCAACTCCCGATACCACCGCTAAAGCGTATAGGGAGTTGATTCTTTATTTTTTCATATTATTTTCTTCATAAACTgggg    <  1:2675411/100‑3 (MQ=255)                                                               tataAGAGGATATATGTTCATCAACTCCCGATACCACCGCTAAAGCGTATAGGGAGTTGATTCTTTATTTTTTCATATTATTTTCTTCATAAACTGGtta   <  1:3773963/100‑2 (MQ=255)                                                                ataAGAGGATATATGTTCATCAACTCCCGATACCACCGCTAAAGCGTATAGGGAGTTGATTCTTTATTTTTTCATATTATTTTCTTCATAAACTGGtttg  <  1:324834/100‑2 (MQ=255)                                                                                                   |                                                                 TTTTTTC‑TTTTTTTCATGTATATCCTCCTTTACTTTTAAAAAGGTTATTTATAAAATAAGTTATAAGAGGATATATGTTCATCAACTCCCGATACCATCGCTAAAGCGTATAGGGAGTTGATTCTTTATTTTTTCATATTATTTTCTTCATAAACTGGTTTA  >  NC_007793/1631946‑1632107

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