Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	2,904,103	Δ61 bp	intergenic (+729/+584)	ygcE → / ← queE	putative sugar kinase YgcE/putative 7‑carboxy‑7‑deazaguanine synthase QueE

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	2904103	2904163	61	21 [0]	[5] 26	ygcE/queE	putative sugar kinase YgcE/putative 7‑carboxy‑7‑deazaguanine synthase QueE

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2904102	0 (0.000)	21 (0.410)	21/432	1.1	100%	intergenic (+728/+645)	ygcE/queE	putative sugar kinase YgcE/putative 7‑carboxy‑7‑deazaguanine synthase QueE
	?	NC_000913	2904164 =	0 (0.000)					intergenic (+790/+583)	ygcE/queE	putative sugar kinase YgcE/putative 7‑carboxy‑7‑deazaguanine synthase QueE

TAAATTTAATTTTTTGTATCGATAATTGTGAACCTCTCTGGCATGGAGAACTATTTTGAACATGAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAG  >  NC_000913/2903872‑2904319                                                                                                                                                                                                                                       |                                                                                                                                                                                                                          tAAATTTAATTTTTTGTATCGATAATTGTGAACCTCTCTGGCATGGAGAACTATTTTGAACATGAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaatacaaag                                                                                                                                                                                                        <  1:239171‑M1/250‑20 (MQ=255)        aatttTTTGTATCGATAATTGTGAACCTCTCTGGCATGGAGAACTATTTTGAACATGAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCTTCAGGCTCCGGATTCTCCCCGCAGGCGCGGGGAACACgggcgcacgtaatacaaagcctcgtat                                                                                                                                                                                                >  2:99829‑M1/1‑224 (MQ=255)            ttttGTATCGATAATTGTGAACCTCTCTGGCATGGAGAACTATTTTGAACATGAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaatacaaagccgtgtatctgc                                                                                                                                                                                            >  2:319029‑M1/1‑220 (MQ=255)                                      ctGGCATGGAGAACTATTTTGAACATGAGGTGTTACGTGTATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaatacaaagccgtgtatctgctcggtttatccccgctggcgcgg                                                                                                                                                                     <  2:43307‑M1/248‑55 (MQ=255)                                          cATGGAGAACTATTTTGAACATGAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaatacaaagccgtgtatctgctcggtttatccccgctggcgcggggaacac                                                                                                                                                              <  1:117348‑M1/251‑62 (MQ=255)                                             ggAGAACTATTTTGAACATGAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaatacaaagccgtgtatctgctcggtttatccccgctggcgcggggaacactg                                                                                                                                                            <  1:371952‑M1/250‑64 (MQ=255)                                                     aTTTTGAACATGAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaatacaaagccgtgtatctg                                                                                                                                                                                             >  1:396111‑M1/1‑179 (MQ=255)                                                     aTTTTGAACATGAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGGCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcatggaatacaaagccgtgtatctg                                                                                                                                                                                             <  2:396111‑M1/209‑31 (MQ=255)                                                              aTGAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaatacaaagccgtgtatctgctcggtttatccccgctggcgcggggaacactggctctgcaacagcagcac                                                                                                                                          >  1:169776‑M1/1‑170 (MQ=255)                                                                            tgaTATGTTGCTTATTACAAGTTCTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATGCCCGCGGGCGCGTGGAACTCGACAGAACGGCCTCAGTAGGCTCGTCAGGCGCCGGTTTATCCCCGCGGGCGCGGTGAACACgggcgcacggaatacaaagccgtgtgtctgctcggtttatcccctcgggcgcggggaacactggctcttgaactgcagcatccagggccacgtc                                                                                                                             <  2:296400‑M1/249‑95 (MQ=255)                                                                                                     gCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGCGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaatacaaagccgtgtatctgctcggtttatccccgctggcgcggggcacactggctctgcaacagcagcacccatgaccacgtcggtttatccccgctggcgcggggaac                                                                                                   >  1:188911‑M1/1‑131 (MQ=255)                                                                                                                 aCTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaatacaaagccgtgtatctgctcggtttatccccgctggcgcggggaacactggctctgcaacagcagcacccatgaccacgtcggtttatccccgctggcgcggggaacacgaaatgctgg                                                                                       <  1:88458‑M1/251‑133 (MQ=255)                                                                                                                             gaTAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaatacaaagccgtgtatctgctcggtttatccccgctggcgcggggaacactggctctgcaacagcagcacccatgaccacgtcggtttatccccgctggcgcggggaacacgaaatgctggtgagcgttaatg                                                                           >  1:300498‑M1/1‑107 (MQ=255)                                                                                                                                  cGGGTTATATGGTGGTTTATCCCCGCTGGACCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaatacaaagccgtgtatctgctcggtttatccccgctggcgcggggaacactggctctgcaacagcagcaaccatgaccacgtcggtttatccccgctggcgcggggaacacgaaatgctggtgagcgttaatgccgc                                                                       >  2:94651‑M1/1‑102 (MQ=255)                                                                                                                                        atatGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaatacaaagccgtgtatctgctcggtttatccccgctggcgcggggaacactggctctgcaacagcagcacccatgaccacgtcggtttatccccgctggcgcggggaacacgaaatgctggtgagcgttaatgccgcaaacac                                                                 <  2:402641‑M1/250‑155 (MQ=255)                                                                                                                                          atGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaatacaaagccgtgtatctgctcggtttatccccgctggcgcggggaacactggctctgcaacagcagcacccatgaccacgtcggtttatccccgctggcgcggggaacacgaaatgctggtgagcgttaatgccgcaaacacag                                                               <  1:54514‑M1/250‑157 (MQ=255)                                                                                                                                                       cccGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaatacaaagccgtgtatctgctcggtttatccccgctggcgcggggaacactggctctgcaacagcagcacccatgaccacgtcggtttatccccgctggcgcggggaacacgaaatgctggtgagcgttaatgccgcaaacacaggtttatccccgctg                                                 >  1:208309‑M1/1‑81 (MQ=255)                                                                                                                                                          gCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaatacaaagccgtgtatctgctcggtttatccccgctggcgcggggaacactggctctgcaacagcagcacccatgaccacgtcggtttatccccgctggcgcggggaacacgaaatgctggtgagcgttaatgccgcaaacccaggtttatccccgctggcg                                              >  1:424634‑M1/1‑78 (MQ=255)                                                                                                                                                                       aCTCGACAGAACGGCCGTAGTAGTCGCTTCAGGTTCCGTTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaatacaaagccgtgtatctgctcggtttatccccgctggcgcggggaacactggcgctgcaacagcagcacccatgaccacgtcggtttatccccgctggcgctgggaacacgaaatgctggggagcgttaatgccgcaaacacaggtgtatccccgctggcgctgggaacacatt                                 <  2:289591‑M1/251‑187 (MQ=255)                                                                                                                                                                                            agtCTCGTCAGGCTCTGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaatacaaagccgtgtatctgctcggtttatccccgctggcgcggggaacactggctctgcaacagcagcacacatgaccacgtcggtttatccccgctggcgcggggaacacgaaatgctggtgagcgttaatgccgcaaacacaggtttatccccgctggcgcggggaacacattacgcctttttgcgattgccc             <  1:525206‑M1/250‑207 (MQ=255)                                                                                                                                                                                                      ggCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaatacaaagccgtgtatctgctcggtttatccccgctggcgcggggaacactggctctgcaacagcagcacccatgaccacgtcggtttatccccgctggcgcggggaacacgaaatgctggtgagcgttaatgccgcaaacaaaggtttatccccgctggcgcggggaacacattacgcctttttgcgattgcccggtttttgccg  >  2:245353‑M1/1‑34 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                          TAAATTTAATTTTTTGTATCGATAATTGTGAACCTCTCTGGCATGGAGAACTATTTTGAACATGAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAG  >  NC_000913/2903872‑2904319

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