Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,791,162	G→A	G555E (GGG→GAG)	SAUSA300_RS14330 →	PTS mannose transporter subunit IIABC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,791,162	0	G	A	100.0%	68.7 / NA	27	G555E (GGG→GAG)	SAUSA300_RS14330	PTS mannose transporter subunit IIABC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/2);  total (25/2)

ATTCATAAACTTAAAATTTGCCGTTATGTAGAGCACACATCGCATTTAAAAAATGCAATACTAGAACGTGAAATGGGGTCTACGACAGCTATTGGCATGAACGTGGCCATTCCACATGCAAAGTCTGATGTTGTAAAACAACCGATTGTAGCTGTTATGAAAAATAATC  >  NC_007793/2791086‑2791254                                                                             |                                                                                             aTTCATAAACTTAAAATTTGCCGTTATGTAGAGCACACATCGCATTTAAAAAATGCAATACTAGAACGTGAAATGGAGTCTACGACAGCTATTGGCATgg                                                                       <  1:5709410/100‑2 (MQ=255)       aaaCTTAAAATTTGCCGTTATGTAGAGCACACATCGCATTTAAAAAATGCAATACTAGAACGTGAAATGGAGTCTACGACAGCTATTGGCATGAACGTgg                                                                 >  1:2587189/1‑100 (MQ=255)                  ttGCCGTTATGTAGAGCACACATCGCATTTAAAAAATGCAATACTAGAACGTGAAATGGAGTCTACGACAGCTATTGGCATGAACGTGGCCATTCCACAt                                                      >  1:5353376/1‑100 (MQ=255)                  ttGCCGTTATGTAGAGCACACATCGCATTTAAAAAATGCAATACTAGAACGTGAAATGGAGTCTACGACAGCTATTGGCATGAACGTGGCCATTCCACAt                                                      >  1:1686895/1‑100 (MQ=255)                  ttGCCGTTATGTAGAGCACACATCGCATTTAAAAAATGCAATACTAGAACGTGAAATGGAGTCTACGACAGCTATTGGCATGAACGTGGCCATTCCACAt                                                      >  1:2078406/1‑100 (MQ=255)                  ttGCCGTTATGTAGAGCACACATCGCATTTAAAAAATGCAATACTAGAACGTGAAATGGAGTCTACGACAGCTATTGGCATGAACGTGGCCATTCCACAt                                                      >  1:3475921/1‑100 (MQ=255)                    gCCGTTATGTAGAGCACACATCGCATTTAAAAAATGCAATACTAGAACGTGAAATGGAGTCTACGACAGCTATTGGCATGAACGTAGCCATTCCACATGc                                                    >  1:4867830/1‑100 (MQ=255)                      cGTTATGTAGAGCACACATCGCATTTAAAAAATGCAATACTAGAACGTGAAATGGAGTCTACGACAGCTATTGGCATGAACGTGGCCATTCCACATGCaa                                                  <  1:1118749/100‑1 (MQ=255)                       cctATGTAGAGCACACATCGCATTTAAAAAATGCAATACTAGAACGTGAAATGGAGTCTACGACAGCTATTGGCATGAACGTGGCCATTCCACATGCaaa                                                 >  1:5477432/3‑100 (MQ=255)                                 ccgcacaTCGCATTTAAAAAATGCAATACTAGAACGTGAAATGGAGTCTACGACAGCTATTGGCATGAACGTGGCCATTCCACATGCAAAGTCTGAtgct                                       >  1:2424570/4‑98 (MQ=255)                                 acccacaTCGCATTTAAAAAATGCAATACTAGAACGTGACATGGAGTCTACGACAGCTATTGGCATTAACATGGCCATTCCACATGCAAAGTCTGAtgtt                                       >  1:4768186/4‑100 (MQ=255)                                  cccacaTCGCATTTAAAAAATGCAATACTAGAACGTGAAATGGAGTCTACGACAGCTATTGGCATGAACGTGGCCATTCCACATGCAAAGTCTGAtgttg                                      >  1:1377455/3‑100 (MQ=255)                                  cccacaTCGCATTTAAAAAATGCAATACTAGAACGTGAAATGGAGTCTACGACAGCTATTGGCATGAACGTGGCCATTCCACATGCAAAGTCTGAtgttg                                      >  1:2853809/3‑100 (MQ=255)                                  cccacaTCGCATTTAAAAAATGCAATACTAGAACGTGAAATGGAGTCTACGACAGCTATTGGCATGAACGTGGCCATTCCACATGCAAAGTCTGAtgttg                                      >  1:1763855/3‑100 (MQ=255)                                  cccacaTCGCATTTAAAAAATGCAATACTAGAACGTGAAATGGAGTCTACGACAGCTATTGGCATGAACGTGGCCATTCCACATGCAAAGCCTGAtgttg                                      >  1:731727/3‑100 (MQ=255)                                   gcacaTCGCATTTAAAAAATGCAATACTAGAACGTGAAATGGAGTCTACGACAGCTATTGGCATGAACGTGGCCATTCCACATGCAAAGTCTGAtgttgt                                     >  1:4773802/2‑100 (MQ=255)                                   gcacaTCGCATTTAAAAAATGCAATACTAGAACGTGAAATGGAGTCTACGACAGCTATTGGCATGAACGTGGCCATTCCACATGCAAAGTCTGAtgttgt                                     >  1:1855847/2‑100 (MQ=255)                                   ccacaTCGCATTTAAAAAATGCAATACTAGAACGTGAAATGGAGTCTACGACAGCTATTGGCATGAACGTGGCCATTCCACATGCAAAGTCTGAtgttgt                                     >  1:4408154/2‑100 (MQ=255)                                     acaTCGCATTTAAAAAATGCAATACTAGAACGTGAAATGGAGTCTACGACAGCTATTGGCATGAACGTGGCCATTCCACATGCAAAGTCTGATGTTGTaa                                   >  1:2790882/1‑100 (MQ=255)                                                    taTGCAATACTAGAACGTGAAATGGAGTCTACGACAGCTATTGGCATGAACGTGGCCATTCCACATGCAAAGTCTGATGTTGTAAAACAACCGATTGTAg                    >  1:4378655/2‑100 (MQ=255)                                                    aaTGCAATACTAGAACGTGAAATGGAGTCTACGACAGCTATTGGCATGAACGTGGCCATTCCACATGCAAAGTCTGATGTTGTAAAACAACCGATTGTAg                    >  1:957773/1‑100 (MQ=255)                                                       gCAATACTAGAACGTGAAATGGAGTCTACGACAGCTATTGGCATGAACGTGGCCATTCCACATGCAAAGTCTGATGTTGTAAAACAACCGATTGTAGCTg                 >  1:2809602/1‑100 (MQ=255)                                                            ccTAGAACGTGAAATGGAGTCTACGACAGCTATTGGCATGAACGTGGCCATTCCACATGCAAAGTCTGATGTTGTAAAACAACCGATTGTAGCTGTTATg            >  1:6127116/2‑100 (MQ=255)                                                               cgAACGAGAAATGGAGTCTACGACAGCTATTGGCATGAACGTGGCCATTCCACATGCAAAGTCTGATGTTGTAAAACACCCGACTGAAGCTGTTATGaaa         >  1:6434781/2‑100 (MQ=255)                                                               aGAACGTGAAATGGAGTCTACGACAGCTATTGGCATGAACGTGGCCATTCCACATGCAAAGTCTGATGTTGTAAAACAACCGATTGTAGCTGTTATGaaa         >  1:3925061/1‑100 (MQ=255)                                                                     tGAAATGGAGTCTACGACAGCTATTGGCATGAACGTGGCCATTCCACATGCAAAGTCTGATGTTGTAAAACAACCGATTGTAGCTGTTATGAAAaataat   >  1:5192138/1‑100 (MQ=255)                                                                      gAAATGGAGTCTACGACAGCTATTGGCATGAACGTGGCCATTCCACATGCAAAGTCTGATGTTGTAAAACAACCGATTGTAGCTGTTATGAAAAATAATc  >  1:1643132/1‑100 (MQ=255)                                                                             |                                                                                             ATTCATAAACTTAAAATTTGCCGTTATGTAGAGCACACATCGCATTTAAAAAATGCAATACTAGAACGTGAAATGGGGTCTACGACAGCTATTGGCATGAACGTGGCCATTCCACATGCAAAGTCTGATGTTGTAAAACAACCGATTGTAGCTGTTATGAAAAATAATC  >  NC_007793/2791086‑2791254

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