Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,579,936	T→C	E44E (GAA→GAG)	SAUSA300_RS07735 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,579,936	0	T	C	100.0%	44.3 / NA	21	E44E (GAA→GAG)	SAUSA300_RS07735	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (2/19);  total (2/19)

TATAGACCACTTGCCCCGTCATAATCCCTACTGCTTCATCAAGACCAATATCTTCTTTGAGTGCATCTTGCATAGCATTAGGTAAACCTTCAAGTATTTCATCAAACGCTTGCGCTTTCTTATACACGTCTTCAATCTCTTTTAGTAATCCCTCTGTGTCATTGCCGTTATACGCACTAGCA  >  NC_007793/1579848‑1580029                                                                                         |                                                                                              tcttccgatcTTGCCCGGTCATAATCCCTACTGCTTCATCAAGACCAATATCTTCTTTGAGTGCATCTTGCATAGCATTAGGTAAACCCTCAAGTATTtc                                                                                    <  1:2909184/91‑1 (MQ=255) tataGACAACTTGACCCGTCATAATCCCTACTGCTTCATCAAGTTCAATACCTTCTTTAACTGAATGTTGAATAGCATTTGTCATTCCCTCAAGTATTtc                                                                                    <  1:2839188/100‑1 (MQ=255)   taGACAACTTGACCCGTCATAATCCCTACTGCTTCATCAAGTTCAATACCTTCTTTAACTGAATGTTGAATAGCATTTGTCATTCCCTCAAGTATTtcat                                                                                  <  1:3954566/100‑1 (MQ=255)     gACCACTTGCCCCGTCATAATCCCTACTGCTTCATCAAGACCAATATCTTCTTTGAGTGCATCTTGCATAGCATTAGGTAAACCCTCAAGTATTtcatta                                                                                <  1:2429135/100‑3 (MQ=255)     gACAACTTGACCCGTCATAATCCCTACTGCTTCATCAAGTTCAATACCTTCTTTAACTGAATGTTGAATAGCATTTGTCATTCCCTCAAGTATTtcatca                                                                                <  1:2314003/100‑1 (MQ=255)              acccGTCATAATCCCTACTGCTTCAACAAGTTCAATACCTTCTTTAACTGAATGTTGAATAGCATTTGTCATTCCCTCAAGTATTTCATCAAACGCTTgc                                                                       <  1:6138619/99‑1 (MQ=255)                  gTCATAATCCCTACTGCTTCATCAAGTTCAATACCTTCTTTAACTGAATGTTGAATAGCATTTGTCATTCCCTCAAGTATTTCATCAAACGCTTGCGCtg                                                                   <  1:3466162/100‑2 (MQ=255)                             tACTGCTTCATCAAGACCAATATCTTCTTTGAGTGCATCTTGCATAGCATTAGGTAAACCCTCAAGTATTTCATCAAACGCTTGTGCTTTCTTATacacc                                                        <  1:4086671/100‑2 (MQ=255)                                  cTTCATCAAGACCAATATCTTCTTTGAGTGCATCTTGCATAGCATTAGGTAAACCCTCAAGTATTTCATCAAACGCTTGTGCTTTCTTATACACGTCCTc                                                   <  1:3080428/100‑1 (MQ=255)                                            aCCAATATCTTCTTTGAGTGCATCTTGCATAGCATTAGGTAAACCCTCAAGTATTTCATCAAACGCTTGTGCTTTCTTATACACGTCCTCAATCTCtttt                                         <  1:5115415/100‑1 (MQ=255)                                              cAATACCTTCTTTAACTGTATGTTGAATAGCATTTGTCATTCCCTCCAGTATTTCATCAAACGCTTGCGCTTTCTTATATACGTCCTCAATCTCTTTTAg                                       >  1:5519834/1‑100 (MQ=9)                                                atacCTTCTTTAACTGAATGTTGAATAGCATTTGTCATTCCCTCAAGTATTTCATCAAACGCTTGCGCTTTCTTATATACGTCCTCAATCTCTTTTAGTa                                     <  1:640791/100‑1 (MQ=9)                                                      tcttTGAGTGCATCTTGCATAGCATTAGGTAAACCCGCAAGTATTTCATCAAACGCTTGTGCTTTCTTATACACGGCCTCAATCGCTTTTAGCAATCCct                               <  1:3203830/100‑1 (MQ=21)                                                             ctgaatgttGAATAGCATTTGTCATTCCCTCAAGTATTTCATCAAACGCTTGCGCTTTCTTATATACGTCCTCAATCTCTTTTAGTAATCCCTCTGTGTc                        <  1:5710706/93‑1 (MQ=14)                                                             ctgaatgttGAATAGCATTTGTCATTCCCTCAAGTATTTCATCAAACGCTTGCGCTTTCTTATATACGTCCTCAATCTCTTTTAGTAATCCCTCTGTGTc                        <  1:148918/93‑1 (MQ=14)                                                               gaatgttGAATAGCATTTGTCATTCCCTCAAGTATTTCATCAAACGCTTGCGCTTTCTTATATACGTCCTCAATCTCTTTTAGTAATCCCTCTGTGTCag                      <  1:1708752/95‑2 (MQ=14)                                                               gaatgttGAATAGCATTTGTCATTCCCTCAAGTATTTCATCAAACGCTTGCGCTTTCTTATATACGTCCTCAATCTCTTTTAGTAATCCCTCTGTGTCAt                      <  1:5183183/95‑1 (MQ=14)                                                                 atgttGAATAGCATTTGTCATTCCCTCAAGTATTTCATCAAACGCTTGCGCTTTCTTATATACGTCCTCAATCTCTTTTAGTAATCCCTCTGTGTCAtga                    <  1:352948/97‑3 (MQ=14)                                                                     tGCATAGCATTAGGTAAACCCTCAAGTATTTCATCAAACGCTTGTGCTTTCTTATACACGTCCTCAATCTCTTTTAGCAATCCCTCCGCGTCagatcgga                >  1:3569134/1‑93 (MQ=17)                                                                        aTAGCATTAGGTAAACCCTCAAGTATTTCATCAAACGCTTGTGCTTTCTTATACACGTCCTCAATCTCTTTTAGCAATCCCTCTGTGTCATTACCGTtat             <  1:3264323/100‑1 (MQ=17)                                                                                   tAAACCCTCAAGTATTTCATCAAACGCTTGTGCTTTCTTATACACGTCCTCAATCTCTTTTAGCAATCCCTCTGTGTCATTACCGTTATACGCACTAggc  <  1:2068406/100‑3 (MQ=14)                                                                                         |                                                                                              TATAGACCACTTGCCCCGTCATAATCCCTACTGCTTCATCAAGACCAATATCTTCTTTGAGTGCATCTTGCATAGCATTAGGTAAACCTTCAAGTATTTCATCAAACGCTTGCGCTTTCTTATACACGTCTTCAATCTCTTTTAGTAATCCCTCTGTGTCATTGCCGTTATACGCACTAGCA  >  NC_007793/1579848‑1580029

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