Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	1	91,658	Δ4 bp	coding (750‑753/831 nt)	unknown ←	hypothetical protein

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	1	91658	91661	4	48 [0]	[0] 48	unknown	hypothetical protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	1	= 91657	0 (0.000)	48 (0.540)	26/430	2.0	100%	coding (754/831 nt)	unknown	hypothetical protein
	?	1	91662 =	0 (0.000)					coding (749/831 nt)	unknown	hypothetical protein

TCCCTCTTGCTTACGTGTCTAATGTTATATTATATAATATAACCGCGCAACCGCTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTGTCAAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTTTTGGTTGACGAGAAGAAGATATACGGCGGGCGGGTGATGTTGACTAAGCGCAGGAAG  >  1/91462‑91781                                                                                                                                                                                                    |                                                                                                                             tCCCTCTTGCTTACGTGTCTAATATTATATTATATAATATAACTGCGCAACCGGTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaataaagcgcacaa                                                                       >  1:519533‑M1/1‑196 (MQ=255)      ctTGCTTACGTGTCTAATATTATATTATATAATATAACTGCGCAACCGGTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaataaagcgcacaaactcg                                                                  >  1:807424‑M1/1‑191 (MQ=255)              cGTGTCTAATATTATATTATATAATATAACTGCGCAACCGGTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaataaagcgcacaaactcgctc                                                               <  1:218500‑M1/246‑64 (MQ=255)                 gtCTAATATTATATTATATAATATAACTGCGCAACCGGTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaataaagcgcacaaactcgctcttagttga                                                       >  1:248287‑M1/1‑180 (MQ=255)                      atattatattatataaTATAACTGCGCAACCGGTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaataaagcgcacaaactcgctcttagttgacgaga                                                  >  1:117604‑M1/4‑175 (MQ=255)                      atattatattatataaTATAACTGCGCAACCGGTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaataaagcgcacaaactcgctcttagttgacgaga                                                  >  1:73114‑M1/4‑175 (MQ=255)                      atattatattatataaTATAACTGCGCAACCGGTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaataaagcgcacaaactcgctcttagttgacgaga                                                  >  1:504031‑M1/4‑175 (MQ=255)                      atattatattatataaTATAACTGCGCAACCGGTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaat                                                                                  <  1:112826‑M1/216‑45 (MQ=255)                                atataatataACTGCGCAACCGGTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaataaagcgcacaaactcgctcttagttgacgagaagaatatata                                        >  1:822762‑M1/1‑165 (MQ=255)                                atataatataACTGCGCAACCGGTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaataaagcgcacaaactcgctcttagttgacgagaagaatatata                                        >  1:808279‑M1/1‑165 (MQ=255)                                atataatataACTGCGCAACCGGTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaataaagcgcacaaactcgctcttagttgacgagaagaatatata                                        >  1:767160‑M1/1‑165 (MQ=255)                                atataatataACTGCGCAACCGGTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaataaagcgcacaaactcgctcttagttgacgagaagaatat                                           >  1:278575‑M1/1‑165 (MQ=255)                                atataatataACTGCGCAACCGGTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaataa                                                                                >  1:10627‑M1/1‑165 (MQ=255)                                atataatataACTGCGCAACCGGTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatac                                                                                       <  1:812665‑M1/204‑40 (MQ=255)                                atataatataACTGCGCAACCGGTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctac                                                                                                          <  1:617446‑M1/185‑21 (MQ=255)                                 tataatataaCTGCGCAACCGGTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaataaagcgcacaaactcgctcttagttgacgagaagaatatatac                                       <  1:140861‑M1/251‑88 (MQ=255)                                   taatataaCTGCGCAACCGGTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaataaagcgcacaaactcgctcttagttgacgagaagaatatatacgg                                     >  1:895333‑M1/1‑162 (MQ=255)                                   taatataaCTGCGCAACCGGTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaataaagcgcacaaactcgctcttagttgacgagaagaatatata                                        >  1:41569‑M1/1‑162 (MQ=255)                                   taatataaCTGCGCAACCGGTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaataaagcgcacaaactcgctcttagttgacgagaagaatatata                                        <  1:628600‑M1/248‑87 (MQ=255)                                   taatataaCTGCGCAACCGGTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccg                                                                                                  >  1:267486‑M1/1‑162 (MQ=255)                                   taatataaCTGCGCAACCGGTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgcc                                                                                                                    <  1:601220‑M1/172‑11 (MQ=255)                                   taatataaCTGCGCAACCGGTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgc                                                                                                                         <  1:178371‑M1/167‑6 (MQ=255)                                     atataaCTGCGCAACCGGTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaataaagcgcacaaactcgctcttagttgacgagaagaatatatacggcg                                   >  1:809348‑M1/1‑160 (MQ=255)                                     atataaCTGCGCAACCGGTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaataaagcgcacaaactcgctcttagttgacgagaagaatatata                                        <  1:553053‑M1/246‑87 (MQ=255)                                     atataaCTGCGCAACCGGTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatac                                                                                       >  1:622383‑M1/1‑160 (MQ=255)                                        taaCTGCGCAACCGGTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaataaagcgcacaaactcgctcttagttgacgagaagaatatatac                                       >  1:507963‑M1/1‑157 (MQ=255)                                           ctgcgcaaCCGGTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaataaagcgcacaaactcgctcttagttgacgagaagaatatatac                                       <  1:516350‑M1/239‑88 (MQ=255)                                                 aaCCGGTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaataaagcgcacaaactcgctcttagttgacgagaagaatatatacggcgggcgggtgatgt                       >  1:763292‑M1/1‑148 (MQ=255)                                                  aCCGGTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccg                                                                                                  <  1:498922‑M1/175‑29 (MQ=255)                                                   ccggtatttTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatact                                                                                      <  1:746228‑M1/182‑41 (MQ=255)                                                      gtaTTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgc                                                                                                                         <  1:373734‑M1/147‑6 (MQ=255)                                                          tttGTAAGGTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatac                                                                                       <  1:452135‑M1/178‑40 (MQ=255)                                                                        gtAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaataaagcgcacaaactcgctcttagttgacgagaagaatatata                                        >  1:6909‑M1/1‑125 (MQ=255)                                                                        gtAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactt                                                                                                        <  1:845001‑M1/147‑23 (MQ=255)                                                                                       ccccTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaataaagcgcacaaactcgctcttagttgacgagaagaatatatacggcgggcgggtgatgttgactagcaattctttgacta  >  1:730979‑M1/1‑110 (MQ=255)                                                                                       ccccTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatac                                                                                       <  1:649944‑M1/149‑40 (MQ=255)                                                                                       ccccTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtat                                                                                           <  1:344845‑M1/145‑36 (MQ=255)                                                                                        cccTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaataaagcgcacaaactcgctcttagttgacgagaagaatatatacggcgggcgggtgatgttgact                  <  1:372950‑M1/217‑109 (MQ=255)                                                                                        cccTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaataaagcgcacaaactcgctcttagttgacgagaagaatatatac                                       <  1:261590‑M1/196‑88 (MQ=255)                                                                                        cccTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactc                                                                                     >  1:497122‑M1/1‑109 (MQ=255)                                                                                         ccTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaataaagcgcacaaactc                                                                   <  1:245777‑M1/167‑60 (MQ=255)                                                                                          cTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatac                                                                                       >  1:77777‑M1/1‑107 (MQ=255)                                                                                                                                       aGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaataaagcgcacaaactcgctcttagttgacgagaagaatatatac                                       <  1:881290‑M1/149‑88 (MQ=255)                                                                                                                                             tCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaataaagcgcacaaactcgctcttagttgacgagaagaatatatacg                                      >  1:172254‑M1/1‑56 (MQ=255)                                                                                                                                             tCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaataaagcgcacaaactcgctcttagttgacgagaagaatatata                                        >  1:146535‑M1/1‑56 (MQ=255)                                                                                                                                             tCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaataa                                                                                >  1:569545‑M1/1‑56 (MQ=255)                                                                                                                                                         ccATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaataaagcgcacaaactcgctcttagttgacgagaagaatatatac                                       <  1:290215‑M1/131‑88 (MQ=255)                                                                                                                                                               tAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTaacgcttgccagttatctactttatccgtttgtatatactcaataaagcgcacaaactcgctcttagttgacgagaagaatatatac                                       <  1:896805‑M1/125‑88 (MQ=255)                                                                                                                                                                                                    |                                                                                                                             TCCCTCTTGCTTACGTGTCTAATGTTATATTATATAATATAACCGCGCAACCGCTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTGTCAAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTTTTGGTTGACGAGAAGAAGATATACGGCGGGCGGGTGATGTTGACTAAGCGCAGGAAG  >  1/91462‑91781

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