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	23 [0]	[1] 24	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)	23 (0.520)	23/426	0.7	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

TTGAACATGAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACG‑GGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACC‑ACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACATTACGCCTTTTTGCGATTGCCCGGTTTTTGCCGGTTTATCCCCGCTGGCGCGGGGAA  >  NC_000913/2903927‑2904404                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   tttgattaaatttaattttttgtatcgataattgtgaacctctctggcatggagaactattttgaacatgaggtgttacgtggatatgttgcttattacaagtactgctaatataaaaacttgagaaagagataacgggttatatggtggtttatccccgctggcgcggggaactcgacagaacggcctcagtagtctcgtcaggctccggtttatccccgctggcgcggggaacacG‑GGCGCACGGAAt                                                                                                                                                                                                                                       >  1:424715‑M2/238‑250 (MQ=255)  tagattaaatttaattttttgtatcgaaaattgtgaacctctctggcatggagaactatttggaacatgaggtgttacgtggatatgttgcttattacaattactgctaatataaaaacttgagaaagagataacgggttatatggtggtttatccccgctggcgcggggaactcgacagaacggcctcagtagtctcgtcaggctccggtttatccccgctggcgcggggaacacG‑GGCGCACGGAATa                                                                                                                                                                                                                                      <  2:482333‑M2/14‑1 (MQ=255)               aattttttgtatcgataattgtgaacctcactggcatggaggactattttgaacatgaggtgttacgtggatatgttgattattactagtactgctaatataaaaacttgagaaagagataacgggttatatggtggtttatccacgctggcgcggggaactcgaaagaaaggcctcagtagtcccgtcaggctccagtttatccccgctggagcggggacaaGAGGCgcacggaa‑acaaaagcgttgtat                                                                                                                                                                                                                        >  2:549154‑M2/224‑241 (MQ=255)                ttttttgtatcgataattgtgaacctctctggcatggagaactattttgaacatgaggtgttacgtggctatgttgcttattacaagtactgctaatataaaaacttgcgaaagagataacgggttatatggtggtttatccccgctggcgaggggaactcgacagaacggcctcagtagtctcgtcaggctccgggttatccccgctggcgcggggaacacg‑ggcgcacggagtaaaaaaccgagaatct                                                                                                                                                                                                                       >  2:441558‑M2/223‑233 (MQ=255)                       tatcgataattgtgaacctctctggcatggagaactattttgaacatgaggtgttacgtggatatgttgcttattacaagtactgctaatatcaaaacttgagaaagagataacgggttatatggtggtttatccccgctggcgcggggaactcgacagaacggcctcagtagtcgcgtcaggctccggtttatccccgctggcgcggggaacacG‑GGAGCACGGAATACAAAGCCTAGTATCTGCAcggg                                                                                                                                                                                                                >  2:58755‑M2/216‑250 (MQ=255)                        atcgataattgtgaacctctctggcatggagaactattttgaacatgaggtgttacgtggatatgttgcttattacaagtactgctaatataaatacttgagaaagagataacgggttatatggtggtttatacccgctggcgcggggaactcgacagaacggcctcagtagtctcgtcaggctccggggtatccccgctgacgcggggaacacG‑GGCGCCCGGAATACAAAGCCgt                                                                                                                                                                                                                             >  2:154664‑M2/215‑237 (MQ=255)                        atcgataattgtgaacctctctggcatggagaactattttgaacatgaggtgttacgtggatatgttgcttattacaagtactgctaatataaaaacttgagaaagagataacgggttatatggtggtttatccccgctggcgcggggaactcgacagaacggcctcagtagtctcgtcaggctccggtttatccccgctggcgcggggaacacG‑GGCGCACGGAATACAAAGCCgt                                                                                                                                                                                                                             <  1:154664‑M2/23‑1 (MQ=255)                          cgataattgtgaacctctctggcatggagaactattttgaacatgaggtgttacgtggatatgttgcttattacaagtactgctaatataaaaacttgagaaagagataacgggttatatggtggtttatccccgctggcgcggggaactcgacagaacggcctcagtagtctcgtcaggctccggtttatccccgctggcgcggggaacacG‑GGCGCACGGAATACAAAGCCGTGTATCTGCTCGGttt                                                                                                                                                                                                              >  1:198460‑M2/213‑250 (MQ=255)                                     aacctctctggcatggagaactattttgaacaagaggtgttacgtggatatgttgcttattacaagtactgctaatataaaaacttgagaaagagataacgggttatatggtggtttatcaccgctggcgcggggaactcgacagaacggcctcagtagtctcgtcaggctccggtttatccccgctggcgcggggaacacG‑GGCGCACGGAATACACAGCCGTGGATCTGCTCGGTTTATCcccgccgcc                                                                                                                                                                                                  >  2:270805‑M2/202‑247 (MQ=255)                                                 atggagaactattttgaacatgaggtgttacgtggatatgttgcttattacaagtactgctaatataaaaacttgagaaagagataacgggttatatggtggtttatccccgctggcgcggggaactcgacagaacggcctcagtagtctcgtcaggctccggtttatccccgctggcgcggggaacacG‑GGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACt                                                                                                                                                                                      >  1:104866‑M2/190‑251 (MQ=255)                                                                                         ttgcttattacaagtactgctaatataaaaacttgagaaagagataacgggttatatggtggtttatccccgctggcgcggggaactcgacagaacggcctcagtagtctcgtcaggctccggtttatccccgctggcgcggggaacacG‑GGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACC‑ACGTCGGTTTATc                                                                                                                                             >  1:153506‑M2/150‑251 (MQ=255)                                                                                                ttacaagtactgctaatataaaaacttgagaaagagataacgggttatatggtggtttatccccgctggcgcggggaactcgacagaacggcctcagtagtctcgtcaggctccggtttatccccgctggcgcggggaacacG‑GGCGCCCGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGTCTCTGCAACAGCAGCACCCATGACC‑ACGTCGGGTTATCACCGCTg                                                                                                                                      >  2:223380‑M2/143‑251 (MQ=255)                                                                                                         ctgctaatataaaaacttgagaaagagataacgggttatatggtggtttatcccagctggcgcggggaactcgacagaacggcctcagtagtctcgtcaggctccggtttatccccgctggcgcggggaacacG‑GGCGCACGGAATACAAAGCCGTGTATCTTCTCGGTTTATCCCAGCTGGCGCGGGGAACACTGACTCTGCAACAACAGCACCCCT‑ACCAACGTCGGTTTATCCCCACAGGCGCgggg                                                                                                                              >  2:431487‑M2/134‑250 (MQ=255)                                                                                                            ctaatataaaaacttgagaaagagataacgggttatatggtggtttatccccgctggcgcggggaactcgacagaacggcctcagtagtctcgtcaggctccggtttatccccgctggcgcggggaacacG‑GGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACC‑ACGTCGGTTTATCCCCGCTGGCGCGGGGAac                                                                                                                           >  1:471171‑M2/131‑250 (MQ=255)                                                                                                              aatataaaaacttgagaaagagataacgggttatatggtggtttatccccgctggcgcggggaactcgacagaacggcctcagtagtctcgtcaggctccggtttatccccgctggcgcggggaacacG‑GGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACC‑ACGTCGGTTTATCCCCGCTGGCGCGGGGAaca                                                                                                                          <  1:149239‑M2/121‑1 (MQ=255)                                                                                                                    aaaacttgagaaagagataacgggttatatggtggtttatccccgctggcgcggggaactcgacagaacggcctcagtagtctcgtcaggctccggtttatccccgctggcgcggggaacacG‑GGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCGCCCATGACC‑ACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGc                                                                                                                  <  1:441558‑M2/129‑1 (MQ=255)                                                                                                                             gaaagagataacgggttatatggtggtttatccccgctggcgcggggaactcgacagaacggcctcagtagtctcgtcaggctccggtttatccccgctggcgcggggaacacG‑GGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACC‑ACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCg                                                                                                         <  1:508984‑M2/138‑1 (MQ=255)                                                                                                                                         gggttatatggtggtttatccccgctggcgcggggaactcgacagaacggcctcagtagtctcgtcaggctccggtttatccccgctggcgcggggaacacG‑GGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACC‑ACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCa                                                                                              >  1:546867‑M2/102‑250 (MQ=255)                                                                                                                                                              ccgctggcgcggggaactcgacagaacggcctcagtagtctcgtcaggctccggtttatccccgctggcgcggggaacacG‑GGCGCACGGCATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCACACATGACC‑ACGTCGGTTTATCCCCGCTGGCGCGGGGCACACGAAATGCTAGTGAGCGATAATGCCGCAAACACAGGGTTATACCCGCTgg                                                                        >  2:363732‑M2/81‑251 (MQ=255)                                                                                                                                                                   ggcgcggggaactcgacagaacggcctcagtagtctcgtcaggctccggtttatccccgctggcgcggggaacacG‑GGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACC‑ACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCgg                                                                   >  1:532766‑M2/76‑251 (MQ=255)                                                                                                                                                                   ggcgcggggaactcgacagaacggccgcagtagtctcgtcaggctccggtttatccccgctggcgcggggaacacG‑GGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACC‑ACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGgcgc                                                                     <  1:270805‑M2/174‑1 (MQ=255)                                                                                                                                                                        ggggaacgcgacagaacggccacagtagtgtcgtcaggcgccgggttatccccgctgtcgcggggaacacG‑AGCGCACGGAATTCAAAGACGTTAATCTGCTCGGCTTATCCCCACTGGCGCGGGGAACACTGGCTCTTCAACAGCAGCACCCATGACC‑ACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTTGCGCAGGGAac                                                              <  1:165186‑M2/181‑1 (MQ=255)                                                                                                                                                                                                              ctccggtttatccccgctggcgcggggaacacG‑GGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACC‑ACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGATCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACATTACGCCTTTTTGCGATTGCCCGGTTTTTGCCGGt                        >  2:111855‑M2/33‑251 (MQ=255)                                                                                                                                                                                                                                    gggggAACACG‑GGCGCACGTAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACC‑ACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACATTACGCCTTTTTGCGATTGCCCGGTTTTTGCCGGTTTATCCCCGCTGGCGCGGGGaa  <  2:546867/250‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   TTGAACATGAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACG‑GGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACC‑ACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACATTACGCCTTTTTGCGATTGCCCGGTTTTTGCCGGTTTATCCCCGCTGGCGCGGGGAA  >  NC_000913/2903927‑2904404

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