Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,983,085	(A)6→5	coding (174/762 nt)	kduD ←	putative 2‑keto‑3‑deoxy‑D‑gluconate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,983,080	0	A	.	100.0%	133.4 / NA	32	coding (179/762 nt)	kduD	putative 2‑keto‑3‑deoxy‑D‑gluconate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (21/11);  total (21/11)

AATATGACCAAACTCCGCTACCGCGCGATCCAGCAGTGCTGGAATACCATCAATCTTTCGCAGATCGGCGGTCAGGCTTAAAAAACGACGCCCCAGCGCTGTGACCTGCTCGATGGTTTCAGTCGGTTCAACGATGTTAATGCCAACAATGTCACA  >  NC_000913/2983001‑2983156                                                                                |                                                                             aaTATGACCAAACTCCGCTACCGCGCGATCCAGCAGTGCTGGAATACCATCAATCTTTCGCAGATCGGCGGTCAGGCTTaa                                                                             <  2:488381/81‑1 (MQ=255)     tGACCAAACTCCGCTACCGCGCGATCCAGCAGTGCTGGAATACCATCAATCTTTCGCAGATCGGCGGTCAGGCTTaaaaac                                                                         <  1:488383/81‑2 (MQ=255)             cTCCGCTACCGCGCGATCCAGCAGTGCTGGAATACCATCAATCTTTCGCAGATCGGCGGTCAGGCTT‑AAAAACGACGccc                                                                 >  1:488410/1‑80 (MQ=255)                    aCCGCGCGATCCAGCAGTGCTGGAATACCATCAATCTTTCGCAGATCGGCGGTCAGGCTT‑AAAAACGACGCCCCAGCGCtg                                                         >  2:488423/1‑81 (MQ=255)                     ccGCGCGATCCAGCAGTGCTGGAATACCATCAATCTTTCGCAGATCGGCGGTCAGGCTT‑AAAAACGACGCCCCAGCGCtgt                                                        >  1:488413/1‑81 (MQ=255)                       gcgcgATCCAGCAGTGCTGGAATACCATCAATCTTTCGCAGATCGGCGGTCAGGCTT‑AAAAACGACGCCCCAGCGCTGTGa                                                      >  1:488439/1‑81 (MQ=255)                         gcgATCCAGCAGTGCTGGAATACCATCAATCTTTCGCAGATCGGCGGTCAGGCTT‑AAAAACGACGCCCCAGCGCTGTGAcc                                                    >  1:488414/1‑81 (MQ=255)                         gcgATCCAGCAGTGCTGGAATACCATCAATCTTTCGCAGATCGGCGGTCAGGCTT‑AAAAACGACGCCCCAGCGCTGTGAcc                                                    >  1:488418/1‑81 (MQ=255)                           gATCCAGCAGTGCTGGAATACCATCAATCTTTCGCAGATCGGCGGTCAGGCTT‑AAAAACGACGCCCCAGCGCTGTGACCTg                                                  >  2:488430/1‑81 (MQ=255)                            aTCCAGCAGTGCTGGAATACCATCAATCTTTCGCAGATCGGCGGTCAGGCTT‑AAAAACGACGCCCCAGCGCTGTGACCTGc                                                 >  1:488406/1‑81 (MQ=255)                              ccAGCAGTGCTGGAATACCATCAATCTTTCGCAGATCGGCGGTCAGGCTT‑AAAAACGACGCCCCAGCGCTGTGACCTGCTc                                               >  2:488416/1‑81 (MQ=255)                                 gcagTGCTGGAATACCATCAATCTTTCGCAGATCGGCGGTCAGGCTT‑AAAAACGACGCCCCAGCGCTGTGACCTGCTCGAt                                            >  2:488415/1‑81 (MQ=255)                                 gcagTGCTGGAATACCATCAATCTTTCGCAGATCGGCGGTCAGGCTT‑AAAAACGACGCCCCAGCGCTGTGACCTGCTCGAt                                            <  2:488383/81‑1 (MQ=255)                                  cagTGCTGGAATACCATCAATCTTTCGCAGATCGGCGGTCAGGCTT‑AAAAACGACGCCCCAGCGCTGTGACCTGCTCGATg                                           <  2:488384/81‑1 (MQ=255)                                   agTGCTGGAATACCATCAATCTTTCGCAGATCGGCGGTCAGGCTT‑AAAAACGACGCCCCAGCGCTGTGACCTGCTCGATgg                                          >  2:488434/1‑81 (MQ=255)                                     tGCTGGAATACCATCAATCTTTCGCAGATCGGCGGTCAGGCTT‑AAAAACGACGCCCCAGCGCTGTGACCTGCTCGATGGtt                                        <  2:488385/81‑1 (MQ=255)                                      gCTGGAATACCATCAATCTTTCGCAGATCGGCGGTCAGGCTT‑AAAAACGACGCCCCAGCGCTGTGACCTGCTCGATGGttt                                       >  2:488427/1‑81 (MQ=255)                                       cTGGAATACCATCAATCTTTCGCAGATCGGCGGTCAGGCTT‑AAAAACGACGCCCCAGCGCTGTGACCTGCTCGATGGTTTc                                      <  1:488387/81‑1 (MQ=255)                                           aaTACCATCAATCTTTCGCAGATCGGCGGTCAGGCTT‑AAAAACGACGCCCCAGCGCTGTGACCTGCTCGATGGTTTCAGTc                                  >  1:488407/1‑81 (MQ=255)                                                 aTCAATCTTTCGCAGATCGGCGGTCAGGCTT‑AAAAACGACGCCCCAGCGCTGTGACCTGCTCGATGGTTTCAGTCGGTTCa                            <  1:488388/81‑1 (MQ=255)                                                    aaTCTTTCGCAGATCGGCGGTCAGGCTT‑AAAAACGACGCCCCAGCGCTGTGACCTGCTCGATGGTTTCAGTCGGTTCAAc                          >  1:488420/1‑80 (MQ=255)                                                      tCTTTCGCAGATCGGCGGTCAGGCTT‑AAAAACGACGCCCCAGCGCTGTGACCTGCTCGATGGTTTCAGTCGGTTCAACGAt                       >  2:488429/1‑81 (MQ=255)                                                         ttCGCAGATCGGCGGTCAGGCTT‑AAAAACGACGCCCCAGCGCTGTGACCTGCTCGATGGTTTCAGTCGGTTCAACGATGt                     >  2:488426/1‑80 (MQ=255)                                                          tCGCAGATCGGCGGTCAGGCTT‑AAAAACGACGCCCCAGCGCTGTGACCTGCTCGATGGTTTCAGTCGGTTCAACGAt                       <  2:488388/77‑1 (MQ=255)                                                             cAGATCGGCGGTCAGGCTT‑AAAAACGACGCCCCAGCGCTGTGACCTGCTCGATGGTTTCAGTCGGTTCAACGATGTTAATg                >  1:488422/1‑81 (MQ=255)                                                                aTCGGCGGTCAGGCTT‑AAAAACGACGCCCCAGCGCTGTGACCTGCTCGATGGTTTCAGTCGGTTCAACGATGTTAATGcc              >  1:488421/1‑80 (MQ=255)                                                                   ggcggTCAGGCTT‑AAAAACGACGCCCCAGCGCTGTGACCTGCTCGATGGTTTCAGTCGGTTCAACGATGTTAATGCcaac           <  2:488389/80‑1 (MQ=255)                                                                   ggcggTCAGGCTT‑AAAAACGACGCCCCAGCGCTGTGACCTGCTCGATGGTTTCAGTCGGTTCAACGATGTTAATGCcaac           >  2:488425/1‑80 (MQ=255)                                                                     cggTCAGGCTT‑AAAAACGACGCCCCAGCGCTGTGACCTGCTCGATGGTTTCAGTCGGTTCAACGATGTTAATGCCAACAAt        >  2:488432/1‑81 (MQ=255)                                                                       gTCAGGCTT‑AAAAACGACGCCCCAGCGCTGTGACCTGCTCGATGGTTTCAGTCGGTTCAACGATGTTAATGCCAACAATGt      <  1:725804/81‑1 (MQ=255)                                                                        tCAGGCTT‑AAAAACGACGCCCCAGCGCTGTGACCTGCTCGATGGTTTCAGTCGGTTCAACGATGTTAATGCCAACAATGTc     <  1:488391/81‑1 (MQ=255)                                                                          aGGCTT‑AAAAACGACGCCCCAGCGCTGTGACCTGCTCGATGGTTTCAGTCGGTTCAACGATGTTAATGCCAACAATGTcac   <  1:488392/81‑1 (MQ=255)                                                                           ggCTT‑AAAAACGACGCCCCAGCGCTGTGACCTGCTCGATGGTTTCAGTCGGTTCAACGATGTTAATGCCAACAATGTcaca  <  2:488392/81‑1 (MQ=255)                                                                           ggCTT‑AAAAACGACGCCCCAGCGCTGTGACCTGCTCGATGGTTTCAGTCGGTTCAACGATGTTAATGCCAACAATGTca    >  1:488423/1‑79 (MQ=255)                                                                                |                                                                             AATATGACCAAACTCCGCTACCGCGCGATCCAGCAGTGCTGGAATACCATCAATCTTTCGCAGATCGGCGGTCAGGCTTAAAAAACGACGCCCCAGCGCTGTGACCTGCTCGATGGTTTCAGTCGGTTCAACGATGTTAATGCCAACAATGTCACA  >  NC_000913/2983001‑2983156

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