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–258675	258683	9–776	35 [1]	[0] 34	insB9–[crl]	insB9,insA9,[crl]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 257907	1 (0.010)	34 (0.380)	21/572	2.9	98.6%	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

ACTGGAAGAACGATCCGTTCAGCCCGGAAGAGGACTCACGCTTAATCACCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCACCACGTCAAAGGGTGACAGCAGGCTCATAAGACGCCCCAGCGTCGCCA  >  NC_000913/257611‑258188                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                          aCTGGAAGAACGATCCGTTCAGCCCGGAAGAGGACTCACGCTTAATCACCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCACGTATATTCGTGAAGGTAagt                                                                                                                                                                                                                                                                                        >  1:70829‑M1/1‑297 (MQ=255) aCTGGAAGAACGATCCGTTCAGCCCGGAAGAGGACTCACGCTTAATCACCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAAAAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagt                                                                                                                                                                                                                                                                                        >  1:476646‑M1/1‑297 (MQ=255) aCTGGAAGAACGATCCGTTCAGCCCGGAAGAGGACTCACGCTTAATCACCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCAAAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagt                                                                                                                                                                                                                                                                                        >  1:476487‑M1/1‑297 (MQ=255) aCTGGAAGAACAATCCGTTCAGCCCGGAAGAGGACTCACGCTTAATCACCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCAAGTATATTCGCGAAGGTAagt                                                                                                                                                                                                                                                                                        >  1:475838‑M1/1‑297 (MQ=255)        gaaCGATCCGTTCAGCCCGGAAGAGGACTCACGCTTAATCACCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaaga                                                                                                                                                                                                                                                                                 <  1:577228‑M1/300‑11 (MQ=255)                 gTTCAGCCCGGAAGAGGACTCACGCTTAATCACCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACNCCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgatt                                                                                                                                                                                                                                                                        >  1:154479‑M1/1‑281 (MQ=255)                      gCCCGGAAGAGGACTCNCGCTTAATCACCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGACTAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATNTTCGTGAAGGTAagtgcanagataatcgattctttt                                                                                                                                                                                                                                                                   >  1:621555‑M1/1‑276 (MQ=255)                           gAAGAGGACTCACGCTTAATCACCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgat                                                                                                                                                                                                                                                              >  1:569818‑M1/1‑271 (MQ=255)                                      aCGCTTAATCACCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctg                                                                                                                                                                                                                                                    <  1:279664‑M1/299‑40 (MQ=255)                                                                                      aTTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACGGTTGCATCAAAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggct                                                                                                                                                                                                    >  1:7596‑M1/1‑212 (MQ=255)                                                                                      aTTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggct                                                                                                                                                                                                    >  1:7610‑M1/1‑212 (MQ=255)                                                                                                 tgtATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagc                                                                                                                                                                                        <  1:78759‑M1/300‑100 (MQ=255)                                                                                                        gAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCTAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagc                                                                                                                                                                                 <  2:446984‑M1/300‑107 (MQ=255)                                                                                                                                        aCTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGTAGATAGCAATGACGTTACCGAGTGGACACCCGAAGACTAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccgtaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagt                                                                                                                                                 <  2:621555‑M1/300‑139 (MQ=255)                                                                                                                                              tCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccgtaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtc                                                                                                                                           <  2:258825‑M1/300‑145 (MQ=255)                                                                                                                                                gAAAAACGCTTGTTTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtc                                                                                                                                           <  2:264252‑M1/298‑145 (MQ=255)                                                                                                                                                gAAAAACCCTTGTTTTAAAAATCTGCTAAATTTTGCCAATTTGGTAAACCAGTTGCATCACAACAGGAGATAGCAAGGACGTTACCGAGTGTACACCATAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCTTATGTTCTTGAAGGTAaatgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtc                                                                                                                                           <  2:261740‑M1/298‑145 (MQ=255)                                                                                                                                                         ttgttgtaaacatTTGCTAAATTTTGCCAATTTGGTAAAACATTTGCATCACAACAGTAGATAGCAAGGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTGACCGCACTAGTCCCGTATATTCGTTAAGGTAaatgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagttcgtgaattcttgggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataa                                                                                                                                 <  2:103269‑M1/293‑155 (MQ=255)                                                                                                                                                         ttGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGGTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTATGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcttcaacgtgaaaccggcaccggaagttcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataa                                                                                                                                 <  2:102248‑M1/299‑155 (MQ=255)                                                                                                                                                                aaaaaTTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcnattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcg                                                                                                                         <  2:556349‑M1/300‑163 (MQ=255)                                                                                                                                                                                                              acaGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctaaagttaacagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagaaactgaagtggttgaacgact                                                                           >  2:295681‑M1/1‑92 (MQ=255)                                                                                                                                                                                                                          aaTGACGTTACCGAGTGTACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactgtaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccc                                                                <  2:317126‑M1/299‑220 (MQ=255)                                                                                                                                                                                                                          aaTGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccc                                                                <  1:278296‑M1/299‑220 (MQ=255)                                                                                                                                                                                                                          aaTGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccc                                                                <  2:284911‑M1/299‑220 (MQ=255)                                                                                                                                                                                                                          aaTGACGTGACCGAGTGGACACCCTAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggngctggtggatggagctggaagcgcaggaatccngttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccc                                                                <  2:569818‑M1/299‑220 (MQ=255)                                                                                                                                                                                                                                gTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcac                                                                    <  1:132420‑M1/289‑216 (MQ=255)                                                                                                                                                                                                                                gTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcac                                                                    >  2:132420‑M1/1‑74 (MQ=255)                                                                                                                                                                                                                                   aCCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGACGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacatcctgcgtgagtt                                                      >  2:310381‑M1/1‑71 (MQ=255)                                                                                                                                                                                                                                               cccGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgtgagtttcacgagaagct                                          >  1:128047‑M1/1‑59 (MQ=255)                                                                                                                                                                                                                                                               aTCAAAAAATTTACCGTACTAGGCTCGTATATTTGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgttaacttgaaaccggcaccgtaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccaggttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgtgagtttcacgagaagctgtgtgaactgttgacg                          <  2:475838‑M1/300‑258 (MQ=255)                                                                                                                                                                                                                                                               aTCAAAAAATTTACCGCACTATGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacttgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgtgagtttcacgagaagctgcgtgaactgctgacg                          <  2:70829‑M1/300‑258 (MQ=255)                                                                                                                                                                                                                                                               aTCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgtgagtttcacgagaagctgcgtgaactgctgacg                          <  2:476487‑M1/300‑258 (MQ=255)                                                                                                                                                                                                                                                               aTCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgtgaatttcacgagaagctgcgtgaactgctgacg                          <  2:476646‑M1/300‑258 (MQ=255)                                                                                                                                                                                                                                                                               cACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgtgagtttcacgagaagctgcgtgaactgctgacgacgctgaatctgaagc          >  2:415908‑M1/1‑27 (MQ=255)                                                                                                                                                                                                                                                                                       gtattaagtaacgttgtaaggtaATGACTCAAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCACCAAGTCAAAGGGTGAAAGCAGGCTCATACGACGCCCCAGCGTCGCCa  >  2:350845/19‑300 (MQ=11)                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                          ACTGGAAGAACGATCCGTTCAGCCCGGAAGAGGACTCACGCTTAATCACCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCACCACGTCAAAGGGTGACAGCAGGCTCATAAGACGCCCCAGCGTCGCCA  >  NC_000913/257611‑258188

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑