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

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	257908–258675	258683	9–776	64 [0]	[3] 67	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)	64 (0.610)	57/480	0.5	100%	intergenic (+8/+16)	crl/insB9	RNA polymerase holoenzyme assembly factor Crl/IS1 family transposase B
	?	NC_000913	258684 =	0 (0.000)					pseudogene (9/331 nt)	crl	RNA polymerase holoenzyme assembly factor Crl

GCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACT  >  NC_000913/258442‑258933                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                          atctgctgacggtaaattattagagatcccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCaa                                                                                                                                                                                                                                                     >  1:566538‑M2/243‑249 (MQ=255)  tctgctgacggtaaattattagagatcccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCaa                                                                                                                                                                                                                                                     >  1:547682‑M2/242‑248 (MQ=255)     gcggacggtaaattattagagatcccaattaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcgtgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGAt                                                                                                                                                                                                                                                 <  2:870398‑M2/11‑1 (MQ=255)        gacggtaaattattagagatcccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATc                                                                                                                                                                                                                                             <  1:746995‑M2/15‑1 (MQ=255)                       gagatcccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaagttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtgAGTGCAAAGATAATCGATTCTTTTTCGATTg                                                                                                                                                                                                                             <  2:9581‑M2/31‑1 (MQ=255)                             ccatttaacccggtgtatcggaattttgacaaaggtcttcagtaaatcaccgactggatcgaaaaactcttgtgttaaaaatttgctaaattttgccaattgtgtaaaacagttgcatcacaacaggagatagcaatgtctttaccgagtagacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatgttcgtgaaggtaAGTACAAAGATAATCGATTCTTTTTCGATTGTCtgga                                                                                                                                                                                                                       <  2:272866‑M2/37‑2 (MQ=255)                               atttaacacggtgtatcggaattttgacaacggtcttcaggaaatcacagactggagcgaaaaacgcttgtgttaaaaatttgctaacttttgccaatttggaaaaacagttgcatcaaaacaggagatagcaatgacgttaccgcgtggacacccgaagagaagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATGGTCTGGCTg                                                                                                                                                                                                                     >  1:704103‑M2/213‑251 (MQ=255)                                tttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGt                                                                                                                                                                                                                    >  1:34043‑M2/212‑251 (MQ=255)                                 ttagcccggtgtatcggaattttgacaaaggtctgcatgaaatcaccgactgtatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatgtggtaaaacagttgcatcacaacaggagctagcaatgtcgttaccgagtggacacccgaagagcagattgatcaaaaaatataccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAAACGATTCATTTTCGATTGTCTGGCTGt                                                                                                                                                                                                                    <  2:422381‑M2/40‑1 (MQ=255)                                    acccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGc                                                                                                                                                                                                                >  1:68928‑M2/208‑251 (MQ=255)                                     cccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCg                                                                                                                                                                                                               >  2:158151‑M2/207‑251 (MQ=255)                                      ccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGt                                                                                                                                                                                                              >  1:78175‑M2/206‑251 (MQ=255)                                          tgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAAc                                                                                                                                                                                                          >  1:265866‑M2/202‑251 (MQ=255)                                           gtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACg                                                                                                                                                                                                         <  1:450965‑M2/51‑1 (MQ=255)                                                ggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAGAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCTTCAACGTGaa                                                                                                                                                                                                     >  1:33370‑M2/196‑250 (MQ=255)                                                ggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGaa                                                                                                                                                                                                     >  1:14759‑M2/196‑250 (MQ=255)                                                        gacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCAcc                                                                                                                                                                                            <  2:699775‑M2/64‑1 (MQ=255)                                                                   tcaggaaatccccgactggatcgaaaaacgcttgtgttaaaaatttgctaaatgttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagcggactcccgaagagcagattgagcaaaaaatttaccgcactaggcccgtatagtcgtgaaggtaAGTGCAAAGATAATCGAATCTTGGGCGATTGTCTGGCTGTATTTGTCAACGTTAAACCGGCACGGGAAGTGCGt                                                                                                                                                                                  <  2:219736‑M2/74‑1 (MQ=255)                                                                   tcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCg                                                                                                                                                                                   <  2:34043‑M2/73‑1 (MQ=255)                                                                         aatcaccgacttgatcgaaaaacgcttgtgttaaaaatttgctaaattttcccaatttggtaaaacagttgcatcaaaacaggagatagcaatgacgttaccgagtttacacccgaagagcagattggtcaaaaaattttccgcacttggcccgtatattcgtgaaagtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTTGCTGTATGCGTGAACGTGAAACCTGCACCTGAAGTGCGTGAATtcg                                                                                                                                                                           <  2:709935‑M2/81‑2 (MQ=255)                                                                         aatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCt                                                                                                                                                                           <  1:716005‑M2/81‑1 (MQ=255)                                                                                gactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttacagagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCtg                                                                                                                                                                    >  1:472114‑M2/164‑251 (MQ=255)                                                                                gactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcaaaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCtg                                                                                                                                                                    >  2:841792‑M2/164‑251 (MQ=255)                                                                                      atcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGtggat                                                                                                                                                              <  2:572675‑M2/94‑1 (MQ=255)                                                                                        cgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGtggatgg                                                                                                                                                            >  2:70252‑M2/156‑251 (MQ=255)                                                                                          aaaaacgcttgtgttaaaaatttgctaaattttaccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcacaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGtggatgg                                                                                                                                                            >  1:682136‑M2/154‑249 (MQ=255)                                                                                            aatcgcttgtgttaaaaattcgctaagttttgccaatttggtaaaacagttgaatcacaacaggagaaagcaatggcgttaccgagtggacacccgaagagcaggttgatcaaaaaatttaccgcacaagacccgtatattcgtgaaggtaAGAGAAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCGACGGGAAACAGGCACCGGAAGTGCGAGAATTGTGGGGATGGTGGTTGGAGCt                                                                                                                                                        >  2:652229‑M2/152‑251 (MQ=255)                                                                                              acgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTgg                                                                                                                                                      >  1:76345‑M2/150‑251 (MQ=255)                                                                                                gcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCt                                                                                                                                                                                                                                       <  2:466334‑M2/21‑1 (MQ=255)                                                                                                gcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCt                                                                                                                                                                                                                                       >  1:466334‑M2/148‑168 (MQ=255)                                                                                                        taaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAgg                                                                                                                                             >  1:646806‑M2/140‑250 (MQ=255)                                                                                                                 gctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGAGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGttt                                                                                                                                   >  1:48189‑M2/131‑251 (MQ=255)                                                                                                                    aaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTAAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTAc                                                                                                                                >  1:93697‑M2/128‑251 (MQ=255)                                                                                                                        tttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTAc                                                                                                                            >  1:766330‑M2/124‑251 (MQ=255)                                                                                                                         ttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTAc                                                                                                                            >  1:764599‑M2/123‑250 (MQ=255)                                                                                                                           gccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGt                                                                                                                         >  2:423882‑M2/121‑251 (MQ=255)                                                                                                                                  tggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGt                                                                                                                  >  2:582135‑M2/114‑251 (MQ=255)                                                                                                                                   ggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGtt                                                                                                                 >  2:143607‑M2/113‑251 (MQ=255)                                                                                                                                      aaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTgg                                                                                                              >  2:424062‑M2/110‑251 (MQ=255)                                                                                                                                      aaaacagttgcagcacaacaggagatagcaatgacggtaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTg                                                                                                               <  2:93697‑M2/141‑1 (MQ=255)                                                                                                                                              tgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGCGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGTTCTGTTCg                                                                                                      >  2:207977‑M2/102‑251 (MQ=255)                                                                                                                                               gcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGa                                                                                                     <  2:171819‑M2/151‑1 (MQ=255)                                                                                                                                                catcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCCGTTTGGTCTGTTCGAt                                                                                                    >  2:818960‑M2/100‑251 (MQ=255)                                                                                                                                                 atcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGAt                                                                                                    >  1:712300‑M2/99‑250 (MQ=255)                                                                                                                                                     caacaggagatagcaatgacgttaccgcgtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGc                                                                                               >  1:200307‑M2/95‑251 (MQ=255)                                                                                                                                                          ggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCg                                                                                          <  1:537289‑M2/162‑1 (MQ=255)                                                                                                                                                           gagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgaatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGGGAATTCTGGGGATGGTGGATGGAGCTGGCCGCGCAGGAATCCCGTTTTACCTACAGTTACCCGTTTGGACTGTTCGTTAAAGCAGGCg                                                                                          >  2:562663‑M2/89‑250 (MQ=255)                                                                                                                                                            agatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGAc                                                                                        <  2:78175‑M2/164‑1 (MQ=255)                                                                                                                                                                   aatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCCCCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGCGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGa                                                                                    >  2:834297‑M2/81‑248 (MQ=255)                                                                                                                                                                           taccgagtggacacccgaagagcagattgatcaaaaaattttccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTAGATAAAGCAGGCGACTGGAAGAGTGTtct                                                                          <  2:33370‑M2/178‑2 (MQ=255)                                                                                                                                                                              cgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTa                                                                      <  1:349153‑M2/182‑1 (MQ=255)                                                                                                                                                                               gagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtataatcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGCGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTGGGTCTGTTCGATAAAGCAGGAGACTGGAAGAGAGTTccggga                                                                      >  2:94731‑M2/69‑248 (MQ=255)                                                                                                                                                                                agtggacacccgaagagcagattgttcaaaaactttaccgcactaggcacgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGTAACCGGCACCGGAACTGCGTGAATTCTGGGGCTGGTGGATGGCGCTGGAAGCGCAGGACTCCCGTTTTACCTACAGTTACCAGTTTGGCCTGTTCGATAATGCAGGCGACGGGAAGAGTGTTCCTGTaa                                                                     >  2:99252‑M2/68‑250 (MQ=255)                                                                                                                                                                                            aagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGt                                                        >  1:602729‑M2/56‑251 (MQ=255)                                                                                                                                                                                              gagcagattgatcaaaacatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTgg                                                      >  1:164011‑M2/54‑251 (MQ=255)                                                                                                                                                                                                   gattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGa                                                  >  2:860956‑M2/49‑250 (MQ=255)                                                                                                                                                                                                       gatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGAc                                             >  1:453094‑M2/45‑251 (MQ=255)                                                                                                                                                                                                       gatcaaaaaaattaccgcacgaggcccgtatattcgagaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCGACAGTTACCAGTTGGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTgaaccac                                             <  2:48189‑M2/207‑4 (MQ=255)                                                                                                                                                                                                          caaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTgg                                          >  1:570012‑M2/42‑251 (MQ=255)                                                                                                                                                                                                                atttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGcac                                     >  2:622478‑M2/36‑250 (MQ=255)                                                                                                                                                                                                                     ccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTg                               <  1:70252‑M2/221‑1 (MQ=255)                                                                                                                                                                                                                       ccactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCg                             <  1:53190‑M2/223‑1 (MQ=255)                                                                                                                                                                                                                                  gtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACg                  <  1:143607‑M2/234‑1 (MQ=255)                                                                                                                                                                                                                                   tatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACg                  >  1:562882‑M2/17‑250 (MQ=255)                                                                                                                                                                                                                                      attggttgtggTAAGTGCAAAGATAATCGATTCTTTTTCGGTTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGGAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTTCGTTAGTTTCACGAGaa              <  2:845813/242‑1 (MQ=255)                                                                                                                                                                                                                                             aaGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCgg       >  1:385403/1‑250 (MQ=255)                                                                                                                                                                                                                                                  aaGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAAct  >  2:204627/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                          GCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACT  >  NC_000913/258442‑258933

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: ‑