Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,127,574	C→T	S321S (TCC→TCT)	SAUSA300_RS10835 →	integrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,127,574	0	C	T	100.0%	55.1 / NA	20	S321S (TCC→TCT)	SAUSA300_RS10835	integrase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/2);  total (18/2)

AGACACACTCATTGTTCTTATTTACTAGCAAAAGGTGTATCTATACATTACATTTCTAAAAGATTAGGTCATAAAAATATAGCAATAACTACATCCGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGA  >  NC_007793/2127479‑2127668                                                                                                |                                                                                               aGACACACTCATTGTTCTTATTTACTAGCAAAAGGTGTATCTATACATTACATTTCTAAAAGATTAGGTCATAAAAATATAGCAATAACTACATCTgtgt                                                                                            <  1:4178676/100‑1 (MQ=255)          cATTGTTCTTATTTACTAGCAAAAGGTGTATCTATACATTACATTTCTAAAAGATTAGGTCATAAAAATATAGCAATAACTACATCTGTGTATTCTCAtt                                                                                   >  1:4569120/1‑100 (MQ=255)              gTTCTTATTTACTAGCAAAAGGTGTATCTATACATTACATTTCTAAAAGATTAGGTCATAAAAATATAGCAATAACTACATCTGTGTATTCTCATTTgat                                                                               >  1:152405/1‑98 (MQ=255)                       tACTAGCAAAAGGTGTATCTATACATTACATTTCTAAAAGATTAGGTCATAAAAATATAGCAATAACTACATCTGTGTATTCTCATTTGTTAGAAGaaaa                                                                      >  1:2331241/1‑100 (MQ=255)                       tACTAGCAAAAGGTGTATCTATACATTACATTTCTAAAAGATTAGGTCATAAAAATATAGCAATAACTACATCTGTGTATTCTCATTTGTTAGAAGaaaa                                                                      >  1:4004592/1‑100 (MQ=255)                       tACTAGCAAAAGGTGTATCTATACATTACATTTCTAAAAGATTAGGTCATAAAAATATAGCAATAACTACATCTGTGTATTCTCATTTGTTAGAAGaaaa                                                                      >  1:4147525/1‑100 (MQ=255)                       tACTAGCAAAAGATGTATCTATACATTACATTTCTAAAAGATTAGGTCATAAAAATATAGCAATAACTACATCTGTGTATTCTCATTTGTTAGAAGaaaa                                                                      >  1:3560636/1‑100 (MQ=255)                        aCTAGCAAAAGGTGTATCTATACATTACATTTCTAAAAGATTAGGTCATAAAAATATAGCAATAACTACATCTGTGTATTCTCATTTGTTAGAAGaaaaa                                                                     >  1:2287476/1‑100 (MQ=255)                          aaGCAAAAGGTGTATCTATACATTACATTTCTAAAAGATTAGGTCATAAAAATATAGCAATAACTACATCTGTGTATTCTCATTTGTTAGAAGAAAAAtt                                                                   >  1:3703768/2‑100 (MQ=255)                           aGCAAAAGGTGTATCTATACATTACATTTCTAAAAGATTAGGTCATAAAAATATAGCAATAACTACATCTGTGTATTCTCATTTGTTagaagatcggaag                                                                  >  1:2853535/1‑93 (MQ=255)                           aGCAAAAGGTGTATCTATACATTACATTTCTAAAAGATTAGGTCATAAAAATATAGCAATAACTACATCTGTGTATTCTCATTTGTTAGAAGAAAAAttt                                                                  >  1:2046497/1‑100 (MQ=255)                           aGCAAAAGGTGTATCTATACATTACATTTCTAAAAGATTAGGTCATAAAAATATAGCAATAACTACATCTGTGTATTCTCATTTGTTAGAAGAAAAAttt                                                                  >  1:4328097/1‑100 (MQ=255)                           aGCAAAAGGTGTATCTATACATTACATTTCTAAAAGATTAGGTCATAAAAATATAGCAATAACTACATCTGTGTATTCTCATTTGTTAGAAGAAAAAttc                                                                  >  1:593749/1‑99 (MQ=255)                                           ataCATTACATTTCTAAAAGATTAGGTCATAAAAATATAGCAATAACTACATCTGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACaaaa                                                  >  1:4559185/1‑100 (MQ=255)                                           ataCATTACATTTCTAAAAGATTAGGTCATAAAAATATAGCAATAACTACATCTGTGTATTCTCATTTGTTAGAAGAAAAATTTAACGAagagatcggaa                                                  >  1:1309261/1‑93 (MQ=255)                                                                                 aGCAATAACTACATCTGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGAc            >  1:4643272/1‑100 (MQ=255)                                                                                    aaTAACTACATCTGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCa         >  1:2683774/1‑100 (MQ=255)                                                                                      tcttccgaTCTGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCAtt       <  1:1611322/93‑1 (MQ=255)                                                                                       aaCTACATCTGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAACATTTTAGAAAGTATGTAATTTAGGGACCCATTa      >  1:3819596/1‑100 (MQ=255)                                                                                           aCATCTGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGa  >  1:4204568/1‑100 (MQ=255)                                                                                                |                                                                                               AGACACACTCATTGTTCTTATTTACTAGCAAAAGGTGTATCTATACATTACATTTCTAAAAGATTAGGTCATAAAAATATAGCAATAACTACATCCGTGTATTCTCATTTGTTAGAAGAAAAATTTAATGAAGAGGACAAAAAAACAACTAAAATTTTAGAAAGTATGTAATTTAGGGACCCATTAGGGA  >  NC_007793/2127479‑2127668

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