BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	158,711	C→T	R128H (CGC→CAC)	pcnB ←	poly(A) polymerase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	158,711	0	C	T	99.6%	911.4 / NA	258	R128H (CGC→CAC)	pcnB	poly(A) polymerase I
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (142/115); total (143/115)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 266 total aligned reads displayed.

TACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATGA  >  NC_000913/158621‑158800

tACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCACgtgt                                                                                         >  2:4055873/1‑93 (MQ=255)

 aCAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCgtgtg                                                                                        >  1:4219633/1‑92 (MQ=255)

 aCAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGtt                                                                                        <  2:533786/93‑1 (MQ=255)

   aGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTgg                                                                                      <  1:4686700/93‑1 (MQ=255)

   aGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTgg                                                                                      >  1:4098858/1‑93 (MQ=255)

   aGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTgg                                                                                      <  2:7050475/93‑1 (MQ=255)

    ggCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTggg                                                                                     >  2:4906099/1‑93 (MQ=255)

    ggCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTggg                                                                                     >  1:2287138/1‑93 (MQ=255)

     gCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTgggg                                                                                    >  2:5936028/1‑92 (MQ=255)

       tGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAgg                                                                                  >  1:390857/1‑93 (MQ=255)

         ttGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTc                                                                                >  1:319771/1‑93 (MQ=255)

          tGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCg                                                                               >  1:5060470/1‑93 (MQ=255)

          tGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCg                                                                               >  2:2343951/1‑93 (MQ=255)

           gATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGt                                                                              <  2:5183831/93‑1 (MQ=255)

           gATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGt                                                                              >  1:2984242/1‑93 (MQ=255)

             tAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGc                                                                            >  1:7541099/1‑93 (MQ=255)

             tAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGc                                                                            <  2:1093729/93‑1 (MQ=255)

             tAGTGAAATCGCGGCGCTGGGCGTCTTATTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGc                                                                            >  1:579769/1‑93 (MQ=255)

               gTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCgg                                                                          >  1:1784479/1‑93 (MQ=255)

               gTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCgg                                                                          <  2:5584785/93‑1 (MQ=255)

                tGAACTCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGAAACATTCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGt                                                                         <  2:6920731/93‑1 (MQ=255)

                tGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGt                                                                         <  2:3029476/93‑1 (MQ=255)

                tGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGt                                                                         >  1:3580053/1‑93 (MQ=255)

                tGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGt                                                                         <  2:3364404/93‑1 (MQ=255)

                tGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGt                                                                         <  1:1296526/93‑1 (MQ=255)

                      cgcgGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGa                                                                   <  1:6926568/93‑1 (MQ=255)

                      cgcgGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGa                                                                   <  2:6745412/93‑1 (MQ=255)

                       gcggcgCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGAc                                                                  >  2:2615295/1‑93 (MQ=255)

                        cggcgCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACg                                                                 >  2:4626362/1‑93 (MQ=255)

                         ggcgCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGt                                                                <  2:5469441/93‑1 (MQ=255)

                          gcgcTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGCNGtt                                                               >  2:6308434/1‑93 (MQ=255)

                          gcgcTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGtt                                                               >  2:4469515/1‑93 (MQ=255)

                          gcgcTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGtt                                                               >  2:5682117/1‑93 (MQ=255)

                             cTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTAcc                                                            <  2:338193/93‑1 (MQ=255)

                              tGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCt                                                           >  1:2438407/1‑93 (MQ=255)

                              tGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCt                                                           <  2:5099452/93‑1 (MQ=255)

                               gggCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCtt                                                          <  1:1587528/93‑1 (MQ=255)

                                ggCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTc                                                         <  2:77399/93‑1 (MQ=255)

                                ggCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTc                                                         >  2:318964/1‑93 (MQ=255)

                                 gCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCg                                                        <  2:1543285/93‑1 (MQ=255)

                                 gCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCg                                                        >  2:134905/1‑93 (MQ=255)

                                 gCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTAGTGCGGTCGCTGACGTTACATTCg                                                        >  1:6751115/1‑93 (MQ=255)

                                  ggTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCgt                                                       <  1:5985307/92‑1 (MQ=255)

                                   gTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCgtg                                                      <  1:6269914/93‑1 (MQ=255)

                                    tcttctTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCgtgg                                                     >  1:5531739/1‑93 (MQ=255)

                                    tcttctTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCgtgg                                                     >  1:5758102/1‑93 (MQ=255)

                                    tcttctTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCgtgg                                                     >  2:5813617/1‑93 (MQ=255)

                                      ttcttcGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGACGTCGTGCGGTCGTTGACGTTACCTTCgtggtg                                                   >  1:1003547/1‑93 (MQ=255)

                                       tcttcGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGgtgt                                                  <  1:2851596/93‑1 (MQ=255)

                                       tcttcGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGgtgt                                                  >  2:3267240/1‑93 (MQ=255)

                                         ttcGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTcc                                                >  2:4141168/1‑93 (MQ=255)

                                         ttcGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTcc                                                >  1:3229241/1‑93 (MQ=255)

                                          tcGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCa                                               <  1:3656191/93‑1 (MQ=255)

                                           cGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCAc                                              <  1:6384130/93‑1 (MQ=255)

                                           cGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCAc                                              <  2:492137/93‑1 (MQ=255)

                                           cGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCAc                                              >  1:5633433/1‑93 (MQ=255)

                                              tGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGa                                           <  1:7686032/93‑1 (MQ=255)

                                              tGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGa                                           >  1:2006850/1‑93 (MQ=255)

                                              tGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGa                                           >  1:261163/1‑93 (MQ=255)

                                              tGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGa                                           >  2:2423070/1‑93 (MQ=255)

                                               ggAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGaa                                          <  1:1708673/93‑1 (MQ=255)

                                               ggAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGaa                                          <  1:2913838/93‑1 (MQ=255)

                                                gAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAg                                         >  2:3384611/1‑93 (MQ=255)

                                                  gCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGt                                       >  1:4481251/1‑93 (MQ=255)

                                                   ccGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTc                                      <  2:2251574/93‑1 (MQ=255)

                                                      aaaaTGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCa                                   >  2:2683847/1‑93 (MQ=255)

                                                      aaaaTGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCa                                   >  1:2244461/1‑93 (MQ=255)

                                                      aaaaTGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCa                                   <  2:3799577/93‑1 (MQ=255)

                                                        aaTGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAAc                                 <  1:7753974/93‑1 (MQ=255)

                                                        aaTGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAAc                                 >  1:419560/1‑93 (MQ=255)

                                                         aTGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACt                                >  2:622269/1‑93 (MQ=255)

                                                         aTGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACt                                >  2:1736521/1‑93 (MQ=255)

                                                          tgttgtCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACtt                               >  1:2923669/1‑93 (MQ=255)

                                                          tgttgtCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACtt                               >  1:7313627/1‑93 (MQ=255)

                                                          tgttgtCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACtt                               <  1:386476/93‑1 (MQ=255)

                                                           gttgtCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTc                              <  1:5163032/93‑1 (MQ=255)

                                                           gttgtCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTc                              <  1:3097037/93‑1 (MQ=255)

                                                           gttgtCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTc                              <  1:1367649/93‑1 (MQ=255)

                                                           gttgtCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTc                              <  2:1709570/93‑1 (MQ=255)

                                                            ttgtCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCg                             >  1:1133357/1‑93 (MQ=255)

                                                              gtCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGAt                           <  2:7906087/93‑1 (MQ=255)

                                                                  cgcAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAAtc                       <  1:5872918/93‑1 (MQ=255)

                                                                   gcagcaACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAAtct                      >  1:6575429/1‑93 (MQ=255)

                                                                      gcaACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCg                   <  2:1606247/93‑1 (MQ=255)

                                                                      gcaACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCg                   <  2:1267874/93‑1 (MQ=255)

                                                                         aCATGCCGTTTTGCCCGTGTTGGGAGGCCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGc                >  2:2025493/1‑93 (MQ=255)

                                                                          cATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGcc               >  1:7611076/1‑93 (MQ=255)

                                                                               cGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACa          <  2:434787/93‑1 (MQ=255)

                                                                               cGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACa          >  1:343541/1‑93 (MQ=255)

                                                                                gTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACAt         >  1:4934696/1‑93 (MQ=255)

                                                                                 ttttGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACAtt        >  1:7907106/1‑93 (MQ=255)

                                                                                 ttttGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACAtt        >  1:5284956/1‑93 (MQ=255)

                                                                                  tttGCCCGTGTTGGGAGGTCGTGCGGTCGCTGAAGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTa       <  2:1040313/93‑1 (MQ=255)

                                                                                   ttccccGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTAc      <  1:5832310/90‑1 (MQ=255)

                                                                                    tGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACa     <  2:4244396/93‑1 (MQ=255)

                                                                                      cccGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATg   <  2:6596222/93‑1 (MQ=255)

                                                                                       ccGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATGa  >  2:7364942/1‑93 (MQ=255)

                                                                                       ccGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATGa  >  2:7643640/1‑93 (MQ=255)

                                                                                       ccGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATGa  >  1:2225768/1‑93 (MQ=255)

                                                                                       ccGTGTTGGGAGGACGCGCGGCCGCTGACGTTACCTTCGTGGAGTCCACGGAAGGACGCAACTTCGATAATGTCTGCGACAAACAATACACGa  >  2:7236179/1‑93 (MQ=255)

TACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATGA  >  NC_000913/158621‑158800

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