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	98.3%	620.0 / ‑2.4	178	R128H (CGC→CAC)	pcnB	poly(A) polymerase I
Reads supporting (aligned to +/- strand): ref base C (2/1); new base T (98/77); total (100/78)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

Only 100 of 190 total aligned reads displayed.

TACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATGAG  >  NC_000913/158621‑158801

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

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

  cAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTg                                                                                        >  1:2902912/1‑93 (MQ=255)

  cAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTg                                                                                        >  2:3989499/1‑93 (MQ=255)

  cAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTg                                                                                        <  2:1820070/93‑1 (MQ=255)

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

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

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

     gCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGa                                                                                     <  2:1647483/93‑1 (MQ=255)

     gCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGa                                                                                     <  1:1620490/93‑1 (MQ=255)

      gtGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAg                                                                                    <  1:3775514/92‑1 (MQ=255)

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

        gTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGt                                                                                  >  2:3147885/1‑93 (MQ=255)

        gTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGt                                                                                  <  1:3042142/93‑1 (MQ=255)

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

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

            aTAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTg                                                                              >  2:1351608/1‑93 (MQ=255)

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

              aGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTCTGCCCGTGTTGGGAGGTCGGGCg                                                                            >  2:1901431/1‑93 (MQ=255)

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

                   aaTCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGc                                                                       <  2:3205253/93‑1 (MQ=255)

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

                       gcggcgCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTTGTGCAGCAACATGCCGCTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGAc                                                                   >  2:1520891/1‑93 (MQ=255)

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

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

                          tcgcTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGtt                                                                <  2:507347/92‑1 (MQ=255)

                            gcTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTAc                                                              <  1:2040128/93‑1 (MQ=255)

                              tGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCACGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCt                                                            <  2:2521328/93‑1 (MQ=255)

                               gggCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCtt                                                           >  1:3393188/1‑93 (MQ=255)

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

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

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

                               gggCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCtt                                                           >  1:2081281/1‑93 (MQ=255)

                                ggCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTc                                                          >  2:4500866/1‑93 (MQ=255)

                                   gTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCgtg                                                       >  2:440120/1‑93 (MQ=255)

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

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

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

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

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

                                      ttcttcGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCgtggtg                                                    >  2:446577/1‑93 (MQ=255)

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

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

                                        cttcGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTc                                                  <  1:538745/93‑1 (MQ=255)

                                         ttcGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTAACTTCGTGGTGTcc                                                 >  1:4469761/1‑93 (MQ=255)

                                           tgATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCAc                                               <  1:3816443/92‑1 (MQ=255)

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

                                            gATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCCCg                                              >  1:1262321/1‑93 (MQ=255)

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

                                               ggAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGaa                                           >  1:3894711/1‑93 (MQ=255)

                                                  gCCGAAAATGTTGTCGCGCAGCAACATGTCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTCCCTTCGTGGTGTCCACGGAAGGt                                        <  2:2741790/93‑1 (MQ=255)

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

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

                                                  gCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGt                                        <  2:4421068/93‑1 (MQ=255)

                                                  gCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGt                                        <  1:4082275/93‑1 (MQ=255)

                                                  gCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGCAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGt                                        <  2:2810403/93‑1 (MQ=255)

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

                                                    cGAAATTGTTGTCGCGCAGTAACATGCCGGTTCGCCCGTGTTGGGAGGTCGTGCGGTCGTTGACGTTACCTTCGTGGTGTCCACGGAAGGTCg                                      >  1:4199607/1‑93 (MQ=255)

                                                     gAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGc                                     >  1:457801/1‑93 (MQ=255)

                                                     gAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGc                                     >  2:1915988/1‑93 (MQ=255)

                                                     gAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGACGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGc                                     >  2:846828/1‑93 (MQ=255)

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

                                                      aaaaTGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGCGTCCACGGAAGGTCGCa                                    >  2:1699276/1‑93 (MQ=255)

                                                       aaaTGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCaa                                   >  2:1047633/1‑93 (MQ=255)

                                                       aaaTGTTGTCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCaa                                   >  1:3162445/1‑93 (MQ=255)

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

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

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

                                                           gttgtCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGCCCACGGAAGGTCGCAACttg                               >  1:355654/1‑92 (MQ=255)

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

                                                               tCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATa                           <  1:2589275/93‑1 (MQ=255)

                                                               tCGCGCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATa                           <  2:1277646/93‑1 (MQ=255)

                                                                cgcgCAGCAACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATaa                          <  1:1072103/93‑1 (MQ=255)

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

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

                                                                   gcagcaACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAAtct                       <  1:2248765/93‑1 (MQ=255)

                                                                     agcaACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTcc                     <  1:1097746/93‑1 (MQ=255)

                                                                     agcaACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTcc                     <  2:3339169/93‑1 (MQ=255)

                                                                     agcaACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTcc                     >  1:1969356/1‑93 (MQ=255)

                                                                     agcaACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTcc                     >  2:2383831/1‑93 (MQ=255)

                                                                     agcaACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTcc                     <  1:4329173/93‑1 (MQ=255)

                                                                     agcaACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTcc                     <  1:3317459/93‑1 (MQ=255)

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

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

                                                                       caACATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCgg                   <  2:3275555/93‑1 (MQ=255)

                                                                         aCATGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGCCGCAACTTCGATAATCTCCGGGc                 >  1:2315482/1‑93 (MQ=255)

                                                                           aTGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCa               >  2:2307961/1‑93 (MQ=255)

                                                                           aTGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCa               >  2:4627350/1‑93 (MQ=255)

                                                                            tGCCGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCaa              >  1:1658064/1‑93 (MQ=255)

                                                                              ccGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAAc            >  2:3050545/1‑93 (MQ=255)

                                                                              ccGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAAc            >  1:2643426/1‑93 (MQ=255)

                                                                              ccGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAAc            <  1:403348/93‑1 (MQ=255)

                                                                              ccGTTTTGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAAc            <  2:3614001/93‑1 (MQ=255)

                                                                                 ttttGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACAtt         <  2:1103833/93‑1 (MQ=255)

                                                                                 ttttGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACAtt         <  1:2373851/93‑1 (MQ=255)

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

                                                                                 ttttGCCCGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACAtt         <  2:724589/93‑1 (MQ=255)

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

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

                                                                                      cccGTGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGCAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATg    >  2:4443103/1‑93 (MQ=255)

TACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCCGCGTTGGGAGGTCGTGCGGTCGCTGACGTTACCTTCGTGGTGTCCACGGAAGGTCGCAACTTCGATAATCTCCGGGCCAAACATTACATGAG  >  NC_000913/158621‑158801

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