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

TTATTTTTACTATCAAAAATAGCTCATTGCCTTACTATAATTAATATATAAATTCGGTAATTAATTCTTAACATGCTTTTACTAATAATCTCAATTGCTGGCCCTATAATATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGAACTGAAAGCCCGTCCTATAATGGTAGGTTTTTTCTTCACCGGTAGCGCTATTATGGGAATTCTGGTCAGTCAATTTCTGGCAAGGCACTCCGATAAACAAGGCGA  >  NC_000913/3836723‑3837188                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                         ttATTTTTACTATCAAAAATAGCGCATTGCCTTACTATAATTAATATATAAATTCGGTAATTAATTCTTAACATGCTTTTACTAATAATCTCAATTGCTGGCCCTATAATATTGCGCTAGCATTGATTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccgg                                                                                                                                                                                                                         >  1:5194‑M1/1‑234 (MQ=255)            aTCAAAAATAGCTCATTGCCTTACTATAATTAATATATAAATTCGGTAATTAATTCTTAACATGCTTTTACTAATAATCTCAATTGCTGGCCCTATAATATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgg                                                                                                                                                                                                              >  1:410044‑M1/1‑223 (MQ=255)                   aTAGCTCATTGCCTTACTATAATTAATATATAAATTCGGTAATTAATTCTTAACATGCTTTTACTAATAATCTCAATTGCTGGCCCTATAATATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaa                                                                                                                                                                                                         >  2:527952‑M1/1‑216 (MQ=255)                                   cTATAATTAATATATAAATTCGGTAATTAATTCTTAACATGCTTTTACTAATAATCTCAATTGCTGGCCCTATAATATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgttt                                                                                                                                                                                        >  1:58943‑M1/1‑200 (MQ=255)                                       aattaatATATAAATTCGGTAATTAATTCTTAACATGCTTTTACTAATAATCTCAATTGCTGGCCCTATAATATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgtttatgcc                                                                                                                                                                                   >  1:75891‑M1/1‑196 (MQ=255)                                       aattaatATATAAATTCGGTAATTAATTCTTAACATGCTTTTACTAATAATCTCAATTGCTGGCCCTATAATATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaaggcaccagttacgtttatgcc                                                                                                                                                                                   <  1:571752‑M1/251‑56 (MQ=255)                                        atgaatataTAAATTCGGTAATTAATTCTTAACATGCTTTTACTAATAACCTCAATGGCTGGCCCTATAATATTGCGCTAGCATGGCTTGCTGGTGGTTTCAGCGATCCTCAAAAATTCTTAATACAATACTCATTCGACCGGTACTTATTTGTAACTCAGTTATATTTTTTCGTCCGGTGATTCATAAGACTGCtaataataaaagcccggcgtcatgccgggaaaacgtcacctgttgcctttatgcca                                                                                                                                                                                  <  2:45103‑M1/248‑57 (MQ=255)                                                ataAATTTGGTAATTAATTCTTAAAATGCTTTTACTAATAATCTCAATTGCTGGCCCTATAATATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaatagcccggcgtcatgccgggcaaaagacaccagttacgtttaagccactgtcaac                                                                                                                                                                          >  2:196968‑M1/1‑187 (MQ=255)                                                  aaaTTCGGTAATTAATTCTTAACATGCTTTTACTAATAATCTCAATTGCTGGCCCTATAATATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgtttatgccactgtcaactg                                                                                                                                                                        <  1:118704‑M1/251‑67 (MQ=255)                                                                    ttAACATGCTTTTACTAATAATCTCAATTGCTGGCCCTATAATATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgtttatgccactgtcaactgctgaatttttttctcgcg                                                                                                                                                      >  1:222244‑M1/1‑167 (MQ=255)                                                                                          ctcAATTGCTGGCCCTATAATATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgtttatgccactgtcaactgctgaatttttttctcgcggcggattttgcgttcttccatt                                                                                                                                >  2:434061‑M1/1‑145 (MQ=255)                                                                                               ttGCTGGCCCTATAATATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgtttatgccactgtcaactgctgaatttttttctcgcggcggattttgcgttcttccattaccgc                                                                                                                           >  1:23901‑M1/1‑140 (MQ=255)                                                                                                 gCTGGACCTATAATATTGCGTTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgtttatgccactgtcaactgctgaatttttttctcgcggcggattttgcgttcttccattaccgcc                                                                                                                          <  1:222048‑M1/250‑113 (MQ=255)                                                                                                               atTGCGCTAGCATAACTTTCTGGATGAATAAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTCTTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAAAAGAATGCtaatacaaaaagcactgcgtcctgccgggcaaaagtcacccgttacgtttatgccactgtcaactgctgaatttttttctcgcggcggattttacgttcttacattacccccactatccccatccgg                                                                                                           >  2:373357‑M1/1‑124 (MQ=255)                                                                                                                   cgcTAGCATTGCTTTCTGGTTGTATAAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagaccggcgtcatgccgggcaacagtcaccagttaagtttatgccactgtcaactgctgaatttttttctcgcggcggattttgcgatcttccattaccgccactatcgccatcaggc                                                                                                          >  1:378377‑M1/1‑120 (MQ=255)                                                                                                                    gcTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgtttatgccactgtcaactgctgaatttttttctcgcggcggatttttcgttcttccattaccgccactatcgccatcaggccgat                                                                                                      >  2:92783‑M1/1‑119 (MQ=255)                                                                                                                              cTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagacaccagttacgtttatgccactgtcaactgctgaatttttttctcgcggcggattttgccatcttccattaccgccactatcgccatcaggcagatcaaaccaatc                                                                                            >  1:15111‑M1/1‑109 (MQ=255)                                                                                                                                         tATCAGCGATACTCAAAAATTCTTAATACACTACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgtttatgccactgtcaactgctgaatttttttctcgcggcggattttgcgttcttccattaccgccactatcgccatcaggcagaaacacccaccagcggcgatatc                                                                                 >  2:412803‑M1/1‑98 (MQ=255)                                                                                                                                            cAGCGATACGCAAAAATTCTTAATACAATACTCATTCGACTGGAACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgtttatgccactgtcaactgctgaatttttttctcgcggcggattttgcgttcttccattaccgccactatcgccatcaggcagatacaaccaatcgcggcgatatcc                                                                                >  1:44764‑M1/1‑95 (MQ=255)                                                                                                                                                   aCTCAAAAATTCTTAATACAATACTCATCCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagaccggcgtcatgccgggcaaaagtcaccagttacgtttatgccactgtcaactgctgaatttttttctcgcggcgggttttgcattcttccattaccgccactatcgccatcaggcagatacaaacaatcgcgactatatccagcgcggcg                                                                       >  1:144999‑M1/1‑88 (MQ=255)                                                                                                                                                       aaaaaTTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgtttatgccactgtcaactgctgaatttttttctcgcggcggattttgcgttcttccattaccgccactatcgccatcaggcagatacaaccaatcgcggcgatatc                                                                                 >  1:411709‑M1/1‑84 (MQ=255)                                                                                                                                                         aaaTTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgtttatgccactgtcaactgctgaatttttttctcgcggcggattttgcgttcttccattaccgccactatcgccatcaggcagatacaaccaatcgcggcgatatccagcgcggcgaaggtg                                                                 >  1:410993‑M1/1‑82 (MQ=255)                                                                                                                                                          aaTTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGCAGAATGCtaataaaaaaagccaggcgtcatgccgggcaaaagtcaccagttacgtttatgccactgtcaactgctgaatttttttctcgcggcggattttgcgttcttccattaccgccactatcgccatcaggcagataccaccaatcccgtccatatcccgcgcggcgaaggtgc                                                                >  2:357632‑M1/1‑81 (MQ=255)                                                                                                                                                                tAATACAATACTCATTCGACTGGTCCTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgtttatgccactgtcaactgctgaatttttttctcgcggcggattttgcgttcttccattaccgccactatcgccatcaggcagatacaaccaatcgcggcgatatccagcgcggcgcaggggcctgccc                                                          >  2:459970‑M1/1‑75 (MQ=255)                                                                                                                                                                         acacaatcgactgGTACTTAATTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGCAGAATGCtaataaaaaaagcccggcgtcaagccgggcaaaagtcaccatttacgtttatgccactatcaactgctgaatttttttctcgcggcggattttgcgttcttccatgaccgccactaacgccatcaggcagatacaaccaatctcggcgatatccagcgcggcgaaggtgcctgcacagccggtc                                                  >  2:464021‑M1/7‑66 (MQ=255)                                                                                                                                                                                      gTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgtttatgccactgtcaactgctgaatttttttctcgcggcggattttgcgttcttccattaccgccactatcgccatcaggcagatacaaccaatcgcggcgatatccagcgcggcgaaggtgcctgcccagccggtaaggccgaataccgg                                    >  1:91890‑M1/1‑53 (MQ=255)                                                                                                                                                                                                            tatTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgtttatgccactgtcaactgctgaatttttttctcgcggcggattttgcgttcttccattaccgccactatcgccatcaggcagatacaaccaatcgcggcgatatccagcgcggcgaaggtgcctgcccagccggtcaggccgaataccggcgtcccatcggcaatcattccc              >  1:232634‑M1/1‑31 (MQ=255)                                                                                                                                                                                                                 tttCGCCCGGTGATTCAGAAGAATGCtactaaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgtttatgccactgtcaactgctgaatttttttctcgcggcggattttgcgttcttccattaccgcccctatcgccatcaggcagatacaaccaatcgcggcgctatccagcgcggcgaaggtgcctgcccagccggtaaggccgaataccggcgtaccatcggcaatcattcccagacc         >  2:222992‑M1/1‑26 (MQ=255)                                                                                                                                                                                                                 tttCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgtttatgccactgtcaactgctgaatttttttctcgcggcggattttgcgttcttccattaccgccactatcgccatcaggcagatacaaccaatcgcggcgatatccagcgcggcgaaggtgcctgcccagccggtaaggccgaataccggcgtcccatcggcaatcattcccagacc         >  1:300222‑M1/1‑26 (MQ=255)                                                                                                                                                                                                                        cGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgtttatgccactgtcaactgctgaatttttttctcgcggcggattttgcgttcttccattaccgccactatcgccatcaggcagatacaaccaatcgcggcgatatccagcgcggcgaaggtgcctgcccagccggtaaggccgaataccggcgtcccatcggccatcattcccagacctaacttg  >  2:236353‑M1/1‑19 (MQ=255)                                                                                                                                                                                                                         ggTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgtttatgccactgtcaactgctgaatttttttctcgcggcggattttgcgttcttccattaccgccactatcgccatcaggcagatacaaccaatcgcggcgatatccagcgcggcgaaggtgcctgcccagccggt                                                   >  2:493233‑M1/1‑18 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                         TTATTTTTACTATCAAAAATAGCTCATTGCCTTACTATAATTAATATATAAATTCGGTAATTAATTCTTAACATGCTTTTACTAATAATCTCAATTGCTGGCCCTATAATATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGAACTGAAAGCCCGTCCTATAATGGTAGGTTTTTTCTTCACCGGTAGCGCTATTATGGGAATTCTGGTCAGTCAATTTCTGGCAAGGCACTCCGATAAACAAGGCGA  >  NC_000913/3836723‑3837188

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