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	59 [0]	[3] 62	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)	59 (0.820)	56/472	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

GGCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACT  >  NC_000913/258441‑258933                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                          catctgctgacggtaaattattagagatcccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaccagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCaa                                                                                                                                                                                                                                                     >  1:26305‑M2/244‑250 (MQ=255)         gacggtaaattattagagatcccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAAt                                                                                                                                                                                                                                              <  2:178407‑M2/14‑1 (MQ=255)                             cccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTgg                                                                                                                                                                                                                        >  1:113003‑M2/216‑251 (MQ=255)                             cacattgaacccggtgtatcggaattttgacaaagggcttcaagaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgaggggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTgg                                                                                                                                                                                                                        <  2:760‑M2/36‑1 (MQ=255)                                  gtaacccggtgtatcggaattttgacaaaggtcgtcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaagggaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGt                                                                                                                                                                                                                    <  2:168261‑M2/40‑1 (MQ=255)                                    aacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATg                                                                                                                                                                                                                 <  2:59204‑M2/43‑1 (MQ=255)                                                   aattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAAcc                                                                                                                                                                                                  >  1:420043‑M2/194‑251 (MQ=255)                                                       ttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacccttgtattacaaatttgctaaattttgcccatttggtaaaaaagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGgcc                                                                                                                                                                                              >  1:319301‑M2/190‑250 (MQ=255)                                                         gacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCAcc                                                                                                                                                                                            >  2:453293‑M2/188‑251 (MQ=255)                                                                gtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTg                                                                                                                                                                                     <  2:263380‑M2/71‑1 (MQ=255)                                                                    tcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaattttgtaaaacagttgcatcacaacaggagatagcaatggcgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccggatattcgtgaaggtaAGGGCAAAGATAATCGATTCTTGGTCGAGTGTCTGGCTGTATGCGTCAACGTGAGTCCGGCACTGGAAGTGCGt                                                                                                                                                                                  <  2:598971‑M2/74‑1 (MQ=255)                                                                    tcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTg                                                                                                                                                                                 >  1:379552‑M2/177‑251 (MQ=255)                                                                       ggaaatcaccgactggatcgacaaacgcttgtgttaaaaatttgctaaattttgcaaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttactgagtggaaccccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAGCGTGAAAACGGCACCGGAAGTGCGTGAAt                                                                                                                                                                              >  2:343254‑M2/174‑251 (MQ=255)                                                                       ggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAAt                                                                                                                                                                              >  1:571937‑M2/174‑251 (MQ=255)                                                                                       atcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGtggat                                                                                                                                                              <  2:652441‑M2/94‑1 (MQ=255)                                                                                          gaaaaacgcttatgttaaaaatttgctaaattttccccatttggtaaaacagtcgcatcacaacaggagatagcaatgacgataccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatcttcgagaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAATGGCGTGAATTCTGGGGCTGGtggatgg                                                                                                                                                            >  2:291683‑M2/155‑250 (MQ=255)                                                                                             aaacgcttgtgttaaaaatttgctaaattttgccaatttggtaacacagttgcataaaaaaaggagataacaatgacgttaccgaggggacaccagaagagcagattgatcaaaaaatttaccgcactaggcccgtatatccgtgagggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGACGTGCGTGAATTCTGGGGCTGGCGGATGGAgcg                                                                                                                                                        >  2:324466‑M2/152‑250 (MQ=255)                                                                                                  cttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGa                                                                                                                                                     >  2:525217‑M2/147‑249 (MQ=255)                                                                                                           aaaatttgctaaattttgccaatttggtaaaacagttgcataacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggaaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGaa                                                                                                                                           >  2:446599‑M2/138‑250 (MQ=255)                                                                                                               tttgctaaattttgccaatttggtaaaacagttgcatcacaaacggagatagcacaaacgttaccgagtggacaccagaagagcagattgatcaacaaatttaccgcactaggcccgtatattcgtgaagataAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGAGTCAACGGGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGTATGGAGCTGGCCGCgcaggaaacacg                                                                                                                                      >  2:586916‑M2/134‑246 (MQ=255)                                                                                                                     aaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGtttt                                                                                                                                  <  2:376778‑M2/122‑1 (MQ=255)                                                                                                                      aattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGACTTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTAcc                                                                                                                               >  2:135921‑M2/127‑251 (MQ=255)                                                                                                                            gccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAg                                                                                                                          <  1:453293‑M2/130‑1 (MQ=255)                                                                                                                             ccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCTTTATACCTACAGtt                                                                                                                        <  2:66362‑M2/132‑1 (MQ=255)                                                                                                                                atttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTAcc                                                                                                                     <  2:113003‑M2/135‑1 (MQ=255)                                                                                                                                   tggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGt                                                                                                                  >  1:454378‑M2/114‑251 (MQ=255)                                                                                                                                     gtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGttt                                                                                                                >  1:226474‑M2/112‑251 (MQ=255)                                                                                                                                       aaaacagttgcatcacaccaggagatagcaatgacgttaccgagtggacacccgaagcgcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgacggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGACTTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTgg                                                                                                              >  2:673728‑M2/110‑251 (MQ=255)                                                                                                                                           cagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTg                                                                                                          >  2:177947‑M2/106‑251 (MQ=255)                                                                                                                                            agttgcatcacaacaggagctagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGt                                                                                                         >  2:75954‑M2/105‑251 (MQ=255)                                                                                                                                            agttgcatcacaacaggagatagcaatgacgttaccgcgtggacacccgaagcgcagattgctcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGt                                                                                                         >  1:638440‑M2/105‑251 (MQ=255)                                                                                                                                                 catcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGAt                                                                                                    >  2:230650‑M2/100‑251 (MQ=255)                                                                                                                                                     acaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCAGTTTTACCTACAGTTACCAGTTTGGTCTGTTAGATaa                                                                                                  >  1:108225‑M2/96‑249 (MQ=255)                                                                                                                                                        acaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCATGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAg                                                                                             <  1:673728‑M2/159‑1 (MQ=255)                                                                                                                                                        acaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGc                                                                                               >  2:531383‑M2/93‑249 (MQ=255)                                                                                                                                                             agatagcaatgacgataccgggtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaagggaAGTGCAAAGATACTCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGCGGCTGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTGCAGTTACCAGCTTGGTCTGTTCGATAAAGCAGTCGAc                                                                                        >  2:62963‑M2/88‑251 (MQ=255)                                                                                                                                                              gatagcaatgacggcaccgagtggacacccgaagcgcagattgatcaaaaaatttaccgcactaggccagcatattcgtgaaggtaAGTGCAACGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGGAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTCCGATAAAGCAGGcgatt                                                                                       >  2:642475‑M2/87‑249 (MQ=255)                                                                                                                                                                         cgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGtt                                                                            <  1:177947‑M2/176‑1 (MQ=255)                                                                                                                                                                          gttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTc                                                                           >  1:177783‑M2/75‑251 (MQ=255)                                                                                                                                                                            taccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCTGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTcc                                                                          <  1:291683‑M2/178‑1 (MQ=255)                                                                                                                                                                             accgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcatgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCgg                                                                        <  1:518687‑M2/180‑1 (MQ=255)                                                                                                                                                                                       cacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCGGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGacac                                                              <  2:453123‑M2/190‑1 (MQ=255)                                                                                                                                                                                        acccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACt                                                             >  1:593597‑M2/61‑251 (MQ=255)                                                                                                                                                                                         cccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTg                                                            >  2:564560‑M2/60‑251 (MQ=255)                                                                                                                                                                                          ccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTg                                                            >  1:245182‑M2/59‑250 (MQ=255)                                                                                                                                                                                          ccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggaaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGGGGATGGAGCTGGAAGCGCAGGAATCCCGGTTTACCTACAGTTACCAGTTTGGTCTGTTCGAtaaagcaggcgactggaagagtaatcgcggtaagctgttcct                                                           >  1:248374‑M2/59‑231 (MQ=255)                                                                                                                                                                                            gaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGa                                                           >  1:652572‑M2/57‑249 (MQ=255)                                                                                                                                                                                             aagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTGCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGt                                                        <  1:164091‑M2/196‑1 (MQ=255)                                                                                                                                                                                              agagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGCTTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCACGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTg                                                       >  1:657893‑M2/55‑251 (MQ=255)                                                                                                                                                                                                   agattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTg                                                   >  2:513512‑M2/50‑250 (MQ=255)                                                                                                                                                                                                    gattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAg                                                                  >  1:434558‑M2/49‑234 (MQ=255)                                                                                                                                                                                                    gattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAAGCCCGTTTTACCGACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCTACTGGAAGAGTGTTCCGGTAAAAg                                                                  <  2:434558‑M2/186‑1 (MQ=255)                                                                                                                                                                                                     attgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCGGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAAc                                                >  2:210878‑M2/48‑251 (MQ=255)                                                                                                                                                                                                         atcaaaaaacttaccgcactaggcccgtatagtcgtgaatgtaAGTGCAAAGATAATCGATTCTTATTCGATTGTCTGGATGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAAGTCGGGGGCTGGTGGATGGAGCTGGAAGCGCAGCACTCCCGTTGTACCTACAGTTGCCAGTTTGGTCTGTGCGATAAAGCAGGCGACTGGAAGAGGGATCCGGTAAAAGACACTGAAGTTGTTTAACGACt                                            <  2:56897‑M2/208‑1 (MQ=255)                                                                                                                                                                                                             aaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTg                                                                                                                                                       <  2:97380‑M2/101‑1 (MQ=255)                                                                                                                                                                                                             aaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTg                                                                                                                                                       >  1:97380‑M2/40‑140 (MQ=255)                                                                                                                                                                                                              aaaatttaccgcactaggcccgtatattcgtgtaggtaAGTGCAATGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAGGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGCAACGCAGGCATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGTCGACTGGAAGAGTGTTCCGTTAAAAGACACTGAAGAGGTTGAACGACTGGAGc                                       >  2:29599‑M2/39‑251 (MQ=255)                                                                                                                                                                                                                         cactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGCACAAGACACTGAAGTGGTTGAACGACTGGAGCACATCCTGCGt                            >  2:566592‑M2/28‑251 (MQ=255)                                                                                                                                                                                                                                    tatattcgtgaaggtaCGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACg                  >  2:573103‑M2/17‑250 (MQ=255)                                                                                                                                                                                                                                             gAAGGCAAGTGCAAAGATAATCGCTTCTTTTCCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACCGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAGCGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCg        >  2:222541/1‑251 (MQ=255)                                                                                                                                                                                                                                                ggTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTATACCTGCATGCGTCAACGTGAAAACGGCACCGGAAGTGCGTGAATTCAGGGGCTGGTGGATGGAGATGGAAGCGAAGGCATCCCGTTTTACCTACAGTTACCAGTTTGGCCTGTTCGATAAAGCAGGCGACTGGAAGCGTGTTCCGGTAAAAGACACTGAAGTGGTAGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTg      >  1:74285/1‑250 (MQ=255)                                                                                                                                                                                                                                                   aaGTGCAACGATCATCGATTCTTTTTCGATTCTCTGGCTGTATGCGTCCACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAAct  >  2:675683/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                          GGCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACT  >  NC_000913/258441‑258933

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