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.2%	407.4 / NA	124	A162G (GCA→GGA)	btuB	BtuB
Reads supporting (aligned to +/- strand): ref base C (0/1); new base G (71/52); total (71/53)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 138 total aligned reads displayed.

GGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATT  >  NC_000913/4161612‑4161791

ggCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTcaggag                                                                                         <  1:3329798/93‑1 (MQ=255)

ggCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTcaggag                                                                                         <  1:457714/93‑1 (MQ=255)

 gCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCaggagg                                                                                        <  2:3527856/93‑1 (MQ=255)

 gCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCaggagg                                                                                        >  2:4459456/1‑93 (MQ=255)

 gCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCCCGATGAACCCGGAACGGAAATTTCaggagg                                                                                        >  2:4042789/1‑93 (MQ=255)

  ccGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAggg                                                                                       <  1:3435815/93‑1 (MQ=255)

   cgcgCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGt                                                                                      >  1:1266860/1‑93 (MQ=255)

   cgcgCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGt                                                                                      <  1:1872263/93‑1 (MQ=255)

    gcgcTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGTGATGAACCCGGAACGGAAATTTCAGGAGGGTg                                                                                     <  2:1417288/93‑1 (MQ=255)

     cgcTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTgg                                                                                    >  1:1208887/1‑93 (MQ=255)

        tCCGCTGTTTATGGTTCCTATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGa                                                                                 >  1:2385063/1‑93 (MQ=255)

        tCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGa                                                                                 <  1:374696/93‑1 (MQ=255)

         ccGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGaa                                                                                <  1:4130223/93‑1 (MQ=255)

         ccGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGaa                                                                                >  2:2564272/1‑93 (MQ=255)

         ccGCTGTTTATGGTTCCGATGCAATACGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGaa                                                                                >  1:4726386/1‑93 (MQ=255)

          cGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAg                                                                               >  2:3452210/1‑93 (MQ=255)

          cGCTGTTTATGGGTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAg                                                                               >  2:4026769/1‑93 (MQ=255)

            cTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCa                                                                             >  1:4401328/1‑93 (MQ=255)

            cTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCa                                                                             >  2:151689/1‑93 (MQ=255)

             tGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCaa                                                                            >  2:4291733/1‑93 (MQ=255)

               tttATGGTTCCGATGCAATAGTCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATa                                                                          >  1:834945/1‑93 (MQ=255)

               tttATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATa                                                                          >  2:3716063/1‑93 (MQ=255)

               tttATGGTTCCGATGCAATAGGCGGGGCGGTGAATATCATCACGACCCGCGATGAACCCGGAACGGAAATTTCAGGAGGGCGGGGAAGCAAtg                                                                          >  1:271924/1‑92 (MQ=255)

                ttATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAg                                                                         <  1:2918584/93‑1 (MQ=255)

                ttATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAg                                                                         <  2:4382103/93‑1 (MQ=255)

                 tATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGt                                                                        >  2:289214/1‑93 (MQ=255)

                 tATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGt                                                                        <  1:3911622/93‑1 (MQ=255)

                 tATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGt                                                                        >  1:844277/1‑93 (MQ=255)

                  aTGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGtt                                                                       >  1:2083885/1‑93 (MQ=255)

                  aTGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGtt                                                                       >  2:1267253/1‑93 (MQ=255)

                   tGGTTGCGATTCAATAGGCGGGGTGGGTAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTa                                                                      <  1:1422160/93‑1 (MQ=255)

                   tGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTa                                                                      <  1:385570/93‑1 (MQ=255)

                     gTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATc                                                                    >  1:2064161/1‑93 (MQ=255)

                     gTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATc                                                                    >  2:4055095/1‑93 (MQ=255)

                      ttCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCa                                                                   >  1:903246/1‑93 (MQ=255)

                       tCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAg                                                                  >  2:2660095/1‑93 (MQ=255)

                          gATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAAc                                                               >  1:541639/1‑93 (MQ=255)

                           aTGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACt                                                              >  2:1455964/1‑93 (MQ=255)

                            tGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTa                                                             >  2:1282922/1‑93 (MQ=255)

                            tGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTa                                                             >  2:139552/1‑93 (MQ=255)

                            tGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTa                                                             <  1:3633605/93‑1 (MQ=255)

                             gCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTat                                                            <  1:2825158/93‑1 (MQ=255)

                             gCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTat                                                            >  1:2632241/1‑93 (MQ=255)

                             gCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTat                                                            <  2:415145/93‑1 (MQ=255)

                              cAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatg                                                           >  1:4001115/1‑93 (MQ=255)

                               aaTAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatga                                                          <  2:336841/93‑1 (MQ=255)

                               aaTAGGCGGGGTGGTGAAAATCATCATGACGCGCGATGAACCCGGAACGGACATTTCAGGAGGGTGGGGAAACCATAGTTATAAGAACTatga                                                          >  2:3869327/1‑93 (MQ=255)

                                aTAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAATTatgag                                                         >  1:201604/1‑92 (MQ=255)

                                aTAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGGGGGGAGGCAATAGTTATCAGAACTatgat                                                         >  2:1260502/1‑93 (MQ=255)

                                 tAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatgatg                                                        <  2:830484/93‑1 (MQ=255)

                                 tAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatgatg                                                        <  1:436648/93‑1 (MQ=255)

                                   ggCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGtc                                                      <  1:1296270/93‑1 (MQ=255)

                                   ggCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGtc                                                      <  2:3179852/93‑1 (MQ=255)

                                   ggCGGGGTGCTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGtc                                                      <  2:2253809/93‑1 (MQ=255)

                                      ggggTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTctct                                                   <  1:4506155/93‑1 (MQ=255)

                                       gggTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTa                                                  <  2:1084315/93‑1 (MQ=255)

                                        ggtggtGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTAc                                                 >  1:4381188/1‑93 (MQ=255)

                                        ggtggtGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTAc                                                 <  1:633382/93‑1 (MQ=255)

                                          tggtgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgc                                               >  2:125369/1‑93 (MQ=255)

                                          tggtgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgc                                               >  1:927174/1‑93 (MQ=255)

                                               aaTATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaa                                          >  2:3231410/1‑93 (MQ=255)

                                               aaTATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaa                                          >  1:638093/1‑93 (MQ=255)

                                                 tatCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaaca                                        >  2:3894274/1‑93 (MQ=255)

                                                 tatCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaaca                                        >  1:450838/1‑93 (MQ=255)

                                                 tatCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGCGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaaca                                        >  2:3609612/1‑93 (MQ=255)

                                                  atcatcACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaacaa                                       <  2:2327337/93‑1 (MQ=255)

                                                  atcatcACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaacaa                                       >  2:3554359/1‑93 (MQ=255)

                                                  atcatcACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaacaa                                       >  1:2271745/1‑93 (MQ=255)

                                                   tcatcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCaacaac                                      <  1:2202912/93‑1 (MQ=255)

                                                   gcatcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCaacaac                                      <  2:2499392/92‑1 (MQ=255)

                                                      tcaCGACGCGCGCTGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAGGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTgg                                   <  2:4213213/93‑1 (MQ=255)

                                                      tcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTgg                                   <  1:4168992/93‑1 (MQ=255)

                                                      tcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTgg                                   >  1:976275/1‑93 (MQ=255)

                                                        acgtcgCGCGATGAACCCGGAACGGAAATTTCAGGAGGGGGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTgggg                                 >  2:3980502/1‑93 (MQ=255)

                                                        acgacgCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTgggg                                 >  1:720963/1‑93 (MQ=255)

                                                             gcgcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATaa                            <  2:505273/93‑1 (MQ=255)

                                                             gcgcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATaa                            <  1:1455062/93‑1 (MQ=255)

                                                               gcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGa                          >  1:2094398/1‑93 (MQ=255)

                                                                cgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGac                         <  2:2895205/93‑1 (MQ=255)

                                                                  aTGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGacac                       >  1:1511546/1‑93 (MQ=255)

                                                                    gAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACgg                     <  2:4603688/93‑1 (MQ=255)

                                                                      aCCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGt                   >  1:28945/1‑93 (MQ=255)

                                                                        ccGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTaa                 >  2:3836624/1‑93 (MQ=255)

                                                                        ccGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTaa                 >  1:2790513/1‑93 (MQ=255)

                                                                         cggaacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTGTCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAAc                <  2:1937871/93‑1 (MQ=255)

                                                                         cggaacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAAc                >  2:2943915/1‑93 (MQ=255)

                                                                           gaactgaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGAGAAGACACGGGTAACGc              <  2:1898643/93‑1 (MQ=255)

                                                                           gaacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGc              >  2:2445879/1‑93 (MQ=255)

                                                                           gaacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGc              <  1:3623686/93‑1 (MQ=255)

                                                                            aacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCt             <  2:370453/93‑1 (MQ=255)

                                                                            aacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCt             <  2:2758909/93‑1 (MQ=255)

                                                                              cggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGt           >  2:1989316/1‑93 (MQ=255)

                                                                              cggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGt           >  2:3505314/1‑93 (MQ=255)

                                                                              cggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGt           >  2:3341454/1‑93 (MQ=255)

                                                                               ggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGtt          <  2:2308638/93‑1 (MQ=255)

                                                                                   aTTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGc      <  1:2754331/93‑1 (MQ=255)

                                                                                    tttCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCg     >  2:457636/1‑93 (MQ=255)

                                                                                       caggagGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTACCGCTGTTGGGCGAtt  <  2:3545568/93‑1 (MQ=255)

                                                                                       caggagGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGAtt  <  1:1716638/93‑1 (MQ=255)

                                                                                       caggagGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGAtt  <  1:900446/93‑1 (MQ=255)

GGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATT  >  NC_000913/4161612‑4161791

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