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)

CGCAATTTACCGATAAAAACGGCAAAGCAAAAGCGGTGCGCGTGGCTGATTTGGGACGAATGGCAGATGAAATCCATGAACCAGCGGCAAGTTATTTGTTTGGCGAAATTTTTTATGATTCAGAAATTGAACATCAAAACATCGCCAGCGAAGAAATTGAAGGTGTAACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTACCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGG  >  1/1010510‑1010927                                                                                                                                                                        |                                                                                                                                                                                                                                                           tagtattaccaataaattggcggataatcagctgattgatgataacgatcatttcgtcagcagtgccgcttatcaacaggtgaaacagctgttcccgaatgcctttgccgcagtcttgcataaaggcaaaatcactgatgcgaaaaccaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTatat                                                                                                                                                                          <  1:490172‑M2/83‑1 (MQ=255) tagtattaccaataaattggcggataatcagctgattgatgataacgatcatttcgtcagcagtgccgcttatcaacaggtgaaacagctgttcccgaatgcctttgccgcagtcttgcataaaggcaaaatcactgatgcgaaaaccaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATc                                                                                                                                                                         <  1:272779‑M2/84‑1 (MQ=255) tagtattaccaataaattggcggataatcagctgattgatgataacgatcatttcgtcagcagtgccgcttatcaacaggtgaaacagctgttcccgaatgcctttgccgcagtcttgcataaaggcaaaatcactgatgcgaaaaccaaaggcgatgttattcctaACCGTGTTTACTAAAATCACGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATc                                                                                                                                                                         <  1:371442‑M2/84‑1 (MQ=255)         ccaataaattggcggataatcagctgattgatgataacgatcatttcgtcagcagtgccgcttatcaacaggtgaaacagctgttcccgaatgcctttgccgcagtcttgcataaaggcaaaatcactgatgcgaaaaccaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTcaaaa                                                                                                                                                                  >  1:403156‑M2/160‑250 (MQ=255)              aaattggcggataatcagctgattgatgataacgatcatttcgtcagcagtgccgcttatcaacaggtgaaacagctgttcccgaatgcctttgccgcagtcttgcataaaggcaaaatcactgatgcgaaaaccaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTcaaaaccaaaa                                                                                                                                                            >  1:186513‑M2/155‑251 (MQ=255)               aattggcggataatcagctgattgatgataacgatcatttcgtcagcagtgccgcttatcaacaggtgaaacagctgttcccgaatgcctttgccgcagtcttgcataaaggcaaaatcactgatgcgaaaaccaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTcaaaaccaaaac                                                                                                                                                           <  1:735720‑M2/98‑1 (MQ=255)                       gataatcagctgattgatgataacgatcatttcgtcagcagtgccgcttatcaacaggtgaaacagctgttcccgaatgcctttgccgcagtcttgcataaaggcaaaatcactgatgcgaaaaccaaaggcgatgttattcctaACCGTGTTTACTaa                                                                                                                                                                                                                                               >  1:595683‑M2/146‑159 (MQ=255)                       gataatcagctgattgatgataacgatcatttcgtcagcagtgccgcttatcaacaggtgaaacagctgttcccgaatgcctttgccgcagtcttgcataaaggcaaaatcactgatgcgaaaaccaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCgaga                                                                                                                                                   >  1:7521‑M2/146‑251 (MQ=255)                                                             cagtgccgcttatcaacaggtgaaacagctgttcccgaatgcctttgccgcagtcttgcataaaggcaaaatcactgatgcgaaaaccaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGAt                                                                                                                <  1:663044‑M2/141‑1 (MQ=255)                                                             cagtgccgcttatcaacaggtgaaacagctgttcccgaatgcctttgccgcagtcttgcataaaggcaaaatcactgatgcgaaaaccaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATg                                                                                                               <  1:531151‑M2/142‑1 (MQ=255)                                                                        atcaacaggtgaaacagctgttcccgaatgcctttgccgcagtcttgcataaaggcaaaatcactgatgcgaaaaccaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACtt                                                                                                  >  1:37742‑M2/97‑251 (MQ=255)                                                                         tcaacaggtgaaacagctgttcccgaatgcctttgccgcagtcttgcataaaggcaaaatcactgatgcgaaaaccaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTa                                                                                                 <  1:901634‑M2/156‑1 (MQ=255)                                                                              aggtgaaacagctgttcccgaatgcctttgccgcagtcttgcataaaggcaaaatcactgatgcgaaaaccaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGaaa                                                                                            >  1:364251‑M2/91‑251 (MQ=255)                                                                                 tgaaacagctgttcccgaatgcctttgccgcagtcttgcataaaggcaaaatcactgatgcgaaaaccaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACGTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAgg                                                                                         <  1:71499‑M2/164‑1 (MQ=255)                                                                                              cccgaatgcctttgccgcagtcttgcataaaggcaaaatcactgatgcgaaaaccaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCaaa                                                                            >  1:813174‑M2/75‑251 (MQ=255)                                                                                                      cctttgccgcagtcttgcataaaggcaaaatcactgatgcgaaaaccaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAaca                                                                    >  1:786107‑M2/67‑251 (MQ=255)                                                                                                                 gtcttgcataaaggcaaaatcactgatgcgaaaaccaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCaggagg                                                                                                                                                                                                        >  1:640665‑M2/56‑108 (MQ=255)                                                                                                                 gtcttgcataaaggcaaaatcactgatgcgaaaaccaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGaa                                                                                             >  1:185152‑M2/56‑215 (MQ=255)                                                                                                                 gtcttgcataaaggcaaaatcactgatgcgaaaaccaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGtt                                                         >  1:176467‑M2/56‑251 (MQ=255)                                                                                                                 gtcttgcataaaggcaaaatcactgatgcgaaaaccaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGtt                                                         >  1:117197‑M2/56‑251 (MQ=255)                                                                                                                      gcataaaggcaaaatcactgatgcgaaaaccaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATc                                                                         >  1:626861‑M2/51‑230 (MQ=255)                                                                                                                          aaaggcaaaatcactgatgcgaaaaccaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGccc                                                <  1:733705‑M2/205‑1 (MQ=255)                                                                                                                             ggcaaaatcactgatgcgaaaaccaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTatat                                                                                                                                                                          >  1:900655‑M2/44‑126 (MQ=255)                                                                                                                                      actgatgcgaaaaccaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAag                                    <  1:162789‑M2/217‑1 (MQ=255)                                                                                                                                         gatgcgaaaaccaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGaaa                                                                                            >  1:620831‑M2/32‑192 (MQ=255)                                                                                                                                               aaaaccaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAac                                                                     >  1:629056‑M2/26‑209 (MQ=255)                                                                                                                                                 aaccaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGt                          >  1:75313‑M2/24‑250 (MQ=255)                                                                                                                                                 aaccaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGt                          >  1:571232‑M2/24‑250 (MQ=255)                                                                                                                                                  accaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCg                        >  1:605010‑M2/23‑251 (MQ=255)                                                                                                                                                  accaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCg                        >  1:230783‑M2/23‑251 (MQ=255)                                                                                                                                                  accaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCg                        >  1:94995‑M2/23‑251 (MQ=255)                                                                                                                                                            atgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTc                                                                                                           <  1:355201‑M2/146‑1 (MQ=255)                                                                                                                                                            atgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACtt                                                                                                  <  1:410335‑M2/155‑1 (MQ=255)                                                                                                                                                            atgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCACAAGAGGTGAAAGTCGTGTTTaaaa               >  1:858713‑M2/13‑250 (MQ=255)                                                                                                                                                                 attcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAAt                                                                                                                                           <  1:896995‑M2/114‑1 (MQ=255)                                                                                                                                                                 attcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGc                                                                 <  1:612904‑M2/188‑1 (MQ=255)                                                                                                                                                                 attcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGtt                                                         <  1:212968‑M2/196‑1 (MQ=255)                                                                                                                                                                 attcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGt         >  1:865805‑M2/8‑251 (MQ=255)                                                                                                                                                                  ttcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGACGAAGAACGCAAAATCAAacac                                                                   >  1:446118‑M2/7‑192 (MQ=255)                                                                                                                                                                   tcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGttt       >  1:492655‑M2/6‑251 (MQ=255)                                                                                                                                                                   tcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGttt       <  1:485454‑M2/246‑1 (MQ=255)                                                                                                                                                                        aCCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAggg  <  1:554768/251‑1 (MQ=255)                                                                                                                                                                        aCCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAggg  <  1:482110/251‑1 (MQ=255)                                                                                                                                                                        aCCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAggg  <  1:253462/251‑1 (MQ=255)                                                                                                                                                                        |                                                                                                                                                                                                                                                           CGCAATTTACCGATAAAAACGGCAAAGCAAAAGCGGTGCGCGTGGCTGATTTGGGACGAATGGCAGATGAAATCCATGAACCAGCGGCAAGTTATTTGTTTGGCGAAATTTTTTATGATTCAGAAATTGAACATCAAAACATCGCCAGCGAAGAAATTGAAGGTGTAACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTACCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGG  >  1/1010510‑1010927

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