Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	1	1,010,036	Δ641 bp	coding (293‑933/1236 nt)	unknown →	Type III restriction‑modification system restriction subunit (EC 3.1.21.5)

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	1	1010036	1010676	641	41 [0]	[0] 44	unknown	Type III restriction‑modification system restriction subunit (EC 3.1.21.5)

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	1	= 1010035	0 (0.000)	41 (0.490)	28/426	1.8	100%	coding (292/1236 nt)	unknown	Type III restriction‑modification system restriction subunit (EC 3.1.21.5)
	?	1	1010677 =	0 (0.000)					coding (934/1236 nt)	unknown	Type III restriction‑modification system restriction subunit (EC 3.1.21.5)

GAAAAATATCGCATTGCGTTTGATTAATCAGCTGATTGATGATAACGATCATTTCGTCAGCAGTGCCGCTTATCAACAGGTGAAACAGCTGTTCCCGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTATGCGTGAAGGCAAATATGAAGAGTTAAAAGCCTTGTGGGAGCTGATTCATCACAAAGCGATTTTGCAGTATAAAATTGATTCCGAAGCCGACTTTTTAGCCTTATTTGTCGTTTATTTGCAGGAAAATGCAGATAAATTCAAACAAGCAGGCATTCGCACAAAAGTGTCGGAAACCTACATAAAAAATAACGTGATGCTGACCCGAGACAAGGAAAATCTAGACGATACGCCGATTCGTTTTAACA  >  1/1009869‑1010281                                                                                                                                                                       |                                                                                                                                                                                                                                                       tagtattaccaataaattggcggataatcagctgattgatgataacgatCATTTCGTCAGCAGTGCCGCTTATCAACAGGTGAAACAGCTGTTCCCGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatc                                                                                                                                                                    <  1:272779‑M1/227‑85 (MQ=255) tagtattaccaataaattggcggataatcagctgattgatgataacgatCATTTCGTCAGCAGTGCCGCTTATCAACAGGTGAAACAGCTGTTCCCGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatat                                                                                                                                                                     <  1:490172‑M1/226‑84 (MQ=255) tagtattaccaataaattggcggataatcagctgattgatgataacgatCATTTCGTCAGCAGTGCCGCTTATCAACAGGTGAAACAGCTGTTCCCGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcacgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatc                                                                                                                                                                    <  1:371442‑M1/227‑85 (MQ=255)         ccaataaattggcggataatcagctgattgatgATAACGATCATTTCGTCAGCAGTGCCGCTTATCAACAGGTGAAACAGCTGTTCCCGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaa                                                                                                                                                             >  1:403156‑M1/17‑159 (MQ=255)              aaattggcggataatcagctgatTGATGATAACGATCATTTCGTCAGCAGTGCCGCTTATCAACAGGTGAAACAGCTGTTCCCGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaa                                                                                                                                                       >  1:186513‑M1/12‑154 (MQ=255)               aattggcggataatcagctgaTTGATGATAACGATCATTTCGTCAGCAGTGCCGCTTATCAACAGGTGAAACAGCTGTTCCCGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaac                                                                                                                                                      <  1:735720‑M1/241‑99 (MQ=255)                       gataaTCAGCTGATTGATGATAACGATCATTTCGTCAGCAGTGCCGCTTATCAACAGGTGAAACAGCTGTTCCCGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgaga                                                                                                                                              >  1:7521‑M1/3‑145 (MQ=255)                       gataaTCAGCTGATTGATGATAACGATCATTTCGTCAGCAGTGCCGCTTATCAACAGGTGAAACAGCTGTTCCCGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaa                                                                                                                                                                                                                                          >  1:595683‑M1/3‑145 (MQ=255)                                                             cagTGCCGCTTATCAACAGGTGAAACAGCTGTTCCCGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatg                                                                                                          <  1:531151‑M1/249‑143 (MQ=255)                                                             cagTGCCGCTTATCAACAGGTGAAACAGCTGTTCCCGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggat                                                                                                           <  1:663044‑M1/248‑142 (MQ=255)                                                                        aTCAACAGGTGAAACAGCTGTTCCCGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataactt                                                                                             >  1:37742‑M1/1‑96 (MQ=255)                                                                         tCAACAGGTGAAACAGCTGTTCCCGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataactta                                                                                            <  1:901634‑M1/251‑157 (MQ=255)                                                                              aGGTGAAACAGCTGTTCCCGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaa                                                                                       >  1:364251‑M1/1‑90 (MQ=255)                                                                                 tGAAACAGCTGTTCCCGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacgtattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaagg                                                                                    <  1:71499‑M1/251‑165 (MQ=255)                                                                                              cccGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaa                                                                       >  1:813174‑M1/1‑74 (MQ=255)                                                                                                      ccTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaaca                                                               >  1:786107‑M1/1‑66 (MQ=255)                                                                                                                 gTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaacacgccgaacgtt                                                    >  1:117197‑M1/1‑55 (MQ=255)                                                                                                                 gTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaacacgccgaacgtt                                                    >  1:176467‑M1/1‑55 (MQ=255)                                                                                                                 gTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaa                                                                                        >  1:185152‑M1/1‑55 (MQ=255)                                                                                                                 gTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggagg                                                                                                                                                                                                   >  1:640665‑M1/1‑55 (MQ=255)                                                                                                                      gCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatc                                                                    >  1:626861‑M1/1‑50 (MQ=255)                                                                                                                          aaaGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaacacgccgaacgtttatttgccc                                           <  1:733705‑M1/251‑206 (MQ=255)                                                                                                                             ggCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatat                                                                                                                                                                     >  1:900655‑M1/1‑43 (MQ=255)                                                                                                                                      aCTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaacacgccgaacgtttatttgcccaaatctcaaaag                               <  1:162789‑M1/251‑218 (MQ=255)                                                                                                                                         gATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaa                                                                                       >  1:620831‑M1/1‑31 (MQ=255)                                                                                                                                               aaaaCCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaac                                                                >  1:629056‑M1/1‑25 (MQ=255)                                                                                                                                                 aaCCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaacacgccgaacgtttatttgcccaaatctcaaaagaggtgaaagt                     >  1:75313‑M1/1‑23 (MQ=255)                                                                                                                                                 aaCCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaacacgccgaacgtttatttgcccaaatctcaaaagaggtgaaagt                     >  1:571232‑M1/1‑23 (MQ=255)                                                                                                                                                  aCCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaacacgccgaacgtttatttgcccaaatctcaaaagaggtgaaagtcg                   >  1:605010‑M1/1‑22 (MQ=255)                                                                                                                                                  aCCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaacacgccgaacgtttatttgcccaaatctcaaaagaggtgaaagtcg                   >  1:230783‑M1/1‑22 (MQ=255)                                                                                                                                                  aCCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaacacgccgaacgtttatttgcccaaatctcaaaagaggtgaaagtcg                   >  1:94995‑M1/1‑22 (MQ=255)                                                                                                                                                            aTGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaacacgccgaacgtttatttgcccaaatctcacaagaggtgaaagtcgtgtttaaaa          >  1:858713‑M1/1‑12 (MQ=255)                                                                                                                                                            aTGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataactt                                                                                             <  1:410335‑M1/167‑156 (MQ=255)                                                                                                                                                            aTGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttc                                                                                                      <  1:355201‑M1/158‑147 (MQ=255)                                                                                                                                                                 aTTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaacacgccgaacgtttatttgcccaaatctcaaaagaggtgaaagtcgtgtttaaaactcagt    >  1:865805‑M1/1‑7 (MQ=255)                                                                                                                                                                 aTTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaacacgccgaacgtt                                                    <  1:212968‑M1/203‑197 (MQ=255)                                                                                                                                                                 aTTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaacacgc                                                            <  1:612904‑M1/195‑189 (MQ=255)                                                                                                                                                                 aTTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaat                                                                                                                                      <  1:896995‑M1/121‑115 (MQ=255)                                                                                                                                                                  ttCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaagacgaagaacgcaaaatcaaacac                                                              >  1:446118‑M1/1‑6 (MQ=255)                                                                                                                                                                   tCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaacacgccgaacgtttatttgcccaaatctcaaaagaggtgaaagtcgtgtttaaaactcagttt  >  1:492655‑M1/1‑5 (MQ=255)                                                                                                                                                                   tCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaacacgccgaacgtttatttgcccaaatctcaaaagaggtgaaagtcgtgtttaaaactcagttt  <  1:485454‑M1/251‑247 (MQ=255)                                                                                                                                                                       |                                                                                                                                                                                                                                                       GAAAAATATCGCATTGCGTTTGATTAATCAGCTGATTGATGATAACGATCATTTCGTCAGCAGTGCCGCTTATCAACAGGTGAAACAGCTGTTCCCGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTATGCGTGAAGGCAAATATGAAGAGTTAAAAGCCTTGTGGGAGCTGATTCATCACAAAGCGATTTTGCAGTATAAAATTGATTCCGAAGCCGACTTTTTAGCCTTATTTGTCGTTTATTTGCAGGAAAATGCAGATAAATTCAAACAAGCAGGCATTCGCACAAAAGTGTCGGAAACCTACATAAAAAATAACGTGATGCTGACCCGAGACAAGGAAAATCTAGACGATACGCCGATTCGTTTTAACA  >  1/1009869‑1010281

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