Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,125,650	T→C	intergenic (+20/‑58)	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,650	0	T	C	100.0%	52.1 / NA	23	intergenic (+20/‑58)	SAUSA300_RS10825/SAUSA300_RS10830	hypothetical protein/toxin MazF
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (23/0);  total (23/0)

ATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCG  >  NC_007793/2125554‑2125690                                                                                                 |                                         ttAATGCGCGTAAATACTTGGAATTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGCCTa                                       >  1:5261562/2‑100 (MQ=255)   aaTGCGCGTAAATACTTGGAATTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGCCTACa                                     >  1:3567417/1‑100 (MQ=255)   aaTGCGCGTAAATACTTGGAATTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGCCTACa                                     >  1:42298/1‑100 (MQ=255)   aaTGCGCGTAAATACTTGGAATTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGCCTACa                                     >  1:5077637/1‑100 (MQ=255)      gcgcGTAAATACTTGGAATTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGCCTACAAgt                                  >  1:104729/1‑100 (MQ=255)            aaaTACTAGGAATTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGCCTACAAGTGTAGAc                            >  1:1830505/1‑100 (MQ=255)               tACTTGGAATTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGCCTACAAGTGTAGACGtt                         >  1:860886/1‑100 (MQ=255)                aCTTGGAATTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGCCTACAAGTGTAGACGTTg                        >  1:1714362/1‑100 (MQ=255)                aCTTGGAATTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGCCTACAAGTGTAGACGTTg                        >  1:2812171/1‑100 (MQ=255)                aCTTGGAATTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGCCTACAAGTGTAGACGTTg                        >  1:5019572/1‑100 (MQ=255)                aCTTGGAATTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGCCTACAAGTGTAGACGTTg                        >  1:4895812/1‑100 (MQ=255)                   tGGAATTTATTTCAC‑AAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGCCTACAAGTGTAGACGTTGAAtg                    >  1:1480/1‑100 (MQ=255)                             ttCACAAAAAATAGTTGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGCCTACAAGTGTAGACGTTGAATGGTGGTgaga           >  1:6057431/1‑100 (MQ=255)                             ttCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGCCTACAAGTGTAGACGTTGAATGGTGGTgaga           >  1:114355/1‑100 (MQ=255)                             ttCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGCCTACAAGTGTAGACGTTGAATGGTGGTgaga           >  1:6543318/1‑100 (MQ=255)                             ttCACAAAAAATAGATGGAGGTAAGTTAACCAAAGAAGATTCATTATGATTTGACTAATCAAAAAACGCCTACAAGTGTAGACGTTGAACGGTGGTgaga           >  1:3100975/1‑100 (MQ=37)                              tCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGCCTACAAGTGTAGACGTTGAATGGTGGTgagaa          >  1:4310520/1‑99 (MQ=255)                              acacaAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGCCTACAAGTGTAGACGTTGAATGGTGGTgagaa          >  1:4205764/2‑99 (MQ=255)                               gacaAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGCCTACAAGTGTAGACGTTGAATGGTGGTgagatc         >  1:6408828/2‑98 (MQ=255)                               cacaAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGCCTACAAGTGTAGACGTTGAATGGTGGTgagaat         >  1:3921438/1‑98 (MQ=255)                                gcaAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGCCTACAAGTGTAGACGTTGAATGGTGGTgagaatt        >  1:4499636/2‑97 (MQ=255)                                      aaTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGCCTACAAGTGTAGACGTTGAATGGTGGTgagaattttatgg  >  1:2831543/1‑91 (MQ=25)                                      aaTAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGCCTACAAGTGAAGACGTTGAATGGTGGTgagaattttatgg  >  1:948533/1‑91 (MQ=25)                                                                                                 |                                         ATAATGCGCGTAAATACTTGGAATTTATTTCACAAAAAATAGATGGCGATAAGTTAACCAAAGAAGATTCATTATGATTTTACTAATCAAAAAACGTCTACAAGTGTAGACGTTGAATGGTGGTGAGAGTGTGAGCG  >  NC_007793/2125554‑2125690

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