Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,127,387	C→G	A259G (GCC→GGC)	SAUSA300_RS10835 →	integrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,127,387	0	C	G	100.0%	69.7 / NA	26	A259G (GCC→GGC)	SAUSA300_RS10835	integrase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (20/6);  total (20/6)

CCTAGATATATCAGTATCGCTAAATCTGATATGAAACACATTATGGACGTCATAAGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGCCGGATCTATAATTAACCTTCATGCTATCAATAATGCTTTGAAATCAGCCTGTAGAGTCAATAATATATCAATTATTACATCGCACGCATTAAG  >  NC_007793/2127290‑2127480                                                                                                  |                                                                                              ccTAGATATATCAGTATCGCTAAATCTGATATGAAACACATTATGGACGTCATAAGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGGCg                                                                                             >  1:5533561/1‑100 (MQ=255) ccTAGATATATCAGAATCGCTAAATCTGATATGAAACACATTATGGACGTCATAAGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGGCg                                                                                             >  1:4362962/1‑100 (MQ=255)  cTAGATATATCAGTATCGCTAAATCTGATATGAAACACATTATGGACGTCATAAGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGGCgg                                                                                            >  1:3322134/1‑100 (MQ=255)    aGATATATCAGTATCGCTAAATCTGATATAAAACACATTATGGACGTCATAAGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGGCGGAt                                                                                          >  1:1691543/1‑100 (MQ=255)                     tAAATCTGATATGAAACACATTATGGACGTCATAAGTACATTTGCAATTAGCTATGATGGTTACTTTTTCAAAGAATGCGGATCTATAATTAACCTTCAt                                                                         >  1:2508500/1‑100 (MQ=255)                                  caaCACATTATGGACGTCATAAGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGGCGGATCTATAATTAACCTTCATGCTATCAATAATg                                                            >  1:3692473/2‑100 (MQ=255)                                   tacacATTATGGACGTCATAAGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGGCGGATCTATAATTAACCTTCATGCTATCAATAATGc                                                           >  1:1656310/2‑100 (MQ=255)                                          gaTGGACGTCATAAGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGGCGGATCTATAATTAACCTTCATGCTATCAATAATGCTTTGaaa                                                    >  1:2263203/2‑100 (MQ=255)                                          atgggACGTCATAAGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGGCGGATCTATAATTAACCTTCATGCTATCAATAATGCTTTGaaa                                                    >  1:3892923/4‑100 (MQ=255)                                           aTGGACGTCATAAGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGGCGGATCTATAATTAACCTTCATGCTATCAATAATGCTTTGAAAt                                                   >  1:4068777/1‑100 (MQ=255)                                            tGGACGTCATAAGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGGCGGATCTATAATTAACCTTCATGCTATCAATAATGCTTTGAAATc                                                  <  1:4484410/100‑1 (MQ=255)                                                   cATAAGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGGCGGATCTATAATTAACCTTCATGCTATCAATAATGCTTTGAAATCAGCCTgg                                           >  1:3279219/1‑99 (MQ=255)                                                   cATAAGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGGCGGATCTATAATTAACCTTCATGCTATCAATAATGCTTTGAAATCAGCCTGt                                           >  1:1800360/1‑100 (MQ=255)                                                   cATAAGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGGCGGATCTATAATTAACCTTCATGCTATCAATAATGCTTTGAAATCAGCCTGt                                           <  1:2207911/100‑1 (MQ=255)                                                   cATAAGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGGCGGATCTATAATTAACCTTCATGCTATCAATAATGCTTTGAAATCAGCCTGt                                           <  1:2221825/100‑1 (MQ=255)                                                    aTAAGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGGCGGATCTATAATTAACCTTCATGCTATCAATAATGCTTTGAAATCAGCCTGTa                                          <  1:1686693/100‑1 (MQ=255)                                                     tAAGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGGCGGATCTATAATTAACCTTCATGCTATCAATAATGCTTTGAAATCAGCCTGTag                                         <  1:756150/100‑1 (MQ=255)                                                     tAAGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGGCGGATCTATAATTAACCTTCATGCTATCAATAATGCTTTGAAATCAGCCTGTag                                         <  1:3650639/100‑1 (MQ=255)                                                      aaGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGGCGGATCTATAATTAACCTTCATGCTATCAATAATGCTTTGAAATCAGCCTGTaga                                        >  1:3512401/1‑100 (MQ=255)                                                              ttGCAATTAGCTATGATGGTTACATTTTCAAAGAAGGCGGATCTATAATTAACCTTCCTGCTATCAATAATGCTTTGAAATCAGCCTGTAGAGTCagtaa                                >  1:1806139/1‑100 (MQ=255)                                                              ttGCAATTAGCTATGATGGTTACATTTTCAAAGAAGGCGGATCTATAATTAACCTTCATGCTATCAATAATGCTTTGAAATCAGCCTGTAGAGTCaataa                                >  1:545681/1‑100 (MQ=255)                                                                           tgatgGTTACATTTTCAAAGAAGGCGGATCTATAATTAACCTTCATGCTATCAATAATGCTTTGAAATCAGCCTGTAGAGTCAATAATATACCAattatt                   >  1:3242615/1‑100 (MQ=255)                                                                            gctgGTTACATTTTCAAAGAAGGCGGATCTATAATTAACCTTCATGCTATCAATAATGCTTTGAAATCAGCCTGTAGAGTCAATAATATACCAAttatta                  >  1:1410035/3‑100 (MQ=255)                                                                            gatgGTTACATTTTCAAAGAAGGCGGATCTATAATTAACCTTCATGCTATCAATAATGCTTTGAAATCAGCCTGTAGAGTCAATAATATACCAAttatta                  >  1:1393412/1‑100 (MQ=255)                                                                                  tACATTTTCAAAGAAGGCGGATCTATAATTAACCTTCATGCTATCAATAATGCTTTGAAATCAGCCTGTAGAGTCAATAATATACCAATTATTACATcgc            >  1:5360360/1‑100 (MQ=255)                                                                                            aagaagGCGGATCTATAATTAACCTTCATGCTATCAATAATGCTTTGAAATCAGCCTGTAGAGTCAATAATATACCAATTATTACATcgcacacagatcg  >  1:135394/1‑95 (MQ=255)                                                                                                  |                                                                                              CCTAGATATATCAGTATCGCTAAATCTGATATGAAACACATTATGGACGTCATAAGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGCCGGATCTATAATTAACCTTCATGCTATCAATAATGCTTTGAAATCAGCCTGTAGAGTCAATAATATATCAATTATTACATCGCACGCATTAAG  >  NC_007793/2127290‑2127480

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