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	41 [0]	[1] 43	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)	41 (0.740)	37/474	0.2	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

TACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTG  >  NC_000913/258467‑258934                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                           atcccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatatccgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCt                                                                                                                                                                                                                           <  2:360362‑M2/34‑1 (MQ=255) atcccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagagtgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCt                                                                                                                                                                                                                           <  2:406814‑M2/34‑1 (MQ=255)   cccatttaacccggagtatcggcattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagctagcaatgacgttaccgagtggacaaccgaatagcagattgatcaaaaaatttaccgcactaggcccgtctattcgtgaaggaaAGTGCAAAGATAATCGTATCTTTTTCGATTGTCTgg                                                                                                                                                                                                                         >  2:249846‑M2/216‑251 (MQ=255)    ccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGCTTCTTTTTCGCTTGTCTGGc                                                                                                                                                                                                                        >  2:73486‑M2/215‑251 (MQ=255)         taacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTa                                                                                                                                                                                                                    <  1:284382‑M2/41‑1 (MQ=255)          aaccaggttcatcggaattttgacaaaggtcttcaggaaatcaccgaatggatcgaacaacgcttgtgttaaaaatttgctaaattttggcaatttggaaaaacagtttcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgccctcagcccgtatattcgtgaaggtcAGAGCAAAGATACTCGCTTCTTTTTCGATTGTCTGGctgtctg                                                                                                                                                                                                                  >  2:68364‑M2/209‑248 (MQ=255)               ggtgtatcggaattttgacaacggtcttcaggacatcaccgcctggatcgaaaaacgctcgtgttaaaaatttgctaaattttgctaatttggaaaaaaagttgcatcacaaaaggagatagcaatgacgttaccgagtggacacccgacgatcagattgatcaaaaaatttaccgcactaggcccgtatattagtgacggtaAGTGCAAGGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGGGTc                                                                                                                                                                                                              >  2:313951‑M2/204‑250 (MQ=255)               ggtgtatcggaattttgacaaaggtattcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCa                                                                                                                                                                                                             <  2:4763‑M2/48‑1 (MQ=255)                    atcggaattttgacaaaggtcttcaggacatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagaagattgatcaaaaaatttaccgcactaggcccgtatattcgtgatggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTg                                                                                                                                                                                                        >  2:114376‑M2/199‑251 (MQ=255)                      cggaatttagacaaaggtcttcaggaaatcaccgactagatcgaaaaacgcttgtgttaaaaatgtgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGAAAAGATCATCGATTCTTTTTCGCTTTTCTGGCTGTATGCGACaacgaga                                                                                                                                                                                                       >  2:315983‑M2/197‑247 (MQ=255)                       ggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGaaa                                                                                                                                                                                                     <  1:381822‑M2/56‑1 (MQ=255)                         aattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAAcc                                                                                                                                                                                                   >  1:392368‑M2/194‑251 (MQ=255)                               gacaaaggtcgtcaggatatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcagcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCAcc                                                                                                                                                                                             <  2:240316‑M2/64‑1 (MQ=255)                                     ggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgcgaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGt                                                                                                                                                                                       <  2:421290‑M2/70‑1 (MQ=255)                                            aggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGCGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGa                                                                                                                                                                                 >  2:35200‑M2/175‑250 (MQ=255)                                               aaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgcccatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTc                                                                                                                                                                             >  1:219263‑M2/172‑251 (MQ=255)                                                                      cgcttgtgttaaaaatttgctaaattttgccaatttggtcaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaatgtaAGCGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGAGTCAACGTGAAACCGACACCTGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTgg                                                                                                                                                       >  2:185566‑M2/149‑250 (MQ=255)                                                                       gcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGaa                                                                                                                                                     <  1:209410‑M2/104‑1 (MQ=255)                                                                         ttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcgcaaccggagatcgcaatgacgttaccgagtggacacccgacgagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGCGCGTGAATTCAGGTGCTGGATGATGGAGCtggacgc                                                                                                                                                   >  1:37544‑M2/146‑248 (MQ=255)                                                                                        gctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGttt                                                                                                                                    >  2:163976‑M2/131‑251 (MQ=255)                                                                                              ttttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCt                                                                                                                               >  1:188084‑M2/125‑250 (MQ=255)                                                                                               tttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAACGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTAc                                                                                                                             >  2:374067‑M2/124‑251 (MQ=255)                                                                                               tttgccaatttggtaaaacagttgaatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATCGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTAc                                                                                                                             >  1:150329‑M2/124‑251 (MQ=255)                                                                                                ttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTAc                                                                                                                             <  1:208040‑M2/128‑1 (MQ=255)                                                                                                         tggtaaaaaagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGt                                                                                                                   >  2:346478‑M2/114‑251 (MQ=255)                                                                                                          ggtaaaacagttgcatcacaacaggagatagcactgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGtt                                                                                                                  >  1:127492‑M2/113‑251 (MQ=255)                                                                                                           gtaaaacagttgcattacaacaggcgatagcattgacgttaccgagtggacacctgaagaccagattgatcaaaaaatttaccgcactaggcccgcatattcgggaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGAGTCAACGTGAAACAGGCACTGGAAGTGCGAGAATTCTGGGGCTGGTGGATGGTGCTTGATGCGCAGAAATCCGGTTTTACCTATAGTTACCCGttt                                                                                                                 >  2:178874‑M2/112‑251 (MQ=255)                                                                                                                accgttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGGGCTGGAAGCGCAGGACTCCCGTTTTACCTACAGTTACCAGTTTggtgt                                                                                                            >  2:111131‑M2/107‑249 (MQ=255)                                                                                                                    ttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGCAGCGCAGGAATCCCNTTTTACCTACAGTTACCAGTTTGGTCTGTTc                                                                                                        >  2:225731‑M2/103‑251 (MQ=255)                                                                                                                     tgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCACAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCg                                                                                                       >  1:292206‑M2/102‑251 (MQ=255)                                                                                                                       catcacaacaggagatagcaatgacgttaccnagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTGTGGTCTGTTCGAt                                                                                                     <  2:213278‑M2/152‑1 (MQ=255)                                                                                                                                aggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgccctaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGTTCTGTTCGATAAAGCAGGc                                                                                            <  1:159163‑M2/161‑1 (MQ=255)                                                                                                                                    gatagctatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgaactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCGGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCTACt                                                                                        <  1:220880‑M2/165‑1 (MQ=255)                                                                                                                                                   accgagtggacacccgaagagcagattgatcaaaaaatttaccgcactatgcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCgg                                                                         >  1:259590‑M2/72‑251 (MQ=255)                                                                                                                                                                         agattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTg                                                    >  1:108994‑M2/50‑250 (MQ=255)                                                                                                                                                                            ttgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACg                                                >  1:199859‑M2/47‑251 (MQ=255)                                                                                                                                                                                   aaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTgg                                           >  1:447221‑M2/40‑249 (MQ=255)                                                                                                                                                                                    aaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGc                                        <  2:108994‑M2/213‑1 (MQ=255)                                                                                                                                                                                       atttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACAGGCACCGGAAGGGCGTGAATTCTGGGGCTGGAGGGTGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGAGGTTGAACGACTGGAgcacc                                     >  1:289305‑M2/36‑250 (MQ=255)                                                                                                                                                                                                  taggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTg                            >  2:97628‑M2/25‑249 (MQ=255)                                                                                                                                                                                                                cgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGCATTCTGGGGCTGGTGGATGGAGCTGGAAACGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTACGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCt            >  2:442040‑M2/11‑251 (MQ=255)                                                                                                                                                                                                                         aaGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAAc    <  1:163976/250‑1 (MQ=255)                                                                                                                                                                                                                          aGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTCTGAGTCAACGTGAAACCGGCACCGGCAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATAACGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGCGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActg  >  1:168289/1‑251 (MQ=255)                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                           TACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTG  >  NC_000913/258467‑258934

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