Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	996,519	T→G	H248Q (CAT→CAG)	SAUSA300_RS04890 →	RNA pseudouridine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	996,519	0	T	G	100.0%	70.7 / NA	26	H248Q (CAT→CAG)	SAUSA300_RS04890	RNA pseudouridine synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (25/1);  total (25/1)

CGGGACGTACACATCAAATTCGTGTACATTTTCAACATATTGGGCATCCAATTGTGGGAGATTCTTTGTATGATGGTTTTCATGACAAAATTCATGGTCAAGTACTGCAATGTACGCAAATATATTTTGTTCATCCAATCAATAAGAACAATATTTATATTACAATTGATTATAAGCAATTACTTAAATTA  >  NC_007793/996425‑996615                                                                                               |                                                                                                 cGGGACGTACACATCAAATTCGTGTACATTTTCAACATATTGGGCATCCAATTGTGGGAGATTCTTTGTATGATGGTTTTCATGACAAAATTCAGGGTCa                                                                                             >  1:950698/1‑100 (MQ=255)                        gtACATTTTCAACATATTGGGCATCCAATTGTGGGAGATTCTTTGTATGATGGTTTTCATGACAAAATTCAGGGTCAAGTACTGCAATGTACGCAAatat                                                                      >  1:4004790/1‑100 (MQ=255)                                tCAACATATTGGGCATCCAATTGTGGGAGATTCTTTGTATGATGGTTTTCATGACAAAATTCAGGGTCAAGTACTGCAATGTACGCAAATATATTTTGtt                                                              >  1:3779706/1‑100 (MQ=255)                                gcAACATATTGGGCATCCAATTGTGGGAGATTCTTTGTATGATGGTTTTCATGACAAAATTCAGGGTCAAGTACTGCAATGTACGCAAATATATTTTGtt                                                              >  1:3267458/2‑100 (MQ=255)                                   atctcgtGGGCATCCAATTGTGGGAGATTCTTTGTATGATGGTTTTCATGACAAAATTCAGGGTCAAGTACTGCAATGTACGCAAATATATTTTGTTCAt                                                           >  1:2626192/7‑100 (MQ=255)                                        ctGGGCATCCAATTGTGGGAGATTCTTTGTATGATGGTTTTCATGACAAAATTCAGGGTCAAGTACTGCAATGTACGCAAATATATTTTGTTCATCcaat                                                      >  1:3570995/2‑100 (MQ=255)                                         ttcgCATCCAATTGTGGGAGATTCTTTGTATGATGGTTTTCATGACAAAATTCAGGGTCAAGTACTGCAATGTACGCAAATATATTTTGTTCATCcaatc                                                     >  1:1006796/4‑100 (MQ=255)                                         tGGGCATCCAATTGTGGGAGATTCTTTGTATGATGGTTTTCATGACAAAATTCAGGGTCAAGTACTGCAATGTACGCAAATATATTTTGTCCATCcaatc                                                     >  1:520161/1‑100 (MQ=255)                                          gggCATCCAATTGTGGGAGATTCTTTGTATGATGGTTTTCATGACAAAATTCAGGGTCAAGTACTGCAATGTACGCAAATATATTTTGTTCATCcaatca                                                    >  1:2656494/1‑100 (MQ=255)                                             ttgtgcaTTGTGGGAGATTCTTTGTATGATGGTTTTCATGACAAAATTCAGGGTCAAGTACTGCAATGTACGCAAATATATTTTGTTCATCCAATCAATa                                                 >  1:5472203/7‑100 (MQ=255)                                               tCCAATTGTGGGAGATTCTTTGTATGATGGTTTTCATGACAAAATTCAGGGTCAAGTACTGCAATGTACGCAAATATATTTTGTTCATCCAATCAATAAg                                               >  1:1037880/1‑100 (MQ=255)                                                ccAATTGTGGGAGATTCTTTGTATGATGGTTTTCATGACAAAATTCAGGGTCAAGTACTGCAATGTACGCAAATATATTTTGTTCATCCAATCAATAAGa                                              >  1:5371619/1‑100 (MQ=255)                                                    ctgtgGGAGATTCTTTGTATGATGGTTTTCATGACAAAATTCAGGGTCAAGTACTGCAATGTACGCAAATATATTTTGTTCATCCAATCAATAAGAACaa                                          >  1:584910/2‑100 (MQ=255)                                                     tgtgGGAGATTCTTTGTATGATGGTTTTCATGACAAAATTCAGGGTCAAGTACTGCAATGTACGCAAATATATTTTGTTCATCCAATCAATAAGAACAat                                         >  1:1853424/1‑100 (MQ=255)                                                     tgtgGGAGATTCTTTGTATGATGGTTTTCATGACAAAATTCAGGGTCAAGTACTGCAATGTACGCAAATATATTTTGTTCATCCAATCAATAAGAACAat                                         >  1:1459083/1‑100 (MQ=255)                                                     tgtgGGAGATTCTTTGTATGATGGTTTTCATGACAAAATTCAGGGTCAAGTACTGCAATGTACGCAAATATATTTTGTTCATCCAATCAATAAGAACAat                                         >  1:3960786/1‑100 (MQ=255)                                                     tgtgGGAGATTCTTTGTATGATGGTTTTCATGACAAAATTCAGGGTCAAGTACTGCAATGTACGCAAATATATTTTGTTCATCCAATCAATAAGAACAat                                         >  1:1156466/1‑100 (MQ=255)                                                      ttgGGAGATTCTTTGTATGATGGTTTTCATGACAAAATTCAGGGTCAAGTACTGCAATGTACGCAAATATATTTTGTTCATCCAATCAATAAGAACAata                                        >  1:1224531/2‑100 (MQ=255)                                                      gtgGGAGATTCTTTGTATGATGGTTTTCATGACAAAATTCAGGGTCAAGTACTGCAATGTACGCAAATATATTTTGTTCATCCAATCAATAAGAACAata                                        >  1:1241356/1‑100 (MQ=255)                                                                   tGTATGATGGTTTTCATGACAAAATTCAGGGTCAAGTACTGCAATGTACGCAAATATATTTTGTTCATCCAATCAATAAGAACAATATTTATATTACAAt                           >  1:1676133/1‑100 (MQ=255)                                                                              tttCATGACAAAATTCAGGGTCAAGTACTGCAATGTACGCAAATATATTTTGTTCATCCAATCAATAAGAACAATATTTATATTACAATTGATTATAAGc                >  1:3834124/1‑100 (MQ=255)                                                                                tCATGACAAAATTCAGGGTCAAGTACTGCAATGTACGCAAATATATTTTGTTCATCCAATCAATAAGAACAATATTTATATTACAATTGATTATAAGCag              <  1:2715670/100‑2 (MQ=255)                                                                                caaTGACAAAATTCAGGGTCAAGTACTGCAATGTACGCAAATATATTTTGTTCATCCAATCAATAAGAACAATATTTATATTACAATTGATTATAAGCaa              >  1:1404646/3‑100 (MQ=255)                                                                                    gACAAAATTCAGGGTCAAGTACTGCAATGTACGCAAATATATTTTGTTCATCCAATCAATAAGAACAATATTTATATTACAATTGATTATAAGCAATTAc          >  1:309642/1‑100 (MQ=255)                                                                                            tCAGGGTCAAGTACTGCAATGTACGCAAATATATTTTGTTCATCCAATCAATAAGAACAATATTTATATTACAATTGATTATAAGCAATTACTTAAatta  >  1:4750727/1‑100 (MQ=255)                                                                                            tCAGGGTCAAGTACTGCAATGTACGCAAATATATTTTGTTCATCCAATCAATAAGAACAATATTTATATTACAATTGATTATAAGCAATTACTTAAatta  >  1:1233091/1‑100 (MQ=255)                                                                                               |                                                                                                 CGGGACGTACACATCAAATTCGTGTACATTTTCAACATATTGGGCATCCAATTGTGGGAGATTCTTTGTATGATGGTTTTCATGACAAAATTCATGGTCAAGTACTGCAATGTACGCAAATATATTTTGTTCATCCAATCAATAAGAACAATATTTATATTACAATTGATTATAAGCAATTACTTAAATTA  >  NC_007793/996425‑996615

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