Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,631,272	A→G	S8S (AGT→AGC)	SAUSA300_RS08060 ←	integrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,631,272	0	A	G	100.0%	128.9 / NA	51	S8S (AGT→AGC)	SAUSA300_RS08060	integrase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (0/51);  total (0/51)

GATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTCTTCAAATAGAGCAAACTCCTTACCTCTATCCATAGTTATACTTTTGGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTATTTAATCTACTAGCATATAAATATCTTGTC  >  NC_007793/1631175‑1631368                                                                                                  |                                                                                                 gATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGTTATGCt                                                                                                <  1:4567181/100‑1 (MQ=255) gATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGTTATGCt                                                                                                <  1:5382323/100‑1 (MQ=255)  atatGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGTTATGCtt                                                                                               <  1:2832346/100‑1 (MQ=255)  atatGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGTTATGCtt                                                                                               <  1:6331595/100‑1 (MQ=255)  atatGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGTTATGCtt                                                                                               <  1:5789476/100‑1 (MQ=255)   tatGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGTTATGCttt                                                                                              <  1:1663822/100‑1 (MQ=255)   tatGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGTTATGCttg                                                                                              <  1:4442829/100‑2 (MQ=255)    atGGACAACCCAAATCGCTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGTTATGCtttg                                                                                             <  1:994071/100‑2 (MQ=255)                 aTCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGTTATGCTTTTTGGTAATAATGGt                                                                                <  1:5810732/100‑1 (MQ=255)                   cACTAAAATAAATACTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTGCCTCTATCCATAGTTATGCTTTTTGGTAATAATGGttt                                                                              <  1:6484515/100‑1 (MQ=255)                         aataaataCTACAATTCAATTCATCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGTTATGCTTTTTGGTAATAATGGTTTCataat                                                                        <  1:2311478/100‑1 (MQ=255)                                           attcaTCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACa                                                      <  1:6299408/100‑1 (MQ=255)                                           attcaGCTTCTATTTGTGCAAATAGAGCGAACTCCTTTCCGCTACCCATATTTATGCTTTGTGGTAAGAATGGTTTCATAATATTTATAATCTCATTACa                                                      <  1:3063454/100‑1 (MQ=255)                                               aTCTTCTATTTGTTCAAATAGAGCGAACTCCTTACCTTTATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACt                                                  <  1:363823/100‑1 (MQ=255)                                                   tctATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACag                                              <  1:3705594/100‑2 (MQ=255)                                                    ctATTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACAtt                                             <  1:5732396/100‑1 (MQ=255)                                                      aTTTGTTCAAATAGAGCGAACTCCTTACCTCTATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTgg                                           <  1:3757010/100‑1 (MQ=255)                                                                aTAGAGCGAACTCCTTACCTCTATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACgtgt                                 <  1:3023433/100‑1 (MQ=255)                                                                aTAGAGCGAACTCCCTACCTCTATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACgtgg                                 <  1:4991132/100‑2 (MQ=255)                                                                   gagCGAACTCCTTACCTCTATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTAtt                              <  1:6476554/100‑1 (MQ=255)                                                                    tccgatcTCCTTGCCTCTATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTAttt                             <  1:5178461/94‑1 (MQ=255)                                                                            tccgatctctTTCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTATTTAATctact                     <  1:1712382/94‑1 (MQ=255)                                                                            tccgatctctATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAAGCTCATTACAAACTACATTGGATTTACGTGTATTTAATctact                     <  1:6240186/94‑1 (MQ=255)                                                                            cttccgatctATCCCTAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTATTTAATctact                     <  1:121991/93‑1 (MQ=255)                                                                            ccTTACCTCTATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTATTTAATctact                     <  1:1856023/100‑1 (MQ=255)                                                                            ccTTACCTCTATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTATTTAATctact                     <  1:6317888/100‑1 (MQ=255)                                                                              ttACCTCTATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTATTTAATctactgt                   <  1:2041683/100‑3 (MQ=255)                                                                              ttACCTCTATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTATTTAATCTACTAg                   <  1:6129942/100‑1 (MQ=255)                                                                               tACCTCTATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTATTTAATCTACTAGc                  <  1:1798159/100‑1 (MQ=255)                                                                               ccgatctGTCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTATTTAATCTACTAGc                  <  1:2510701/96‑1 (MQ=255)                                                                                cgatctATCCCTAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAAACTCATTACAGACTACATTGGATTTACGTGTATTTAATCTACTAGCa                 <  1:6525467/97‑1 (MQ=255)                                                                                cgatctATCCACAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTATTTAATCTATTAGCa                 <  1:5234330/97‑1 (MQ=255)                                                                                aCCTCTATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTATTTAATCTACTAGCa                 <  1:5340125/100‑1 (MQ=255)                                                                                aCCTCTATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTATTTAATCTACTAGCa                 <  1:4360567/100‑1 (MQ=255)                                                                                 tctctATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTATTTAATCTACTAGCat                <  1:3755751/99‑1 (MQ=255)                                                                                 gatctATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTATTTAATCTACTAGCat                <  1:5844628/98‑1 (MQ=255)                                                                                 gatctATCCACAATTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTATTTAATCTACTAGCat                <  1:5502988/98‑1 (MQ=255)                                                                                 cgatctttcctAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTATTTAATCTACTAGCat                <  1:1710451/90‑1 (MQ=255)                                                                                 cgatcttCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTATTTAATCTACTAGCat                <  1:4656884/94‑1 (MQ=255)                                                                                 ccTCTATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTATTTAATCTACTAGCat                <  1:6105518/100‑1 (MQ=255)                                                                                 ccTCTATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTATTTAATCTACTAGCat                <  1:3981388/100‑1 (MQ=255)                                                                                 ccTCTATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTATTTAACCTACTAGCat                <  1:3120551/100‑1 (MQ=255)                                                                                  ctctATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTATTTAATCTACTAGCata               <  1:5834488/100‑1 (MQ=255)                                                                                  ctctATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTATTTAATCTACTAGCata               <  1:5310432/100‑1 (MQ=255)                                                                                   tctATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTATTTAATCTACTAGCatca              <  1:6652753/100‑3 (MQ=255)                                                                                   tctATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTATTTAATCTACTAGCatat              <  1:3050920/100‑1 (MQ=255)                                                                                   tctATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTATTTAATCTACTAGCaaaa              <  1:6105601/100‑4 (MQ=255)                                                                                   tctATCCATAGTTACGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTATTTAATCTACTAGCatag              <  1:4912448/100‑2 (MQ=255)                                                                                   ctcctCCGTAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTATTTAATCTACTAGCatat              <  1:563100/96‑1 (MQ=255)                                                                                    ctATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTGTAATCTCATTACAAACTACATTGGATTTACGTGTATTTAATCTACTAGCAtata             <  1:6432582/100‑1 (MQ=255)                                                                                    ctATCCATAGTTATGCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTATTTAATCTACTAGCAtata             <  1:771994/100‑1 (MQ=255)                                                                                               tatgCTTTTTGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTATTTAATCTACTAGCATATAAATATCTTGtt  <  1:1554155/100‑2 (MQ=255)                                                                                                  |                                                                                                 GATATGGACAACCCAAATCACTAAAATAAATACTACAATTCAATTCATCTTCTATTTCTTCAAATAGAGCAAACTCCTTACCTCTATCCATAGTTATACTTTTGGGTAATAATGGTTTCATAATATTTATAATCTCATTACAAACTACATTGGATTTACGTGTATTTAATCTACTAGCATATAAATATCTTGTC  >  NC_007793/1631175‑1631368

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