BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	2,910,754	Δ300 bp		[mazF]	[mazF]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	2910754	2911053	300	77 [1]	[2] 79	[mazF]	[mazF]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2910753	1 (0.010)	76 (0.860)	70/488	0.1	99.4%	intergenic (‑68/+3)	mazG/mazF	nucleoside triphosphate pyrophosphohydrolase/endoribonuclease toxin MazF
*	?	NC_000913	2911054 =	0 (0.000)	76 (0.860)	70/488	0.1	99.4%	coding (38/336 nt)	mazF	endoribonuclease toxin MazF

GGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTTCATCATCAATATTCAGATTGAGCGCCTGCATTAACGTAGCCGGGATCCGCACCGCCGGT  >  NC_000913/2910806‑2911304

gcccagttcaccgcgaagatcgtcaaaatcttcacgggcgatggcgtccagcacttcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacntgtgaccagaatagaagtgagttagtaaATc                                                                                                                                                                                                                                                          >  2:355934‑M2/249‑251 (MQ=255)

  ccagttcaccgcgaagatcgtcaaaatcttcacgggcgatggcgtccagcacttcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAg                                                                                                                                                                                                                                                        <  2:606926‑M2/5‑1 (MQ=255)

  ccagttcaccgcgaagatcgtcaaaatcttcacgggcgatggcgtccagcacttcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAg                                                                                                                                                                                                                                                        <  2:382719‑M2/5‑1 (MQ=255)

     gttcaccgcgaagatcgtcaaaatcttcacgggcgatggcgtccagcacttcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATc                                                                                                                                                                                                                                                          >  1:637393‑M2/244‑246 (MQ=255)

      ttcaccgcgaagatcgtcaaaatcttcacgggcgatggcgtccagcacttcgtatgtttcttcaagggtgtaaggcgcaattgtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcagtcgctgcataatagtgggcaaacggtcgatttgattcattgaattgtcctgaaaagtgcggttctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATc                                                                                                                                                                                                                                                     <  2:86441‑M2/8‑1 (MQ=255)

       tcaccgcgaagatcgtcaaaatcttcacgggcgatggcgtccagcacttcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCg                                                                                                                                                                                                                                                    <  1:256074‑M2/9‑1 (MQ=255)

         accgcgaagatcgtcaaaatcttcacgggcgatggcgtccagcacttcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGccc                                                                                                                                                                                                                                                 <  1:59239‑M2/12‑1 (MQ=255)

           cgcgaagatcgtcaaaatcttcacgggcgatggcgtccagcacttcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaccctgtgcccagaatagaagtgagttagtaaATCAGATCGCCCat                                                                                                                                                                                                                                               >  2:274107‑M2/238‑251 (MQ=255)

           cgcgaagatcgtcaaaatcttcacgggcgatggcgtccagcacttcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCat                                                                                                                                                                                                                                               <  1:392371‑M2/14‑1 (MQ=255)

            gcgaagatcgtcaaaatcttcacgggcgatggcgtccagcacttcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCat                                                                                                                                                                                                                                               >  1:679223‑M2/237‑250 (MQ=255)

              gaagatcgtcaaaatcttcacgggcgatggcgtccagcacttcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATc                                                                                                                                                                                                                                            <  1:50620‑M2/17‑1 (MQ=255)

               aagatcgtcaaaatctgcacgggcgatggcgtccagcacttcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCg                                                                                                                                                                                                                                           <  2:75744‑M2/18‑1 (MQ=255)

                     gtcaaaatcttcacgggcgatggcgtccagcacttcgtaggtttcttcaagggagtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtactgaaaattgcgggtctgtcaggtggaaacatgtgaccagaatagcagtgagttagtaaATCAGATCGCCCATATCGGGTACg                                                                                                                                                                                                                                     >  2:698303‑M2/228‑251 (MQ=255)

                     gtcaaaatctgcacgggcgatggcgtccagcacttcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgatgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACg                                                                                                                                                                                                                                     <  2:677937‑M2/24‑1 (MQ=255)

                         aaatcttcacgggcgatggcgtccagcacttcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctcttgatcccacgggcagccgttttccggatcgcgcaggcgctgcataagagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATc                                                                                                                                                                                                                                 <  1:82368‑M2/28‑1 (MQ=255)

                               tcacgggcgatggcgtccagcacttcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgctttgattcattgaattgtcctgcaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCtt                                                                                                                                                                                                                            >  2:177681‑M2/218‑250 (MQ=255)

                                     gcgatggcgtccagcacttcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgatttagtaaATCAGATCGCCCATATCTGGTACGTATCGGCttaccatta                                                                                                                                                                                                                     <  1:559396‑M2/40‑1 (MQ=255)

                                       gatggcgtccagcacttcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCttaccattacc                                                                                                                                                                                                                   <  1:344428‑M2/42‑1 (MQ=255)

                                            cgtccagcacttcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACt                                                                                                                                                                                                              <  2:248276‑M2/47‑1 (MQ=255)

                                              tccagcacttcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTc                                                                                                                                                                                                            >  1:29593‑M2/203‑251 (MQ=255)

                                                  gcacttcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccactggcagccgttttccggatcgcgcaggcgctgcataatagttagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTa                                                                                                                                                                                                        <  1:460318‑M2/53‑1 (MQ=255)

                                                       tcgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCtt                                                                                                                                                                                                    <  1:478397‑M2/57‑1 (MQ=255)

                                                        cgtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTc                                                                                                                                                                                                  <  1:605964‑M2/59‑1 (MQ=255)

                                                         gtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCg                                                                                                                                                                                                 <  1:564759‑M2/60‑1 (MQ=255)

                                                         gtaggtttcttcaagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCg                                                                                                                                                                                                 >  1:499236‑M2/192‑251 (MQ=255)

                                                                     aagggtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGt                                                                                                                                                                                     >  1:474200‑M2/180‑251 (MQ=255)

                                                                         gtgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggagctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGat                                                                                                                                                                                 >  1:51797‑M2/176‑251 (MQ=255)

                                                                          tgtaaggcgcaatggtggcaaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGat                                                                                                                                                                                 >  1:587172‑M2/175‑250 (MQ=255)

                                                                                          ggcaaatgtctgctctttatccctcgggcagccgttttccggagcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGtttt                                                                                                                                                                <  2:101059‑M2/93‑1 (MQ=255)

                                                                                          ggcaaatgtctgctctttatcacacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgagaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGCTATTCTCGAGGAGGtttt                                                                                                                                                                >  2:27963‑M2/159‑251 (MQ=255)

                                                                                            caaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTcc                                                                                                                                                              <  1:364377‑M2/95‑1 (MQ=255)

                                                                                            caaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcggatctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACcagaatt                                                                                                                                                                                                             >  2:698790‑M2/157‑201 (MQ=255)

                                                                                            caaatgtctgctctttagcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACtt                                                                                                                                                                                                             <  1:698790‑M2/48‑1 (MQ=255)

                                                                                             aaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTAt                                                                                                                                                                                                       <  2:340489‑M2/54‑1 (MQ=255)

                                                                                             aaatgtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTAt                                                                                                                                                                                                       >  1:340489‑M2/156‑209 (MQ=255)

                                                                                                 gtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTTCCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGt                                                                                                                                                         <  1:635425‑M2/100‑1 (MQ=255)

                                                                                                 gtctgctctttatcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGt                                                                                                                                                         <  2:635206‑M2/100‑1 (MQ=255)

                                                                                                        ctttatcccacgggcagccgttttccggattgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgacccgaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCAGGCGTGATGTCg                                                                                                                                                  >  2:101025‑M2/145‑251 (MQ=255)

                                                                                                            atcccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTg                                                                                                                                               <  1:466985‑M2/110‑1 (MQ=255)

                                                                                                              cccacgggcagccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGAcc                                                                                                                                            <  1:177681‑M2/113‑1 (MQ=255)

                                                                                                                 acgggcagccgttttccgtatcgcgcaggcgctgcataatagtgagcaaccggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGCGGTTTTCCGGCGTGATGTCGTTGACCAGt                                                                                                                                         >  2:611645‑M2/136‑251 (MQ=255)

                                                                                                                 acgggcagccgttttccggatcgcgcaggcgctgcataatagggagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggaggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACAAGAATTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGCGTTTCCGGCGTGATGTCGTTGACCAGt                                                                                                                                         >  2:80528‑M2/136‑251 (MQ=255)

                                                                                                                       agccgttgtccgtatcgcgcaggcgctgcataatagtgagcaaacggtcgatttggttcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGTTGTCGTTGACCAGTTCAGCa                                                                                                                                   <  2:307825‑M2/122‑1 (MQ=255)

                                                                                                                        gccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCaa                                                                                                                                  <  1:81220‑M2/123‑1 (MQ=255)

                                                                                                                         ccgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAg                                                                                                                                 <  1:44612‑M2/124‑1 (MQ=255)

                                                                                                                          cgttttccggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGc                                                                                                                                <  1:602119‑M2/125‑1 (MQ=255)

                                                                                                                                 cggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAAt                                                                                                                          >  1:31613‑M2/120‑250 (MQ=255)

                                                                                                                                 cggatcgcgcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAAt                                                                                                                          >  1:477507‑M2/120‑250 (MQ=255)

                                                                                                                                         gcaggcgctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGctct                                                                                                                 <  1:477628‑M2/140‑1 (MQ=255)

                                                                                                                                          caggcgctgcatcatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgacaagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTACGGCGAGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCtt                                                                                                                >  2:620175‑M2/111‑251 (MQ=255)

                                                                                                                                               gctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGTCTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGc                                                                                                           <  1:27963‑M2/146‑1 (MQ=255)

                                                                                                                                                ctgcataatagtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGc                                                                                                           >  2:211432‑M2/105‑250 (MQ=255)

                                                                                                                                                       atagtgagcaaacggtcgatttgattcattgaattgtcctgaaacttgcggttctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCCTATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCa                                                                                                                                                                                       >  1:252491‑M2/98‑167 (MQ=255)

                                                                                                                                                         agtgagcaaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaaactgtgaccagaatagaagagagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTTCCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCa                                                                                                  >  1:329037‑M2/96‑250 (MQ=255)

                                                                                                                                                                aaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTc                                                                                                                                                   <  2:391065‑M2/106‑1 (MQ=255)

                                                                                                                                                                aaacggtcgatttgattcattgaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTc                                                                                                                                                   >  1:391065‑M2/89‑194 (MQ=255)

                                                                                                                                                                                     gaattgtcctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAAt                                                                     <  1:355934‑M2/184‑1 (MQ=255)

                                                                                                                                                                                            cctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATcttcact                                                              <  1:515001‑M2/191‑1 (MQ=255)

                                                                                                                                                                                            cctgaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATcttcac                                                               <  2:125148‑M2/190‑1 (MQ=255)

                                                                                                                                                                                               gaaaattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATcttcacttc                                                            <  1:19320‑M2/193‑1 (MQ=255)

                                                                                                                                                                                                  aattgcgggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGCTGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACAAGGTCAATCTTCACTtcatc                                                         >  1:707003‑M2/55‑250 (MQ=255)

                                                                                                                                                                                                        gggtctgtcaggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTTCATCATCAata                                                  <  1:274107‑M2/203‑1 (MQ=255)

                                                                                                                                                                                                                caggtggaaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCAATCCACCAGGTCAATCTTCACTTCATCATCAATATTCAGAtt                                          >  2:661087‑M2/41‑251 (MQ=255)

                                                                                                                                                                                                                       aaacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTTCATCATCAATATTCAGATTGAGCGcc                                   >  1:608307‑M2/34‑251 (MQ=255)

                                                                                                                                                                                                                        aacctgtgaccagaatagcagagagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTACGGCGTGATGTCGTTGACCAGTTCCGCACGCGTAAATACGGGATCTTTCCGCACTGGCACAATAATTAATTTGCCATCCCCCAGGCCAATCTTCACTTCATCATCAATATTCAGATTGAACGCCt                                  >  2:479040‑M2/33‑251 (MQ=255)

                                                                                                                                                                                                                        aacctgtgaccagaatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTTCATCATCAATATTCAGATtgagcgact                                  <  2:637393‑M2/219‑4 (MQ=255)

                                                                                                                                                                                                                              tgaccagaatagaagtgagttaggaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTAGCTTGCGGCTCTCCCCAGTCGATAGTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTTCATCATCAATAGTCAGATTGAGCGCCTGCAtt                             <  2:641697‑M2/224‑1 (MQ=255)

                                                                                                                                                                                                                                    gaatagaagtgatttagtaaACCCGATCGCCAATAGCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTACGTCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGTCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGGCAATCTTCACTTAATCATCAATATTCAGATTGAGCGCCTGCTTTAACGTAg                      <  2:558464‑M2/231‑1 (MQ=255)

                                                                                                                                                                                                                                     aatagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTTCATCATCAATATTCAGATTGAGCGCCTGCATTACCGTAGc                     >  1:450624‑M2/20‑251 (MQ=255)

                                                                                                                                                                                                                                      atagaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTTCATCATCAATATTCAGATTGAGCGCCTGCATTAACGTAGcc                    <  1:185703‑M2/233‑1 (MQ=255)

                                                                                                                                                                                                                                         gaagtgagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTTCATCATCAATATTCAGATTGAGCGCCTGCATTAACGTAGCCggg                 >  1:513518‑M2/16‑251 (MQ=255)

                                                                                                                                                                                                                                         gaagcgagtaagaaaaTCAGATCGCCCATCTCAGGTACGTATCGGCTAACCATTACAAGACTTCCTTATCTTCCGGCTCTCCCCAGTCGATATTCCCGTGGAGGGTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAAAACGGGCTCTTTAAGCACTGGCTCTATAATTAATTTGCCACCCCCCCGGTCAATCTCCACTTCCTCATCAATCTTCTGATTGAGCGCCTGGATTCACGTAGCCggg                 >  2:414953/14‑251 (MQ=255)

                                                                                                                                                                                                                                              gagttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCTCTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTTCATCATCAATATTCAGATTGAGCGCCTGCATTAACGTAGCCGGGATcc             <  1:234433‑M2/240‑1 (MQ=255)

                                                                                                                                                                                                                                               agttagtaaATCAGCTCCCCCATATCGGGGACGTATCGGCTTACCATTACCAGACTTACTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCTGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGTGCTCGTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTTCATCATCAATATTCAGATTGAGCGCCTGCATTAACTTAGCCGGGatcccc           <  2:474200‑M2/242‑3 (MQ=255)

                                                                                                                                                                                                                                                gttagtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTTCATCATCAATATTCAGATTGAGCGCCTGCATTAACGTAGCCGGGATCCGc           >  1:475628‑M2/9‑250 (MQ=255)

                                                                                                                                                                                                                                                  tagtaaATCAGATCGCCAATATCGGGTACGTATCGGCTTACCATTACCAGACTTACTTATCTTTCGGCTCTCCCCAGTCGATATTATCGTGGAGGTTTTACGGCGTGCTGTTGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCAGCCACCAGGTAAATCTTCACTTCATCATCAAAATTCAGATTGCGCGCCTGCATTAACGTATCCGGGCTccgctcc        >  2:670402‑M2/7‑248 (MQ=255)

                                                                                                                                                                                                                                                   agtaaATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTTCATCATCAATATTCAGATTGAGCGCCTGCATTAACGTAGCCGGGATCCGCAccg       >  1:56787‑M2/6‑251 (MQ=255)

                                                                                                                                                                                                                                                       aaTCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTTCATCATCAATATTCAGATTGAGCGCCTGCATTAACGTAGCCGGGATCCGCACCGCCgg   >  1:330018/1‑251 (MQ=255)

                                                                                                                                                                                                                                                        aTCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTACCATCCACCAGGTCAATCTTCACTTCATCATCAATATTCAGATTGAGCGCCTGCATTAACGTATCCGGGATCCGCAccgccgtt  >  1:258315/1‑249 (MQ=255)

GGGCAACTGTTCCTTTCTTCGTTGCTCCTCTTGCCCGCCAGGCGATACTTTTTACCTGATCAGCTAACGCTACGCCATCACGTTCCTGACCGGATAAAACAACTTCGAACGGATATCCTTTTGATTGCGTTGTACAAGGAACACACAGACACATACCTGTTTTGTTGTTGTACATGAAAGGACTCAGGACAACAGCTGGACGATGTCCAGCTTGCTCGCTACCTTTTGTCGGGTCAAAATCAACCCAAATCAGATCGCCCATATCGGGTACGTATCGGCTTACCATTACCAGACTTCCTTATCTTTCGGCTCTCCCCAGTCGATATTCTCGTGGAGGTTTTCCGGCGTGATGTCGTTGACCAGTTCAGCAAGCGTAAATACGGGCTCTTTACGCACTGGCTCAATAATTAATTTGCCATCCACCAGGTCAATCTTCACTTCATCATCAATATTCAGATTGAGCGCCTGCATTAACGTAGCCGGGATCCGCACCGCCGGT  >  NC_000913/2910806‑2911304

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