Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	4,101,787	Δ13 bp	coding (98‑110/984 nt)	kdgT →	2‑dehydro‑3‑deoxy‑D‑gluconate:H(+) symporter

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	4101787	4101799	13	53 [0]	[0] 53	kdgT	2‑dehydro‑3‑deoxy‑D‑gluconate:H(+) symporter

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4101786	0 (0.000)	53 (0.610)	38/274	1.1	100%	coding (97/984 nt)	kdgT	2‑dehydro‑3‑deoxy‑D‑gluconate:H(+) symporter
	?	NC_000913	4101800 =	0 (0.000)					coding (111/984 nt)	kdgT	2‑dehydro‑3‑deoxy‑D‑gluconate:H(+) symporter

TTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGTTTTTTTGCATGGGGGCGTCAATAAAATTAAGCGCGACGGGAACGGTACTGCGTAAATCCGGTACGCT  >  NC_000913/4101649‑4101919                                                                                                                                          |                                                                                                                                      ttACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGat                                                                                                                                     >  2:794833‑M1/1‑138 (MQ=255) ttACACATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGat                                                                                                                                     >  1:839944‑M1/1‑138 (MQ=255)   aCGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatca                                                                                                                                   <  1:794833‑M1/140‑5 (MQ=255)     gCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGTGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcatt                                                                                                                                 <  1:299527‑M1/140‑7 (MQ=255)         cAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGAGTGAGAAAATCCCGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcacc                                                                                                                             >  2:1265234‑M1/1‑130 (MQ=255)           gCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaa                                                                                                                           >  2:1445423‑M1/1‑128 (MQ=255)           gCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaa                                                                                                                           >  2:695958‑M1/1‑128 (MQ=255)               ttAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacggg                                                                                                                       <  1:695958‑M1/140‑17 (MQ=255)                    aaTATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGAGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgat                                                                                                                  <  1:1265234‑M1/140‑22 (MQ=255)                      tatGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgatta                                                                                                                >  2:1148205‑M1/1‑117 (MQ=255)                             aCAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtac                                                                                                         >  1:867445‑M1/1‑110 (MQ=255)                             aCAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtac                                                                                                         >  1:205509‑M1/1‑110 (MQ=255)                                 gTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtg                                                                                                     <  1:1148205‑M1/140‑35 (MQ=255)                                  tGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgc                                                                                                    >  1:456379‑M1/1‑105 (MQ=255)                                      ggAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccat                                                                                                >  1:1108413‑M1/1‑101 (MQ=255)                                        aaaTGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattc                                                                                              >  1:295329‑M1/1‑99 (MQ=255)                                        aaaTGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccagtacggtgcccattc                                                                                              <  2:839944‑M1/140‑42 (MQ=255)                                         aaTGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattct                                                                                             >  2:709310‑M1/1‑98 (MQ=255)                                           tGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctgg                                                                                           >  2:1328712‑M1/1‑96 (MQ=255)                                             cAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcg                                                                                         >  1:151479‑M1/1‑94 (MQ=255)                                                  aaaaCGCTCGATTGAGAAAATCCCGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcggtgtg                                                                                    >  2:114370‑M1/1‑89 (MQ=255)                                                   aaaCGCTCGATTGAGAAAATCCCGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcggtgtgg                                                                                   >  2:1362235‑M1/1‑88 (MQ=255)                                                   aaaCGCTCGATTGAGAAAATCCCGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcggtgtgg                                                                                   <  2:331570‑M1/140‑53 (MQ=255)                                                        cTCGATTGAGAAAATCCCGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcggtgtggttttt                                                                              <  2:205509‑M1/140‑58 (MQ=255)                                                        cTCGATTGAGAAAATCCCGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcggtgtggttttt                                                                              <  2:867445‑M1/140‑58 (MQ=255)                                                        cTCGATTGAGAAAATCCCGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcggtgtggttttt                                                                              <  2:441090‑M1/140‑58 (MQ=255)                                                         tCGATTGAGAAAATCCCGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcggtgtggtttttt                                                                             >  1:743156‑M1/1‑82 (MQ=255)                                                                      tCCCGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcgccaacgggatgattaccggtacggtgcccattctggcggtgtggtttttttgcatgggggcgt                                                                <  1:1445423‑M1/140‑72 (MQ=255)                                                                       cccGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcggtgtggtttttttgcatgggggcgtc                                                               >  2:860945‑M1/1‑68 (MQ=255)                                                                          gggggggATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcggtgtggtttttttgcatgggggcgtcaat                                                            >  1:790507‑M1/1‑65 (MQ=255)                                                                                       cTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcggtgtggtttttttgcatgggggcgtcaataaaattaagcgcg                                               >  2:249608‑M1/1‑52 (MQ=255)                                                                                         cGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcggtgtggtttttttgcatgggggcgtcaataaaattaagcgcgac                                             <  1:860945‑M1/140‑91 (MQ=255)                                                                                         cGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcggtgtggtttttttgcatgggggcgtcaataaaattaagcgcgac                                             <  2:521197‑M1/140‑91 (MQ=255)                                                                                         cGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcggtgtggtttttttgcatgggggcgtcaataaaattaagcgcgac                                             <  1:709310‑M1/140‑91 (MQ=255)                                                                                         cGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcggtgtggtttttttgcatgggggcgtcaataaaattaagcgcgac                                             <  1:647953‑M1/140‑91 (MQ=255)                                                                                         cGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcggtgtggtttttttgcatgggggcgtcaataaaattaagcgcgac                                             <  1:158875‑M1/140‑91 (MQ=255)                                                                                         cGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcggtgtggtttttttgcatgggggcgtcaataaaattaagcgcgac                                             <  2:295329‑M1/140‑91 (MQ=255)                                                                                          gTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcggtgtggtttttttgcatgggggcgtcaataaaattaagcgcgacg                                            >  2:427017‑M1/1‑49 (MQ=255)                                                                                                cTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcggtgtggtttttttgcatgggggcgtcaataaaattaagcgcgacgggaacg                                      <  2:376414‑M1/140‑98 (MQ=255)                                                                                                   ttCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcggtgtggtttttttgcatgggggcgtcaataaaattaagcgcgacgggaacggta                                   <  1:114370‑M1/140‑101 (MQ=255)                                                                                                   ttCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcggtgtggtttttttgcatgggggcgtcaataaaattaagcgcgacgggaacggta                                   <  1:427017‑M1/140‑101 (MQ=255)                                                                                                             cACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcggtgtggtttttttgcatgggggcgtcaataaaattaagcgcgacgggaacggtactgcgtaaat                         >  1:895347‑M1/1‑30 (MQ=255)                                                                                                             cACTGTGCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcggtgtggtttttttgcatgggggcgtcaataaaattaagcgcgacgggaacggtactgcgtaaat                         >  2:659390‑M1/1‑30 (MQ=255)                                                                                                                tgtgCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcggtgtggtttttttgcatgggggcgtcaataaaattaagcgcgacgggaacggtactgcgtaaatccg                      >  1:223733‑M1/1‑27 (MQ=255)                                                                                                                tgtgCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcggtgtggtttttttgcatgggggcgtcaataaaattaagcgcgacgggaacggtactgcgtaaatccg                      >  1:1156951‑M1/1‑27 (MQ=255)                                                                                                                  tgCCACACCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcggtgtggtttttttgcatgggggcgtcaataaaattaagcgcgacgggaacggtactgcgtaaatccggt                    >  2:1305954‑M1/1‑25 (MQ=255)                                                                                                                     cacaCCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcggtgtggtttttttgcatgggggcgtcaataaaattaagcgcgacgggaacggtactgcgtaaatccggtacg                 <  2:895347‑M1/140‑119 (MQ=255)                                                                                                                        acCTTCTCGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcggtgtggtttttttgcatgggggcgtcaataaaattaagcgcgacgggaacggtactgcgtaaatccggtacgctg              <  1:907556‑M1/140‑122 (MQ=255)                                                                                                                            tctcGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcggtgtggtttttttgcatgggggcgtcaataaaattaagcgcgacgggaacggtactgcgtaaatccggtacgctggtgg          >  1:719528‑M1/1‑15 (MQ=255)                                                                                                                            tctcGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcggtgtggtttttttgcatgggggcgtcaataaaattaagcgcgacgggaacggtactgcgtaaatccggtacgctggtgg          >  1:379495‑M1/1‑15 (MQ=255)                                                                                                                            tctcGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcggtgtggtttttttgcatgggggcgtcaataaaattaagcgcgacgggaacggtactgcgtaaatccggtacgctggtgg          >  1:1476256‑M1/1‑15 (MQ=255)                                                                                                                             ctcGCCGGGGGCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcggtgtggtttttttgcatgggggcgtcaataaaattaagcgcgacgggaacggtactgcgtaaatccggtacgctggtggt         >  1:791219‑M1/1‑14 (MQ=255)                                                                                                                                    ggggCGGatcattcaccaacgggatgattaccggtacggtgcccattctggcggtgtggtttttttgcatgggggcgtcaataaaattaagcgcgacgggaacggtactgcgtaaatccggtacgctggtggtaactaaa  >  1:933187‑M1/1‑7 (MQ=255)                                                                                                                                          |                                                                                                                                      TTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGTTTTTTTGCATGGGGGCGTCAATAAAATTAAGCGCGACGGGAACGGTACTGCGTAAATCCGGTACGCT  >  NC_000913/4101649‑4101919

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: ‑