Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,551,712	G→A	G26D (GGC→GAC)	uxuA →	D‑mannonate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,551,712	0	G	A	100.0%	169.5 / NA	46	G26D (GGC→GAC)	uxuA	D‑mannonate dehydratase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/17);  total (29/17)

ATGGAACAGACCTGGCGCTGGTACGGCCCAAACGATCCGGTTTCTTTAGCTGATGTCCGTCAGGCGGGCGCAACTGGCGTGGTTACCGCGCTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCA  >  NC_000913/4551636‑4551789                                                                             |                                                                              aTGGAACAGACCTGGCGCTGGTACGGCCCAAACGATCCGGTTTCTTTAGCTGATGTCCGTCAGGCGGGCGCAACTGACGt                                                                            <  1:915495/80‑1 (MQ=255)    gAACAGACCTGGCGCTGGTACGGCCCAAACGATCCGGTTTCTTTAGCTGATGTCCGTCAGGCGGGCGCAACTGACGTGGtt                                                                        <  2:915497/81‑1 (MQ=255)    gAACAGACCTGGCGCTGGTACGGCCCAAACGATCCGGTTTATTTAGCTGATGTCCGTCAGGCGGGCGAAACTGAAGTGGtt                                                                        <  2:915496/81‑1 (MQ=255)       cAGACCTGGCGCTGGTACGGCCCAAACGATCCGGTTTCTTTAGCTGATGTCCGTCAGGCGGGCGCAACTGACGTGGTTAcc                                                                     <  1:915498/81‑1 (MQ=255)       cAGACCTGGCGCTGGTACGGCCCAAACGATCCGGTTTCTTTAGCTGATGTCCGTCAGGCGGGCGCAACTGACGTGGTTAcc                                                                     >  1:915500/1‑81 (MQ=255)       cAGACCTGGCGCTGGTACGGCCCAAACGATCCGGTTTCTTTAGCTGATGTCCGTCAGGCGGGCGCAACTGACGTGGTTAcc                                                                     <  2:915499/81‑1 (MQ=255)       cAGACCTGGCGCTGGTACGGCCCAAACGATCCGGTTTCTTTAGCTGATGTCCGTCAGGCGGGCGCAACTGACGTGGTTAcc                                                                     >  1:915529/1‑81 (MQ=255)        aGACCTGGCGCTGGTACGGCCCAAACGATCCGGTTTCTTTAGCTGATGTCCGTCAGGCGGGCGCAACTGACGTGGTTACcg                                                                    >  1:915547/1‑81 (MQ=255)             tGGCGCTGGTACGGCCCAAACGATCCGGTTTCTTTAGCTGATGTCCGTCAGGCGGGCGCAACTGACGTGGTTACCGCGCTg                                                               >  2:915555/1‑81 (MQ=255)                cgcTGGTACGGCCCAAACGATCCGGTTTCTTTAGCTGATGTCCGTCAGGCGGGCGCAACTGACGTGGTTACCGCGCTGCAc                                                            <  2:915501/81‑1 (MQ=255)                      tACGGCCCAAACGATCCGGTTTCTTTAGCTGATGTCCGTCAGGCGGGCGCAACTAACGTGGTTACCGCGCTGCACCATATc                                                      <  1:915502/81‑1 (MQ=255)                       aCGGCCCAAACGATCCGGTTTCTTTAGCTGATGTCCGTCAGGCGGGCGCAACTGACGTGGTTACCGCGCTGCACCATATcc                                                     >  1:915540/1‑81 (MQ=255)                        cGGCCCAAACGATCCGGTTTCTTTAGCTGATGTCCGTCAGGCGGGCGCAACTGACGTGGTTACCGCGCTGCACCATATccc                                                    <  1:915503/81‑1 (MQ=255)                               aaCGATCCGGTTTCTTTAGCTGATGTCCGTCAGGCGGGCGCAACTGACGTGGTTACCGCGCTGCACCATATCCCGAACGGc                                             >  1:915549/1‑81 (MQ=255)                                 cGATCCGGTTTCTTTAGCTGATGTCCGTCAGGCGGGCGCAACTGACGTGGTTACCGCGCTGCACCATATCCCGAACGGCg                                            >  2:915568/1‑80 (MQ=255)                                  gATCCGGTTTCTTTAGCTGATGTCCGTCAGGCGGGCGCAACTGACGTGGTTACCGCGCTGCACCATATCCCGAACGGCg                                            >  1:915556/1‑79 (MQ=255)                                  gATCCGGTTTCTTTAGCTGATGTCCGTCAGGCGGGCGCAACTGACGTGGTTACCGCGCTGCACCATATCCCGAACGGCg                                            >  1:915531/1‑79 (MQ=255)                                  gATCCGGTTTCTTTAGCTGATGTCCGTCAGGCGGGCGCAACTGACGTGGTTACCGCGCTGCACCATATCCCGAACGGCg                                            >  2:915536/1‑79 (MQ=255)                                    tCCGGTTTCTTTAGCTGATGTCCGTCAGGCGGGCGCAACTGACGTGGTTACCGCGCTGCACCATATCCCGAACGGCGAAGt                                        <  1:915504/81‑1 (MQ=255)                                       ggTTTCTTTAGCTGATGTCCGTCAGGCGGGCGCAACTGACGTGGTTACCGCGCTGCACCATATCCCGAACGGCGAAGTATg                                     >  2:915543/1‑81 (MQ=255)                                        gTTTCTTTAGCTGATGTCCGTCAGGCGGGCGCAACTGACGTGGTTACCGCGCTGCACCATATCCCGAACGGCGAAGTATgg                                    >  1:915548/1‑81 (MQ=255)                                        gTTTCTTTAGCTGATGTCCGTCAGGCGGGCGCAACTGACGTGGTTACCGCGCTGCACCATATCCCGAACGGCGAAGTATgg                                    >  1:915539/1‑81 (MQ=255)                                            cTTTAGCTGATGTCCGTCAGGCGGGCGCAACTGACGTGGTTACCGCGCTGCACCATATCCCGAACGGCGAAGTATGGTCCg                                >  2:915535/1‑81 (MQ=255)                                               tAGCTGATGTCCGTCAGGCGGGCGCAACTGACGTGGTTACCGCGCTGCACCATATCCCGAACGGCGAAGTATGGTCCGtt                              <  2:915505/80‑1 (MQ=255)                                               tAGCTGATGTCCGTCAGGCGGGCGCAACTGACGTGGTTACCGCGCTGCACCATATCCCGAACGGCGAAGTATGGTCCGtt                              <  1:915506/80‑1 (MQ=255)                                               tAGCTGATGTCCGTCAGGCGGGCGCAACTGACGTGGTTACCGCGCTGCACCATATCCCGAACGGCGAAGTATGGTCCGt                               <  2:915507/79‑1 (MQ=255)                                               tAGCTGATGTCCGTCAGGCGGGCGCAACTGACGTGGTTACCGCGCTGCACCATATCCCGAACGGCGAAGTATGGTCCGt                               <  1:915508/79‑1 (MQ=255)                                                aGCTGATGTCCGTCAGGCGGGCGCAACTGACGTGGTTACCGCGCTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTg                             >  1:915545/1‑80 (MQ=255)                                                  cTGATGTCCGTCAGGCGGGCGCAACTGACGTGGTTACCGCGCTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAag                          <  2:915509/81‑1 (MQ=255)                                                    gATGTCCGTCAGGCGGGCGCAACTGACGTGGTTACCGCGCTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGaa                           >  1:915553/1‑78 (MQ=255)                                                     aTGTCCGTCAGGCGGGCGCAACTGACGTGGTTACCGCGCTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAagag                        <  1:915510/80‑1 (MQ=255)                                                        tCCGTCAGGCGGGCGCAACTGACGTGGTTACCGCGCTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATcc                    >  1:915565/1‑81 (MQ=255)                                                         ccGTCAGGCGGGCGCAACTGACGTGGTTACCGCGCTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCt                   >  2:915558/1‑81 (MQ=255)                                                           gTCAGGCGGGCGCAACTGACGTGGTTACCGCGCTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTc                  >  2:915557/1‑80 (MQ=255)                                                           gTCAGGCGGGCGCAACTGACGTGGTTACCGCGCTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTc                  >  1:915564/1‑80 (MQ=255)                                                           gTCAGGCGGGCGCAACTGACGTGGTTACCGCGCTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTc                  >  1:915551/1‑80 (MQ=255)                                                             cAGGCGGGCGCAACTGACGTGGTTACCGCGCTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCaaa               >  2:915571/1‑81 (MQ=255)                                                               ggcgggcgCAACTGACGTGGTTACCGCGCTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACg             <  2:915511/81‑1 (MQ=255)                                                                   ggcgCAACTGACGTGGTTACCGCGCTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAg         <  2:915512/81‑1 (MQ=255)                                                                     cgcAACTGACGTGGTTACCGCGCTGCACCATATCCCGAACCGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGc       >  1:915570/1‑81 (MQ=255)                                                                      gcAACTGACGTGGTTACCGCGCTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGc       >  1:915550/1‑80 (MQ=255)                                                                       cAACTGACGTGGTTACCGCGCTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGa     >  1:915554/1‑81 (MQ=255)                                                                        aaCTGACGTGGTTACCGCGCTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGAt    >  2:915544/1‑81 (MQ=255)                                                                        aaCTGACGTGGTTACCGCGCTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGAt    >  1:915546/1‑81 (MQ=255)                                                                        aaCTGACGTGGTTACCGCGCTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGAt    >  2:915576/1‑81 (MQ=255)                                                                          cTGACGTGGTTACCGCGCTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCa  >  1:915584/1‑81 (MQ=255)                                                                             |                                                                              ATGGAACAGACCTGGCGCTGGTACGGCCCAAACGATCCGGTTTCTTTAGCTGATGTCCGTCAGGCGGGCGCAACTGGCGTGGTTACCGCGCTGCACCATATCCCGAACGGCGAAGTATGGTCCGTTGAAGAGATCCTCAAACGCAAGGCGATCA  >  NC_000913/4551636‑4551789

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