Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,836,957	Δ8,771 bp		[setC]–[adeD]	[setC], yicL, yicU, nlpA, yicS, nepI, yicN, adeQ, [adeD]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3836957	3845727	8771	31 [0]	[0] 31	[setC]–[adeD]	[setC],yicL,yicU,nlpA,yicS,nepI,yicN,adeQ,[adeD]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3836956	0 (0.000)	31 (0.430)	29/492	0.9	100%	coding (4/1185 nt)	setC	putative arabinose exporter
	?	NC_000913	3845728 =	0 (0.000)					coding (1765/1767 nt)	adeD	adenine deaminase

CAAAGTCATCTGCCGTTACCCATTGCCGGGCTGATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTAACTTG  >  NC_000913/3845494‑3845959                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                        ttatttttactatcaaaaatagcgcattgccttactataattaatatataaattcggtaattaattcttaacatgcttttactaataatctcaattgctggccctataatattgcgctagcattgatttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCgg                                                                                                                                                                                                                         >  1:5194‑M2/235‑251 (MQ=255)            atcaaaaatagctcattgccttactataattaatatataaattcggtaattaattcttaacatgcttttactaataatctcaattgctggccctataatattgcgctagcattgctttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCgg                                                                                                                                                                                                              >  1:410044‑M2/224‑251 (MQ=255)                   atagctcattgccttactataattaatatataaattcggtaattaattcttaacatgcttttactaataatctcaattgctggccctataatattgcgctagcattgctttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCaaa                                                                                                                                                                                                         >  2:527952‑M2/217‑249 (MQ=255)                                   ctataattaatatataaattcggtaattaattcttaacatgcttttactaataatctcaattgctggccctataatattgcgctagcattgctttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGttt                                                                                                                                                                                        >  1:58943‑M2/201‑250 (MQ=255)                                       aattaatatataaattcggtaattaattcttaacatgcttttactaataatctcaattgctggccctataatattgcgctagcattgctttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGcc                                                                                                                                                                                   >  1:75891‑M2/197‑251 (MQ=255)                                       aattaatatataaattcggtaattaattcttaacatgcttttactaataatctcaattgctggccctataatattgcgctagcattgctttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGGCACCAGTTACGTTTATGcc                                                                                                                                                                                   <  1:571752‑M2/55‑1 (MQ=255)                                        atgaatatataaattcggtaattaattcttaacatgcttttactaataacctcaatggctggccctataatattgcgctagcatggcttgctggtggtttcagcgatcctcaaaaattcttaatacaatactcattcgaccggtacttatttgtaactcagttatattttttcgtccggtgattcataagactgcTAATAATAAAAGCCCGGCGTCATGCCGGGAAAACGTCACCTGTTGCCTTTATGCCa                                                                                                                                                                                  <  2:45103‑M2/56‑1 (MQ=255)                                                ataaatttggtaattaattcttaaaatgcttttactaataatctcaattgctggccctataatattgcgctagcattgctttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAATAGCCCGGCGTCATGCCGGGCAAAAGACACCAGTTACGTTTAAGCCACTGTCAAc                                                                                                                                                                          >  2:196968‑M2/188‑251 (MQ=255)                                                  aaattcggtaattaattcttaacatgcttttactaataatctcaattgctggccctataatattgcgctagcattgctttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAActg                                                                                                                                                                        <  1:118704‑M2/66‑1 (MQ=255)                                                                    ttaacatgcttttactaataatctcaattgctggccctataatattgcgctagcattgctttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTcgcg                                                                                                                                                      >  1:222244‑M2/168‑251 (MQ=255)                                                                                          ctcaattgctggccctataatattgcgctagcattgctttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCAtt                                                                                                                                >  2:434061‑M2/146‑251 (MQ=255)                                                                                               ttgctggccctataatattgcgctagcattgctttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGc                                                                                                                           >  1:23901‑M2/141‑251 (MQ=255)                                                                                                 gctggacctataatattgcgttagcattgctttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGcc                                                                                                                          <  1:222048‑M2/112‑1 (MQ=255)                                                                                                               attgcgctagcataactttctggatgaataagcgatactcaaaaattcttaatacaatactcattcgactggtacttctttgtaactcagttatattttttcgcccggtgattcaaaagaatgcTAATACAAAAAGCACTGCGTCCTGCCGGGCAAAAGTCACCCGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTACGTTCTTACATTACCCCCACTATCCCcatccgg                                                                                                           >  2:373357‑M2/125‑248 (MQ=255)                                                                                                                   cgctagcattgctttctggttgtataagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGACCGGCGTCATGCCGGGCAACAGTCACCAGTTAAGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGATCTTCCATTACCGCCACTATCGCCATCAGGc                                                                                                          >  1:378377‑M2/121‑248 (MQ=255)                                                                                                                    gctagcattgctttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTTCGTTCTTCCATTACCGCCACTATCGCCATCAGGCCGAt                                                                                                      >  2:92783‑M2/120‑251 (MQ=255)                                                                                                                              ctttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGACACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCCATCTTCCATTACCGCCACTATCGCCATCAGGCAGATCAAACCAATc                                                                                            >  1:15111‑M2/110‑251 (MQ=255)                                                                                                                                         tatcagcgatactcaaaaattcttaatacactactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGAAACACCCACCAGCGGCGATATc                                                                                 >  2:412803‑M2/99‑251 (MQ=255)                                                                                                                                            cagcgatacgcaaaaattcttaatacaatactcattcgactggaacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATcc                                                                                >  1:44764‑M2/96‑249 (MQ=255)                                                                                                                                                   actcaaaaattcttaatacaatactcatccgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGACCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGGTTTTGCATTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAAACAATCGCGACTATATCCAGCgcggcg                                                                       >  1:144999‑M2/89‑251 (MQ=255)                                                                                                                                                       aaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATc                                                                                 >  1:411709‑M2/85‑237 (MQ=255)                                                                                                                                                         aaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGtg                                                                 >  1:410993‑M2/83‑251 (MQ=255)                                                                                                                                                          aattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagcagaatgcTAATAAAAAAAGCCAGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACCACCAATCCCGTCCATATCCCGCGCGGCGAAGGtgc                                                                >  2:357632‑M2/82‑251 (MQ=255)                                                                                                                                                                taatacaatactcattcgactggtccttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGCAGGGGCCTGccc                                                          >  2:459970‑M2/76‑251 (MQ=255)                                                                                                                                                                         acacaatcgactggtacttaattgtaactcagttatattttttcgcccggtgattcagcagaatgcTAATAAAAAAAGCCCGGCGTCAAGCCGGGCAAAAGTCACCATTTACGTTTATGCCACTATCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATGACCGCCACTAACGCCATCAGGCAGATACAACCAATCTCGGCGATATCCAGCGCGGCGAAGGTGCCTGCACAGCCGgtc                                                  >  2:464021‑M2/67‑249 (MQ=255)                                                                                                                                                                                      gtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCgg                                    >  1:91890‑M2/54‑251 (MQ=255)                                                                                                                                                                                                            tattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTCAGGCCGAATACCGGCGTCCCATCGGCAATCATTccc              >  1:232634‑M2/32‑251 (MQ=255)                                                                                                                                                                                                                 tttcgcccggtgattcagaagaatgcTACTAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCCCTATCGCCATCAGGCAGATACAACCAATCGCGGCGCTATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTACCATCGGCAATCATTCCCAGAcc         >  2:222992‑M2/27‑251 (MQ=255)                                                                                                                                                                                                                 tttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGAcc         >  1:300222‑M2/27‑251 (MQ=255)                                                                                                                                                                                                                        cggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCCATCATTCCCAGACCTAACTTg  >  2:236353‑M2/20‑251 (MQ=255)                                                                                                                                                                                                                         ggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGt                                                   >  2:493233‑M2/19‑201 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                        CAAAGTCATCTGCCGTTACCCATTGCCGGGCTGATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTAACTTG  >  NC_000913/3845494‑3845959

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