Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,125,640	T→A	intergenic (+10/‑68)	SAUSA300_RS10825 → / → SAUSA300_RS10830	hypothetical protein/toxin MazF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,125,640	0	T	A	100.0%	47.4 / NA	19	intergenic (+10/‑68)	SAUSA300_RS10825/SAUSA300_RS10830	hypothetical protein/toxin MazF
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (18/1);  total (18/1)

TTATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATC‑AAAAAACGTCTACAAGTGTAGACGTTGA  >  NC_007793/2125544‑2125669                                                                                                 |                               ttATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAAC‑aa                            >  1:5279555/1‑100 (MQ=255) ttATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAAC‑aa                            >  1:4003007/1‑100 (MQ=255)  gcggggcacaTAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGCTTTTACTAAAC‑aaa                           >  1:2717669/7‑100 (MQ=255)    acGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAAC‑aaaaa                         >  1:2536742/2‑100 (MQ=255)     cGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAAC‑aaaaaa                        >  1:897557/1‑100 (MQ=255)        cacaTAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAAC‑aaaaaaacg                     >  1:2280214/1‑97 (MQ=255)        cacaTAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAAC‑aaaaaaacg                     >  1:4551210/1‑97 (MQ=255)           gtAATGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAAC‑aaaaaaacgcct                  >  1:3240433/2‑94 (MQ=255)             aaTGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAAC‑aaaaaaacgcctac                >  1:3889151/1‑92 (MQ=255)             aaTGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAAC‑aaaaaaacgcctac                >  1:412968/1‑92 (MQ=255)             aaTGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAAC‑aaaaaaacgcccag                >  1:4519552/1‑92 (MQ=255)             aaTGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAAC‑AAAAAACGCCTAct                >  1:3470874/1‑99 (MQ=255)              ctGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAAC‑aaaaaaacgcctagg               <  1:4191690/99‑10 (MQ=255)               tGCGCGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAAC‑aaaaaaacgcctacta              >  1:780784/1‑90 (MQ=255)                gcgcGTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAAC‑AAAAAACGCCTACTAgt             >  1:1873586/1‑100 (MQ=255)                    gTAAATACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAACAAAAAAACGCCTACTAgtgtg          >  1:1192966/1‑99 (MQ=255)                      caaTACTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAACAAAAAAACGCCTACTAgtgtgaa        >  1:5339638/2‑97 (MQ=255)                         caCTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAAC‑AAAAAACGCCTACTAgtgtgaaaacg    >  1:5324164/2‑93 (MQ=255)                           cTTGGAATTTATTTCACAAAAATTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAAAC‑AAAAAACGCCTACTAgtgtgaaaacgta  >  1:225185/1‑91 (MQ=255)                                                                                                 |                               TTATCGTCACATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATC‑AAAAAACGTCTACAAGTGTAGACGTTGA  >  NC_007793/2125544‑2125669

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