BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
MC JC	NC_000913	1,978,503	Δ776 bp	100%		insB1–insA	insB1, insA

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1978503–1979270	1979278	9–776	166 [0]	[1] 164	insB1–insA	insB1, insA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1978502	0 (0.000)	132 (1.040)	80/162	NT	100%	intergenic (‑305/+16)	flhD/insB1	flagellar class II regulon transcriptional activator, with FlhC/IS1 transposase B
*	?	NC_000913	1979279 =	0 (0.000)	132 (1.040)	80/162	NT	100%	intergenic (‑64/‑474)	insA/uspC	IS1 repressor TnpA/universal stress protein

Only 100 of 164 total aligned reads displayed.

GAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGT  >  NC_000913/1979181‑1979378

gtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgta                                                                                                    >  1:3119919‑M2/99‑100 (MQ=255)

  actcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtata                                                                                                  >  2:256938‑M2/97‑100 (MQ=255)

       tcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCa                                                                                             >  2:1958968‑M2/92‑100 (MQ=255)

        cactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCAt                                                                                            <  1:3863503‑M2/10‑1 (MQ=255)

        cactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCAt                                                                                            <  1:3588727‑M2/10‑1 (MQ=255)

         actacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATa                                                                                           <  2:2624270‑M2/11‑1 (MQ=255)

         actacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATa                                                                                           >  1:1320445‑M2/90‑100 (MQ=255)

          ctacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATAc                                                                                          >  2:356429‑M2/89‑100 (MQ=255)

             cacgcacatacaacggaggggggctgcgattatcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATAc                                                                                          >  2:1478067‑M2/86‑97 (MQ=255)

               cgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttg                                                                                     >  1:2811425‑M2/84‑100 (MQ=255)

               cgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttg                                                                                     >  1:832237‑M2/84‑100 (MQ=255)

                   catacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtttt                                                                                 >  1:2498706‑M2/80‑100 (MQ=255)

                    atacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttttg                                                                                <  2:2545275‑M2/22‑1 (MQ=255)

                    atacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttttg                                                                                >  1:616696‑M2/79‑100 (MQ=255)

                    atacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtttt                                                                                 >  1:3789340‑M2/79‑99 (MQ=255)

                     tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtttt                                                                                 >  2:2360176‑M2/78‑98 (MQ=255)

                      acaacggagggtggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                                                              <  2:1181678‑M2/24‑1 (MQ=255)

                      acaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                                                              <  1:1128640‑M2/24‑1 (MQ=255)

                      acaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacgtaacataaatgTATAAGTCATACTTTTGTTTTggg                                                                              <  1:3253741‑M2/24‑1 (MQ=255)

                      acaacagaggggggctgcgattttcaataatgcgtgatgcagatcacacataacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                                                              <  2:3032497‑M2/24‑1 (MQ=255)

                       caacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                             >  2:1338298‑M2/76‑100 (MQ=255)

                       caacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                             >  1:3383892‑M2/76‑100 (MQ=255)

                       caacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                             >  2:449711‑M2/76‑100 (MQ=255)

                       caacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                             >  2:183385‑M2/76‑100 (MQ=255)

                        aacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgtg                                                                            >  2:338167‑M2/75‑100 (MQ=255)

                           ggaggggggctgctattttcaatattgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATATGTCATACTTTTGTTTTggg                                                                              <  2:1759923‑M2/24‑1 (MQ=255)

                                 gggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                                                   <  2:544335‑M2/35‑1 (MQ=255)

                                  ggctgcgtttttcaataatgtgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                                                                  <  1:1729510‑M2/36‑1 (MQ=255)

                                  ggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                                                                  <  1:1992634‑M2/36‑1 (MQ=255)

                                   gctgcgattttcaataatgcgtgatgcagatcacacaaaacactcatttacttaacataaatgTATAATTCATACTTTTGTTTTGGGTGTATTTCAATCt                                                                 <  2:3622459‑M2/37‑1 (MQ=255)

                                   gctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCt                                                                 >  1:666837‑M2/64‑100 (MQ=255)

                                      gcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATCCTTTTGTTTTGGGTGTATTTCAATCTGtt                                                              <  2:2733846‑M2/40‑1 (MQ=255)

                                      gcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGtt                                                              <  2:1443089‑M2/40‑1 (MQ=255)

                                        gattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaa                                                            <  1:963543‑M2/42‑1 (MQ=255)

                                        gattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTa                                                             <  1:803059‑M2/41‑1 (MQ=255)

                                          ttttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaa                                                          <  1:2974132‑M2/44‑1 (MQ=255)

                                           tttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaaa                                                         >  1:3623794‑M2/56‑100 (MQ=255)

                                            ttcaataatgcgtgatgcagatcacacaaaactctcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAg                                                        <  1:2331423‑M2/46‑1 (MQ=255)

                                            ttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAg                                                        <  2:2457777‑M2/46‑1 (MQ=255)

                                               aataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGACATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtt                                                      >  1:3648157‑M2/52‑99 (MQ=255)

                                                 taatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAACAAGttttt                                                   >  2:26368‑M2/50‑100 (MQ=255)

                                                  aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTc                                                  <  1:58104‑M2/52‑1 (MQ=255)

                                                  aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTc                                                  >  2:3503302‑M2/49‑100 (MQ=255)

                                                  aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTc                                                  >  1:2053296‑M2/49‑100 (MQ=255)

                                                  aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGATTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTc                                                  >  1:1080356‑M2/49‑100 (MQ=255)

                                                   atgcgttatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcg                                                 <  2:746213‑M2/53‑1 (MQ=255)

                                                   atgcgtgatgcagatcacacaaaacactctattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATATGTTAAAAAGTTTTTcg                                                 <  2:573475‑M2/53‑1 (MQ=255)

                                                   atgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcg                                                 <  1:3503943‑M2/53‑1 (MQ=255)

                                                   atgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcg                                                 >  1:2898557‑M2/48‑100 (MQ=255)

                                                    tgcgtgatgcagatcacacaaatctctcaattatttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcg                                                 <  1:1046646‑M2/53‑1 (MQ=255)

                                                    tgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAagtttttag                                                 >  2:458286‑M2/47‑97 (MQ=255)

                                                    tgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgc                                                <  2:1992644‑M2/54‑1 (MQ=255)

                                                      cgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgcta                                              <  2:2450898‑M2/56‑1 (MQ=255)

                                                      cgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGATTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgcta                                              >  1:3737137‑M2/45‑100 (MQ=255)

                                                       gtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctac                                             <  1:3246232‑M2/57‑1 (MQ=255)

                                                        tgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacg                                            >  2:524584‑M2/43‑100 (MQ=255)

                                                        tgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacg                                            >  1:310577‑M2/43‑100 (MQ=255)

                                                         gatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgc                                           >  2:3790723‑M2/42‑100 (MQ=255)

                                                         gatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGATGTATTTCAATCTGTCAAAAAGTTTTTcgatacgc                                           >  1:2604560‑M2/42‑100 (MQ=255)

                                                           tgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgct                                          >  2:3304775‑M2/40‑99 (MQ=255)

                                                              agatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagca                                      >  1:438385‑M2/37‑100 (MQ=255)

                                                              agatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagca                                      >  2:2908512‑M2/37‑100 (MQ=255)

                                                                 tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaagc                                   >  2:2968275‑M2/34‑100 (MQ=255)

                                                                 tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaagc                                   <  1:449717‑M2/67‑1 (MQ=255)

                                                                    cacaaaacactcaattgcttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaa                                >  1:943637‑M2/31‑100 (MQ=255)

                                                                      caaaacactcaattacttaacataaatgTATAAGTCATACTTTTTTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaaa                              <  2:3787698‑M2/72‑1 (MQ=255)

                                                                      caaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaaa                              >  2:1435605‑M2/29‑100 (MQ=255)

                                                                      caaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTATAATACAATCTGTTAAAAAGTTTTTCGCTACGCTATCAAGCaaaaa                              >  1:1154081‑M2/29‑100 (MQ=255)

                                                                       naaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAAt                             >  2:2204084‑M2/28‑100 (MQ=255)

                                                                        aaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATg                            <  1:1438148‑M2/74‑1 (MQ=255)

                                                                        aaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATg                            <  1:1968741‑M2/74‑1 (MQ=255)

                                                                           cttccgatctcttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaaa                         <  2:2725612‑M2/77‑1 (MQ=255)

                                                                           cactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaaa                         >  2:234340‑M2/24‑100 (MQ=255)

                                                                            ccgctattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTACAATATGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAAc                        <  2:1607796‑M2/78‑1 (MQ=255)

                                                                             ctcaattacttaacataaatgTATAAGTCATTCTTTTGTTTTGTGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAAc                        <  1:190473‑M2/78‑1 (MQ=255)

                                                                             ctcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAAc                        <  1:2796391‑M2/78‑1 (MQ=255)

                                                                             ctcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACa                       >  2:2293985‑M2/22‑100 (MQ=255)

                                                                             ctcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACa                       <  1:183390‑M2/79‑1 (MQ=255)

                                                                              tcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAATGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAg                      <  2:3110853‑M2/80‑1 (MQ=255)

                                                                               caattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAg                      >  1:3289428‑M2/20‑99 (MQ=255)

                                                                               caattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGAAAGAAAAAAAAAATGAAACAgg                     <  2:3920297‑M2/81‑1 (MQ=255)

                                                                                aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGa                    <  1:2499808‑M2/82‑1 (MQ=255)

                                                                                aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGa                    >  1:1901997‑M2/19‑100 (MQ=255)

                                                                                 attacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaa                   >  1:3224466‑M2/18‑100 (MQ=255)

                                                                                   tacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaata                 <  1:3376244‑M2/85‑1 (MQ=255)

                                                                                     cttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAagcaaaaatgaaaagatcggaag               >  1:2725636‑M2/14‑90 (MQ=255)

                                                                                     cttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat               >  2:2857798‑M2/14‑100 (MQ=255)

                                                                                      ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc              <  1:2263342‑M2/88‑1 (MQ=255)

                                                                                         acataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaa           <  1:267836‑M2/91‑1 (MQ=255)

                                                                                         acataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaa           <  2:2960116‑M2/91‑1 (MQ=255)

                                                                                         aaaTAAATGTAAAAGACAAACTTTTGTTTTGGATGTATTTCAATCTGTTAAAAAGTTTTTCGCTAAAAAAAAAAAAAAAAATGAAACAGGAATAATCGaa           <  2:2792906/98‑1 (MQ=255)

                                                                                          cataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa          >  2:2783921‑M2/9‑100 (MQ=255)

                                                                                          cataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa          >  1:911239‑M2/9‑100 (MQ=255)

                                                                                          cataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa          >  2:1061871‑M2/9‑100 (MQ=255)

                                                                                            taaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATg        <  1:1943547‑M2/94‑1 (MQ=255)

                                                                                            taaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATg        >  1:3917582‑M2/7‑100 (MQ=255)

                                                                                              aatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATggg      <  1:2999681‑M2/96‑1 (MQ=255)

                                                                                                 gTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAACCAGGAATAATCGAAATGGGATg   >  2:3196769‑M2/2‑100 (MQ=255)

                                                                                                 gTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg   <  1:3047323‑M2/99‑1 (MQ=255)

                                                                                                  tataAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGt  >  1:3610351/1‑100 (MQ=255)

GAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGT  >  NC_000913/1979181‑1979378

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