Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,640,117	C→T	M59I (ATG→ATA)	SAUSA300_RS08100 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,640,117	0	C	T	100.0%	44.9 / NA	17	M59I (ATG→ATA)	SAUSA300_RS08100	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (5/12);  total (5/12)

TTTCCATCTTGTTTTAACTTAACAATATATTTCGTTATATAATAAGGTGCAACCGTAAAGTAAATAAACATAAGTCCACCAATCGACATAATTACAGGCTGTAATAAATACCAAAATAACAAAAATAATATAAACCACACAATCAAAAATGTAGCCCAAAAATCTAAACGAAAACCCT  >  NC_007793/1640031‑1640208                                                                                       |                                                                                            ttccgaTCTTGTCTTAACTTAACAATATATTTCGTTATATAATAAGGTGCAACCGTAAAGTAAATAAACATAAGTCCACCAATCGATATAATTACAGGCt                                                                                <  1:4301617/95‑1 (MQ=255)  tccgaTCTTGTTCTAACTTAACAATATATTTCGTTATATAATAAGGTGCAACCGTAAAGTAAATAAACATAAGTCCACCAATCGATATAATTACAGGCTg                                                                               <  1:99567/96‑1 (MQ=255)   ccgaTCTTGTTTTAACTTAACAATATATTTCGTTATATAATAAGGTGCAACCGTAAAGTAAATAAACATAAGTCCACCAATCGATATAATTACAGGCTGt                                                                              <  1:2922710/97‑1 (MQ=255)    cgaTCTTGTTTTAACTTAACAATATATTTCGTTATATAATAAGGTGCAACCGTAAAGTAAATAAACATAAGTCCACCAATCGATATAATTACAGGCTGta                                                                             <  1:421036/98‑1 (MQ=255)    cgaTCTTGTTTAAACTTAACAATATATTTCGTTATATAATAAGGTGCAACCGTAAAGTAAATAAACATAAGTCCACCAATCGATATAATTACGGGCTGtt                                                                             <  1:3405925/98‑2 (MQ=255)               taacttaacAATATATTTCGTTATATAATAAGGTGCAACCGTAAAGTAAATAAACATAAGTCCACCAATCGATATAATTACAGGCTGTAATAAATACCat                                                                  <  1:575867/100‑2 (MQ=255)               taacttaacAATATATTTCGTTATATAATAAGGTGCAACCGTAAAGTAAATAAACATAAGTCCACCAATCGATATAATTACAGGCTGTAATAAATACCaa                                                                  >  1:1758924/1‑100 (MQ=255)                                              ggTGCAACCGTAAAGTAAATAAACATAAGTCCACCAATCGATATAATTACAGGCTGTAATAAATACCAAAATAACAAAAATAATATAAACCACACAATCa                                   >  1:2853200/1‑100 (MQ=255)                                                tGCAACCGTAAAGTAAATAAACATAAGTCCACCAATCGATATAATTACAGGCTGTAATAAATACCAAAATAACAAAAATAATATAAACCACACAATCaaa                                 >  1:5695439/1‑100 (MQ=255)                                                 gCAACCGTAAAGTAAATAAACATAAGTCCACCAATCGATATAATTACAGGCTGTAATAAATACCAAAATAACAAAAATAATATAAACcacacagatcgga                                >  1:626261/1‑93 (MQ=255)                                                    aCCGTAAAGTAAATAAACATAAGTCCACCAATCGATATAATTACAGGCTGTAATAAATACCAAAATAACAAAAATAATATAAACCACACAATCAAAAATg                             >  1:3653981/1‑100 (MQ=255)                                                     ccGTAAAGTAAATAAACATAAGTCCACCAATCGATATAATTACAGGCTGTAATAAATACCAAAATAACAAAAATAATATAAACCACACAATCAAAAATgc                            <  1:4461680/100‑2 (MQ=255)                                                     ccGTAAAGTAAATAAACATAAGTCCACCAATCGATATAATTACAGGCTGTAATAAATACCAAAATAACAAAAATAATATAAACCACACAATCAAAAATGt                            <  1:2732063/100‑1 (MQ=255)                                                      cGTAAAGTAAATAAACATAAGTCCACCAATCGATATAATTACAGGCTGTAATAAATACCAAAATAACAAAAATAATATAAACCACACAATCAAAAATGTa                           <  1:1398166/100‑1 (MQ=255)                                                        taaagtaaaTAAACATAAGTCCACCAATCGATATAATTACAGGCTGTAATAAATACCAAAATAACAAAAATAATATAAACCACACAATCAAAAATGTAGc                         <  1:3071742/100‑1 (MQ=255)                                                          aagtaaaTAAACATAAGTCCACCAATCGATATAATTACAGGCTGTAATAAATACCAAAATAACAAAAATAATATAAACCACACAATCAAAAATGTAGccc                       <  1:1651988/100‑1 (MQ=255)                                                                               ccaATCGATATAATTACAGGCTGTAATAAATACCAAAATAACAAAAATAATATAAACCACACAATCAAAAATGTAGCCCAAAAATCTAAACGAAAACCCt  <  1:2579417/100‑1 (MQ=255)                                                                                       |                                                                                            TTTCCATCTTGTTTTAACTTAACAATATATTTCGTTATATAATAAGGTGCAACCGTAAAGTAAATAAACATAAGTCCACCAATCGACATAATTACAGGCTGTAATAAATACCAAAATAACAAAAATAATATAAACCACACAATCAAAAATGTAGCCCAAAAATCTAAACGAAAACCCT  >  NC_007793/1640031‑1640208

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