Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,631,930	G→A	D11D (GAC→GAT)	SAUSA300_RS08065 ←	integrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,631,930	0	G	A	100.0%	109.3 / NA	40	D11D (GAC→GAT)	SAUSA300_RS08065	integrase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (1/39);  total (1/39)

AAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGTTTTTTTC‑TTTTTTTCATGTATATCCTCCTTTACTTTTAAAAAGGTTATTTATAAAATAAGTTATA  >  NC_007793/1631837‑1632010                                                                                              |                                                                                  aaGAACCGCATCATCACGACCTAAAATCTTAGAAAGCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAAt                                                                             <  1:3828313/100‑1 (MQ=255)        gCATCATCACGACCTAAAATCCTAGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTAc                                                                      >  1:2443377/1‑100 (MQ=255)           tcatcaCGACCTAAAATCTTAGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGtt                                                                   <  1:5293796/100‑1 (MQ=255)            catcaCGACCTAAAATCTTAGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGtgg                                                                  <  1:4366518/100‑3 (MQ=255)              tcaCGACCTAAAATCTTAGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGttttt                                                                <  1:2966705/100‑1 (MQ=255)              tcaCGACCTAAAATCTTAGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGttttg                                                                <  1:1986219/100‑2 (MQ=255)               cgcggccTAAAATCTTAGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAAGCTAATTCTTTACGttttct                                                               <  1:5601670/95‑3 (MQ=255)                aCGACCTAAAATCTTAGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTAccgtgtgct                                                              <  1:675358/100‑9 (MQ=255)                            cTTAGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACG‑TTTTTTC‑CTTTTTTCATGt                                                <  1:2038242/100‑1 (MQ=255)                             ttAGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTC‑TTTTTTTCATGt                                                <  1:6161523/100‑1 (MQ=255)                             ttAGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTC‑TTTTTTTCATGt                                                <  1:6039118/100‑1 (MQ=255)                             ttAGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTC‑TTTTTTTCATGt                                                <  1:595323/100‑1 (MQ=255)                             ttAGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTC‑TTTTTTTCATGt                                                <  1:5895464/100‑1 (MQ=255)                             ttAGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTC‑TTTTTTTCATGt                                                <  1:5534897/100‑1 (MQ=255)                             ttAGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATGCTTTACGTTTTTTTC‑TTTTTTTCATGt                                                <  1:42343/100‑1 (MQ=255)                             ttAGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGAAGTGTTTTATAATCTAATTCTTTACGTTTTTTTC‑TTTTTTTCATGt                                                <  1:4976028/100‑1 (MQ=255)                               aGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTC‑TTTTTTTCATGtat                                              <  1:4212023/100‑1 (MQ=255)                               aGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTC‑TTTTTTTCATGtat                                              <  1:6361727/100‑1 (MQ=255)                               aGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTC‑TTTTTTTCATGtat                                              <  1:3695896/100‑1 (MQ=255)                               aGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTC‑TTTTTTTCATGtat                                              <  1:1059394/100‑1 (MQ=255)                               aGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTC‑TTTTTTTCATGtat                                              <  1:5203421/100‑1 (MQ=255)                                 aaaTCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTC‑TTTTTTTCATGtataa                                            <  1:6277892/100‑2 (MQ=255)                                 aaaTCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTC‑TTTTTTTCATGTatat                                            <  1:3344043/100‑1 (MQ=255)                                         atgatgTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTC‑TTTTTTTCATGTATATCCTCCttt                                    <  1:3426637/100‑1 (MQ=255)                                          tgatgTTTAAACCATCATTTAAATACTTTGCGATAGCTAGCCGTGTTTTATAATCTAATTCTTTACGTTTTTTTC‑TTTTTTTCATGTATATCCTCCTTTa                                   <  1:6511169/100‑1 (MQ=255)                                                           tttAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTC‑TTTTTTTCATGTATATCCTCCTTTACTTTTAAAAAGGttaga                  <  1:6076200/100‑3 (MQ=255)                                                                           tAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTC‑TTTTTTTCATGTATATCCTTCTTTACTTTTAAAAAGGTTATTTATAAAATAAGTtata  <  1:5433179/100‑1 (MQ=255)                                                                           tAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTC‑TTTTTTTCATGTATATCCTCCTTTACTTTTAAAAAGGTTATTTATAAAATAAGTtata  <  1:6541255/100‑1 (MQ=255)                                                                           tAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTC‑TTTTTTTCATGTATATCCTCCTTTACTTTTAAAAAGGTTATTTATAAAATAAGTtata  <  1:484311/100‑1 (MQ=255)                                                                           tAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTC‑TTTTTTTCATGTATATCCTCCTTTACTTTTAAAAAGGTTATTTATAAAATAAGTtata  <  1:3562659/100‑1 (MQ=255)                                                                           tAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTC‑TTTTTTTCATGTATATCCTCCTTTACTTTTAAAAAGGTTATTTATAAAATAAGTtata  <  1:602872/100‑1 (MQ=255)                                                                           tAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTC‑TTTTTTTCATGTATATCCTCCTTTACTTTTAAAAAGGTTATTTATAAAATAAGTtata  <  1:3776119/100‑1 (MQ=255)                                                                           tAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTC‑TTTTTTTCATGTATATCCTCCTTTACTTTTAAAAAGGTTATTTATAAAATAAGTtata  <  1:6043858/100‑1 (MQ=255)                                                                           tAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTC‑TTTTTTTCATGTATATCCTCCTTTACTTTTAAAAAGGTTATTTATAAAATAAGTtata  <  1:2175500/100‑1 (MQ=255)                                                                           tAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTC‑TTTTTTTCATGTATATCCTCCTTTACTTTTAAAAAGGTTATTTATAAAATAAGTtata  <  1:6244343/100‑1 (MQ=255)                                                                           tAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTC‑TTTTTTTCATGTATATCCTCCTTTACTTTTAAAAAGGTGATTTATAAAATAAGTtata  <  1:1802074/100‑1 (MQ=255)                                                                           tAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTC‑TTTTTTCCATGTATATCCTCCTTTACTTTTAAAAAGGTTATTTATAAAATAAGTtata  <  1:2666435/100‑1 (MQ=255)                                                                           gatcTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTC‑TTTTTTTCATGTATATCCTCCTTTACTTTTAAAAAGGTTATTTATAAAATAAGTtata  <  1:4339728/97‑1 (MQ=255)                                                                           gatcTAGACGTGTTTTATAATCTAATGCTTTACGTTTTTTTC‑TTTTTTTCATGTATATCCTCCTGTTCGTTTAAAAAGGTTATTTATAAAATAAGTtata  <  1:371196/97‑1 (MQ=255)                                                                            aGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTCTTTTTTTTCATGTATATCCTACTTTACTTTTAAAAAGGTTATGTATAAAATAAGTtata  <  1:5265468/100‑1 (MQ=255)                                                                                              |                                                                                  AAGAACTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGTTTTTTTC‑TTTTTTTCATGTATATCCTCCTTTACTTTTAAAAAGGTTATTTATAAAATAAGTTATA  >  NC_007793/1631837‑1632010

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