Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,631,629	A→G	Y112H (TAT→CAT)	SAUSA300_RS08065 ←	integrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,631,629	0	A	G	100.0%	166.6 / NA	66	Y112H (TAT→CAT)	SAUSA300_RS08065	integrase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (1/65);  total (1/65)

ATTTCTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATACCCTGTATCTTTTTCAATTTGACTTGCTATCTGTCTAGGAGACCATTTTAATGAAAGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACG  >  NC_007793/1631536‑1631725                                                                                              |                                                                                                 tcttccgatcTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGt                                                                                            <  1:1771413/91‑1 (MQ=255) tcttCTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGt                                                                                            <  1:4447104/98‑1 (MQ=255) tcgtccgatcTGCTTGTCTCATATCACGCTTAACATCAACTTATACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGt                                                                                            <  1:5571565/91‑1 (MQ=255) cgatctcGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGt                                                                                            <  1:1786143/97‑1 (MQ=255) ccgatctGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGt                                                                                            <  1:578527/94‑1 (MQ=255) ccgatctGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGt                                                                                            <  1:2249098/94‑1 (MQ=255) atTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGt                                                                                            <  1:1019520/100‑1 (MQ=255) atTTTTTGCATGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGt                                                                                            <  1:3506831/100‑1 (MQ=255)  ttttttGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATGACTTATATAATGCCCTGTa                                                                                           <  1:5538863/100‑1 (MQ=255)  ttttttGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTa                                                                                           <  1:512858/100‑1 (MQ=255)  ttttttGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTa                                                                                           <  1:852201/100‑1 (MQ=255)  tttCTTGCCCGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTa                                                                                           <  1:3408458/100‑1 (MQ=255)  cgatctgtcTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTa                                                                                           <  1:2681495/92‑1 (MQ=255)  cgatctGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTa                                                                                           <  1:3722764/95‑1 (MQ=255)  cgatctGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTa                                                                                           <  1:6372360/95‑1 (MQ=255)     tttGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTAtgg                                                                                        <  1:3186465/99‑3 (MQ=255)     tttGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCt                                                                                        <  1:6229866/99‑1 (MQ=255)     tctGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCt                                                                                        <  1:1924867/98‑1 (MQ=255)      ttGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATcag                                                                                       <  1:5153848/100‑3 (MQ=255)      ttGCCCGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCtt                                                                                       >  1:2909955/1‑100 (MQ=255)       ttccTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTagccag                                                                                      <  1:4695056/98‑6 (MQ=255)       tGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGACCTGTATCttg                                                                                      <  1:5651962/100‑2 (MQ=255)       tGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCAGTGCGAATATATCTATATATAGTAGGATAACTTACATAATGCCCTGTATCtgg                                                                                      <  1:5687213/100‑3 (MQ=255)        gCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTaccccta                                                                                     <  1:6348466/100‑7 (MQ=255)          cgatctGTCTCGTATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTc                                                                                   <  1:2623299/95‑1 (MQ=255)           tGCTTGTCTCATGTCACGCTGATCATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACATATATAATGCCCTGTATCTTTTTCa                                                                                  <  1:5244845/100‑1 (MQ=255)           tGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCa                                                                                  <  1:4668344/100‑1 (MQ=255)                 tctcATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTgc                                                                            <  1:6059699/100‑2 (MQ=255)                  ctcATATCACGCTTAACATCAACTTTTTCCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTggg                                                                           <  1:4303244/100‑3 (MQ=255)                      tatCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGc                                                                       <  1:1379138/100‑1 (MQ=255)                       atCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGCt                                                                      <  1:3896686/100‑1 (MQ=255)                        tCACACTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGCTa                                                                     <  1:5595159/100‑1 (MQ=255)                         taCGCCTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGCTAt                                                                    <  1:3190314/99‑1 (MQ=255)                          aCGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGCTAta                                                                   <  1:134842/100‑2 (MQ=255)                             cTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGCTACCTGt                                                                <  1:5330460/100‑1 (MQ=255)                              ttAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGCTATCTGTc                                                               <  1:3512781/100‑1 (MQ=255)                               tgaCATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGCTATCTGTCt                                                              <  1:5167955/98‑1 (MQ=255)                               tAACATCATCTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGCTATCTGTCt                                                              <  1:5448519/100‑1 (MQ=255)                               tAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGCTATCTGTCt                                                              <  1:3655860/100‑1 (MQ=255)                               tAACATCAACTTTTACCAAGCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGATATCTGTCt                                                              <  1:4720098/100‑1 (MQ=255)                                  cATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGCTATCTGTCTAgg                                                           <  1:919095/100‑1 (MQ=255)                                  cATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGCTATCTGTCTAgg                                                           <  1:5925976/100‑1 (MQ=255)                                  cATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGCTATCTGTCTAgg                                                           <  1:6124628/100‑1 (MQ=255)                                  cATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGCTATCTGTCTAgg                                                           <  1:1316131/100‑1 (MQ=255)                                  cATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGCTATCTGTCTAgg                                                           <  1:5499074/100‑1 (MQ=255)                                  cATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGCTATCTGTCTAgg                                                           <  1:6438223/100‑1 (MQ=255)                                  cATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGCTATCTGTCTAgg                                                           <  1:5080075/100‑1 (MQ=255)                                  cATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGCTATCTGTCTAgg                                                           <  1:3483845/100‑1 (MQ=255)                                  cATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGCTATCTGTCTAgg                                                           <  1:4470524/100‑1 (MQ=255)                                  cATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGCTATCTGTCTAgg                                                           <  1:4211100/100‑1 (MQ=255)                                  cATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGCTATCTGTCTAgg                                                           <  1:3925040/100‑1 (MQ=255)                                  cATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGCTATCTGTCTAgg                                                           <  1:3524777/100‑1 (MQ=255)                                  cATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGCTATCTGTCTAgg                                                           <  1:5352252/100‑1 (MQ=255)                                  cATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGCTATCTGTCTAgg                                                           <  1:3639804/100‑1 (MQ=255)                                  cATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGCTATCTGTCTAgg                                                           <  1:3839286/100‑1 (MQ=255)                                  cATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGCTATCTGTCTAgg                                                           <  1:3829098/100‑1 (MQ=255)                                  cATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGCTATCTGTCTAgg                                                           <  1:3719068/100‑1 (MQ=255)                                  cATCAACGTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGCTATCTGTCTAgg                                                           <  1:5624215/100‑1 (MQ=255)                                    tCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGCTATCTGTCTAGgag                                                         <  1:5828163/100‑1 (MQ=255)                                      aaCTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATC‑TTTTCAATTTGACTTGCTATCTGTCTAGGAGAcc                                                      <  1:3893859/100‑1 (MQ=255)                                       aCTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGCTATCTGTCTAGGAGAcc                                                      <  1:4029900/100‑1 (MQ=255)                                                 aTCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGCTATCTGTCTAGGAGACCATTTTAATGa                                            <  1:3896336/100‑1 (MQ=255)                                                                        ataGTAGGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGCTATCTGTCTAGGAGACCATTTTAATGAAAGTTTTTCTTCAATATATGCtt                     <  1:2676711/100‑1 (MQ=255)                                                                          agtagGATAACTTATATAATGCCCTGTATCTTTTTCAATTTGACTTGCTATCTGTCTAGGAGACCATTTTAATGAAAGTTTTTCTTCAATATATGCtttg                   <  1:933759/100‑2 (MQ=255)                                                                                       atataATGCCCTGTATCTTTTTCAATTTGACTTGCTATCTGTCTAGGAGACCATTTTAATGAAAGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCa      <  1:6109290/100‑1 (MQ=255)                                                                                       atataATGCCCTGTATCTTTTTCAATTTGACTTGCTATCTGTCTAGGAGACCATTTTAATGAAAGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCa      <  1:4646261/100‑1 (MQ=255)                                                                                           aatgCCCTGTATCTTTTTCAATTTGACTTGCTATCTGTCTAGGAGACCATTTTAATGAAAGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACg  >  1:5826687/1‑100 (MQ=255)                                                                                              |                                                                                                 ATTTCTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATACCCTGTATCTTTTTCAATTTGACTTGCTATCTGTCTAGGAGACCATTTTAATGAAAGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACG  >  NC_007793/1631536‑1631725

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: ‑