Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,933,321	G→A	G428G (GGC→GGT)	edd ←	phosphogluconate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,933,321	0	G	A	100.0%	104.8 / NA	31	G428G (GGC→GGT)	edd	phosphogluconate dehydratase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/18);  total (13/18)

ACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCC‑TTCCCGCCAGTCCAGTTCACCATTATTCAGCC  >  NC_000913/1933184‑1933448                                                                                                                                          |                                                                                                                                 acaacaACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGAcc                                                                                                                                <  1:1124605/140‑1 (MQ=255) acaacaACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGAcc                                                                                                                                <  2:352539/140‑1 (MQ=255) acaacaACGACACAATCGCGATCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGAcc                                                                                                                                <  1:535817/140‑1 (MQ=255)                     gTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGACCCAGGTTACCGCTTAACACtt                                                                                                            >  1:1062158/1‑140 (MQ=255)                          gCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGACCCAGGTTACCGCTTAACACTTTTGTc                                                                                                       <  1:529028/140‑1 (MQ=255)                          acAACCCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCCCCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGACCCAGGTTACCGCTTAACACTTTTGTc                                                                                                       <  1:785672/139‑1 (MQ=255)                            aaaCCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGACCCAGGTTACCGCTTAACACTTTTGTccc                                                                                                     <  1:631749/140‑1 (MQ=255)                               cccGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGACCCAGGTTACCGCTTAACACTTTTGTCccacc                                                                                                  <  1:144844/140‑1 (MQ=255)                                ccGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGACCCAGGTTACCGCTTAACACTTTTGTCccacca                                                                                                 <  1:1056014/140‑1 (MQ=255)                                   cTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGACCCAGGTTACCGCTTAACACTTTTGTCCCACCatga                                                                                              >  2:505051/1‑140 (MQ=255)                                    ttCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGACCCAGGTTACCGCTTAACACTTTTGTCCCACCatgat                                                                                             <  2:929460/140‑1 (MQ=255)                                        aaGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGACCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATgaga                                                                                         >  1:240855/1‑140 (MQ=255)                                         aGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGACCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATgagag                                                                                        <  2:831611/140‑1 (MQ=255)                                              ggCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGACCCAGGTTACCGCTTAACACTTTTGTCCCACCatgat                                                                                             <  2:793714/130‑1 (MQ=255)                                              ggCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGACCCAGGTTACCGCTTAACACTTTTGTCCCACCatgat                                                                                             >  1:793714/1‑130 (MQ=255)                                              ggCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGACCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAgg                                                                                   <  2:557654/140‑1 (MQ=255)                                                          tGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGACCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCgaagga                                                                       >  2:799309/1‑140 (MQ=255)                                                             tGGCTTTTAAAAACCACCGCCGGCGCTTCAAACCCCTTGTTCTCAACCCGGCCCGCAGAGGTTTTTATAACCCCCAGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGc                                                                    <  2:857605/140‑1 (MQ=255)                                                               gCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGACCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCgaag                                                                         >  1:277762/1‑133 (MQ=255)                                                               gCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGACCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCgaag                                                                         <  2:277762/133‑1 (MQ=255)                                                               gCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGACCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGa                                                                  >  1:697038/1‑140 (MQ=255)                                                                         acaacCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGACCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCt                                                        <  1:657516/140‑1 (MQ=255)                                                                            acCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGACCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTc                                                     <  2:136084/140‑1 (MQ=255)                                                                               gCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGACCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAg                                                  >  2:497121/1‑140 (MQ=255)                                                                                           aTCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGACCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGc                                      <  1:505051/140‑1 (MQ=255)                                                                                           aTCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGACCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGc                                      <  2:356198/140‑1 (MQ=255)                                                                                           aTCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGACCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGc                                      <  1:497121/140‑1 (MQ=255)                                                                                                    ttCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGACCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCC‑TTCCCg                            >  2:1113686/1‑140 (MQ=255)                                                                                                        cAACCGGCACGGCAGAGGTTTTCATAACCGCACGACCCAGGTTACCGCTTAACGCTGTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCC‑TTCCCGccag                        >  2:689945/1‑140 (MQ=255)                                                                                                         aaCCGGCACGGCAGAGGTTTTCATAACCGCACGACCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCC‑TTCCCGccagt                       >  1:6363/1‑140 (MQ=255)                                                                                                          aCCGGCACGGCAGAGGTTTTCATAACCGCACGACCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCC‑TTCCCGccagtc                      >  1:264212/1‑140 (MQ=255)                                                                                                             ggcacggcaGAGGTTTTCATAACCGCACGACCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCC‑TTCCCGccagtccag                   >  1:313725/1‑140 (MQ=255)                                                                                                               cccggaaGAGTTTTTCATAACCGCACGCCCCAGGTTCCCGTTTAACTCTTTTGTCCCACCATGATGAGAAAAGGGTTGTTCGAAGGAGGCGATCACTTTGCTGTGGATTGTTTTTCCCCCCCC‑TTCCCGCCAGTCCAGtt                 <  2:697038/138‑1 (MQ=255)                                                                                                                        gtttttCAAACCCGCACCCCCTAGGTTACCGCTAAACTCTTTCGCCCCCCCATGATGATAGAATGGTCGTGTAAAGGAGGGCATCCCTTTGCCGTCTAGT‑ATTTTTCCCCCCCTTTCCCGCCAGTCCAGTTCACCattat        <  2:313725/138‑1 (MQ=255)                                                                                                                              agacagcGCACGACCCAGTTTCCCGCTTAACCTTTTTGCCCCCCCATGAGGAGAGAAAGTTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGTTTTTCCCCCCCC‑TCCCCGCCAGTCCAGTTCACCATTATTCAGcc  <  2:1057185/134‑1 (MQ=255)                                                                                                                                          |                                                                                                                                 ACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCC‑TTCCCGCCAGTCCAGTTCACCATTATTCAGCC  >  NC_000913/1933184‑1933448

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