Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,087,424	T→C	intergenic (+79/+16)	SAUSA300_RS10530 → / ← SAUSA300_RS10535	chemotaxis inhibitory protein/peptidoglycan hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,087,424	0	T	C	100.0%	56.1 / NA	22	intergenic (+79/+16)	SAUSA300_RS10530/SAUSA300_RS10535	chemotaxis inhibitory protein/peptidoglycan hydrolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (10/12);  total (10/12)

AATGAATATGCATACTAATAGTAGTTACATAAATTAAAAGGTAGATATTTCTTTTTTATATAAAGGTTTGGCAGACATTTCATAACTTGCCAAACCTTTATATATCTAATTATCAAACTGCACTAAACTTACCAAAACTGCTTATTCTATTACCTGCCTTGTCTACCTCTCCTGTCGCTATATAACGACGTT  >  NC_007793/2087328‑2087519                                                                                                 |                                                                                                aaTGAATATGCATACTAATAGTAGTTACATAAATTAAAAGGTAGATATTTCTTTTTTATATAAAGGTTTGGCAGACATTTCATAACTTGCCAAACCCTta                                                                                              >  1:654465/1‑100 (MQ=255)                          tACATAAATTAAAAGGTAGATATTTCTTTTTTATATAAAGGTTTGGCAGACATTTCATAACTTGCCAAACCCTTATATATCTAATTATCAAACTGCAcag                                                                     >  1:1111957/1‑98 (MQ=255)                           aCATAAATTAAAAGGTAGATATTTCTTTTTTATATAAAGGTTTGGCAGACATTTCATAACTTGCCAAACCCTTATATATCTAATTATCAAACTGCACTaa                                                                    >  1:2493149/1‑100 (MQ=255)                               aacttAAACGGTAGATATTTCTTTTTTATATAAAGGATTGGCAGACATTTCATAACTTGCCAAACCCTTATATATCTAATTATCAAACTGCACTAAACtt                                                                >  1:2692108/4‑100 (MQ=255)                                  ttAAAAGGTAGATATTTCTTTTTTATATAAAGGTTTGGCAGACATTTCATAACTTGCCAAACCCTTATATATCTAATTATCAAACTGCACTAAACTTAcc                                                             <  1:4522646/100‑1 (MQ=255)                                               atTTCTTTTTTATATAAAGGTTTGGCAGACATTTCATAACTTGCCAAACCCTTATATATCTAATTATCAAACTGCACTAAACTTACCAAAACTGCTTAtt                                                <  1:3571470/100‑1 (MQ=255)                                                              aaaGGTTTGGCAGACATTTCATAACTTGCCAAACCCTTATATATCTAATTATCAAACTGCACTAAACTTACCAAAACTGCTTATTCTATTACCTGCCTTg                                 <  1:4394434/100‑1 (MQ=255)                                                              aaaGGTTTGGCAGACATTTCATAACTTGCC‑AACCCTTATATATCTAATTATCAAACTGCACTAAACTTACCAAAACTGCTTATTCTATTACCTGCCTTGt                                <  1:4503473/100‑1 (MQ=255)                                                                   tttGGCAGACATTTCATAACTTGCCAAACCCTTATATATCTAATTATCAAACTGCACTAAACTTACCAAAACTGCTTATTCTATTACCTGCCTTGTCTAc                            <  1:1929045/100‑1 (MQ=255)                                                                               ttccgatcTTGCCAGACCCTTATATATCTAATTATCAAACTGCACTAAACTTACCAAAACTGCTTATTCTATTACCTGCCTTGTCTACCTCTCCTGTCGc                <  1:4781993/93‑1 (MQ=21)                                                                                tCATAACTTGCCAAACCCTTATATATCTAATTATCAAACTGCACTAAACTTACCAAAACTGCTTATTCTATTACCTGCCTTGTCTACCTCTCCTGTCGca               <  1:4334027/100‑2 (MQ=255)                                                                                 ccgatcTTGCAAAACCCTTATATATCTAATTATCAAACTGCACTAAACTTACCAAAACTGCTTATTCTATTACCTGCCTTGTCTACCTCTCCTGTCGCtt              <  1:1889646/95‑2 (MQ=21)                                                                                  aTAACGTGCCAAACCCTTATATATCTAATTATCAAACTGCACTAAACTTACCAAAACTGCTTATTCTATTACCTGCCTTGTCTACCTCTCCTGTCGCtat             <  1:1009049/100‑1 (MQ=255)                                                                                   gatcTTGCCAAACCCTTATATATCTAATTATCAAACTGCACTAAACTTACCAAAACTGCTTATTCTATTACCTGCCTTGTCTACCTCTCCTGTCGCtata            <  1:1812047/97‑1 (MQ=255)                                                                                    aaCTTGCCAAACCCTTATATATCTAATTATCAAACTGCACTAAACTTACCAAAACTGCTTATTCTATTACCTGCCTTGTCTACCTCTCCTGTCGCtatat           <  1:1805443/100‑1 (MQ=255)                                                                                     tcTTGTCAAACCCTTATATATCTAATTATCCAACTGCACTAAACTTACCAAAACTGCTTATTCTATTACCTGCCTTGTCTACCTCTCCTGTCGCtatata          <  1:2741657/99‑1 (MQ=255)                                                                                     aCTTGCCAAACCCTTATATATCTAATTATCAAACTGCACTAAACTTACCAAAACTGCTTATTCTATTACCTGCCTTGTCTACCTCTCCTGTCGCtatata          >  1:2469340/1‑100 (MQ=255)                                                                                      tttGCCAAACCCTTATATATCTAATTATCAAACTGCACTAAACTTACCAAAACTGCTTATTCTATTACCTGCCTTGTCTACCTCTCCTGTCGCTATATaa         >  1:5605760/2‑100 (MQ=255)                                                                                        tGCCAAACCCTTATATATCTAATTATCAAACTGCACTAAACTTACCAAAACTGCTTATTCTATTACCTGCCTTGTCTACCTCTCCTGTCGCTATATAacg       >  1:5786361/1‑100 (MQ=255)                                                                                        tGCCAAACCCTTATATATCTAATTATCAAACTGCACTAAACTTACCAAAACTGCTTATTCTATTACCTGCCTTGTCTACCTCTCCTGTCGCTATATAacg       >  1:607743/1‑100 (MQ=255)                                                                                          ccAAACCCTTATATATCTAATTATCAAACTGCACTAAACTTACCAAAACTGCTTATTCTATTACCTGCCTTGTCTACCTCTCCTGTCGCTATATAacgac     >  1:2342122/1‑100 (MQ=255)                                                                                             aaCCCTTATATATCTAATTATCAAACTGCACTAAACTTACCAAAACTGCTTATTCTATTACCTGCCTTGTCTACCTCTCCTGTCGCTATATAACGACGtt  >  1:5568460/1‑100 (MQ=255)                                                                                                 |                                                                                                AATGAATATGCATACTAATAGTAGTTACATAAATTAAAAGGTAGATATTTCTTTTTTATATAAAGGTTTGGCAGACATTTCATAACTTGCCAAACCTTTATATATCTAATTATCAAACTGCACTAAACTTACCAAAACTGCTTATTCTATTACCTGCCTTGTCTACCTCTCCTGTCGCTATATAACGACGTT  >  NC_007793/2087328‑2087519

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