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

TCTTGCTTACGTGTCTAATGTTATATTATATAATATAACCGCGCAACCGCTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTGTCAAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTTTTGGTTGACGAGAAGAAGATATACGGCGGGCGGGTGATGTTGACTAAGCGCAGGAAGTCGA  >  1/91466‑91785                                                                                                                                                                                                     |                                                                                                                            tccctcttgcttacgtgtctaatattatattatataatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACaa                                                                       >  1:519533‑M2/197‑251 (MQ=255)      cttgcttacgtgtctaatattatattatataatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCg                                                                  >  1:807424‑M2/192‑251 (MQ=255)              cgtgtctaatattatattatataatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGctc                                                               <  1:218500‑M2/63‑2 (MQ=255)                 gtctaatattatattatataatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGa                                                       >  1:248287‑M2/181‑251 (MQ=255)                      atattatattatataatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAAt                                                                                  <  1:112826‑M2/44‑1 (MQ=255)                      atattatattatataatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACgaga                                                  >  1:73114‑M2/176‑251 (MQ=255)                      atattatattatataatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACgaga                                                  >  1:117604‑M2/176‑251 (MQ=255)                      atattatattatataatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACgaga                                                  >  1:504031‑M2/176‑251 (MQ=255)                                atataatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTAc                                                                                                          <  1:617446‑M2/20‑1 (MQ=255)                                atataatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATAc                                                                                       <  1:812665‑M2/39‑1 (MQ=255)                                atataatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATaa                                                                                >  1:10627‑M2/166‑211 (MQ=255)                                atataatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGacgagaagaatat                                           >  1:278575‑M2/166‑245 (MQ=255)                                atataatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAAtatata                                        >  1:822762‑M2/166‑251 (MQ=255)                                atataatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAAtatata                                        >  1:808279‑M2/166‑251 (MQ=255)                                atataatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAAtatata                                        >  1:767160‑M2/166‑251 (MQ=255)                                 tataatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAATATATAc                                       <  1:140861‑M2/87‑1 (MQ=255)                                   taatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGc                                                                                                                         <  1:178371‑M2/5‑1 (MQ=255)                                   taatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGcc                                                                                                                    <  1:601220‑M2/10‑1 (MQ=255)                                   taatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCg                                                                                                  >  1:267486‑M2/163‑190 (MQ=255)                                   taatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAAtatata                                        >  1:41569‑M2/163‑248 (MQ=255)                                   taatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAAtatata                                        <  1:628600‑M2/86‑1 (MQ=255)                                   taatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAATATATAcgg                                     >  1:895333‑M2/163‑251 (MQ=255)                                     atataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATAc                                                                                       >  1:622383‑M2/161‑199 (MQ=255)                                     atataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAAtatata                                        <  1:553053‑M2/86‑1 (MQ=255)                                     atataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAATATATAcggcg                                   >  1:809348‑M2/161‑251 (MQ=255)                                        taactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAATATATAc                                       >  1:507963‑M2/158‑244 (MQ=255)                                           ctgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAATATATAc                                       <  1:516350‑M2/87‑1 (MQ=255)                                                 aaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAATATATACGGCGGGCGGGTGATGt                       >  1:763292‑M2/149‑251 (MQ=255)                                                  accggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCg                                                                                                  <  1:498922‑M2/28‑1 (MQ=255)                                                   ccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACt                                                                                      <  1:746228‑M2/40‑1 (MQ=255)                                                      gtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGc                                                                                                                         <  1:373734‑M2/5‑1 (MQ=255)                                                          tttgtaaggtttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATAc                                                                                       <  1:452135‑M2/39‑1 (MQ=255)                                                                        gtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACtt                                                                                                        <  1:845001‑M2/22‑1 (MQ=255)                                                                        gtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAAtatata                                        >  1:6909‑M2/126‑211 (MQ=255)                                                                                       ccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGtat                                                                                           <  1:344845‑M2/35‑1 (MQ=255)                                                                                       ccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATAc                                                                                       <  1:649944‑M2/39‑1 (MQ=255)                                                                                       ccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAATATATACGGCGGGcgggtgatgttgactagcaattctttgacta  >  1:730979‑M2/111‑219 (MQ=255)                                                                                        cccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTc                                                                                     >  1:497122‑M2/110‑150 (MQ=255)                                                                                        cccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAATATATAc                                       <  1:261590‑M2/87‑1 (MQ=255)                                                                                        cccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAATATATACGGCGGGCGGGTGATGTTGACt                  <  1:372950‑M2/108‑1 (MQ=255)                                                                                         ccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTc                                                                   <  1:245777‑M2/59‑1 (MQ=255)                                                                                          cttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATAc                                                                                       >  1:77777‑M2/108‑146 (MQ=255)                                                                                                                                       agattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAATATATAc                                       <  1:881290‑M2/87‑1 (MQ=255)                                                                                                                                             tcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATaa                                                                                >  1:569545‑M2/57‑102 (MQ=255)                                                                                                                                             tcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAAtatata                                        >  1:146535‑M2/57‑142 (MQ=255)                                                                                                                                             tcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAATATATAcg                                      >  1:172254‑M2/57‑144 (MQ=255)                                                                                                                                                         ccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAATATATAc                                       <  1:290215‑M2/87‑1 (MQ=255)                                                                                                                                                               tagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAATATATAc                                       <  1:896805‑M2/87‑1 (MQ=255)                                                                                                                                                                                                     |                                                                                                                            TCTTGCTTACGTGTCTAATGTTATATTATATAATATAACCGCGCAACCGCTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTGTCAAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTTTTGGTTGACGAGAAGAAGATATACGGCGGGCGGGTGATGTTGACTAAGCGCAGGAAGTCGA  >  1/91466‑91785

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