Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	656,344	2 bp→TT	coding (127‑128/630 nt)	SAUSA300_RS03115 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	656,344	0	A	T	100.0%	60.0 / NA	23	N43Y (AAT→TAT)	SAUSA300_RS03115	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (7/16);  total (7/16)
*	NC_007793	656,345	0	A	T	100.0%	59.2 / NA	23	N43I (AAT→ATT)	SAUSA300_RS03115	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (7/16);  total (7/16)

ATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGAATAGTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATACTTTCTGGCCTGTTGC  >  NC_007793/656256‑656438                                                                                         ||                                                                                              atatAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCaaa                                                                                     <  1:4281159/100‑1 (MQ=255)     aGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGt                                                                                 >  1:947805/1‑100 (MQ=255)          tattaAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTAcgcg                                                                            <  1:1641751/100‑1 (MQ=17)                   cAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATtagt                                                                   <  1:5308104/100‑1 (MQ=17)                   cAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAAAtagt                                                                   <  1:2980800/100‑1 (MQ=17)                     ataatGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATtagtag                                                                 <  1:407524/100‑1 (MQ=255)                        atGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAgaga                                                              <  1:5823200/100‑1 (MQ=18)                        atGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGCAgaga                                                              <  1:4276375/100‑1 (MQ=17)                                        tGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTa                                              >  1:1144135/1‑100 (MQ=255)                                                 aagtaCTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTc                                     <  1:39371/100‑1 (MQ=25)                                                  agtaCTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCa                                    <  1:489252/100‑1 (MQ=37)                                                  agtaCTTGGATTGCTTATTTTTTGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAAAAGAATTTGAAAGTTTAAATATTGTCa                                    <  1:4510615/100‑1 (MQ=255)                                                                      tttGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATACt                >  1:1668312/1‑100 (MQ=255)                                                                       ttGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATACtt               >  1:5423096/1‑100 (MQ=255)                                                                        tGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATACttt              >  1:4313179/1‑100 (MQ=255)                                                                        tGCCCATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATACttt              >  1:1242174/1‑100 (MQ=25)                                                                          accATGATTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATACTTTCt            >  1:3156774/2‑100 (MQ=25)                                                                                tctAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATACTTTCTGGCCTg      <  1:336619/98‑1 (MQ=21)                                                                                aTTAATTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATACTTTCTGGCCta      <  1:153393/100‑2 (MQ=255)                                                                                    aTTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATACTTTCTGGCCTGTTGc  <  1:2339032/100‑1 (MQ=255)                                                                                    aTTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATACTTTCTGGCCTGTTGc  <  1:600546/100‑1 (MQ=255)                                                                                    aTTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATACTTTCTGGCCTGTTGc  <  1:813629/100‑1 (MQ=255)                                                                                    aTTGGTTTATTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATACTTTCTGGCCTGTTGc  <  1:1737984/100‑1 (MQ=255)                                                                                         ||                                                                                              ATATAGAATTATTAAATACAATAATGAATATTTCATGGTTGATTTAGTAAGTACTTGGATTGCTTATTTTTTGCCCATGATTAATTGGAATAGTCCCAAAAAGTACGCGAAAATTAGTAGAGAAGAATTTGAAAGTTTAAATATTGTCAAACCCGCTAAAAATAATACTTTCTGGCCTGTTGC  >  NC_007793/656256‑656438

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