BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,161,702	C→G	A162G (GCA→GGA)	btuB →	BtuB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,161,702	0	C	G	99.8%	1976.4 / ‑4.6	490	A162G (GCA→GGA)	btuB	BtuB
Reads supporting (aligned to +/- strand): ref base C (0/1); new base G (249/240); total (249/241)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.96e-01

Only 100 of 515 total aligned reads displayed.

GTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTAT  >  NC_000913/4161609‑4161793

 tGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCaggagg                                                                                          <  2:6767321/96‑1 (MQ=255)

      cgcgCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTggg                                                                                     <  2:7732992/96‑1 (MQ=255)

        cgcTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGGGGGGa                                                                                   >  1:771273/1‑96 (MQ=255)

         gcTTCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGaa                                                                                  <  1:2498654/96‑1 (MQ=255)

         gcTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGaa                                                                                  <  2:9525429/96‑1 (MQ=255)

           tCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCTCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGc                                                                                <  1:3265920/96‑1 (MQ=255)

              gCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATAACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGAAAt                                                                             >  2:10956056/1‑96 (MQ=255)

              gCTGTTTATGGTTCCGATGCAATAGGCGGGGGGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAAt                                                                             >  1:12673976/1‑96 (MQ=255)

                tGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGACCCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAg                                                                           <  2:3157012/96‑1 (MQ=255)

                tGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAAGGGAAATTTCAGGAGGGAGGGGAAGCAATAg                                                                           >  1:5862257/1‑96 (MQ=255)

                cgTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAg                                                                           >  1:11319546/2‑96 (MQ=255)

                 gTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGt                                                                          <  2:5198808/96‑1 (MQ=255)

                  tttATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGtt                                                                         <  2:12872096/96‑1 (MQ=255)

                   ttATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTa                                                                        <  2:14281506/96‑1 (MQ=255)

                    tATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAGCGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTAt                                                                       <  1:6799194/96‑1 (MQ=255)

                     aTGGTTCCGATGCAATAGGCGGGGGGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATc                                                                      >  1:10014568/1‑96 (MQ=255)

                       ggTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAg                                                                    <  2:7393441/96‑1 (MQ=255)

                       ggTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAg                                                                    >  1:8121129/1‑96 (MQ=255)

                       ggTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAg                                                                    <  1:11006729/96‑1 (MQ=255)

                       ggTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAg                                                                    <  2:7923811/96‑1 (MQ=255)

                       ggTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTAGCAg                                                                    >  1:9534532/1‑96 (MQ=255)

                        gTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGa                                                                   >  1:6987996/1‑96 (MQ=255)

                          tCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGCGGGTGGGGAAGCAATAGTTATCAGAAc                                                                 <  2:12711768/96‑1 (MQ=255)

                             gATGCAATAGGCGGGGGGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCATGAGGGTGGGGAAGCAATAGTAATCAGAACTat                                                              >  1:4841376/1‑96 (MQ=255)

                              aTGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatg                                                             <  1:3438243/96‑1 (MQ=255)

                              aTGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatg                                                             <  1:5977952/96‑1 (MQ=255)

                              aTGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatg                                                             >  2:647482/1‑96 (MQ=255)

                               tGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatga                                                            >  2:10628340/1‑96 (MQ=255)

                                gCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAATCAATAGTTATCAGAACTatgat                                                           <  2:6026455/96‑1 (MQ=255)

                                gCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatgat                                                           <  1:2800543/96‑1 (MQ=255)

                                gCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatgat                                                           <  2:4324287/96‑1 (MQ=255)

                                 cAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatgatg                                                          >  2:10343101/1‑96 (MQ=255)

                                 cAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatgatg                                                          >  2:3780381/1‑96 (MQ=255)

                                  aaTAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGt                                                         >  2:274867/1‑96 (MQ=255)

                                  aaTAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGt                                                         <  2:3370378/96‑1 (MQ=255)

                                  aaTAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGGGGGGAAGCAATAGTTATCAGAACTATGATGt                                                         >  1:2680938/1‑96 (MQ=255)

                                    tAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGtct                                                       <  1:10717072/96‑1 (MQ=255)

                                    tAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGtct                                                       >  2:8455074/1‑96 (MQ=255)

                                    tAGGCGGGGTGGGGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGGGGGGAAGCAATAGTTATCAGAACTATGATGtct                                                       >  1:14046823/1‑96 (MQ=255)

                                     aGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGtctc                                                      <  2:10101547/96‑1 (MQ=255)

                                      ggCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTctct                                                     >  2:4049213/1‑96 (MQ=255)

                                        cGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTAc                                                   <  2:13742455/96‑1 (MQ=255)

                                         ggggTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGTTGTCTCTACg                                                  >  2:13191953/1‑96 (MQ=255)

                                         ggggTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACg                                                  >  1:7781447/1‑96 (MQ=255)

                                         ggggTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACg                                                  >  2:11615086/1‑96 (MQ=255)

                                           ggtggtGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgca                                                >  1:5119829/1‑96 (MQ=255)

                                            gtggtgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgcag                                               <  2:9165167/96‑1 (MQ=255)

                                             tggtgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgcagc                                              >  2:3513356/1‑96 (MQ=255)

                                              ggtgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgcagca                                             <  1:3119933/96‑1 (MQ=255)

                                               gtgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaa                                            >  2:11748820/1‑96 (MQ=255)

                                                 gAATATCATCACGCCGCGCGATGAACTCGGAACGGAAATCTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGAGGTCTCTACGCAGcaaaa                                          <  1:11595979/96‑3 (MQ=255)

                                                 gAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaaca                                          <  2:6908357/96‑1 (MQ=255)

                                                 gAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaaca                                          <  1:9970307/96‑1 (MQ=255)

                                                 gAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaaca                                          >  1:12128385/1‑96 (MQ=255)

                                                 gAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGGGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaaca                                          >  2:7609101/1‑96 (MQ=255)

                                                  aaTATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaacaa                                         >  2:585844/1‑96 (MQ=255)

                                                  aaTATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaa                                            <  1:9839618/93‑1 (MQ=255)

                                                   atatCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCCACGCAGCaacaac                                        <  1:13441766/96‑1 (MQ=255)

                                                   atatCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGGGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCaacaac                                        >  2:13420394/1‑96 (MQ=255)

                                                   atatCATCACGACGCGCGATGAAACCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCaacaac                                        <  1:10212617/96‑1 (MQ=255)

                                                     atcatcACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTg                                      <  2:5406796/96‑1 (MQ=255)

                                                            cgacgCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATa                               >  1:10485766/1‑96 (MQ=255)

                                                            ccacgCGCGATGAACCCGGGACGGAAATTTCAGGAGGGGGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATa                               >  1:9008608/3‑96 (MQ=255)

                                                             gacgCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATaa                              >  2:121981/1‑96 (MQ=255)

                                                             gacgCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATaa                              >  2:5597135/1‑96 (MQ=255)

                                                             gacgCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATaa                              >  1:3722514/1‑96 (MQ=255)

                                                               cgcgcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGa                            <  2:11332719/96‑1 (MQ=255)

                                                                gcgcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTgggg                                   >  2:2644128/1‑88 (MQ=255)

                                                                gcgcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGac                           <  1:1870673/96‑1 (MQ=255)

                                                                 cgcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGaca                          >  2:708603/1‑96 (MQ=255)

                                                                   cgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACg                        <  1:14428456/96‑1 (MQ=255)

                                                                    gATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTCCGCAGCAACAACTGGGGGATAAGACACgg                       <  2:8875117/96‑1 (MQ=255)

                                                                    gATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACgg                       >  1:6608422/1‑96 (MQ=255)

                                                                    gATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACgg                       <  1:8417980/96‑1 (MQ=255)

                                                                     aTGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACggg                      >  2:4831804/1‑96 (MQ=255)

                                                                     aTGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACggg                      >  1:3798807/1‑96 (MQ=255)

                                                                        aaCCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTaa                   >  1:7297313/1‑96 (MQ=255)

                                                                         aCCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAAc                  >  2:4983368/1‑96 (MQ=255)

                                                                         aCCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAAc                  <  1:445613/96‑1 (MQ=255)

                                                                          cccGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTaagg                 <  1:11869601/96‑3 (MQ=255)

                                                                            cggaacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCt               >  1:630919/1‑96 (MQ=255)

                                                                            cggaacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCt               >  2:1582128/1‑96 (MQ=255)

                                                                            cggaacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCt               <  2:1997062/96‑1 (MQ=255)

                                                                             ggaacggaaATTTCAGGAGGGTGGGGAAGCAATAGGTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATATGACACGGGTAACGCTg              >  1:5525943/1‑96 (MQ=255)

                                                                              gaacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCGGCAACAACTGGGGGATAAGACACGGGTAACGCTGt             >  1:4275477/1‑96 (MQ=255)

                                                                              gaacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGt             >  1:13710604/1‑96 (MQ=255)

                                                                               aacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGtt            <  1:8793871/96‑1 (MQ=255)

                                                                               aacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGtt            <  1:13118389/96‑1 (MQ=255)

                                                                               aacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGAGATAAGACACGGGTAACGCTGtt            >  1:12491777/1‑96 (MQ=255)

                                                                                acggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGtag           <  1:7555095/96‑3 (MQ=255)

                                                                                  ggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTggg         >  2:11120849/1‑96 (MQ=255)

                                                                                  ggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTggg         >  2:1749825/1‑96 (MQ=255)

                                                                                   gaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGc        <  2:2462877/96‑1 (MQ=255)

                                                                                   gaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGc        <  1:7192129/96‑1 (MQ=255)

                                                                                    aaaTTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCg       <  1:5322016/96‑1 (MQ=255)

                                                                                     aaTTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGTAAAGCCA‑GGGGAAAGCCGTTggggggc     >  2:11973635/1‑92 (MQ=255)

                                                                                       tttCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGAtt    >  1:1080613/1‑96 (MQ=255)

                                                                                        ttCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTa   >  2:2177324/1‑96 (MQ=255)

                                                                                        ttCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTa   <  1:2228545/96‑1 (MQ=255)

                                                                                         tCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTAt  >  2:8657992/1‑96 (MQ=255)

GTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTAT  >  NC_000913/4161609‑4161793

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 16 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`