Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,819,073	T→C	D116G (GAT→GGT)	SAUSA300_RS14400 ←	flavin reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,819,073	0	T	C	100.0%	71.6 / NA	28	D116G (GAT→GGT)	SAUSA300_RS14400	flavin reductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (1/27);  total (1/27)

TGGAATACGAATCGTTTGAACGTACTCACATTCAAAATGAACTGGCGATTCTTTTACACGATATCCTGGAGCTTCTATACATTTTTCCTTTGTTACATCTGCAAAATTAAATTCATCCTCTTCTGGTGGCAATGCTTTCGATGATAAATTAACTGCTTCTCTTAAATCATACGTTGCCATATTCCACAC  >  NC_007793/2818976‑2819164                                                                                                  |                                                                                            tGGAATACGAATCGTTTGAACGTACTCACATTCAAAATGAACTGGCGATTCTTTTACACGATATCCTGGAGCTTCTATACATTTTTCCTTTGTTACACCt                                                                                           <  1:4690792/100‑1 (MQ=255)  ggAATACGAATCGTTTGAACGTACTCACATTCAAAATGAACTGGCGATTCTTTTACACGATATCCTGGGGCTTCTATACATTTTTCCTTTGTTACACCTg                                                                                          <  1:827726/100‑1 (MQ=255)  ggAATACGAATCGTTTGAACGTACTCACATTCAAAATGAACTGGCGATGCTTTTACACGATATCCTGGAGCTTCTATACATTTTTCCTTTGTTACACCTg                                                                                          <  1:3632476/100‑1 (MQ=255)   gAATACGAATCGTTTGAACGTACTCACATTCAAAATGAACTGGCGATTCTTTTACACGATATCCTGGAGCTTCTATACATTTTTCCTTTGTTACACCTGc                                                                                         <  1:5750416/100‑1 (MQ=255)   gAATACGAATCGTTTGAACGTACTCACATTCAAAATGAACTGGCGATTCTTTTACACGATATCCTGGAGCTTCTATACATTTTTCCTTTGTTACACCTGc                                                                                         <  1:845492/100‑1 (MQ=255)   gAATACGAATCGTTTGAACGTACTCACATTCAAAATGAACTGGCGATTCTTTTACACGATATCCTGGAGCTTCTATACATTTTTCCTTTGTTACACCTGc                                                                                         <  1:1681901/100‑1 (MQ=255)    aaTACGAATCGTTTGAACGTACTCACATTCAAAATGAGCTGGCGATCCTTTTACACGATATCCTGGAGCTTCTATACATTTTTCCTTTGTTACACCTGct                                                                                        <  1:3029811/100‑2 (MQ=255)    aaTACGAATCGTTTGAACGTACTCACATCCAAAATGAACTGGCGATTCTTTTACACGATATCCTGGAGCTTCTATACATTTTTCCTTTGTTACACCTGCa                                                                                        <  1:6266989/100‑1 (MQ=255)    aaTACGAAGCGTTTGAACGTACTCACATTCAAAATGAACTGGCGATTCTTTTACTCGATATGCTGGAGCTTCTATACATTTTTCCTTTGTTACACCTGCa                                                                                        <  1:6422314/100‑1 (MQ=255)     aTACGAATCGTTTGAACGTACTCACATTCAAAATGAACTGGCGATTCTTTTACACGATATCCTGGAGCTTCTATACATTTTTCCTTTGTTACACCTGCac                                                                                       <  1:2251321/100‑2 (MQ=255)                           cacaTTCAAAATGAACTGGCGATTCTTTTACACGATATCCTGGAGCTTCTATACATTTTTCCTTTGTTACACCTGCAAAATTAAATTCATCCTCTTCtgg                                                                 <  1:5269555/100‑1 (MQ=255)                               ttCAAAATGAACTGGCGATTCTTTTACACGATATCCTGGAGCTTCTATACATTTTTCCTTTGTTACTCCTGCAAAATTAAATTCATCCTCTTCTGGTGGc                                                             <  1:2465520/100‑1 (MQ=255)                               ttCAAAATGAACTGGCGATTCTTTTACACGATATCCTGGAGCTTCTATACATTTTTCCTTTGTTACACCTGCAAAATTAAATTCATCCTCTTCtggggca                                                             <  1:4194802/100‑5 (MQ=255)                                tCAAAATGAACTGGCGATTCTTTTACACGATATCCTGGAGCTTCTATACATTTTTCCTTTGTTACACCTGCAAAATTAAATTCATCCTCTTCTGGTGGCa                                                            <  1:2837187/100‑1 (MQ=255)                                tCAAAATGAACTGGCGATTCTTTTACACGATATCCTGGAGCTTCTATACATTTTTCCTTTGTTACACCTGCAAAATTAAATTCATCCTCTTCTGGTGGCa                                                            <  1:5574580/100‑1 (MQ=255)                                tCAAAATGAACTGGCGATTCTTTTACACGATATCCTGGAGCTTCTATACATTTTTCCTTTGTTACACCTGCAAAATTAAATTCATCCTCTTCTGGTGGCa                                                            <  1:6607744/100‑1 (MQ=255)                                tCAAAATGAACTGGCGATTCTTTTACACGATATCCTGGAGCTTCTATACATTTTTCCTTTGTTACACCTGCAAAATTAAATTCATCCTCTGCTGGTGGCa                                                            <  1:3732319/100‑1 (MQ=255)                                tCAAAATGAACTGGCGATTCTTTTACACGATATCCTGGAGCTTCTATAAATTTTTCCTTTGTTACACCTGCAATATTAAATTCATCCTCTTCTGGTGGcc                                                            <  1:3907998/100‑2 (MQ=255)                                 cAAAATGAACTGGCGATTCTTTTACACGATATCCTGGAGCTTCTATACATTTTTCCTTTGTTACACCTGCAAAATTAAATTCATCCTCTTCTGGTGGCaa                                                           <  1:996965/100‑1 (MQ=255)                                 cAAAATGAACTGGCGATTCTTTTACACGATATCCTGGAGCTTCTATACATTATTCCTTTGTTACACCTGCAAAATTAAATTCATCCTCTTCTGGTGGCaa                                                           <  1:5916501/100‑1 (MQ=255)                                  cttccgatcTGGCGGTTCTTTTACACGATATCCTGGAGCTTCTATACATTTTTCCTTTGTTACACCTGCAAAATTAAATTCATCCTCTTCTGGTGGCaag                                                          <  1:3450609/92‑2 (MQ=255)                                  aaaaTGAACTGGCGATTCTTTTACACGATATCCTGGAGCTTCTATACATTTTTCTTTTGTTACACCTGCAAAATTAAATTCATCCTCTTCTGGTGGCAAt                                                          <  1:323550/100‑1 (MQ=255)                                                      ttACACGATATCCTGGAGCTTCTATACATTTTTCCTTTGTTACACCTGCAAAATTAAATTCATCCTCTTCTGGTGGCAATGCTTTCGATGATAAATTAAc                                      >  1:2335952/1‑100 (MQ=255)                                                          cccataaCCTGGAGCCTCCAACCAATTTTCCTTTGTTACACCTGCAAAATTAAATTCATCCTCTTCTGGTGGCAATGCTTTCGATGATAAATTAACTGCt                                  <  1:4326801/97‑1 (MQ=255)                                                                         ttCTATACATTTTTCCTTTGTTACACCTGCAAAATTAAATTCATCCTCTTCTGGTGGCAATGCTTTCGATGATAAATTAACTGCTTCTCTTAAATCATAc                   <  1:6586030/100‑1 (MQ=255)                                                                               tccgatctCCTTTGTTACACCTGCAAAATTAAATTCATCCTCTTCTGGTGGCAATGCTTTCGATGATAAATTAACTGCTTCTCTTAAATCATACGTTGcc             <  1:5403048/93‑1 (MQ=255)                                                                               aCATTTTTCCTTTGTTACACCTGCAAAATTAAATTCATCCTCTTCTGGTGGCAATGCTTTCGATGATAAATTAACTGCTTCTCTTAAATCATACGTTGcc             <  1:100666/100‑1 (MQ=255)                                                                                          ttGTTACACCTGCAAAATTAAATTCATCCTCTTCTGGTGGCAATGCTTTCGATGATAAATTAACTGCTTCTCTTAAATCATACGTTGCCATATTCcacac  <  1:1215398/100‑1 (MQ=255)                                                                                                  |                                                                                            TGGAATACGAATCGTTTGAACGTACTCACATTCAAAATGAACTGGCGATTCTTTTACACGATATCCTGGAGCTTCTATACATTTTTCCTTTGTTACATCTGCAAAATTAAATTCATCCTCTTCTGGTGGCAATGCTTTCGATGATAAATTAACTGCTTCTCTTAAATCATACGTTGCCATATTCCACAC  >  NC_007793/2818976‑2819164

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