Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,567,041	A→C	D331E (GAT→GAG)	glgA ←	glycogen synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,567,041	0	A	C	100.0%	104.2 / NA	27	D331E (GAT→GAG)	glgA	glycogen synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (11/16);  total (11/16)

CTGGGCACCAGAATGACGTCCGCGCCGCCCATAATGCGATGCGAAAATGCTTCGTGATAGCCAATCTGAACGCCCACCTGACCGGGGTATTCCGCTGCCGCCGCAAGGAAACCTTCCTGCAGCACCGGATCGCCCGCGCCGAGTAGCGCCAGCTGCCCGCCCTGCTCCAGAAGACCCGGTAAGGCTTCCAGCACCAGGTCGAGACCTTTCTGGCTGGTCAGACGGCTCACCACTGCAAAAAGCGGCACTTTATCGTC  >  NC_000913/3566913‑3567169                                                                                                                                 |                                                                                                                                 cTGGGCACCAGAATGACGTCCGCGCCGCCCATAATGCGATGCGAAAATGCTTCGTGATAGCCAATCTGAACGCCCACCTGACCGGGGTATTCCGCTGCCGCCGCAAGGAAACCTTCCTGCAGCACCGGCTCGCCCGCGcc                                                                                                                       >  2:140660/1‑140 (MQ=255)         cAGAATGACGTCCGCGCCGCCCATAATGCGATGCGAAAATGCTTCGTGATAGCCAATCTGAACGCCCACCTGACCGGGGTATTCCGCTGCCGCCGCAAGGAAACCTTCCTGCAGCACCGGCTCGCCCGCGCCGAGTAgcg                                                                                                               >  2:741600/1‑140 (MQ=255)                 cGTCCGCGCCGCCCATAATGCGATGCGAAAATGCTTCGTGATAGCCAATCTGAACGCCCACCTGACCGGGGTATTCCGCTGCCGCCGCAAGGAAACCTTCCTGCAGCACCGGCTCGCCCGCGCCGAGTAGCGCCAGCTgc                                                                                                       <  2:790918/140‑1 (MQ=255)                    ccGCGCCGCCCATAATGCGATGCGAAAATGCTTCGTGATAGCCAATCTGAACGCCCACCTGACCGGGGTATTCCGCTGCCGCCGCAAGGAAACCTTCCTGCAGCACCGGCTCGCCCGCGCCGAGTAGCGCCAGCTgcccg                                                                                                    <  1:28652/140‑1 (MQ=255)                      gcgcCGCCCATAATGCGATGCGAAAATGCTTCGTGATAGCCAATCTGAACGCCCACCTGACCGGGGTATTCCGCTGCCGCCGCAAGGAAACCTTCCTGCAGCACCGGCTCGCCCGCGCCGAGTAGCGCCAGCTgcccgcc                                                                                                  <  2:406149/140‑1 (MQ=255)                        gccgccCATAATGCGATGCGAAAATGCTTCGTGATAGCCAATCTGAACGCCCACCTGACCGGGGTATTCCGCTGCCGCCGCAAGGAAACCTTCCTGCAGCACCGGCTCGCCCGCGCCGAGTAGCGCCAGCTGCCCGCCCt                                                                                                >  2:724288/1‑140 (MQ=255)                         ccgccCATAATGCGATGCGAAAATGCTTCGTGATAGCCAATCTGAACGCCCACCTGACCGGGGTATTCCGCTGCCGCCGCAAGGAAACCTTCCTGCAGCACCGGCTCGCCCGCGCCGAGTAGCGCCAGCTGCCCGCCCTg                                                                                               <  1:724288/140‑1 (MQ=255)                                     cgatgcgaAAATGCTTCGTGATAGCCAATCTGAACGCCCACCTGACCGGGGTATTCCGCTGCCGCCGCAAGGAAACCTTCCTGCAGCACCGGCTCGCCCGCGCCGAGTAGCGCCAGCTGCCCGCCCTGCTCCAGAAGAcc                                                                                   <  1:279931/140‑1 (MQ=255)                                                 gCTTCGTGATAGCCAATCTGAACGCCCACCTGACCGGGGTATTCCGCTGCCGCCGCAAGGAAACCTTCCTGCAGCACCGGCTCGCCCGCGCCGAGTAGCGCCAGCTGCCCGCCCTGCTCCAGAAGACCCGGTAAGGCTTc                                                                       <  1:758318/140‑1 (MQ=255)                                                    tCGTGATAGCCAATCTGAACGCCCACCTGACCGGGGTATTCCGCTGCCGCCGCAAGGAAACCTTCCTGCAGCACCGGCTCGCCCGCGCCGAGTAGCGCCAGCTGCCCGCCCTGCTCCAGAAGACCCGGTAAGGCTTCCAg                                                                    >  2:385538/1‑140 (MQ=255)                                                        gATAGCCAATCTGAACGCCCACCTGACCGGGGTATTCCGCTGCCGCCGCAAGGAAACCTTCCTGCAGCACCGGCTCGCCCGCGCCGAGTAGCGCCAGCTGCCCGCCCTGCTCCAGAAGACCCGGTAAGGCTTCCAGCAcc                                                                <  1:875289/140‑1 (MQ=255)                                                                 tCTGAACGCCCACCTGACCGGGGTATTCCGCTGCCGCCGCAAGGAAACCTTCCTGCAGCACCGGCTCGCCCGCGCCGAGTAGCGCCAGCTGCCCGCCCTGCTCCAGAAGACCCGGTAAGGCTTCCAGCACCAGGTCgaga                                                       >  2:147545/1‑140 (MQ=255)                                                                          ccACCTGACCGGGGTATTCCGCTGCCGCCGCAAGGAAACCTTCCTGCAGCACCGGCTCGCCCGCGCCGAGTAGCGCCAGCTGCCCGCCCTGCTCCAGAAGACCCGGTAAGGCTTCCAGCACCAGGTCGAGACCTTTctgg                                              <  2:231690/140‑1 (MQ=255)                                                                           cACCTGACCGGGGTATTCCGCTGCCGCCGCAAGGAAACCTTCCTGCAGCACCGGCTCGCCCGCGCCGAGTAGCGCCAGCTGCCCGCCCTGCTCCAGAAGACCCGGTAAGGCTTCCAGCACCAGGTCGAGACCTTTctggc                                             <  1:416382/140‑1 (MQ=255)                                                                              cTGACCGGGGTATTCCGCTGCCGCCGCAAGGAAACCTTCCTGCAGCACCGGCTCGCCCGCGCCGAGTAGCGCCAGCTGCCCGCCCTGCTCCAGAAGACCCGGTAAGGCTTCCAGCACCAGGTCGAGACCTTTctggctgg                                          <  1:673947/140‑1 (MQ=255)                                                                                   cGGGGTATTCCGCTGCCGCCGCAAGGAAACCTTCCTGCAGCACCGGCTCGCCCGCGCCGAGTAGCGCCAGCTGCCCGCCCTGCTCCAGAAGACCCGGTAAGGCTTCCAGCACCAGGTCGAGACCTTTCTGGCTGGTCAGa                                     <  2:185518/140‑1 (MQ=255)                                                                                   cGGGGTATTCCGCTGCCGCCGCAAGGAAACCTTCCTGCAGCACCGGCTCGCCCGCGCCGAGTAGCGCCAGCTGCCCGCCCTGCTCCAGAAGACCCGGTAAGGCTTCCAGCACCAGGTCGAGACCTTTCTGGCTGGTCAGa                                     <  2:82916/140‑1 (MQ=255)                                                                                      ggTATTCCGCTGCCGCCGCAAGGAAACCTTCCTGCAGCACCGGCTCGCCCGCGCCGAGTAGCGCCAGCTGCCCGCCCTGCTCCAGAAGACCCGGTAAGGCTTCCAGCACCAGGTCGAGACCTTTCTGGCTGGTCAGACgg                                  >  2:843960/1‑140 (MQ=255)                                                                                                tGCCGCCGCAAGGAAACCTTCCTGCAGCACCGGCTCGCCCGCGCCGAGTAGCGCCAGCTGCCCGCCCTGCTCCAGAAGACCCGGTAAGGCTTCCAGCACCAGGTCGAGACCTTTCTGGCTGGTCAGACGGCTCACCACTg                        >  2:627766/1‑140 (MQ=255)                                                                                                tGCCGCCGCAAGGAAACCTTCCTGCAGCACCGGCTCGCCCGCGCCGAGTAGCGCCAGCTGCCCGCCCTGCTCCAGAAGACCCGGTAAGGCTTCCAGCACCAGGTCGAGACCTTTCTGGCTGGTCAGACGGCTCACCACTg                        >  2:180667/1‑140 (MQ=255)                                                                                                tGCCGCCGCAAGGAAACCTTCCTGCAGCACCGACTCGCCCGCGCCGAGTAGCGCCAGCTGCCCGCCCTGCTCCAGAAGACCCGGTAAGGCTTCCAGCACCAGGTCGAGACCTTTCTGGCTGGTCAGACGGCTCACCACTg                        >  2:81178/1‑140 (MQ=255)                                                                                                  ccgccgCAAGGAAACCTTCCTGCAGCACCGGCTCGCCCGCGCCGAGTAGCGCCAGCTGCCCGCCCTGCTCCAGAAGACCCGGTAAGGCTTCCAGCACCAGGTCGAGACCTTTCTGGCTGGTCAGACGGCTCACCACTGCa                      >  2:46882/1‑140 (MQ=255)                                                                                                        cAAGGAAACCTTCCTGCAGCACCGGCTCGCCCGAGCCGAGTAGCGCCAGCTGCCCGCCCTGCTCCAGAAGACCCGGTAAGGCTTCCAGCACCAGGTCGAGACCTTTCTGGCTGGTCAGACGGCTCACCACTGCAAAAAGc                <  2:258819/140‑1 (MQ=255)                                                                                                         aaGGAAACCTTCCTGCAGCACCGGCTCGCCCGCGCCGAGTAGCGCCAGCTGCCCGCCCTGCTCCAGAAGACCCGGTAAGGCTTCCAGCACCAGGTCGAGACCTTTCTGGCTGGTCAGACGGCTCACCACTGCAAAAAGCg               <  1:147545/140‑1 (MQ=255)                                                                                                                    ccTGCAGCACCGGCTCGCCCGCGCCGAGTAGCGCCAGCTGCCCGCCCTGCTCCAGAAGACCCGGTAAGGCTTCCAGCACCAGGTCGAGACCTTTCTGGCTGGTCAGACGGCTCACCACTGCAAAAAGCGGCACTTTATCg    >  2:69995/1‑140 (MQ=255)                                                                                                                      tGCAGCACCGGCTCGCCCGCGCCGAGTAGCGCCAGCTGCCCGCCCTGCTCCAGAAGACCCGGTAAGGCTTCCAGCACCAGGTCGAGACCTTTCTGGCTGGTCAGACGGCTCACCACTGCAAAAAGCGGCACTTTATCGTc  <  1:734382/140‑1 (MQ=255)                                                                                                                      tGCAGCACCGGCTCGCCCGCGCCGAGTAGCGCCAGCTGCCCGCCCTGCTCCAGAAGACCCGGTAAGGCTTCCAGCACCAGGTCGAGACCTTTCTGGCTGGTCAGACGGCTCACCACTGCAAAAAGCGGCACTTTATCGTc  <  1:627766/140‑1 (MQ=255)                                                                                                                                 |                                                                                                                                 CTGGGCACCAGAATGACGTCCGCGCCGCCCATAATGCGATGCGAAAATGCTTCGTGATAGCCAATCTGAACGCCCACCTGACCGGGGTATTCCGCTGCCGCCGCAAGGAAACCTTCCTGCAGCACCGGATCGCCCGCGCCGAGTAGCGCCAGCTGCCCGCCCTGCTCCAGAAGACCCGGTAAGGCTTCCAGCACCAGGTCGAGACCTTTCTGGCTGGTCAGACGGCTCACCACTGCAAAAAGCGGCACTTTATCGTC  >  NC_000913/3566913‑3567169

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