Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,590,291	2 bp→TC	coding (316‑317/615 nt)	SAUSA300_RS07850 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,590,291	0	A	T	100.0%	73.2 / NA	28	I106K (ATA→AAA)	SAUSA300_RS07850	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (8/20);  total (8/20)
*	NC_007793	1,590,292	0	T	C	100.0%	73.5 / NA	28	I106V (ATA→GTA)	SAUSA300_RS07850	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (8/20);  total (8/20)

ATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCTTATACAATCCAATTTCATAAACTTTCTGGTTATATTAATTAACCAGATAATAAATTTAACAATTAAAGGATTAAATACAAACACTGTTAAAACAAAAATAAATAGAAACAAAATATTT  >  NC_007793/1590204‑1590377                                                                                        ||                                                                                      ataAATACTGGTTTACCGTTGTTATTAAACAGCTTGTCTCTTTTGTCTAATAAGCTTATACAATCCAAATTCATAAACTTTCTTGTTTCATTAATTAAcc                                                                            <  1:2052495/100‑1 (MQ=255) ataAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTGTTTTGTCTAATAAGCTTATACAATCCAAATTCATAAACTGTCTTGTTTCATTAATTAAcc                                                                            <  1:4381741/100‑1 (MQ=255) ataAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCTTATACAATCCAAATTCATAAACTTTCTTGTTTCATTAATTAAcc                                                                            <  1:1218352/100‑1 (MQ=255) ataAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCTTATACAATCCAAATTCATAAACTTTCTTGTTTCATTAATTAAcc                                                                            <  1:4371282/100‑1 (MQ=255) ataAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCTTATACAATCCAAATTCATAAACTTTCTTGTTTCATTAATTAAcc                                                                            <  1:1932330/100‑1 (MQ=255) ataAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCTTATACAATCCAAATTCATAAACTTTCTTGTTTCATTAATTAAcc                                                                            <  1:2212743/100‑1 (MQ=255) ataAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCTTATACAATCCAAATTCATAAACTTTCTTGTTTCATTAATTAAcc                                                                            <  1:5418031/100‑1 (MQ=255) ataAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCTTATACAATCCAAATTCATAAACTTTCTTGTTTCATTAATTAAcc                                                                            <  1:2346551/100‑1 (MQ=255)   aaaTACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCTTATACAATCCAAATTCATAAACTTTCTTGTTTCATTAATTAACCAg                                                                          >  1:5220307/1‑100 (MQ=255)   aaaTACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTCGTCTAATAAGCTTATACAATCCAAATTCATAAACTTTCTTGTTTCATTAATTAACCAg                                                                          <  1:1973361/100‑1 (MQ=255)    aaTACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCTTATACAATCCAAATTCATAAACTTTCTTGTTTCATTAATTAACCAgc                                                                         <  1:1048710/100‑2 (MQ=255)     aTACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCTTATACAATCCAAATTCATAAACTTTCTTGTTTCATTAATTAACCAGat                                                                        <  1:1118515/100‑1 (MQ=255)     aTACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCTTATACAATCCAAATTCATAAACTTTCTTGTTTCATTAATTAACCAGat                                                                        <  1:2344385/100‑1 (MQ=255)          ggTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCTTATACAATCCAAATTCATAAACTTTCTTGTTTCATTAATTAACCAGataatga                                                                   >  1:2045981/1‑98 (MQ=255)          ggTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCTTATACAATCCAAATTCATAAACTTTCTTGTTTCATTAATTAACCAGataatga                                                                   >  1:6271959/1‑98 (MQ=255)                             acaacTTGTCTCTTTTGTCTAATAAGCTTATACAATCCAAATTCATAAACTTTCTTGTTTCATTAATTAACCAGATAATGAATTTAACAATTAAAGGAtt                                                <  1:3670082/100‑1 (MQ=255)                                acTTGTCTCTTTTGTCTAATAAGCTTATACAATCCAAATTCATAAACTTTCTTGTTTCATTAATTAACCAGATAATGAATTTAACAATTAAAGGATTaaa                                             >  1:4597020/1‑100 (MQ=255)                                       tctTTTGTCTAATAAGCTTATACAATCCAAATTCATAAACTTTCTTGTTTCATTAATTAACCAGATAATGAATTTAACAATTAAAGGATTAAATACaagc                                      <  1:2585690/100‑3 (MQ=255)                                        ctTTTGTCTAATAAGCTTATACAATCCAAATTCATAAACTTTCTTGTTTCATTAATTAACCAGATAATGGATTTAACAATTAAAGGATTAAATACAAgca                                     >  1:5106136/1‑100 (MQ=255)                                              tCTAATAAGCTTATACAATCCAAATTCATAAACTTTCTTGTTTCATTAATTAACCAGATAATGAATTTAACAATTAAAGGATTAAATACAAGCACTGTTa                               <  1:2840696/100‑1 (MQ=255)                                              tCTAATAAGCTTATACAATCCAAATTCATAAACTTTCTTGTTTCATTAATTAACCAGATAATGAATTTAACAATTAAAGGATTAAATACAAGCACTGTTa                               <  1:5150138/100‑1 (MQ=255)                                                taataaGCTTATACAATCCAAATTCATAAACTTTCTTGTTTCATTAATTAACCAGATAATGAATTTAACAATTAAAGGATTAAATACAAGCACTGTTaaa                             <  1:2341228/100‑1 (MQ=255)                                                 aataaGCTTATACAATCCAAATTCATAAACTTTCTTGTTTCATTAATTAACCAGATAATGAATTTAACAATTAAAGGATTAAATACAAGCACTGTTaaaa                            >  1:7581/1‑100 (MQ=255)                                                          ataCAATCCAAATTCATAAACTTTCTTGTTTCATTAATTAACCAGATAATGAATTTAACAATTAAAGGATTAAATACAAGCACTGTTAAAACAAaaataa                   >  1:1553138/1‑100 (MQ=255)                                                                 ccAAATTCATAAACTTTCTTGTTTCATTAATTAACCAGATAATGAATTTAACAATTAAAGGATTAAATACAAGCACTGTTAAAACAAAAATAATTAGagg            <  1:6264852/100‑3 (MQ=255)                                                                    aattCATAAACTTTTTTGTTTCATTAATTAACCAGATAATGAATTTAACAATTAAAGGATTAAATACAAGCACTGTTAAAACAAAAATAATTAGAAAcca         <  1:4927965/98‑3 (MQ=255)                                                                    aattCATAAACTTTCTTGTTTCATTAATTAACCAGATAATGAATTTAACAATTAAAGGATTAAATACAAGCACTGTTAAAACAAAAATAATTAGAAAcca         <  1:2649252/98‑3 (MQ=255)                                                                           aaaCTTTCTTGTTTCATTAATTAACCAGATAATGAATTTAACAATTAAAGGATTAAATACAAGCACTGTTAAAACAAAAATAATTAGAAACCAAATAttt  >  1:753239/1‑100 (MQ=255)                                                                                        ||                                                                                      ATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCTTATACAATCCAATTTCATAAACTTTCTGGTTATATTAATTAACCAGATAATAAATTTAACAATTAAAGGATTAAATACAAACACTGTTAAAACAAAAATAAATAGAAACAAAATATTT  >  NC_007793/1590204‑1590377

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