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	20 [0]	[0] 20	[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)	20 (0.310)	18/478	1.5	100%	coding (4/1185 nt)	setC	putative arabinose exporter
	?	NC_000913	3845728 =	0 (0.000)					coding (1765/1767 nt)	adeD	adenine deaminase

AACTTGTTTTATTTTTACTATCAAAAATAGCTCATTGCCTTACTATAATTAATATATAAATTCGGTAATTAATTCTTAACATGCTTTTACTAATAATCTCAATTGCTGGCCCTATAATATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGAACTGAAAGCCCGTCCTATAATGGTAGGTTTTTTCTTCACCGGTAGCGCTATTATGGGAATTCTGGTCAGTCAATTTCTGGCAAG  >  NC_000913/3836715‑3837167                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                    aaCTTGTTTTATTTTTACTATCAAAAATAGCTCATTGCCTTACTATAATTAATATATAAATTCGGTAATTAATTCTTAACATGCTTTTACTAATAATCTCAATTGCTGGCCCTATAATATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaa                                                                                                                                                                                                            <  2:89048‑M1/251‑10 (MQ=255)                  cTATCAAAAATAGCTCATTGCCTTACTATAATTAATATATAAATTCGGTAATTAATTCTTAACATGCTTTTACTAATAATCTCAATTGCTGGCCCTATAATATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgcc                                                                                                                                                                                           >  1:226174‑M1/1‑225 (MQ=255)                                  aTTGCCTTACTATAATTAATATATAAATTCGGTAATTAATTCTTAACATACTTTTACTAATAATCTCAATTGCTGGCCCTATAATATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcacc                                                                                                                                                                             >  1:271963‑M1/1‑209 (MQ=255)                                    tGCCTTACTATAATTAATATATAAATTCGGTAATTAATTCTTAACATGCTTTTACTAATAATATCAATTGCTGGCCCTATAATATTGCGCTAGCCTTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAAAACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCACCCGGTGTTTCGGAAGAATTCtaataaaaaaagcccggcgtcatgccgggcaaaagacaccagtt                                                                                                                                                                         >  2:399714‑M1/1‑207 (MQ=255)                                              taattaatATATAAATTCGGTAATTAATTCTTAACATGCTTTTACTAATAATCTCAATTGCTGGCCCTATAATATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgtttatgc                                                                                                                                                               >  1:326326‑M1/1‑197 (MQ=255)                                                attaatATATAAATTCGGTAATTAATTCTTAACATGCTTTTACTAATAATCTCAATTGCTGGCCCTATAATATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgtttatgccc                                                                                                                                                             >  2:263008‑M1/1‑195 (MQ=255)                                                        ataAATTCGGTAATTAATTCTTAACATGCTTTTACTAATAATCTCAATTGCTGGCCCTATAATATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgtttatgccactgtcaac                                                                                                                                                     >  2:352630‑M1/1‑187 (MQ=255)                                                          aaaTTCGGTAATTAATTCTTAACATGCTTTTACTAATAATCTCAATTGCTGGCCCTATAATATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgtttatgccactgtcaact                                                                                                                                                    <  1:9470‑M1/250‑66 (MQ=255)                                                           aaTTCGGTAATTAATTCTTAACATGCTTTTACTAATAATCTCAATTGCTGGCCCTATAATATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgtttatgccactgtcaactgc                                                                                                                                                  >  1:163586‑M1/1‑184 (MQ=255)                                                                                   gCTTTTACTAATAATCTCAATTGCTGGCCCTATAATAGTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgtttatgccactgtcaactgctgaatttttttctcgcggcggatt                                                                                                                          <  1:482717‑M1/251‑92 (MQ=255)                                                                                            aataatCTCAATTGCTGGCCCTATAATATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgtttatgccactgtcaactgctgaatttttttctcgcggcggattttgcgttct                                                                                                                 >  1:117735‑M1/1‑151 (MQ=255)                                                                                                   tcAATTGGTGGGCCTGTAATATTGAGCTAGCATTGCTTTCTGGATGTATCAGCGACACTAAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGTTTCAGAAGAATGCtaataaaaaaagcccggcgccatgcagggcaaaagtcacaagttacgtttatgccactgtcaactgctgaatttttttctcgaggcggattttgagttcttccatt                                                                                                           >  2:429461‑M1/1‑144 (MQ=255)                                                                                                                                                         aTACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgtttatgccactgtcaactgctgaatttttttctcgcggcggattttgcgttcttccattaccgccactatcgccatcaggcagatacaaccaatcgcggcgatatccagcgcgg                                                    >  1:65059‑M1/1‑90 (MQ=255)                                                                                                                                                         aTACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgtttatgccactgtcaactgctgaatttttttctcgcggcggattttgcgttcttccattaccgccactatcgccatcaggcagatacaaccaatcgcggcgatatccagcgcgg                                                    >  1:461264‑M1/1‑90 (MQ=255)                                                                                                                                                                                         gACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgtttatgccactgtcaactgctgaatttttttctcgcggcggattttgcgttcttccattaccgccactatcgccatcaggcagatacaaccaatcgcggcgatatccagcgcggcgaaggtgcctgcccagccggtaaggccgaat                    >  1:340186‑M1/1‑58 (MQ=255)                                                                                                                                                                                         gACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgtttatgccactgtcaactgctgaatttttttctcgcggcggattttgcgttcttccattaccgccactatcgccatcaggcagatacaaccaatcgcggcgatatccagcgcggcgaaggtgcctgcccagccggtaaggccgaa                     >  1:98566‑M1/1‑58 (MQ=255)                                                                                                                                                                                                        gTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGACTGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgtttatgccactgtcaactgctgaatttttttctcgcggcggat                                                                                                                           >  2:75511‑M1/1‑43 (MQ=255)                                                                                                                                                                                                        gTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgcctggcaaaagtcaccagttacgtttatgccactgtcaactgctgaatttttttctcgcgtcggat                                                                                                                           <  1:75511‑M1/133‑91 (MQ=255)                                                                                                                                                                                                           aCTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgtttatgccactgtcaactgctgaatttttttctcgcggcggattttgcgttcttccattaccgccactatcgccatcaggcagatacaaccaatcgcggcgatatccagcgcggcgaaggtgcctgcccagccggtaaggccgaataccggcgtcccatcggc   >  1:317011‑M1/1‑40 (MQ=255)                                                                                                                                                                                                            cTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCtaataaaaaaagcccggcgtcatgccgggcaaaagtcaccagttacgtttatgccactgtcaactgctgaatttttttctcgcggcggattttgcgttcttccattaccgccactatcgccatcaggcagatacaaccaatcgcggcgatatccagcgcggcgaaggtgcctgcccagccggtaaggccgaataccggcgtcccatcggca  >  1:223664‑M1/1‑39 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                    AACTTGTTTTATTTTTACTATCAAAAATAGCTCATTGCCTTACTATAATTAATATATAAATTCGGTAATTAATTCTTAACATGCTTTTACTAATAATCTCAATTGCTGGCCCTATAATATTGCGCTAGCATTGCTTTCTGGTTGTATCAGCGATACTCAAAAATTCTTAATACAATACTCATTCGACTGGTACTTATTTGTAACTCAGTTATATTTTTTCGCCCGGTGATTCAGAAGAATGCAAAAAACGGCTACCACTCCATCAAAAATACTTGATCTCACTGCCGCGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGAACTGAAAGCCCGTCCTATAATGGTAGGTTTTTTCTTCACCGGTAGCGCTATTATGGGAATTCTGGTCAGTCAATTTCTGGCAAG  >  NC_000913/3836715‑3837167

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