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

CTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGTTTTTTTGCATGGGGGCGTCAATAAAATTAAGCGCGACGGGAACGGTACTGCGTAAATCCGGTACGCTGGTGGTAACTAAA  >  NC_000913/4101662‑4101932                                                                                                                                           |                                                                                                                                     ttacgcatcagcccttaaaaatatgccgacaggtgatggaaatgcagataaaacgctcgattgagaaaatcccgggggggatgatgctcgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggAt                                                                                                                                     >  2:794833‑M2/139‑140 (MQ=255) ttacacatcagcccttaaaaatatgccgacaggtgatggaaatgcagataaaacgctcgattgagaaaatcccgggggggatgatgctcgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggAt                                                                                                                                     >  1:839944‑M2/139‑140 (MQ=255)   acgcatcagcccttaaaaatatgccgacaggtgatggaaatgcagataaaacgctcgattgagaaaatcccgggggggatgatgctcgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCa                                                                                                                                   <  1:794833‑M2/4‑1 (MQ=255)     gcatcagcccttaaaaatatgccgacaggtgatggaaatgcagataaaacgctcgattgagaaaatcccggggtggatgatgctcgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCAtt                                                                                                                                 <  1:299527‑M2/6‑1 (MQ=255)         cagcccttaaaaatatgccgacaggtgatggaaatgcagataaaacgctcgagtgagaaaatcccgggggggatgatgctcgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCAcc                                                                                                                             >  2:1265234‑M2/131‑140 (MQ=255)           gcccttaaaaatatgccgacaggtgatggaaatgcagataaaacgctcgattgagaaaatcccgggggggatgatgctcgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCaa                                                                                                                           >  2:1445423‑M2/129‑140 (MQ=255)           gcccttaaaaatatgccgacaggtgatggaaatgcagataaaacgctcgattgagaaaatcccgggggggatgatgctcgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCaa                                                                                                                           >  2:695958‑M2/129‑140 (MQ=255)               ttaaaaatatgccgacaggtgatggaaatgcagataaaacgctcgattgagaaaatcccgggggggatgatgctcgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACggg                                                                                                                       <  1:695958‑M2/16‑1 (MQ=255)                    aatatgccgacaggtgatggaaatgcagataaaacgctcgattgagaaaatcccggggaggatgatgctcgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACGGgatgat                                                                                                                  <  1:1265234‑M2/21‑1 (MQ=255)                      tatgccgacaggtgatggaaatgcagataaaacgctcgattgagaaaatcccgggggggatgatgctcgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTa                                                                                                                >  2:1148205‑M2/118‑140 (MQ=255)                             acaggtgatggaaatgcagataaaacgctcgattgagaaaatcccgggggggatgatgctcgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTAc                                                                                                         >  1:867445‑M2/111‑140 (MQ=255)                             acaggtgatggaaatgcagataaaacgctcgattgagaaaatcccgggggggatgatgctcgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTAc                                                                                                         >  1:205509‑M2/111‑140 (MQ=255)                                 gtgatggaaatgcagataaaacgctcgattgagaaaatcccgggggggatgatgctcgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTg                                                                                                     <  1:1148205‑M2/34‑1 (MQ=255)                                  tgatggaaatgcagataaaacgctcgattgagaaaatcccgggggggatgatgctcgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGc                                                                                                    >  1:456379‑M2/106‑140 (MQ=255)                                      ggaaatgcagataaaacgctcgattgagaaaatcccgggggggatgatgctcgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCAt                                                                                                >  1:1108413‑M2/102‑140 (MQ=255)                                        aaatgcagataaaacgctcgattgagaaaatcccgggggggatgatgctcgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTc                                                                                              >  1:295329‑M2/100‑140 (MQ=255)                                        aaatgcagataaaacgctcgattgagaaaatcccgggggggatgatgctcgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCAGTACGGTGCCCATTc                                                                                              <  2:839944‑M2/41‑1 (MQ=255)                                         aatgcagataaaacgctcgattgagaaaatcccgggggggatgatgctcgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCt                                                                                             >  2:709310‑M2/99‑140 (MQ=255)                                           tgcagataaaacgctcgattgagaaaatcccgggggggatgatgctcgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTgg                                                                                           >  2:1328712‑M2/97‑140 (MQ=255)                                             cagataaaacgctcgattgagaaaatcccgggggggatgatgctcgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCg                                                                                         >  1:151479‑M2/95‑140 (MQ=255)                                                  aaaacgctcgattgagaaaatcccgggggggatgatgctcgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGtgtg                                                                                    >  2:114370‑M2/90‑140 (MQ=255)                                                   aaacgctcgattgagaaaatcccgggggggatgatgctcgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTgg                                                                                   >  2:1362235‑M2/89‑140 (MQ=255)                                                   aaacgctcgattgagaaaatcccgggggggatgatgctcgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTgg                                                                                   <  2:331570‑M2/52‑1 (MQ=255)                                                        ctcgattgagaaaatcccgggggggatgatgctcgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGttttt                                                                              <  2:205509‑M2/57‑1 (MQ=255)                                                        ctcgattgagaaaatcccgggggggatgatgctcgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGttttt                                                                              <  2:867445‑M2/57‑1 (MQ=255)                                                        ctcgattgagaaaatcccgggggggatgatgctcgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGttttt                                                                              <  2:441090‑M2/57‑1 (MQ=255)                                                         tcgattgagaaaatcccgggggggatgatgctcgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGtttttt                                                                             >  1:743156‑M2/83‑140 (MQ=255)                                                                      tcccgggggggatgatgctcgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCGCCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGTTTTTTTGCATGGGGGCGt                                                                <  1:1445423‑M2/71‑1 (MQ=255)                                                                       cccgggggggatgatgctcgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGTTTTTTTGCATGGGGGCGTc                                                               >  2:860945‑M2/69‑140 (MQ=255)                                                                          gggggggatgatgctcgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGTTTTTTTGCATGGGGGCGTCAAt                                                            >  1:790507‑M2/66‑140 (MQ=255)                                                                                       ctcgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGTTTTTTTGCATGGGGGCGTCAATAAAATTAAGcgcg                                               >  2:249608‑M2/53‑140 (MQ=255)                                                                                         cgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGTTTTTTTGCATGGGGGCGTCAATAAAATTAAGCGCGAc                                             <  1:860945‑M2/90‑1 (MQ=255)                                                                                         cgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGTTTTTTTGCATGGGGGCGTCAATAAAATTAAGCGCGAc                                             <  2:521197‑M2/90‑1 (MQ=255)                                                                                         cgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGTTTTTTTGCATGGGGGCGTCAATAAAATTAAGCGCGAc                                             <  1:709310‑M2/90‑1 (MQ=255)                                                                                         cgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGTTTTTTTGCATGGGGGCGTCAATAAAATTAAGCGCGAc                                             <  1:647953‑M2/90‑1 (MQ=255)                                                                                         cgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGTTTTTTTGCATGGGGGCGTCAATAAAATTAAGCGCGAc                                             <  1:158875‑M2/90‑1 (MQ=255)                                                                                         cgtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGTTTTTTTGCATGGGGGCGTCAATAAAATTAAGCGCGAc                                             <  2:295329‑M2/90‑1 (MQ=255)                                                                                          gtcccgctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGTTTTTTTGCATGGGGGCGTCAATAAAATTAAGCGCGACg                                            >  2:427017‑M2/50‑140 (MQ=255)                                                                                                ctattccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGTTTTTTTGCATGGGGGCGTCAATAAAATTAAGCGCGACGGGAACg                                      <  2:376414‑M2/97‑1 (MQ=255)                                                                                                   ttccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGTTTTTTTGCATGGGGGCGTCAATAAAATTAAGCGCGACGGGAACGGTa                                   <  1:114370‑M2/100‑1 (MQ=255)                                                                                                   ttccttggcgcactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGTTTTTTTGCATGGGGGCGTCAATAAAATTAAGCGCGACGGGAACGGTa                                   <  1:427017‑M2/100‑1 (MQ=255)                                                                                                             cactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGTTTTTTTGCATGGGGGCGTCAATAAAATTAAGCGCGACGGGAACGGTACTGCGTAAAt                         >  1:895347‑M2/31‑140 (MQ=255)                                                                                                             cactgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGTTTTTTTGCATGGGGGCGTCAATAAAATTAAGCGCGACGGGAACGGTACTGCGTAAAt                         >  2:659390‑M2/31‑140 (MQ=255)                                                                                                                tgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGTTTTTTTGCATGGGGGCGTCAATAAAATTAAGCGCGACGGGAACGGTACTGCGTAAATCCg                      >  1:223733‑M2/28‑140 (MQ=255)                                                                                                                tgtgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGTTTTTTTGCATGGGGGCGTCAATAAAATTAAGCGCGACGGGAACGGTACTGCGTAAATCCg                      >  1:1156951‑M2/28‑140 (MQ=255)                                                                                                                  tgccacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGTTTTTTTGCATGGGGGCGTCAATAAAATTAAGCGCGACGGGAACGGTACTGCGTAAATCCGGt                    >  2:1305954‑M2/26‑140 (MQ=255)                                                                                                                     cacaccttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGTTTTTTTGCATGGGGGCGTCAATAAAATTAAGCGCGACGGGAACGGTACTGCGTAAATCCGGTACg                 <  2:895347‑M2/118‑1 (MQ=255)                                                                                                                        accttctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGTTTTTTTGCATGGGGGCGTCAATAAAATTAAGCGCGACGGGAACGGTACTGCGTAAATCCGGTACGCtg              <  1:907556‑M2/121‑1 (MQ=255)                                                                                                                            tctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGTTTTTTTGCATGGGGGCGTCAATAAAATTAAGCGCGACGGGAACGGTACTGCGTAAATCCGGTACGCtggtgg          >  1:719528‑M2/16‑140 (MQ=255)                                                                                                                            tctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGTTTTTTTGCATGGGGGCGTCAATAAAATTAAGCGCGACGGGAACGGTACTGCGTAAATCCGGTACGCtggtgg          >  1:379495‑M2/16‑140 (MQ=255)                                                                                                                            tctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGTTTTTTTGCATGGGGGCGTCAATAAAATTAAGCGCGACGGGAACGGTACTGCGTAAATCCGGTACGCtggtgg          >  1:1476256‑M2/16‑140 (MQ=255)                                                                                                                             ctcgccgggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGTTTTTTTGCATGGGGGCGTCAATAAAATTAAGCGCGACGGGAACGGTACTGCGTAAATCCGGTACGCTggtggt         >  1:791219‑M2/15‑140 (MQ=255)                                                                                                                                    ggggcggATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGTTTTTTTGCATGGGGGCGTCAATAAAATTAAGCGCGACGGGAACGGTACTGCGTAAATCCGGTACGCTGGTGGTAACTaaa  >  1:933187‑M2/8‑140 (MQ=255)                                                                                                                                           |                                                                                                                                     CTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCGGTGTGGTTTTTTTGCATGGGGGCGTCAATAAAATTAAGCGCGACGGGAACGGTACTGCGTAAATCCGGTACGCTGGTGGTAACTAAA  >  NC_000913/4101662‑4101932

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