Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,131,773	A→G	S160G (AGC→GGC)	SAUSA300_RS10855 →	succinyl‑diaminopimelate desuccinylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,131,773	0	A	G	100.0%	155.0 / NA	57	S160G (AGC→GGC)	SAUSA300_RS10855	succinyl‑diaminopimelate desuccinylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (45/12);  total (45/12)

CATGGAACGATTAGATTACTGGCTACTGCTGGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATAGCTTAATTATTGCTGAACCAACTGGATCTGGAATTTATTATGCACATAAGGGGTCTATGTCATGTAAAGTAACTGCAACTGGTAAAGC  >  NC_007793/2131677‑2131861                                                                                                 |                                                                                         cATGGAACGATTAGATTACTGGCTACTGCTGGCGAAGAGCAAGAACAAGAAGGTGCCAAATTATTAGCTGTTAAAGGCTATTTAGACGATGTCGATGGCt                                                                                       >  1:1384667/1‑100 (MQ=255) cATGGAACGATTAGATTACTGGCTACTGCTGGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCt                                                                                       >  1:4634642/1‑100 (MQ=255) cATGGAACGATTAGATTACTGGCTACTGCTGGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCt                                                                                       >  1:2362999/1‑100 (MQ=255)     aaggccctagattaCTGGCTACTGCTGGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAat                                                                                   >  1:4480988/8‑100 (MQ=255)                  aCTGGCTACTGCTGGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCa                                                                      >  1:4700773/1‑100 (MQ=255)                  aCTGGCTACTGCTGGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCa                                                                      >  1:3745986/1‑100 (MQ=255)                     ggCTACTGCTGGCGCAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAGCt                                                                   >  1:3241433/1‑100 (MQ=255)                     ggCTACTGCTGGCGAAGAGAAAGAACAAGCAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAaacaacc                                                                   >  1:299582/1‑99 (MQ=255)                     ggCTACTGCTGGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGTCGATGTCGATGGCTTAATTATTGCTGAACCacat                                                                   >  1:2640049/1‑97 (MQ=255)                     ggCTACTGCTGGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTCTTGCTGAACCAACt                                                                   >  1:4881378/1‑100 (MQ=255)                     ggCTACTGCTGGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTgaacgaacg                                                                   >  1:2449126/1‑99 (MQ=255)                     ggCTACTGCTGGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCacta                                                                   >  1:1634198/1‑97 (MQ=255)                     ggCTACTGCTGGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAAcc                                                                   >  1:5170687/1‑99 (MQ=255)                     ggCTACTGCTGGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAAca                                                                   >  1:1354741/1‑99 (MQ=255)                     ggCTACTGCTGGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAACt                                                                   <  1:3166514/100‑1 (MQ=255)                     ggCTACTGCTGGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAACt                                                                   >  1:3648179/1‑100 (MQ=255)                     ggCTACTGCTGGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAACt                                                                   >  1:4201088/1‑100 (MQ=255)                     ggCTACTGCTGGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAACt                                                                   >  1:5348614/1‑100 (MQ=255)                     ggCTACTGCTGGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAACt                                                                   >  1:1293319/1‑100 (MQ=255)                     ccgatctgctgGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAACt                                                                   <  1:2563010/95‑1 (MQ=255)                     cccTACTGCTGGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCactg                                                                   >  1:3754811/3‑97 (MQ=255)                      caaccccctgGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAACTg                                                                  >  1:6339292/8‑100 (MQ=255)                       gatctgctgGTGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAACTgg                                                                 <  1:725418/97‑1 (MQ=255)                       gatctgctgGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAACTgg                                                                 <  1:2464427/97‑1 (MQ=255)                       atACTGCTGGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAATCAACtag                                                                 >  1:1217183/2‑98 (MQ=255)                        tACTGCTGGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAACTgag                                                                >  1:1671715/1‑98 (MQ=255)                        tACTGCTGGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACAAAcagag                                                                >  1:2594339/1‑96 (MQ=255)                        gaCTGCTGGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAACTgag                                                                >  1:655289/2‑98 (MQ=255)                         gctgctgGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAACTGGAt                                                               >  1:893989/2‑100 (MQ=255)                         aCTGCTGGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAATCaatcagat                                                               >  1:559422/1‑94 (MQ=255)                         aCTGCTGGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAAcagatc                                                               >  1:3882491/1‑95 (MQ=255)                         aCTGCTGGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAACTAGAt                                                               >  1:5988223/1‑100 (MQ=255)                           tgctgGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCACCTAGATcg                                                             >  1:3808763/1‑99 (MQ=255)                           atctgGCGAAGAGAAAGAACAAGAAGGTGCCAAACTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCCAcagatcgg                                                             >  1:2673306/3‑93 (MQ=255)                            gctgGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAAGCAACCAGATcgg                                                            >  1:3108833/1‑98 (MQ=255)                            gctgGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAACTAGATcgg                                                            >  1:1031176/1‑98 (MQ=255)                            gctgGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAACTAGATcgg                                                            >  1:3135283/1‑98 (MQ=255)                            gctgGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAACTAGATcgg                                                            >  1:1224262/1‑98 (MQ=255)                            gctgGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAACTAGATcgg                                                            >  1:2878359/1‑98 (MQ=255)                            gctgGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAACCAGATcgg                                                            >  1:3502437/1‑98 (MQ=255)                            gctgGCGAAGAGAAAGAACAAGAAGGGGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCactagatcgga                                                            >  1:4109760/1‑90 (MQ=255)                             ttgGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAACAAGATcgga                                                           >  1:1651706/2‑97 (MQ=255)                               ggCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAACTAGATcggaag                                                         >  1:6481440/1‑95 (MQ=255)                               ggCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACAAACTAGATcggaag                                                         >  1:5402335/1‑95 (MQ=255)                               ggCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAAACAACTAGATcggaag                                                         >  1:2556572/1‑95 (MQ=255)                               ggCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGGCGATGGCTTAATTATTGCTGAAACAACAAGATcggaag                                                         >  1:4781883/1‑95 (MQ=255)                                                     gTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAACTGGATCTGGAATTTATTATGCACATAAGGGGtt                                   <  1:4519392/100‑2 (MQ=255)                                                     gTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAACTGGATCTGGAATTTATTATGCACATAAGGGGta                                   <  1:2200678/100‑2 (MQ=255)                                                     gTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAACTGGATCTGGAATTTATTATGCACATAAGGGGTc                                   <  1:1319818/100‑1 (MQ=255)                                                      tGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAACTGGATCTGGAATTTATTATGCACATAAGGGGTcc                                  <  1:6335489/100‑2 (MQ=255)                                                      tGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAACTGGATCTGGAATTTATTATGCACATAAGGGGTca                                  >  1:4684042/1‑99 (MQ=255)                                                      tGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAACTGGATCTGGAATTTATTATGCACATAAGGGGTca                                  <  1:2075940/100‑2 (MQ=255)                                                      tGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAACTGGATCTGGAAGTTATTATGCACATAAGGGGTCt                                  <  1:1388533/100‑1 (MQ=255)                                                               attaGCTGATAAAGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAACTGGATCTGGAATTTATTATGCACATAAGGGGTCTATGTCATGt                         <  1:4841919/100‑1 (MQ=255)                                                                          caGGCTATTTAGACGATGTCGATGGCTTAATTATTGCTGAACCAACTGGATCTGGAATTTATTATGCACATAAGGGGTCTATGTCATGTAAAGTAACTGc              >  1:5406935/2‑100 (MQ=255)                                                                                   tggACGATGTCGATGGCTTAATTATTGCTGAACCAACTGGATCTGGAATTTATTATGCACATAAGGGGTCTATGTCATGTAAAGTAACTGCAACTGGTaa     >  1:3068207/3‑100 (MQ=255)                                                                                      aCGATGTCGATGGCTTAATTATTGCTGAACCAACTGGATCTGGAATTTATTATGCACATAAGGGGTCTATGTCATGTAAAGTAACTGCAACTGGTAAAgg  <  1:2458577/100‑2 (MQ=255)                                                                                                 |                                                                                         CATGGAACGATTAGATTACTGGCTACTGCTGGCGAAGAGAAAGAACAAGAAGGTGCCAAATTATTAGCTGATAAAGGCTATTTAGACGATGTCGATAGCTTAATTATTGCTGAACCAACTGGATCTGGAATTTATTATGCACATAAGGGGTCTATGTCATGTAAAGTAACTGCAACTGGTAAAGC  >  NC_007793/2131677‑2131861

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