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%	64.2 / NA	23	D11D (GAC→GAT)	SAUSA300_RS08065	integrase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (1/22);  total (1/22)

CTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGTTTTTTTCTTTTTTTCATGTATATCCTCCTTTACTTTTAAAAAGGTTA  >  NC_007793/1631842‑1631992                                                                                         |                                                               ccgCATCATCACGACCTAAAATCTTAGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCttt                                                     <  1:5287472/98‑1 (MQ=255)             cgatctatgtCTTAGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGtttttttca                                         <  1:5501888/91‑2 (MQ=255)                    aaTCTTAGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACG‑TTTTTTCTTTTTTTCa                                 <  1:3039888/100‑1 (MQ=255)                        ttAGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTCTTTTTTTCATGt                              <  1:816265/100‑1 (MQ=255)                        ttAGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTCTTTTTTTCATGt                              <  1:5501407/100‑1 (MQ=255)                        ttAGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTCTTTTTTTCATGt                              <  1:3933637/100‑1 (MQ=255)                        ttAGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTCTTTTTTTCATGt                              <  1:3998311/100‑1 (MQ=255)                         tAGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTCTTTTTTTCATGta                             >  1:176850/1‑100 (MQ=255)                          aGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTCTTTTTTTCATGtat                            <  1:5073321/100‑1 (MQ=255)                          aGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTCTTTTTTTCATGtat                            <  1:267192/100‑1 (MQ=255)                           gAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTCTTTTTTTCATGtata                           <  1:4393267/100‑1 (MQ=255)                              aTCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTCTTTTTTTCATGTATAtcc                        <  1:2937585/100‑1 (MQ=255)                               tCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTCTTTTTTTCATGTATAtcct                       <  1:1076392/100‑1 (MQ=255)                                 tGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTCTTTTTTTCATGTATAtccaca                     <  1:3427980/100‑4 (MQ=255)                                   aaTGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTCTTTTTTTCATGTATATCCTCCtt                   <  1:5156565/100‑1 (MQ=255)                                    atgatgTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTCTTTTTTTCATGTATATCCTCCttt                  <  1:2990362/100‑1 (MQ=255)                                    atgatgTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTCTTTTTTTCATGTATATCCTCCttt                  <  1:2298754/100‑1 (MQ=255)                                       atgTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTCTTTTTTTCATGTATATCCTCCTTTACt               <  1:4235613/100‑1 (MQ=255)                                        tgTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTCTTTTTTTCATGTATATCCTCCTTTACtg              <  1:746496/100‑2 (MQ=255)                                        tgTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACG‑TTTTTTCTTTTTTTCATGTATATCCTCCTTTACttt             <  1:2241418/100‑1 (MQ=255)                                          tttAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTCTTTTTTTCATGTATATCCTCCTTTACtttt            <  1:107996/100‑1 (MQ=255)                                            tAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTCTTTTTTTCATGTATATCCTCCTTTACTTTTat          <  1:1842336/100‑2 (MQ=255)                                                    catTTAAATACTTTGCGATAGCTAGACGTGTTTTATAATCTAATTCTTTACGTTTTTTTCTTTTTTTCATGTATATCCTCCTTTACTTTTAAAAAGGttt  <  1:1515370/100‑2 (MQ=255)                                                                                         |                                                               CTGCATCATCACGACCTAAAATCTTAGAAATCTGAATGATGTTTAAACCATCATTTAAATACTTTGCGATAGCTAGACGTGTTTTATAGTCTAATTCTTTACGTTTTTTTCTTTTTTTCATGTATATCCTCCTTTACTTTTAAAAAGGTTA  >  NC_007793/1631842‑1631992

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