Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,127,409	T→G	H266Q (CAT→CAG)	SAUSA300_RS10835 →	integrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,127,409	0	T	G	100.0%	60.3 / NA	24	H266Q (CAT→CAG)	SAUSA300_RS10835	integrase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (16/8);  total (16/8)

GATATGAAACACATTATGGACGTCATAAGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGCCGGATCTATAATTAACCTTCATGCTATCAATAATGCTTTGAAATCAGCCTGTAGAGTCAATAATATATCAATTATTACATCGCACGCATTAAGA  >  NC_007793/2127317‑2127481                                                                                             |                                                                         gATATGAAACACATTATGGACGTCATAAGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGCCGGATCTATAATTAACCTTCAGGCTATCa                                                                   >  1:4200692/1‑100 (MQ=255)      gAAACACATTATGGACGTCATAAGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGCCGGATCTATAATTAACCTTCAGGCTATccggtac                                                              <  1:2645193/100‑7 (MQ=255)      gAAACACATTATGGACGTCATAAGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGCCGGATCTATAATTAACCTTCAGGCTATCaataat                                                              >  1:103380/1‑100 (MQ=255)        aaCACATTATGGACGTCATAAGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGCCGGATCTATAATTAACCTTCAGGCTATCAATAATGc                                                            >  1:4114064/1‑100 (MQ=255)           acaTTATGGACGTCATAAGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGCCGGATCTATAATTAACCTTCAGGCTATCAATAATGCttt                                                         <  1:382983/100‑1 (MQ=255)            caTTATGGACGTCATAAGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGCCAGATCTATAATTAACCTTCAGGCTATCAATAATGCTTTg                                                        <  1:5503826/100‑1 (MQ=255)               caTGGACGTCATAAGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGCCGGATCTATAATTAACCTTCAGGCTATCAATAATGCTTTGaaa                                                     >  1:5423852/2‑100 (MQ=255)                        cATAAGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGCCGGATCTATAATTAACCTTCAGGCTATCAATAATGCTTTGAAATCAGCCTGt                                            <  1:2838839/100‑1 (MQ=255)                         aTAAGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGCCGGATCTATAATTAACCTTCAGGCTATCAATAATGCTTTGAAATCAGCCTGtc                                           <  1:1736876/100‑2 (MQ=255)                          tAAGTACAGTTGCAAGTAGCTTTGATGGTTACATTTTCAAAGAAGCCGGATCTATAATTAACCTTCAGGCTATCAATAATGCTTTGAAATCAGCCTGtgg                                          <  1:5786864/100‑3 (MQ=255)                           aaGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGCCGGATCTATAATTAACCTTCAGGCTATCAATAATGCTTTGAAATCAGCCTGTaga                                         >  1:3506803/1‑100 (MQ=255)                           aaGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGCCGGATCTATAATTAACCTTCAGGCTATCAATAATGCTTTGAAATCAGCCTGTaga                                         >  1:248666/1‑100 (MQ=255)                               aCATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGCCGGATCTATAATTAACCTTCAGGCTATCAATAATGCTTTGAAATCAGCCTGTAGAGTCa                                     >  1:3330039/1‑100 (MQ=255)                                       aaTTAGCTCTGATGGTTACATTTTCAAAGAAGCCGGATCTCTAATTAACCTTCAGGCTATCAATAATGCTTTGAAATCAGCCTGTAGAGTCAATAatata                             >  1:4440949/1‑100 (MQ=255)                                        cctAGCTATGATGGTTACATTTTCAAAGAAGCCGGATCTATAATTAACCTTCAGGCTATCAATAATGCTTTGAAATCAGCCTGTAGAGTCAATAatatac                            >  1:408359/3‑99 (MQ=255)                                               atgatgGTTACATTTTCAAAGAAGCCGGATCTATAATTAACCTTCAGGCTATCAATAATGCTTTGAAATCAGCCTGTAGAGTCAATAATATACCAattat                     >  1:4040699/1‑100 (MQ=255)                                                       tACATTTTCAAAGAAGCCGGATCTATAATTAACCTTCAGGCTATCAATAATGCTTTGAAATCAGCCTGTAGAGTCAATAATATACCAATTATTACATcgc             <  1:167374/100‑1 (MQ=255)                                                               cAAAGAAGCCGGATCTATAATTAACCTTCAGGCTATCAATAATGCTTTGAAATCAGCCTGTAGAGTCAATAATATACCAATTATTACATcacagatcgga     >  1:362592/1‑93 (MQ=255)                                                               cAAAGAAGCCGGATCTATAATTAACCTTCAGGCTATCAATAATGCTTTGAAATCAGCCTGTAGAGTCAATAATATACCAATTATTACATCGCACGCATTa     >  1:6079364/1‑100 (MQ=255)                                                               cAAAGAAGCCGGATCTATAATTAACCTTCAGGCTATCAATAATGCTTTGAAATCAGCCTGTAGAGTCAATAATATACCAATTATTACATCGCACGCATTa     <  1:1103902/100‑1 (MQ=255)                                                               aaaaGAAGCCGGATCTATAATTAACCTTCAGGCTATCAATAATGCTTTGAAATCAGCCTGTAGAGTCAATAATATACCAATTATTACATcgcacgagatc     >  1:4657276/2‑95 (MQ=255)                                                                 aagaagCCGGATCTATAATTAACCTTCAGGCTATCAATAATGCTTTGAAATCAGCCTGTAGAGTCAATAATATACCAATTATTACGTcgcacacagatcg   >  1:4933369/1‑95 (MQ=255)                                                                  cgaagCCGGATCTATAATTAACCTTCAGGCTATCAATAATGCTTTGAAATCAGCCTGTAGAGTCAATAATATACCAATTATTACATCGCACGCATTAAGa  >  1:2461043/2‑100 (MQ=255)                                                                  agaagCCGGATCTATAATTAACCTTCAGGCTATCAATAATGCTTTGAAATCAGCCTGTAGAGTCAATAATATACCAATTATTACATcgcatgagatcgga  >  1:3663415/1‑90 (MQ=255)                                                                                             |                                                                         GATATGAAACACATTATGGACGTCATAAGTACATTTGCAATTAGCTATGATGGTTACATTTTCAAAGAAGCCGGATCTATAATTAACCTTCATGCTATCAATAATGCTTTGAAATCAGCCTGTAGAGTCAATAATATATCAATTATTACATCGCACGCATTAAGA  >  NC_007793/2127317‑2127481

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