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

TTTGATTAAATTTAATTTTTTGTATCGATAATTGTGAACCTCTCTGGCATGGAGAACTATTTTGAACATGAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGT  >  NC_000913/2903866‑2904321                                                                                                                                                                                                                                             |                                                                                                                                                                                                                            tttGATTAAATTTAATTTTTTGTATCGATAATTGTGAACCTCTCTGGCATGGAGAACTATTTTGAACATGAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaat                                                                                                                                                                                                                >  1:424715‑M1/1‑237 (MQ=255)  tagatTAAATTTAATTTTTTGTATCGAAAATTGTGAACCTCTCTGGCATGGAGAACTATTTGGAACATGAGGTGTTACGTGGATATGTTGCTTATTACAATTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaata                                                                                                                                                                                                               <  2:482333‑M1/248‑15 (MQ=255)              aatttTTTGTATCGATAATTGTGAACCTCACTGGCATGGAGGACTATTTTGAACATGAGGTGTTACGTGGATATGTTGATTATTACTAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCACGCTGGCGCGGGGAACTCGAAAGAAAGGCCTCAGTAGTCCCGTCAGGCTCCAGTTTATCCCCGCTGGAGCGGGG‑ACAAgaggcgcacggaaacaaaagcgttgtat                                                                                                                                                                                                 >  2:549154‑M1/1‑223 (MQ=255)                ttttttGTATCGATAATTGTGAACCTCTCTGGCATGGAGAACTATTTTGAACATGAGGTGTTACGTGGCTATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGCGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGAGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGGTTATCCCCGCTGGCGCGGGGAACACgggcgcacggagtaaaaaaccgagaatct                                                                                                                                                                                                >  2:441558‑M1/1‑222 (MQ=255)                       tATCGATAATTGTGAACCTCTCTGGCATGGAGAACTATTTTGAACATGAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATCAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCGCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggagcacggaatacaaagcctagtatctgcacggg                                                                                                                                                                                         >  2:58755‑M1/1‑215 (MQ=255)                        aTCGATAATTGTGAACCTCTCTGGCATGGAGAACTATTTTGAACATGAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAATACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATACCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGGGTATCCCCGCTGACGCGGGGAACACgggcgcccggaatacaaagccgt                                                                                                                                                                                                      >  2:154664‑M1/1‑214 (MQ=255)                        aTCGATAATTGTGAACCTCTCTGGCATGGAGAACTATTTTGAACATGAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaatacaaagccgt                                                                                                                                                                                                      <  1:154664‑M1/237‑24 (MQ=255)                          cGATAATTGTGAACCTCTCTGGCATGGAGAACTATTTTGAACATGAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaatacaaagccgtgtatctgctcggttt                                                                                                                                                                                       >  1:198460‑M1/1‑212 (MQ=255)                                     aaCCTCTCTGGCATGGAGAACTATTTTGAACAAGAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCACCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaatacacagccgtggatctgctcggtttatccccgccgcc                                                                                                                                                                           >  2:270805‑M1/1‑201 (MQ=255)                                                 aTGGAGAACTATTTTGAACATGAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaatacaaagccgtgtatctgctcggtttatccccgctggcgcggggaacact                                                                                                                                                               >  1:104866‑M1/1‑189 (MQ=255)                                                                                         ttGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaatacaaagccgtgtatctgctcggtttatccccgctggcgcggggaacactggctctgcaacagcagcacccatgaccacgtcggtttatc                                                                                                                       >  1:153506‑M1/1‑149 (MQ=255)                                                                                                ttaCAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcccggaatacaaagccgtgtatctgctcggtttatccccgctggcgcggggaacactgtctctgcaacagcagcacccatgaccacgtcgggttatcaccgctg                                                                                                                >  2:223380‑M1/1‑142 (MQ=255)                                                                                                         cTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCAGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaatacaaagccgtgtatcttctcggtttatcccagctggcgcggggaacactgactctgcaacaacagcacccctaccaacgtcggtttatccccacaggcgcgggg                                                                                                        >  2:431487‑M1/1‑133 (MQ=255)                                                                                                            cTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaatacaaagccgtgtatctgctcggtttatccccgctggcgcggggaacactggctctgcaacagcagcacccatgaccacgtcggtttatccccgctggcgcggggaac                                                                                                     >  1:471171‑M1/1‑130 (MQ=255)                                                                                                              aaTATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaatacaaagccgtgtatctgctcggtttatccccgctggcgcggggaacactggctctgcaacagcagcacccatgaccacgtcggtttatccccgctggcgcggggaaca                                                                                                    <  1:149239‑M1/249‑122 (MQ=255)                                                                                                                    aaaaCTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaatacaaagccgtgtatctgctcggtttatccccgctggcgcggggaacactggctctgcaacagcagcgcccatgaccacgtcggtttatccccgctggcgcggggaacacgaaatgc                                                                                            <  1:441558‑M1/251‑130 (MQ=255)                                                                                                                             gaAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaatacaaagccgtgtatctgctcggtttatccccgctggcgcggggaacactggctctgcaacagcagcacccatgaccacgtcggtttatccccgctggcgcggggaacacgaaatgctggtgagcg                                                                                   <  1:508984‑M1/251‑139 (MQ=255)                                                                                                                                         gggTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaatacaaagccgtgtatctgctcggtttatccccgctggcgcggggaacactggctctgcaacagcagcacccatgaccacgtcggtttatccccgctggcgcggggaacacgaaatgctggtgagcgttaatgccgca                                                                        >  1:546867‑M1/1‑101 (MQ=255)                                                                                                                                                              ccGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggcatacaaagccgtgtatctgctcggtttatccccgctggcgcggggaacactggctctgcaacagcagcacacatgaccacgtcggtttatccccgctggcgcggggcacacgaaatgctagtgagcgataatgccgcaaacacagggttatacccgctgg                                                  >  2:363732‑M1/1‑80 (MQ=255)                                                                                                                                                                   ggCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaatacaaagccgtgtatctgctcggtttatccccgctggcgcggggaacactggctctgcaacagcagcacccatgaccacgtcggtttatccccgctggcgcggggaacacgaaatgctggtgagcgttaatgccgcaaacacaggtttatccccgctggcgcgg                                             >  1:532766‑M1/1‑75 (MQ=255)                                                                                                                                                                   ggCGCGGGGAACTCGACAGAACGGCCGCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaatacaaagccgtgtatctgctcggtttatccccgctggcgcggggaacactggctctgcaacagcagcacccatgaccacgtcggtttatccccgctggcgcggggaacacgaaatgctggtgagcgttaatgccgcaaacacaggtttatccccgctggcgc                                               <  1:270805‑M1/249‑175 (MQ=255)                                                                                                                                                                        ggggAACGCGACAGAACGGCCACAGTAGTGTCGTCAGGCGCCGGGTTATCCCCGCTGTCGCGGGGAACACgagcgcacggaattcaaagacgttaatctgctcggcttatccccactggcgcggggaacactggctcttcaacagcagcacccatgaccacgtcggtttatccccgctggcgcggggaacacgaaatgctggtgagcgttaatgccgcaaacacaggtttatccccgcttgcgcagggaac                                        <  1:165186‑M1/251‑182 (MQ=255)                                                                                                                                                                                                              cTCCGGTTTATCCCCGCTGGCGCGGGGAACACgggcgcacggaatacaaagccgtgtatctgctcggtttatccccgctggcgcggggaacactggctctgcaacagcagcacccatgaccacgtcggtttatccccgctggcgcggggaacacgaaatgctggtgatcgttaatgccgcaaacacaggtttatccccgctggcgcggggaacacattacgcctttttgcgattgcccggtttttgccggt  >  2:111855‑M1/1‑32 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                            TTTGATTAAATTTAATTTTTTGTATCGATAATTGTGAACCTCTCTGGCATGGAGAACTATTTTGAACATGAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGT  >  NC_000913/2903866‑2904321

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