Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	323,065	C→A	G241G (GGC→GGA)	SAUSA300_RS01450 →	glycyl‑glycine endopeptidase LytM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	323,065	0	C	A	100.0%	101.0 / NA	36	G241G (GGC→GGA)	SAUSA300_RS01450	glycyl‑glycine endopeptidase LytM
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (34/2);  total (34/2)

CGGTGTCGACTATGCAATGCCTGAAAATTCACCAGTTTACTCATTAACTGATGGTACAGTAGTACAAGCAGGTTGGAGTAACTATGGTGGCGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGTATATGCATAATAATCGTTTAACTGTTTCAGCTGGTGATAAAGTCAA  >  NC_007793/322975‑323160                                                                                           |                                                                                                atctggcGACTATGCAATGCCTGAAAATTCACCAGTTTACTCATTAACTGATGGTACAGTAGTACAAGCAGGTTGGAGTAACTATGGTGGAGGCAATCac                                                                                        <  1:2215149/94‑2 (MQ=255)      ccGACTATGCAATGCCTGAAAATTCACCAGTTTACTCATTAACTGATGGTACAGTAGTACAAGCAGGTTGGAGTAACTATGGTGGAGGCAATCAAGTAAc                                                                                   >  1:1690980/2‑100 (MQ=255)       cGACTATGCAATGCCTGAAAATTCACCAGTTTACTCATTAACTGATGGTACAGTAGTACAAGCAGGTTGGAGTAACTATGGTGGAGGCAATCAAGTAACg                                                                                  >  1:2128247/1‑100 (MQ=255)                aaTGCCTGAAAATTCACCAGTTTACTCATTAACTGATGGTACAGTAGTACAAGCAGGTTGGAGTAACTATGGTGGAGGCAATCAAGTAACGATTAaagaa                                                                         >  1:5464961/1‑100 (MQ=255)                                          cATTAACTGATGGTACAGTAGTACAAGCAGGTTGGAGTAACTATGGTGGAGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGta                                               >  1:5062501/1‑100 (MQ=255)                                                      ttACAGTAGTACAAGCAGGTTGGAGTAACTATGGTGGAGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGTATATGCAtaataa                                   >  1:964894/2‑100 (MQ=255)                                                        ccAGTAGTACAAGCAGGTTGGAGTAACTATGGTGGAGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGTATATGCATAATAATc                                 >  1:2426454/2‑100 (MQ=255)                                                        aCAGTAGTACAAGCAGGTTGGAGTAACTATGGTGGAGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGTATATGCATAATAATc                                 >  1:1110030/1‑100 (MQ=255)                                                        aCAGTAGTACAAGCAGGTTGGAGTAACTATGGTGGAGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGTATATGCATAATAATc                                 >  1:22833/1‑100 (MQ=255)                                                           gtggtaCAAGCAGGTTGGAGTAACTATGGTGGAGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGTATATGCATAATAATCGtt                              >  1:555765/4‑100 (MQ=255)                                                             ggtaCAAGCAGGTTGGAGTAACTATGGTGGAGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGTATATGCATAATAATCGTTTa                            >  1:3081773/2‑100 (MQ=255)                                                              gtaCAAGCAGGTTGGAGTAACTATGGTGGAGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGTATATGCATAATAATCGTTTaa                           >  1:1424836/1‑100 (MQ=255)                                                               gaCAAGCAGGTTGGAGTAACTATGGTGGAGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGTATATGCATAATAATCGTTTAAc                          >  1:3839129/2‑100 (MQ=255)                                                                aCAAGCAGGTTGGAGTAACTATGGTGGAGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGTATATGCATAATAATCGTTTAACt                         >  1:1517713/1‑100 (MQ=255)                                                                 cAAGCAGGTTGGAGTAACTATGGTGGAGGCAATCAAGTAACGATTAACGAAGCGAACAGTAATAACTACCAATGGTATATGCATAATAATCGTTTAACTg                        >  1:3125688/1‑100 (MQ=255)                                                                  taGCAGGTTGGAGTAACTATGGTGGAGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGTATATGCATAATAATCGTTTAACTGt                       >  1:1058797/2‑100 (MQ=255)                                                                  caGCAGGTTGGAGTAACTATGGTGGAGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGTATATGCATAATAATCGTTTAACTGt                       >  1:195387/2‑100 (MQ=255)                                                                  aaGCAGGTTGGAGTAACTATGGTGGAGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGTATATGCATAATAATCGTTTAACTGt                       >  1:4122581/1‑100 (MQ=255)                                                                   tgCAGGTTGGAGTAACTATGGTGGAGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGTATATGCATAATAATCGTTTAACTGtt                      >  1:1966480/2‑100 (MQ=255)                                                                   cccAGGTTGGAGTAACTATGGTGGAGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGTATATGCATAATAATCGTTTAACTGtt                      >  1:2433522/3‑100 (MQ=255)                                                                   aGCAGGTTGGAGTAACTATGGTGGAGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGTATATGCGTAATAATCGTTTAACTGtt                      >  1:410788/1‑100 (MQ=255)                                                                   aGCAGGTTGGAGTAACTATGGTGGAGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGTATATGCATAATAATCGTTTAACTGtt                      >  1:1678011/1‑100 (MQ=255)                                                                    gCAGGTTGGAGTAACTATGGTGGAGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGTATATGCATAATAATCGTTTAACTGttt                     >  1:1981225/1‑100 (MQ=255)                                                                    acAGGTTGGAGTAACTATGGTGGAGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGTATATGCATAATAATCGTTTAACTGttt                     >  1:3835534/2‑100 (MQ=255)                                                                             aGTAACTATGGTGGAGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGTATATGCATAATAATCGTTTAACTGTTTCAGCTGGTg            >  1:1278944/1‑100 (MQ=255)                                                                             aGTAACTATGGTGGAGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGTATATGCATAATAATCGTTTAACTGTTTCAGCTGGTg            >  1:1785265/1‑100 (MQ=255)                                                                              gTAACTATGGTGGAGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGTATATGCATAATAATCGTTTAACTGTTTCAGCTGGTGa           >  1:4034845/1‑100 (MQ=255)                                                                                aaCTATGGTGGAGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGTATATGCATAATAATCGTTTAACTGTTTCAGCTGGTGATa         >  1:2348858/1‑100 (MQ=255)                                                                                 aCTATGGTGGAGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGTATATGCATAATAATCGTTTAACTGTTTCAGCTGGTGATaa        >  1:2455323/1‑100 (MQ=255)                                                                                    aTGGTGGAGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGTATATGCATAATAATCGTTTAACTGTTTCAGCTGGTGATAAAGt     >  1:3351391/1‑100 (MQ=255)                                                                                    aTGGTGGAGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGTATATGCATAATAATCGTTTAACTGTTTCAGCTGGTGATAAAGt     >  1:3223573/1‑100 (MQ=255)                                                                                    aTGGTGGAGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGTATATGCATAATAATCGTTTAACTGTTTCAGCTGGTGATAAAGt     >  1:1316014/1‑100 (MQ=255)                                                                                    aTGGTGGAGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGTATATGCATAATAATCGTTTAACTGTTTCAGCTGGTGATAAAGt     >  1:5437281/1‑100 (MQ=255)                                                                                    aTGGTGGAGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGTATATGCATAATAATCGTTTAACTGTTTCAGCTGGTGATAAAGt     >  1:775553/1‑100 (MQ=255)                                                                                    aTGGTGGAGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGAATATGCATAATAATCGTTTAACTGTTTCAGCTGGTGATAAAGt     >  1:3785844/1‑100 (MQ=255)                                                                                       gtggAGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGTATATGCATAATAATCGTTTAACTGTTTCAGCTGGTGATAAAGTCaa  <  1:1752636/100‑1 (MQ=255)                                                                                           |                                                                                                CGGTGTCGACTATGCAATGCCTGAAAATTCACCAGTTTACTCATTAACTGATGGTACAGTAGTACAAGCAGGTTGGAGTAACTATGGTGGCGGCAATCAAGTAACGATTAAAGAAGCGAACAGTAATAACTACCAATGGTATATGCATAATAATCGTTTAACTGTTTCAGCTGGTGATAAAGTCAA  >  NC_007793/322975‑323160

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: ‑