BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	1,976,527	Δ776 bp		insB5–insA5	insB5, insA5

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1976527–1977294	1977302	9–776	178 [0]	[0] 177	insB5–insA5	insB5,insA5

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1976526	0 (0.000)	157 (0.690)	91/168	2.0	100%	intergenic (‑305/+16)	b1892/insB5	FlhD/IS1 protein InsB
*	?	NC_000913	1977303 =	0 (0.000)	157 (0.690)	91/168	2.0	100%	intergenic (‑64/‑474)	insA5/yecG	IS1 protein InsA/putative regulator

Only 100 of 177 total aligned reads displayed.

GAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAT  >  NC_000913/1977211‑1977395

gtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgt                                                                                              <  1:1951349‑M2/1‑1 (MQ=255)

  cactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtat                                                                                            <  1:4300426‑M2/3‑1 (MQ=255)

  cactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtat                                                                                            <  2:315986‑M2/3‑1 (MQ=255)

   actacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtata                                                                                           <  1:5876005‑M2/4‑1 (MQ=255)

        acgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCa                                                                                      >  2:6370066‑M2/85‑93 (MQ=255)

        acgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCa                                                                                      <  1:6100853‑M2/9‑1 (MQ=255)

         cgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCAt                                                                                     <  2:4507339‑M2/10‑1 (MQ=255)

          gcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATa                                                                                    >  1:2317310‑M2/83‑93 (MQ=255)

          gcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATa                                                                                    <  2:3993861‑M2/11‑1 (MQ=255)

            acatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACt                                                                                  <  1:1611342‑M2/13‑1 (MQ=255)

            acatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACt                                                                                  <  1:3902394‑M2/13‑1 (MQ=255)

            acatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACt                                                                                  >  1:3089876‑M2/81‑93 (MQ=255)

               tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACtttt                                                                               <  1:6095359‑M2/16‑1 (MQ=255)

               tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACtttt                                                                               <  1:214877‑M2/16‑1 (MQ=255)

               tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACtttt                                                                               <  1:538844‑M2/16‑1 (MQ=255)

               tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACtttt                                                                               <  2:1913269‑M2/16‑1 (MQ=255)

               tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACtttt                                                                               >  1:358193‑M2/78‑93 (MQ=255)

                  aacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtt                                                                            >  2:6130116‑M2/75‑93 (MQ=255)

                   acggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttt                                                                           <  2:6439691‑M2/20‑1 (MQ=255)

                   acggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttt                                                                           >  1:4777348‑M2/74‑93 (MQ=255)

                   acggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttt                                                                           <  1:4324120‑M2/20‑1 (MQ=255)

                    cggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtttt                                                                          <  2:2857653‑M2/21‑1 (MQ=255)

                     ggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttttg                                                                         <  2:4171941‑M2/22‑1 (MQ=255)

                      gaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                        <  1:2632048‑M2/23‑1 (MQ=255)

                       aggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                                                       <  2:4266415‑M2/24‑1 (MQ=255)

                       aggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                                                       <  2:3750375‑M2/24‑1 (MQ=255)

                        ggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                      <  1:5864230‑M2/25‑1 (MQ=255)

                          ggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgtgt                                                                    >  2:4790905‑M2/67‑93 (MQ=255)

                          ggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgtgt                                                                    >  2:2466059‑M2/67‑93 (MQ=255)

                          ggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgtgt                                                                    <  1:6118779‑M2/27‑1 (MQ=255)

                           gggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTa                                                                   >  1:1558267‑M2/66‑93 (MQ=255)

                            ggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAt                                                                  >  1:6032863‑M2/65‑93 (MQ=255)

                            ggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAt                                                                  <  2:3211628‑M2/29‑1 (MQ=255)

                               tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTc                                                               <  2:1092015‑M2/32‑1 (MQ=255)

                               tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTc                                                               <  2:301118‑M2/32‑1 (MQ=255)

                                   attttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                                                           <  2:6319778‑M2/36‑1 (MQ=255)

                                   attttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                                                           <  2:3749508‑M2/36‑1 (MQ=255)

                                    ttttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCt                                                          <  1:1313159‑M2/37‑1 (MQ=255)

                                    ttttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCt                                                          <  1:260492‑M2/37‑1 (MQ=255)

                                      ttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGt                                                        <  2:1930566‑M2/39‑1 (MQ=255)

                                      ttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGt                                                        >  2:1958767‑M2/55‑93 (MQ=255)

                                       tcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGACATACTTTTGTTTTGGGTGTATTTCAATCTGtt                                                       <  2:5538894‑M2/40‑1 (MQ=255)

                                         aataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaa                                                     >  1:6041237‑M2/52‑93 (MQ=255)

                                         aataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaa                                                     >  2:6467140‑M2/52‑93 (MQ=255)

                                          ataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaa                                                    >  2:3210602‑M2/51‑93 (MQ=255)

                                          ataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaa                                                    >  2:1623995‑M2/51‑93 (MQ=255)

                                            aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaaa                                                  >  2:3564172‑M2/49‑93 (MQ=255)

                                             atgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAg                                                 >  2:2511126‑M2/48‑93 (MQ=255)

                                             atgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAg                                                 >  2:2849557‑M2/48‑93 (MQ=255)

                                              tgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGt                                                <  2:4877451‑M2/47‑1 (MQ=255)

                                               gcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtt                                               >  1:6204471‑M2/46‑93 (MQ=255)

                                               gcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtt                                               <  2:4921635‑M2/48‑1 (MQ=255)

                                               gcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtt                                               >  2:5459431‑M2/46‑93 (MQ=255)

                                                cgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttt                                              <  1:16971‑M2/49‑1 (MQ=255)

                                                  tgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt                                            <  1:2391833‑M2/51‑1 (MQ=255)

                                                     tgcagatcacacaaaacactcaattacttaacataaatgTATATGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgc                                         <  1:6375122‑M2/54‑1 (MQ=255)

                                                     tgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgc                                         <  2:255116‑M2/54‑1 (MQ=255)

                                                      gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgct                                        >  1:4469915‑M2/39‑93 (MQ=255)

                                                      gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgct                                        >  2:4430959‑M2/39‑93 (MQ=255)

                                                          atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgc                                    <  1:5645212‑M2/59‑1 (MQ=255)

                                                           tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgct                                   >  2:948873‑M2/34‑93 (MQ=255)

                                                            cacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                  <  1:5549741‑M2/61‑1 (MQ=255)

                                                             acacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTag                                 >  1:285488‑M2/32‑93 (MQ=255)

                                                             acacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTag                                 >  1:995566‑M2/32‑93 (MQ=255)

                                                              cacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagc                                >  2:1270710‑M2/31‑93 (MQ=255)

                                                              cacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagc                                >  1:5288384‑M2/31‑93 (MQ=255)

                                                               acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagca                               <  2:5684240‑M2/64‑1 (MQ=255)

                                                                caaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaa                              <  1:1671824‑M2/65‑1 (MQ=255)

                                                                caaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaa                              <  1:731737‑M2/65‑1 (MQ=255)

                                                                caaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaa                              <  2:4979937‑M2/65‑1 (MQ=255)

                                                                 aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaag                             <  1:3397789‑M2/66‑1 (MQ=255)

                                                                  aaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaagc                            <  1:1591202‑M2/67‑1 (MQ=255)

                                                                  aaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGAATTTCGCTACGCTagcaagc                            >  1:3510150‑M2/27‑93 (MQ=255)

                                                                    acactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAgcaagcaa                          >  2:1606058‑M2/25‑93 (MQ=255)

                                                                      actcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaa                        >  2:3475167‑M2/23‑93 (MQ=255)

                                                                          aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGa                    >  2:6174010‑M2/19‑93 (MQ=255)

                                                                           attacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaa                   <  2:1887779‑M2/76‑1 (MQ=255)

                                                                            ttacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaaa                  >  1:6288801‑M2/17‑93 (MQ=255)

                                                                             tacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAAc                 <  2:3328875‑M2/78‑1 (MQ=255)

                                                                             tacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAAc                 >  1:2571835‑M2/16‑93 (MQ=255)

                                                                             tacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAAc                 >  1:4772472‑M2/16‑93 (MQ=255)

                                                                              acttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAGACa                <  1:2119880‑M2/79‑1 (MQ=255)

                                                                               cttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAg               >  1:2313539‑M2/14‑93 (MQ=255)

                                                                                  aacataaatgTATAAGTCATTCTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaa            <  1:408019‑M2/83‑1 (MQ=255)

                                                                                  aacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAATTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaa            >  2:5858937‑M2/11‑93 (MQ=255)

                                                                                  aacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaa            >  2:1678400‑M2/11‑93 (MQ=255)

                                                                                  aacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaa            <  1:3254803‑M2/83‑1 (MQ=255)

                                                                                    cataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaata          <  2:2971595‑M2/85‑1 (MQ=255)

                                                                                    cataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaata          >  1:958842‑M2/9‑93 (MQ=255)

                                                                                     ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa         <  1:5205448‑M2/86‑1 (MQ=255)

                                                                                     ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa         <  1:1579105‑M2/86‑1 (MQ=255)

                                                                                      taaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat        <  2:3082003‑M2/87‑1 (MQ=255)

                                                                                      taaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat        >  1:1403220‑M2/7‑93 (MQ=255)

                                                                                      taaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat        >  1:170584‑M2/7‑93 (MQ=255)

                                                                                       aaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc       <  1:4912944‑M2/88‑1 (MQ=255)

                                                                                         atgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGa     <  1:790744‑M2/90‑1 (MQ=255)

                                                                                         atgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGa     >  2:2352766‑M2/4‑93 (MQ=255)

                                                                                          tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaa    >  1:1065835‑M2/3‑93 (MQ=255)

                                                                                           gTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa   >  1:6143554‑M2/2‑93 (MQ=255)

                                                                                            tataAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt  >  1:51665/1‑93 (MQ=255)

GAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAT  >  NC_000913/1977211‑1977395

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