Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,127,173	C→T	L188L (CTA→TTA)	SAUSA300_RS10835 →	integrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,127,173	0	C	T	100.0%	61.5 / NA	24	L188L (CTA→TTA)	SAUSA300_RS10835	integrase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/5);  total (19/5)

ACGAAGTCAGAACAAGAAAAATTCATCTCGGTAACTGAATATGAAATACTAAAAGATCACGTCAGAAAGAAGAGAACTCGTTCATCATTAGCGCTATTCATAATGATTTGTACGGGTTGTCGTGTCAGTGGTGCAAGAAATATAAAGATTGAGCATATCAACCAA  >  NC_007793/2127080‑2127244                                                                                              |                                                                        aCGAAGTCAGAACAAGAAAAATTCATCTCGGTAACTGAATATGAAATACTAAAAGATCACGTCAGAAAGCAGAGAACTCGTTCATCATTAGCGTTATTCa                                                                   >  1:1701081/1‑100 (MQ=255) aCGAAGTCAGAACAAGAAAAATTCATCTCGGTAACTGAATATGAAATACTAAAAGATCACGTCAGAAAGAAGAGAACTCGTTCATCATTAGCGTTATTCa                                                                   >  1:1773202/1‑100 (MQ=255)       ccAGAACAAGAAAAATTCATCGCGGTAACTGAATATGAAATACTACAAGATCACGTCAGAAAGAAGAGAACTCGTTCATCATTAGCGTTATTCATAATGa                                                             >  1:3598471/2‑100 (MQ=255)            aCAAGAAAAATTCATCTCGGTAACTGAATATGAAATACTAAAAGATCACGTCAGAAAGAAGAGAACTCGTTCATCATTAGCGTTATTCATAATGATTTGt                                                        <  1:1203077/100‑1 (MQ=255)                        cATCTCGGTAACTGAATATGAAATACTAAAAGATCACGTCAGAAAGAAGAGAACTCGTTCATCATTAGCGTTATTCATAATGATTTGTACGGGTTGTCgt                                            >  1:4422277/1‑100 (MQ=255)                             cGGTAACTGAATATGAAATACTAAAAGATCACGTCAGAAAGAAGAGAACTCGTTCATCATTAGCGTTATTCATAATGATTTGTACGGGTTGTCGTGTCAg                                       >  1:5736590/1‑100 (MQ=255)                             cGGTAACTGAATATGAAATACTAAAAGATCACGTCAGAAAGAAGAGAACTCGTTCATCATTAGCGTTATTCATAATGATTTGTACGGGTTGTCGTGTCAg                                       >  1:4400526/1‑100 (MQ=255)                             cGGTAACTGAATATGAAATACTAAAAGATCACGTCAGAAAGAAGAGAACTCGTTCATCATTAGCGTTATTCATAATGATTTGTACGGGTTGTCGTGTCAg                                       >  1:306713/1‑100 (MQ=255)                             cGGTAACTGAATATGAAATACTAAAAGATCACGTCAGAAAGAAGAGAACTCGTTCATCATTAGCGTTATTCATAATGATTTGTACGGGTTGTCGTGTCAg                                       >  1:2769665/1‑100 (MQ=255)                             cGGGAACTGAATATGAAATACTAAAAGATCACGTCAGAAAGAAGAGAACTCGTTCATCATTAGCGTTATTCATAATGATTTGTACGGGTTGTCGCGTCAg                                       >  1:4111814/1‑100 (MQ=255)                              ggTAACTGAATATGAAATACTAAAAGATCACGTCAGAAAGAAGAGAACTCGTTCATCATTAGCGTTATTCATAATGATTTGTACGGGTTGTCGTGTCAgt                                      >  1:2029648/1‑100 (MQ=255)                                    tGAATATGAAATACTAAAAGATCACGTCAGAAAGAAGAGAACTCGTTCATCATTAGCGTTATTCATAATGATTTGTACGGGTTGTCGTGTCAGTGGTGCa                                <  1:5233047/100‑1 (MQ=255)                                    tGAATATGAAATACTAAAAGATCACGTCAGAAAGAAGAGAACTCGTTCATCATTAGCGTTATTCATAATGAGTTGTACGGGTTGTCGTGTCAGTGGTGCa                                <  1:1532033/100‑1 (MQ=255)                                      aaTATGAAATACTAAAAGATCACGTCAGAAAGAAGAGAACTCGTTCATCATTAGCGTTATTCATAATGATTTGTACGGGTTGTCGTGTCAGTGGTGCAAg                              <  1:1976690/100‑1 (MQ=255)                                                       cctCACGTCAGAAAGAAGAGAACTCGTTCATCATTAGCGTTATTCATAATGATTTGTACGGGTTGTCGTGTCAGTGGTGCAAGAAATATAAAGACTGAGc             >  1:1391665/3‑100 (MQ=255)                                                        ctCACGTCAGAAAGAAGAGAACTCGTTCATCATTAGCGTTATTCATAATGATTTGTACGGGTTGTCGTGTCAGTGGTGCAAGAAATATAAAGACTGAGCa            >  1:6632157/2‑100 (MQ=255)                                                        aTCACGTCAGAAAGAAGAGAACTCGTTCATCATTAGCGTTATTCATAATGATTTGTACGGGTTGTCGTGTCAGTGGTGCAAGAAATATAAAGACTGAGCa            >  1:4360769/1‑100 (MQ=255)                                                          cACGTCAGAAAGAAGAGAACTCGTTCATCCTTAGCGTTATTCATAATGATTTGTACGGGTTGTCGTGTCAGTGGTGCAAGAAACATAAAGAtcggaagag          >  1:3446662/1‑92 (MQ=255)                                                          cACGTCAGAAAGAAGAGAACTCGTTCATCATTAGCGTTATTCATAATGATTTGTACGGGTTGTCGTGTCAGTGGTGCAAGAAATATAAAGAtcggaagag          >  1:4535392/1‑92 (MQ=255)                                                          cACGTCAGAAAGAAGAGAACTCGTTCATCATTAGCGTTATTCATAATGATTTGTACGGGTTGTCGTGTCAGTGGTGCAAGAAATATAAAGACTGAGCata          >  1:2818988/1‑100 (MQ=255)                                                          cACGTCAGAAAGAAGAGAACTCGTTCATCATTAGCGTTATTCATAATGATTTGTACGGGTTGTCGCGTCAGTGGTGCAAGAAATATAAAGACTGAGCata          >  1:2602620/1‑100 (MQ=255)                                                               cAGAAAGAAGAGAACTCGTTCATCATTAGCGTTATTCATAATGATTTGTACGGGTTGTCGTGTCAGTGGTGCAAGAAATATAAAGACTGAGCATATCAAc     >  1:6505154/1‑100 (MQ=255)                                                                 gaaagaaGAGAACTCGTTCATCATTAGCGTTATTCATAATGATTTGTACGGGTTGTCGTGTCAGTGGTGCAAGAAATATAAAGACTGAGCATATCAACCa   <  1:699753/100‑1 (MQ=255)                                                                  ctagaagaGAACTCGTTCATCATTAGCGTTATTCATAATGATTTGTACGGGTTGTCGTGTCAGTGGTGCAAGAAATATAAAGACTGAGCATATCAACCaa  >  1:2805157/3‑100 (MQ=255)                                                                                              |                                                                        ACGAAGTCAGAACAAGAAAAATTCATCTCGGTAACTGAATATGAAATACTAAAAGATCACGTCAGAAAGAAGAGAACTCGTTCATCATTAGCGCTATTCATAATGATTTGTACGGGTTGTCGTGTCAGTGGTGCAAGAAATATAAAGATTGAGCATATCAACCAA  >  NC_007793/2127080‑2127244

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