Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,578,590	G→A	D12D (GAC→GAT)	SAUSA300_RS07710 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,578,590	0	G	A	100.0%	43.2 / NA	22	D12D (GAC→GAT)	SAUSA300_RS07710	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (1/21);  total (1/21)

GATAACAAGTAGCAATGCGAATATACCGAATAATGACAGTATTCCAACAATAAAGATGTCGAACCCATCCATATTTTTAAGTTTTTTAATCATCATTGTCATCTCCAGTATCAATTAAACTAGGCATCATTCTTAACATAGCC  >  NC_007793/1578533‑1578675                                                          |                                                                                      gacaaCAAGTAACAATGCGAATATGCCGAATAATGACAGTTCTCCAACAATAAAGATATCGAACCAATCCATATTTTTAAGTTTTTTAATCATCATTGTc                                             <  1:827633/97‑1 (MQ=14) gacaaCAAGTAACAATGCGAATATGCCGAATAATGACAGTACTCCAACAATAAAGATATCGAACCAATCCATATTTTTAAGTTTTTTAATCATCATTGTc                                             <  1:1074854/97‑1 (MQ=14) gacaaCAAGTAACAATGCGAATATGCCGAATAATGACAGTACTCCAACAATAAAGATATCGAACCAATCCATATTTTTAAGTTTTTTAATCATCATTGTc                                             <  1:2655141/97‑1 (MQ=14)   caaCAAGTAACAATGCGAATATGCCGAATAATGACAGTACTCCAACAATAAAGATATCGAACCAATCCATATTTTTAAGTTTTTTAATCATCATTGTCaa                                           <  1:1597321/99‑2 (MQ=14)   caaCAAGTAACAATGCGAATATGCCGAATAATGACAGTACTCCAACAATAAAGATATCGAACCAATCCATATTTTTAAGTTTTTTAATCATCATTGTCAt                                           <  1:2765804/99‑1 (MQ=14)   caaCAAGTAACAATGCGAATATGCCGAATAATGACAGTACTCCAACAATAAAGATATCGAACCAATCCATATTTTTAAGTTTTTTAATCATCATTGTCAt                                           >  1:3260378/2‑100 (MQ=14)   caaCAAGTAACAATGCGAATATGCCGAATAATGACAGTACTCCAACAATAAAGATATCGAACCAATCCATATTTTTAAGTTTTTTAATCATCATTGTCAt                                           <  1:4967706/99‑1 (MQ=14)         gTAACAATGCGAACATGCCGAATAATGACAGTACTCCAACAATAAAGATATCGAACCAATCCATATTTTTAAGTTTTTTAATCATCATTGTCATCTCCAg                                     <  1:6025373/100‑1 (MQ=17)           aacAATGCGAATATGCCGAATAATGACAGTTCTACAACAATAAAGATATCGAACCAATCCATATTTTTAAGTTTTTTAATCATCATTGTCATCTCCAGtc                                   <  1:5448165/98‑2 (MQ=17)           aacAATGCGAATATGCCGAATAATGACAGTACTCCAACAATAAAGATATCGAACCAATCCATATTTTTAAGTTTTTTAATCATCATTGTCATCTCCAGTa                                   <  1:843059/98‑1 (MQ=18)            acAATGCGAATATGCCGAATAATGACAGTACTCCAACAATAAAGATATCGAACCAATCCATATTTTTAAGTTTTTTAATCATCATTGTCATCTCCAGTaa                                  <  1:2325338/99‑2 (MQ=18)             cAATGCGAATATGCCGAATAATGACAGTACTCCAACAATAAAGATATCGAACCAATCCATATTTTTAAGTTTTTTAATCATCATTGTCATCTCCAGTATc                                 <  1:3855946/100‑1 (MQ=18)             cAATGCGAATATGCCGAATAATGACAGTACTCCAACAATAAAGATATCGAACCAATCCATATTTTTAAGTTTTTTAATCATCATTGTCATCTCCAGTATc                                 <  1:5869015/100‑1 (MQ=18)              aaTGCGAATATGCCGAATAATTACAGTACTCCAACAATAAAGATATCGAACCAATCCATATTTTTAAGTTTTTTAATCACCATTGTCATCTCCAGTATCa                                <  1:6053243/100‑1 (MQ=18)              aaTGCGAATATGCCGAATAATGACAGTACTCCAACA‑‑AAAGATATCGAACCAATCCATATTTTTAAGTTTTTTAATCATCATTGTCATCTCCAGTATCaag                              <  1:3162413/100‑2 (MQ=18)               aTGCGAATATGCCGAATAATGACAGTACTCCAACAATAAAGATATCGAACCAATCCATATTTTTAAGTTTTTTAATCATCATTGTCATCTCCAGTATCag                               <  1:6596744/100‑2 (MQ=18)                 gCGAATATGCCGAATAATGACAGTACTCCAACAATAAAGATATCGAACCAATCCATATTTTTAAGTTTTTTAATCATCATTGTCATCTCCAGTATCaact                             <  1:5265046/100‑3 (MQ=18)                          ccGAATAATGACAGTACTCCAACAATAAAGATATCGAACCAATCCATATTTTTAAGTTTTTTAATCATCATTGTCATCTCCAGTATCAATTAAACTAGGc                    <  1:5705217/100‑1 (MQ=255)                          ccGAATAATGACAGTACTCCAACAATAAAGATATCGAACCAATCCATATTTTTAAGTTTTTTAATCATCATTGTCATCTCCAGTATCAATTAAACTAGGc                    <  1:6201397/100‑1 (MQ=255)                             aataatGACAGTACTCCAACAATAAAGATATCGAACCAATCCATATTTTTAAGTTTTTTAATCATCATTGTCATCTCCAGTATCAATTAAACTAGGcatc                 <  1:2894523/100‑1 (MQ=255)                             aataatGACAGTACTCCAACAATAAAGATATCGAACCAATCCATATTTTTAAGTTTTTTAATCATCATTGTCATCTCCAGTATCAATTAAACTAGGcata                 <  1:1300629/100‑2 (MQ=255)                                            ccAACAATAAAGATATCGAACCAATCCATATTTTTAAGTTTTTTAATCATCATTGTCATCTCCAGTATCAATTAAACTAGGCATCATTCTTAACATAGcc  <  1:3859767/100‑1 (MQ=255)                                                          |                                                                                      GATAACAAGTAGCAATGCGAATATACCGAATAATGACAGTATTCCAACAATAAAGATGTCGAACCCATCCATATTTTTAAGTTTTTTAATCATCATTGTCATCTCCAGTATCAATTAAACTAGGCATCATTCTTAACATAGCC  >  NC_007793/1578533‑1578675

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