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	70 [0]	[0] 73	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)	70 (0.840)	57/472	0.3	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

TTCATGGCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTG  >  NC_000913/258436‑258934                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                           ctcatcatctgctgacggtaaattattagagatcccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGt                                                                                                                                                                                                                                                          >  1:71771‑M2/249‑251 (MQ=255)           gctgacggtaaattattagagatcccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatataccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATaa                                                                                                                                                                                                                                                <  2:517291‑M2/13‑1 (MQ=255)           gctgacggtaaattattagagagcccatttaacccggtgtatcggaattttgacaaaggtcttcaggaattcaccgactggatctaaatacgcttgttttaaaaatttgctaaattttgccaatttggtgaaactgttgcatcacaacaggagttagcaaggacgttaccgagtgggcacacgaagaggagattgatcaaaaaatttacctctataggcacgtatattcgtgaaggtaAGTGCAAAGATaa                                                                                                                                                                                                                                                <  1:701767‑M2/13‑1 (MQ=255)                cggtaaattattagagatcccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGAt                                                                                                                                                                                                                                           <  1:564516‑M2/18‑1 (MQ=255)                acggtaaattattagagatcccatttaacccggtgtatcggattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaagtttgccaattttgtaaaacagttggatcacaacaggagatagcaaggacgttaccgagtggacacgcgaagagtagattgatcaaaaaagttacctcactagggccgtatattcgtgaaggtaAGTGCAAAGATAATCGAt                                                                                                                                                                                                                                           <  1:64920‑M2/18‑1 (MQ=255)                 ggtaaattattagagatcccatttaacccggtgtatcggacttttgacaaaggtcttcaggaactcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagcgcagattgatcaaacaatttacagcactaggcccgtatattcgtgaaggcaAGTGCAAAGATAATCGAtt                                                                                                                                                                                                                                          >  1:721506‑M2/233‑251 (MQ=255)                  gtaaattattagagatcccatgtaacccggtttatctgaattttgactaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatgtggtaaaacagttgcatcacaacaggagaaagcaatgacgtttccgagtggtcgctcgatgagcatattgatcaaaaaatttacctcactaggcccgtatattcgtaaaggtaAGTGCAAAGATAATCGATTc                                                                                                                                                                                                                                         <  1:668276‑M2/20‑1 (MQ=255)                            agagatcccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGAtt                                                                                                                                                                                                                               <  1:634394‑M2/30‑1 (MQ=255)                               gatccaatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcacagactggatcgcacaacgcttgtgttaaaaatttgatcaattttgacaatttggtaaaacagttgcatcacaccaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcgaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTc                                                                                                                                                                                                                            >  2:318788‑M2/219‑251 (MQ=255)                                   ccatttaacccggtgtatcggcattttgacaacggtcttaaggaaatcaccgactggttcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcactgacgttaccgagtggacacccgcagagcagattgatcaaaaaatttaccgccctaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGc                                                                                                                                                                                                                        >  2:707056‑M2/215‑251 (MQ=255)                                        taacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTAt                                                                                                                                                                                                                   >  2:251399‑M2/210‑251 (MQ=255)                                        taacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTAt                                                                                                                                                                                                                   >  1:515072‑M2/210‑251 (MQ=255)                                           cccagtgtagcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatgggctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggaccaccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgcgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGAGTGTCTGGCTGTATGCg                                                                                                                                                                                                                <  2:239750‑M2/45‑1 (MQ=255)                                                           tttgacaaaagtcatcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaaattggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGc                                                                                                                                                                                                >  1:178950‑M2/191‑251 (MQ=255)                                                                         tcaggtcttaaacgactggatagaaaaccgcttgtgttaaaaatttgctaaattttgccaatttggtaaaaaagttgaatcaacacaggagatagcaatgacgttaccgagtggacacccgacgagcagattggtcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCCCCGGAAGAGCGTg                                                                                                                                                                                  >  2:645982‑M2/177‑251 (MQ=255)                                                                                atcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTg                                                                                                                                                                           >  1:333815‑M2/170‑251 (MQ=255)                                                                                       cctggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGGGCGTGACTTCTGGGGCtgg                                                                                                                                                                    >  2:28729‑M2/163‑251 (MQ=255)                                                                                       actggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCtgg                                                                                                                                                                    >  1:225788‑M2/163‑251 (MQ=255)                                                                                        ctggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTc                                                                                                                                                                                                                                         <  2:412115‑M2/20‑1 (MQ=255)                                                                                        ctggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTc                                                                                                                                                                                                                                         >  1:412115‑M2/162‑181 (MQ=255)                                                                                        ctggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGctggg                                                                                                                                                                   >  2:254421‑M2/162‑250 (MQ=255)                                                                                          ggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCtggtgg                                                                                                                                                                 >  2:509138‑M2/160‑251 (MQ=255)                                                                                                  aaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCt                                                                                                                                                         >  1:684550‑M2/152‑251 (MQ=255)                                                                                                          gtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAgcg                                                                                                                                                  <  1:691758‑M2/107‑1 (MQ=255)                                                                                                          gtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAgcg                                                                                                                                                  <  1:251399‑M2/107‑1 (MQ=255)                                                                                                              taaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggcaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGAGGGAGCTGGAAGCGCAgg                                                                                                                                              >  2:591269‑M2/140‑250 (MQ=255)                                                                                                                            attttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCt                                                                                                                               >  2:136559‑M2/126‑251 (MQ=255)                                                                                                                              tttgccaatttggtcaaacagttgcatcacaacaggcgatagcaatgacgttacagagtggaaacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAACGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGCGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTCCCTAc                                                                                                                             >  2:270274‑M2/124‑251 (MQ=255)                                                                                                                                tgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGCATCCCGTTTTACCTAc                                                                                                                             >  1:443377‑M2/122‑249 (MQ=255)                                                                                                                                    aatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTAc                                                                                                                       >  1:518656‑M2/118‑251 (MQ=255)                                                                                                                                       ttggtaaaacagttgcatcacaacaggagctagaaataacgttaccgagaggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatcttcgtgaaggtcAGTGCAACGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTaca                                                                                                                      >  2:497863‑M2/115‑248 (MQ=255)                                                                                                                                          gtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGttt                                                                                                                 >  2:324695‑M2/112‑251 (MQ=255)                                                                                                                                          gtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGttt                                                                                                                 >  2:137195‑M2/112‑251 (MQ=255)                                                                                                                                           taaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCTTGAATTCTgggg                                                                                                                                                                        <  1:254717‑M2/85‑1 (MQ=255)                                                                                                                                           taaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTgggg                                                                                                                                                                        >  2:254717‑M2/111‑195 (MQ=255)                                                                                                                                                  gttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaacatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGCTTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGAGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGtt                                                                                                         >  1:628927‑M2/104‑251 (MQ=255)                                                                                                                                                         cacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATaa                                                                                                   >  2:332577‑M2/97‑250 (MQ=255)                                                                                                                                                           caacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGc                                                                                                >  1:162526‑M2/95‑251 (MQ=255)                                                                                                                                                            aacaggagatagcactgccgttaccgagtggccacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggcaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACCGTTACCAGTTTGGTCTGTTCGATAAAGc                                                                                                >  2:671103‑M2/94‑250 (MQ=255)                                                                                                                                                              caggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGTAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTGCGATAAAGCAgg                                                                                             <  2:683749‑M2/160‑1 (MQ=255)                                                                                                                                                                   gatagcaatgacgttaccgcgtggacaaccgaagcgaagattgatcaaaaaatttaccgcactaggcccgtatattcgagcaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCCGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGTGACt                                                                                        >  2:464401‑M2/87‑251 (MQ=255)                                                                                                                                                                   gatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAg                                                                                                                    <  2:568930‑M2/137‑1 (MQ=255)                                                                                                                                                                   gatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAg                                                                                                                    >  1:568930‑M2/87‑223 (MQ=255)                                                                                                                                                                   gatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTGTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACt                                                                                        <  2:684550‑M2/165‑1 (MQ=255)                                                                                                                                                                      agcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGa                                                                                     <  1:509138‑M2/168‑1 (MQ=255)                                                                                                                                                                      agcaatgacgtgaccgagtggacacccgaagagcagattgatcaaaaaagttacctccctaggcccgtatagtcgtgcaggtaAGTGCAAAGATAATCGAGTCTTTTTCTATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTTCGTGAATTCTGTGGCTGGTGGAAGGAGCTGGAAGCGCATGAATCCCGTTTTACCCACAGTTACCAGTTTGGTCGTTTCGATAAAGCAGGCGACTgg                                                                                      <  2:326454‑M2/167‑1 (MQ=255)                                                                                                                                                                             acgttaccgagtggacgacccaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTTAATTCTGGGTCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTGTACCTACAGTGACCAGTTTGGTCTGTTCGATAAAGCAGGCCACTGGAAGAgtgt                                                                              <  2:128650‑M2/175‑1 (MQ=255)                                                                                                                                                                             acgttaccgagtggacacccgaagagcagattgagcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTGTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGGTTTGCCTACAGTGATCAGTTTGGTCTGTTCGATAAAGCAGTCGACTGGAAGAgtgt                                                                              <  2:65140‑M2/175‑1 (MQ=255)                                                                                                                                                                              cgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAGCCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGtt                                                                             <  2:443377‑M2/176‑1 (MQ=255)                                                                                                                                                                                 taccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTGCGATAAAGCAGGCGACTGGAAGAGTGTTcc                                                                           <  2:582030‑M2/178‑1 (MQ=255)                                                                                                                                                                                     gagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTaa                                                                      <  1:671103‑M2/183‑1 (MQ=255)                                                                                                                                                                                      agtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATaaagcaggcgactggagaagtgtccggttaa                                                                      >  1:716312‑M2/68‑235 (MQ=255)                                                                                                                                                                                          gacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGTAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGac                                                                 <  1:497863‑M2/188‑1 (MQ=255)                                                                                                                                                                                             acccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGTGCTGGTGGATGTAGCTGGAAGCTCAGGAATCCCGTTTTACCTTCAGTTACCAGTTGGGTCGGGTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACt                                                              <  1:464401‑M2/191‑1 (MQ=255)                                                                                                                                                                                             acccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTGCAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACt                                                              >  1:149896‑M2/61‑251 (MQ=255)                                                                                                                                                                                              cccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGCCACTg                                                             >  2:209055‑M2/60‑251 (MQ=255)                                                                                                                                                                                                        agattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTg                                                    >  1:567039‑M2/50‑250 (MQ=255)                                                                                                                                                                                                        agattgatcaaaaaatttaccgaactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTg                                                    >  1:198121‑M2/50‑250 (MQ=255)                                                                                                                                                                                                             gatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACTTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGAGTGAGCTGGAAGCGCAGGAATCCCTTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGAc                                              <  1:707056‑M2/207‑1 (MQ=255)                                                                                                                                                                                                                caaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTg                                            <  2:198121‑M2/209‑1 (MQ=255)                                                                                                                                                                                                                  aaaaatttaccgcactaggcacgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGACGTGGTTGAACGACTgg                                           >  2:722666‑M2/40‑249 (MQ=255)                                                                                                                                                                                                                            cgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGc                               <  1:318788‑M2/222‑1 (MQ=255)                                                                                                                                                                                                                             gcactcggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCCCCGGAAGTGCGTGACTTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATCAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCg                              >  2:45605‑M2/29‑251 (MQ=255)                                                                                                                                                                                                                             gcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTGCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGATTTGGTTGAACGACTGGAGCATACCCTGCg                              <  1:254421‑M2/223‑1 (MQ=255)                                                                                                                                                                                                                              cactaggcccgtatatttgagaaggtaAGTGCAAAGATCCTCGATTCTTTTTCGATTGTCTGGCTGTATGAGTCAAAGTGAAACCGGAACCGGAAGTGCGTGAATTATGGGGCTGGTGGATGGCGCTGGAAGCGCAGGAATCCAGTTTTACATACAGTTACCAGTTTGGACTGTTCGATAAAGCAGGCGACTGGAAGCGTGTTCCGGAACATGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGt                             >  2:650155‑M2/28‑251 (MQ=255)                                                                                                                                                                                                                                  aggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGCAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGt                         >  2:240516‑M2/24‑251 (MQ=255)                                                                                                                                                                                                                                  aggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGt                         >  1:443234‑M2/24‑251 (MQ=255)                                                                                                                                                                                                                                    gcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGttt                       >  1:586638‑M2/22‑251 (MQ=255)                                                                                                                                                                                                                                       cgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCAc                    <  1:28729‑M2/233‑1 (MQ=255)                                                                                                                                                                                                                                          atattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGGTTGACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTGCCGGTAAAAGACACTGAAGTGGGTGAACGACTGGAGCACACCCTGCGTGAGTTTCACgag                 <  1:681114‑M2/236‑1 (MQ=255)                                                                                                                                                                                                                                                         aGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAg                                                          <  1:703465/195‑1 (MQ=255)                                                                                                                                                                                                                                                         aGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAg                                                          >  2:703465/1‑195 (MQ=255)                                                                                                                                                                                                                                                         aGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActg  >  1:203230/1‑251 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                           TTCATGGCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTG  >  NC_000913/258436‑258934

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