Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,578,574	T→C	I18V (ATA→GTA)	SAUSA300_RS07710 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,578,574	0	T	C	100.0%	52.8 / NA	26	I18V (ATA→GTA)	SAUSA300_RS07710	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (1/25);  total (1/25)

TGTAATTCTTTGTGTTGGTAACTAGCCACTGTATAGATAGGCAATGTGATAACAAGTAGCAATGCGAATATACCGAATAATGACAGTATTCCAACAATAAAGATGTCGAACCCATCCATATTTTTAAGTTTTTTAATCATCATTGTCATCTCCAGTATCAATTAAACTAGGCATC  >  NC_007793/1578486‑1578660                                                                                         |                                                                                       tgcgtttCTTTGTGTTGGTAACTAGCCACTGTATAGATAGGAAATGTGACAACAAGTAACAATGCGAATATGCCGAATAATGACAGTACTCCAACAATaa                                                                             <  1:2252707/95‑1 (MQ=14) tgcgtttCTTTGTGTTGGTAACTAGCCACTGTATAGATAGGAAATGTGACAACAAGTAACAATGCGAATATGCCGAATAATGACAGTACTCCAACAATaa                                                                             <  1:5817031/95‑1 (MQ=14) tgcgtttCTTTGTGTTGGTAACTAGCCACTGTATAGATAGGAAATGTGACAACAAGTAACAATGCGAATATGCCGAATAATGACAGTACTCCAACAATaa                                                                             <  1:3221420/95‑1 (MQ=14) tgcgtttCTTTGTGTTGGTAACTAGCCACTGTATAGATAGGAAATGTGACAACAAGTAACAATGCGAATATGCCGAATAATGACAGTACTCCAACAATaa                                                                             <  1:4634423/95‑1 (MQ=14)   cgtttCTTTGTGTTGGTAACTAGCCACTGTATAGATAGGAAATGTGACAACAAGTAACAATGCGAATATGCCGAATAATGACAGTACTCCAACAATAAAg                                                                           <  1:5513359/97‑1 (MQ=14)                                                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)                                                                                         |                                                                                       TGTAATTCTTTGTGTTGGTAACTAGCCACTGTATAGATAGGCAATGTGATAACAAGTAGCAATGCGAATATACCGAATAATGACAGTATTCCAACAATAAAGATGTCGAACCCATCCATATTTTTAAGTTTTTTAATCATCATTGTCATCTCCAGTATCAATTAAACTAGGCATC  >  NC_007793/1578486‑1578660

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