Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	976,718	G→A	E216E (GAG→GAA)	mukE →	protein involved in chromosome partitioning

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	976,718	0	G	A	100.0%	136.3 / NA	42	E216E (GAG→GAA)	mukE	protein involved in chromosome partitioning
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/29);  total (13/29)

CCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCT  >  W3110S.gb/976677‑976772                                          |                                                       ccGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAAAATCAGCCAGa                                             <  1:852210/53‑1 (MQ=255) ccGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAAAATCAGCCAGa                                             <  1:1336510/53‑1 (MQ=255) ccGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAAAATCAGCCAGa                                             <  1:297735/53‑1 (MQ=255) ccGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAAAATCAGCCAGa                                             <  1:238524/53‑1 (MQ=255) ccGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAAAATCAGCCAGATAGCGga                                      <  1:3652930/60‑1 (MQ=255) ccGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAAAATCAGCCAGATAGCGga                                      <  1:725048/60‑1 (MQ=255) ccGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAAAATCAGCCAGATAGCGga                                      <  1:2858036/60‑1 (MQ=255) ccGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAAAATCAGCCAGATAGCGga                                      <  1:809336/60‑1 (MQ=255) ccGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAAAATCAGCCAGATAGCGga                                      <  1:2735645/60‑1 (MQ=255) ccGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAAAATCAGCCAGATAGCGga                                      <  1:2617644/60‑1 (MQ=255) ccGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAAAATCAGCCAGATAGCGga                                      <  1:4978/60‑1 (MQ=255) ccGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAAAATCAGCCAGATAGCGga                                      <  1:1234477/60‑1 (MQ=255)  cGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAAAATCAGCCAGATAGCGgaga                                    >  1:1809763/1‑61 (MQ=255)  cGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAAAATCAGCCAGATAGCGgaga                                    >  1:1392960/1‑61 (MQ=255)  cGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAAAATCAGCCAGATAGCGgaga                                    >  1:1011751/1‑61 (MQ=255)     ttGAAAATCATCTGCAACTCAACGATGAAACCGAAGAAAATCAGCCAGATAGCGGAGAGgaa                                <  1:3891940/62‑1 (MQ=255)     ttGAAAATCATCTGCAACTCAACGATGAAACCGAAGAAAATCAGCCAGATAGCGGAGAGgaa                                <  1:3784204/62‑1 (MQ=255)       gAAAATCATCTGCAACTCAACGATGAAACCGAAGAAAATCAGCCAGATAGCGGAGAGgaag                               >  1:2247812/1‑61 (MQ=255)       gAAAATCATCTGCAACTCAACGATGAAACCGAAGAAAATCAGCCAGATAGCGGAGAGgaag                               >  1:257995/1‑61 (MQ=255)       gAAAATCATCTGCAACTCAACGATGAAACCGAAGAAAATCAGCCAGATAGCGGAGAGgaag                               >  1:4524376/1‑61 (MQ=255)       gAAAATCATCTGCAACTCAACGATGAAACCGAAGAAAATCAGCCAGATAGCGGAGAGgaag                               >  1:1423608/1‑61 (MQ=255)       gAAAATCATCTGCAACTCAACGATGAAACCGAAGAAAATCAGCCAGATAGCGGAGAGgaag                               >  1:3712043/1‑61 (MQ=255)        aaaaTCATCTGCAACTCAACGATGAAACCGAAGAAAATCAGCCAGATAGCGGAGAGgaaga                              >  1:3818965/1‑61 (MQ=255)             catcTGCAACTCAACGATGAAACCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATg                        >  1:2880387/1‑62 (MQ=255)                   cAACTCAACGATGAAACCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGaa                  >  1:521199/1‑62 (MQ=255)                   cAACTCAACGATGAAACCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGaa                  >  1:3422100/1‑62 (MQ=255)                   cAACTCAACGATGAAACCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGa                   >  1:3109812/1‑61 (MQ=255)                             aTGAAACCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAAt        <  1:3558539/62‑1 (MQ=255)                             aTGAAACCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAAt        <  1:3446813/62‑1 (MQ=255)                             aTGAAACCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAAt        <  1:4309082/62‑1 (MQ=255)                             aTGAAACCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAAt        <  1:1917392/62‑1 (MQ=255)                                aaaCCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTc     <  1:3009477/62‑1 (MQ=255)                                aaaCCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTc     <  1:526027/62‑1 (MQ=255)                                aaaCCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTc     <  1:4329175/62‑1 (MQ=255)                                aaaCCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTc     <  1:1674783/62‑1 (MQ=255)                                  aCCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTaaa         <  1:723932/56‑1 (MQ=255)                                  aCCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTaaa         <  1:3337075/56‑1 (MQ=255)                                  aCCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTaaa         <  1:421352/56‑1 (MQ=255)                                  aCCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTaaa         <  1:3742892/56‑1 (MQ=255)                                  aCCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTaaa         <  1:3223494/56‑1 (MQ=255)                                  aCCGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTaaa         <  1:2758147/56‑1 (MQ=255)                                   ccGAAGAAAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCt  <  1:1209483/62‑1 (MQ=255)                                          |                                                       CCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCT  >  W3110S.gb/976677‑976772

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