BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	66,533	3540 bp→CCCGAGACCCGCACGCGTCGACCTGCA		[araD]–araB	[araD], araA, araB

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	66533	70072	3540	239 [0]	[0] 244	[araD]–araB	[araD],araA,araB

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 66532	0 (0.000)	178 (1.030) +CCCGAGACCCGCACGCGTCGACCTGCA	87/130	0.3	100%	coding (19/696 nt)	araD	L‑ribulose‑5‑phosphate 4‑epimerase
*	?	NC_000913	70073 =	0 (0.000)	178 (1.030) +CCCGAGACCCGCACGCGTCGACCTGCA	87/130	0.3	100%	intergenic (‑25/‑314)	araB/araC	L‑ribulokinase/transcriptional regulator for ara operon

Only 100 of 244 total aligned reads displayed.

GCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGA  >  NC_000913/69984‑70165

ccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTa                                                                                           <  2:4328286‑M2/4‑1 (MQ=255)

  ccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATg                                                                                         >  1:5660980‑M2/88‑93 (MQ=255)

   catgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATgg                                                                                        >  1:5350161‑M2/87‑93 (MQ=255)

   catgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATgg                                                                                        >  1:2483714‑M2/87‑93 (MQ=255)

   catgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATgg                                                                                        <  1:5458760‑M2/7‑1 (MQ=255)

    atgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGga                                                                                       <  1:4390438‑M2/8‑1 (MQ=255)

    atgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGga                                                                                       >  1:4752025‑M2/86‑93 (MQ=255)

    atgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGga                                                                                       <  1:2043569‑M2/8‑1 (MQ=255)

     tgtgagcgtgaccaggttgtgtttttgcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGGATGgag                                                                                      >  1:2210135‑M2/85‑93 (MQ=255)

     tgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGgag                                                                                      >  2:5668767‑M2/85‑93 (MQ=255)

     tgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGgag                                                                                      >  2:2155769‑M2/85‑93 (MQ=255)

       tgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGaa                                                                                    <  1:560886‑M2/11‑1 (MQ=255)

       tgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGGATGGAGaa                                                                                    >  1:90305‑M2/83‑93 (MQ=255)

       tgagcgtgaccaggttgtgttttggcagcgccaggatggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTAtggagga                                                                                    >  2:4648506‑M2/83‑91 (MQ=255)

        gagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatgcccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGaaa                                                                                   <  2:745705‑M2/12‑1 (MQ=255)

         agcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAAc                                                                                  >  1:3268496‑M2/81‑93 (MQ=255)

           cgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAg                                                                                <  2:2429006‑M2/15‑1 (MQ=255)

           cgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAg                                                                                <  2:3700405‑M2/15‑1 (MQ=255)

           cgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAg                                                                                <  1:1096801‑M2/15‑1 (MQ=255)

             tgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTa                                                                              <  1:3483787‑M2/17‑1 (MQ=255)

             tgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTa                                                                              <  2:4930788‑M2/17‑1 (MQ=255)

             tgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTa                                                                              >  2:2740941‑M2/77‑93 (MQ=255)

              gaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTag                                                                             <  1:1142910‑M2/18‑1 (MQ=255)

              gaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTag                                                                             <  1:2928775‑M2/18‑1 (MQ=255)

              gaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTag                                                                             <  1:4599918‑M2/18‑1 (MQ=255)

               accaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTaga                                                                            >  2:5257755‑M2/75‑93 (MQ=255)

               accaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTaga                                                                            >  1:629832‑M2/75‑93 (MQ=255)

               accaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTaga                                                                            <  1:3033814‑M2/19‑1 (MQ=255)

                ccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTagag                                                                           >  1:1856770‑M2/74‑93 (MQ=255)

                ccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTagag                                                                           >  1:6139271‑M2/74‑93 (MQ=255)

                ccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTagag                                                                           <  1:4869003‑M2/20‑1 (MQ=255)

                ccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTagag                                                                           <  1:5189592‑M2/20‑1 (MQ=255)

                 caggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTagaga                                                                          <  2:2359759‑M2/21‑1 (MQ=255)

                 caggttgtgttgtggcagcgcgaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTagaga                                                                          <  1:4975485‑M2/21‑1 (MQ=255)

                  aggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTagagag                                                                         >  2:416897‑M2/72‑93 (MQ=255)

                   ggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGt                                                                        <  1:2569689‑M2/23‑1 (MQ=255)

                   ggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGt                                                                        <  1:4581845‑M2/23‑1 (MQ=255)

                     ttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTg                                                                      >  1:596519‑M2/69‑93 (MQ=255)

                     ttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGGTg                                                                      >  2:1209093‑M2/69‑93 (MQ=255)

                       gtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCg                                                                    >  2:5190890‑M2/67‑93 (MQ=255)

                         gttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGAt                                                                  <  2:3051320‑M2/29‑1 (MQ=255)

                           tttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATaa                                                                <  1:1957382‑M2/31‑1 (MQ=255)

                             tggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATaaaa                                                              >  2:3338118‑M2/61‑93 (MQ=255)

                              ggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATaaaaa                                                             <  2:1545008‑M2/34‑1 (MQ=255)

                              ggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATaaaaa                                                             <  2:2701590‑M2/34‑1 (MQ=255)

                                 agcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCg                                                          >  2:4286989‑M2/57‑93 (MQ=255)

                                   cgccaggttggcttctaatacccggcgccggagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTc                                                        <  1:5253349‑M2/39‑1 (MQ=255)

                                   cgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTc                                                        <  2:5794850‑M2/39‑1 (MQ=255)

                                    gccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCa                                                       <  2:3962834‑M2/40‑1 (MQ=255)

                                     ccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGGgtccg                                                      >  2:4634742‑M2/53‑91 (MQ=255)

                                     ccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAg                                                      >  1:3420237‑M2/53‑93 (MQ=255)

                                     ccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAg                                                      >  2:6024458‑M2/53‑93 (MQ=255)

                                     ccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAg                                                      <  2:6187971‑M2/41‑1 (MQ=255)

                                       aggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGt                                                    >  2:535654‑M2/51‑93 (MQ=255)

                                        ggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTa                                                   >  2:3603568‑M2/50‑93 (MQ=255)

                                        ggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTa                                                   >  1:4214976‑M2/50‑93 (MQ=255)

                                        ggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTa                                                   >  2:674400‑M2/50‑93 (MQ=255)

                                        ggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTa                                                   >  2:935162‑M2/50‑93 (MQ=255)

                                        ggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTa                                                   >  1:2197756‑M2/50‑93 (MQ=255)

                                        ggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTa                                                   <  1:2666518‑M2/44‑1 (MQ=255)

                                         gttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAg                                                  >  2:249836‑M2/49‑93 (MQ=255)

                                          ttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAgg                                                 <  2:2606729‑M2/46‑1 (MQ=255)

                                           tggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGa                                                <  1:1283068‑M2/47‑1 (MQ=255)

                                            ggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGAt                                               <  2:5152207‑M2/48‑1 (MQ=255)

                                             gcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATc                                              >  2:1949330‑M2/45‑93 (MQ=255)

                                              cttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATcc                                             >  1:5178503‑M2/44‑93 (MQ=255)

                                               ttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCg                                            >  2:5619523‑M2/43‑93 (MQ=255)

                                                tctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGc                                           <  1:4165755‑M2/52‑1 (MQ=255)

                                                  taatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTa                                         <  2:4236335‑M2/54‑1 (MQ=255)

                                                  taatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTa                                         <  2:6230986‑M2/54‑1 (MQ=255)

                                                   aatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTaa                                        >  1:3814549‑M2/39‑93 (MQ=255)

                                                   aatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTaa                                        <  2:3460444‑M2/55‑1 (MQ=255)

                                                    atacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAAt                                       >  2:3524339‑M2/38‑93 (MQ=255)

                                                    atacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGACAGGTAGGATCCGCTAAt                                       >  1:5009309‑M2/38‑93 (MQ=255)

                                                     tacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATc                                      >  1:6215100‑M2/37‑93 (MQ=255)

                                                        ccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTa                                   >  2:2002252‑M2/34‑93 (MQ=255)

                                                         cggcgcgcgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTAt                                  <  2:4843391‑M2/61‑1 (MQ=255)

                                                         cggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTAt                                  >  1:5239797‑M2/33‑93 (MQ=255)

                                                          ggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATg                                 >  2:5696853‑M2/32‑93 (MQ=255)

                                                          cgcgcccgacacccgcacgcgtccacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATg                                 <  2:4468414‑M2/62‑1 (MQ=255)

                                                           gcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATgg                                <  2:1655175‑M2/63‑1 (MQ=255)

                                                           gcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATgg                                <  2:4036707‑M2/63‑1 (MQ=255)

                                                           gcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATgg                                <  2:1480662‑M2/63‑1 (MQ=255)

                                                            cgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGa                               <  1:626138‑M2/64‑1 (MQ=255)

                                                              cccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATa                             >  1:1329907‑M2/28‑93 (MQ=255)

                                                               ccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATaa                            >  1:3567267‑M2/27‑93 (MQ=255)

                                                                     cccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGc                      <  1:5081979‑M2/73‑1 (MQ=255)

                                                                      ccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCt                     <  1:1021125‑M2/74‑1 (MQ=255)

                                                                      ccgcacgcctcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCt                     <  1:3107672‑M2/74‑1 (MQ=255)

                                                                       cgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTa                    <  1:5548538‑M2/75‑1 (MQ=255)

                                                                         cacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATg                  <  2:545625‑M2/77‑1 (MQ=255)

                                                                              gtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATa             <  1:6029118‑M2/82‑1 (MQ=255)

                                                                               tcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAg            >  1:4355781‑M2/11‑93 (MQ=255)

                                                                               tcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAg            >  1:321548‑M2/11‑93 (MQ=255)

                                                                                   cctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCaaa        <  2:1167574‑M2/87‑1 (MQ=255)

                                                                                   cctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCaaa        <  1:2393702‑M2/87‑1 (MQ=255)

                                                                                    ctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAg       <  2:635549‑M2/88‑1 (MQ=255)

                                                                                      gcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAgtg     >  2:3281015‑M2/4‑93 (MQ=255)

                                                                                      gcaGGTATCGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAgtg     <  1:4275276‑M2/90‑1 (MQ=255)

                                                                                       caGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAagtgg    >  2:28986‑M2/3‑92 (MQ=255)

GCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGA  >  NC_000913/69984‑70165

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 19 ≤` `ATCG/ATCG` `< 28 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`