Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	257,908	Δ776 bp		insB9–[crl]	insB9, insA9, [crl]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	257908–258676	258683	8–776	42 [0]	[3] 45	insB9–[crl]	insB9,insA9,[crl]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 257907	0 (0.000)	42 (0.750)	37/474	0.5	100%	intergenic (+8/+16)	crl/insB9	RNA polymerase holoenzyme assembly factor Crl/IS1 family transposase B
	?	NC_000913	258684 =	0 (0.000)					pseudogene (9/331 nt)	crl	RNA polymerase holoenzyme assembly factor Crl

CTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCACCACGTC  >  NC_000913/257668‑258147                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                 ctgctgACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATAGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAACGAGTGGACACCCGAAGGGCAGATTGCTCAAAAAATTTACCGCACTAGGCCCGCATATTCGTGAAGGTAagtgcaa                                                                                                                                                                                                                                           >  1:189893‑M1/1‑240 (MQ=255)      gACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATAACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGAGAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatc                                                                                                                                                                                                                                   <  2:323915‑M1/250‑16 (MQ=255)              attattAGAGATCCCATTTAACCCGGTGTATCGGCATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGTGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGACGGTAagtgcaaagatactcgattctttt                                                                                                                                                                                                                          >  1:226503‑M1/1‑227 (MQ=255)                        aTCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtct                                                                                                                                                                                                                <  1:572442‑M1/251‑35 (MQ=255)                         tCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctg                                                                                                                                                                                                               >  2:404134‑M1/1‑216 (MQ=255)                          cccATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctgg                                                                                                                                                                                                              <  1:501775‑M1/251‑37 (MQ=255)                                    ccGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGAGCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCTCTAGGCCCGTAGATTCGTGAAGGTAagtgcaaagataatcgagtctttttcgattgtctggctgtatgcgt                                                                                                                                                                                                    <  2:35458‑M1/251‑47 (MQ=255)                                                                      aaaTCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattc                                                                                                                                                                  >  1:503093‑M1/1‑171 (MQ=255)                                                                           aCCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCTAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtcgggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctgggg                                                                                                                                                             <  1:578760‑M1/251‑86 (MQ=255)                                                                               aCTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctgg                                                                                                                                                         <  1:208599‑M1/251‑90 (MQ=255)                                                                                       gAAAAACGCTTGTGTTAAAAAGTTGCTAAATTTTGCCAATTTGGTAAAACAGTGGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGTACACCCGAAGAGCAGATTGATCAAAAAATTTGCCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtctacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatgga                                                                                                                                                 <  1:520928‑M1/251‑98 (MQ=255)                                                                                                    gtTAAAAATTTGCTACATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgc                                                                                                                                      >  1:140153‑M1/1‑141 (MQ=255)                                                                                                      tAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtg                                                                                                                                                       >  2:70765‑M1/1‑139 (MQ=255)                                                                                                      tAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtg                                                                                                                                                       <  1:70765‑M1/230‑92 (MQ=255)                                                                                                            tttGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccg                                                                                                                            >  1:423613‑M1/1‑133 (MQ=255)                                                                                                               gCTAAATTTTGCCAATGTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctggttcgattgtctggctgtgggcggcaacgtgaaacgggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttt                                                                                                                         <  2:218732‑M1/251‑122 (MQ=255)                                                                                                                cTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCAGAAGAGCAGATTGATCAAAAAATTTAACGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatagattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatgtagctggacgcgcaggaatcccgtttt                                                                                                                        >  1:573863‑M1/1‑129 (MQ=255)                                                                                                                  aaaTTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgtttt                                                                                                                        <  1:217974‑M1/249‑123 (MQ=255)                                                                                                                        tGCCAATTTGGTAAAACAGCTACGTCACATCAGCAGATAGCAATCACGTTAACGAGTGGACACCCGACGAGCAGATTGATCAACAAATTTACAGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattatttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggggctggaagcgcaggactgccgttttacctacgg                                                                                                                >  2:42824‑M1/1‑121 (MQ=255)                                                                                                                                 ggTAAACCAGTTGCATCACAACAGGAGATAGCAATGACGATGCCGAGTGGACACCCGAAGAGTAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattcttttttgattctctggctgtttgcggcatcgtgaaaccggcgccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagt                                                                                                        <  1:351299‑M1/250‑139 (MQ=255)                                                                                                                                            tGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATGCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttc                                                                                             <  2:536453‑M1/250‑150 (MQ=255)                                                                                                                                                 cacaACAGGAGATAGCAAAGACCGTAACGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTTGTCCCGTATATTCGTGAAGGTAagtgcaaagataagcgattctttgtcgattgtcttgctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtgtatggatctggaagcacaggaatcccgttttaccaacagttaccagtttggtgtgttcgataaa                                                                                       <  2:575237‑M1/251‑156 (MQ=255)                                                                                                                                                     acaGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagc                                                                                     >  1:279564‑M1/1‑92 (MQ=255)                                                                                                                                                     acaGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagc                                                                                     >  1:188419‑M1/1‑92 (MQ=255)                                                                                                                                                      caGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcagg                                                                                  >  1:407128‑M1/1‑91 (MQ=255)                                                                                                                                                        ggAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcg                                                                                <  1:164372‑M1/251‑163 (MQ=255)                                                                                                                                                                 aaTGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaaga                                                                       <  2:108875‑M1/251‑172 (MQ=255)                                                                                                                                                                 aaTGACGTTACCGAGTGGACAACCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtaaaaagataatcgattctttttcgattgtctggctgtatgcgtcaaagagaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactgga                                                                          >  2:562459‑M1/1‑80 (MQ=255)                                                                                                                                                                    gACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtg                                                                    <  2:188419‑M1/251‑175 (MQ=255)                                                                                                                                                                     aCGTTACCGAGTGGACACCAGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGAAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgt                                                                   >  1:513640‑M1/1‑76 (MQ=255)                                                                                                                                                                      cGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaaga                                                                       <  2:465875‑M1/246‑172 (MQ=255)                                                                                                                                                                      cGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaaga                                                                       >  1:465875‑M1/1‑75 (MQ=255)                                                                                                                                                                          aCCGAGTGGACACCCGAAGAGCAGAATGATCAAAAAAGTTACCGCACTAGGCCCGTATATTCGTGAAGGGAagtgcaaagataatcgattcatttgcgattgtctggctgtatgcgtcaacgtgaaaccggcaacggaagtgcgtgaattctggggcgggtggatggagctggaagagcaggaagcccgtggtacctacagttaccagagtggtctgtgcgataaaggaggggagtggaagaatgttccg                                                               <  2:92731‑M1/250‑180 (MQ=255)                                                                                                                                                                            cGAGTGGACACCCGAAGAGAAGCTTGATCAAAAAATTTACCGCCCTAGGCCCGTATATTCGTGACGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccagcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgttaaagcaggcgactggaagagtgttccggt                                                             >  2:372827‑M1/1‑69 (MQ=255)                                                                                                                                                                                              ggagaTTGATCAAAAAAGTTACCGGACTAGACCCGTATATTCGTGAAGGGAagtgcaaagagaagcgattctttttcgattgtgtggctgtatgcgtcaacgtgaaaccggcaccggaagtacgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttgatctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaaggggtt                                          <  2:339927‑M1/249‑201 (MQ=255)                                                                                                                                                                                                    tGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacg                                     <  2:185006‑M1/250‑206 (MQ=255)                                                                                                                                                                                                         aaaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactgga                               <  2:9876‑M1/251‑212 (MQ=255)                                                                                                                                                                                                                             cccGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagcgggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgctgagaga           <  1:245739‑M1/251‑232 (MQ=255)                                                                                                                                                                                                                              ccGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgacttctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggacgagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgctgag              >  2:538921‑M1/1‑19 (MQ=255)                                                                                                                                                                                                                               cGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgctgagagatc         <  1:36350‑M1/251‑234 (MQ=255)                                                                                                                                                                                                                               agtATAGTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgttaaaccggcaccggaagtgcgtgaatgctggggctggtggatggagctggaatcgcaggaatcccgtgttacctacagttaccagtttggactgttcgataaagcaggcgactggaagagtgttcctgtaaaagacgctgaagtgggtgaacgactggatcacaccctgcgctgagagatc         <  2:273576‑M1/250‑234 (MQ=255)                                                                                                                                                                                                                                       cGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgctgagagatcccctcat  >  2:65985‑M1/1‑10 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                 CTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCACCACGTC  >  NC_000913/257668‑258147

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