BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,421,120	Δ1 bp	noncoding (60/121 nt)	rrfF ←	5S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,421,119	0	G	.	94.8%	276.6 / 14.3	115	noncoding (61/121 nt)	rrfF	5S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base G (3/3); new base . (49/60); total (52/63)
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.36e-02

Only 100 of 129 total aligned reads displayed.

CAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAA  >  NC_000913/3421027‑3421207

cAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcg                                                                                          >  1:703768/1‑93 (MQ=255)

cAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcg                                                                                          >  1:3039237/1‑93 (MQ=255)

 aGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgt                                                                                         <  1:4694263/93‑2 (MQ=39)

 aGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgt                                                                                         >  2:5774927/1‑92 (MQ=255)

 aGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgt                                                                                         >  1:2611368/1‑92 (MQ=39)

  gTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgtt                                                                                        >  2:1563974/1‑91 (MQ=255)

  gTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgtt                                                                                        >  2:1402407/1‑91 (MQ=255)

  gTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgtt                                                                                        >  2:1789481/1‑91 (MQ=38)

  gTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGAGCTAcggcgtt                                                                                        >  1:1022242/1‑91 (MQ=255)

   tCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGttt                                                                                       >  1:2113350/1‑93 (MQ=255)

   tCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGttt                                                                                       <  1:5447988/93‑1 (MQ=37)

    cTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGtttc                                                                                      <  2:5339668/93‑2 (MQ=25)

    cTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGtttc                                                                                      <  1:690018/93‑2 (MQ=25)

     ttCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCa                                                                                    >  1:255358/1‑93 (MQ=255)

     ttCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCa                                                                                    >  1:1284828/1‑93 (MQ=255)

      tCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCAc                                                                                   >  2:3208144/1‑93 (MQ=255)

       ccGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTGCTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACt                                                                                  <  1:5935607/93‑1 (MQ=255)

       ccGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACt                                                                                  <  1:6076548/93‑1 (MQ=12)

        cGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACtt                                                                                 <  1:3583228/93‑1 (MQ=9)

        cGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACtt                                                                                 >  1:3098987/1‑93 (MQ=12)

        cGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACtt                                                                                 >  1:2676474/1‑93 (MQ=16)

         gACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTc                                                                                <  2:5317095/93‑1 (MQ=9)

         gACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTc                                                                                <  1:2155447/93‑1 (MQ=9)

         gACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTc                                                                                >  2:1700805/1‑93 (MQ=12)

          aCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCt                                                                               >  2:1106583/1‑93 (MQ=16)

          aCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCt                                                                               >  1:4754346/1‑93 (MQ=12)

          aCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCt                                                                               >  2:744479/1‑93 (MQ=16)

          aCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCt                                                                               <  2:835130/93‑1 (MQ=9)

           cTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTCCTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTg                                                                              <  2:5455358/93‑1 (MQ=9)

           cTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCGTGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTg                                                                              >  2:4926531/1‑93 (MQ=255)

           cTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTg                                                                              >  2:3736350/1‑93 (MQ=16)

           cTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTg                                                                              >  1:5151176/1‑93 (MQ=9)

           cTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTg                                                                              >  1:4887814/1‑93 (MQ=16)

           cTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTg                                                                              >  2:4266416/1‑93 (MQ=16)

            tGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTgg                                                                             >  1:125350/1‑92 (MQ=16)

            tGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGa                                                                             >  1:771839/1‑93 (MQ=12)

             gAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCGCACTACCATCGGCGCTACGGC‑GTTTCACTTTTGAg                                                                            >  2:5153433/1‑93 (MQ=11)

             gAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAg                                                                            >  1:1817180/1‑93 (MQ=9)

             gAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAg                                                                            <  1:2524135/93‑1 (MQ=9)

             gAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAg                                                                            >  1:291461/1‑93 (MQ=16)

             gAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAg                                                                            >  2:3921995/1‑93 (MQ=16)

             gAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAg                                                                            <  2:4273798/93‑1 (MQ=9)

              aGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGt                                                                           >  2:4045441/1‑93 (MQ=16)

              aGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGt                                                                           <  2:1547098/93‑1 (MQ=9)

               gCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGtt                                                                          >  2:1092349/1‑93 (MQ=16)

               gCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGtt                                                                          >  1:2882613/1‑93 (MQ=16)

               gCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGtt                                                                          >  1:5055150/1‑93 (MQ=16)

               gCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGtt                                                                          >  1:3294038/1‑93 (MQ=16)

                ccTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTc                                                                         >  1:5375275/1‑93 (MQ=16)

                 cTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTGCGGC‑GTTTCACTTCTGAGTTCg                                                                        <  2:2869136/93‑1 (MQ=9)

                 cTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCg                                                                        >  2:5074220/1‑93 (MQ=12)

                 cTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCg                                                                        >  2:3833883/1‑93 (MQ=9)

                 cTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCg                                                                        >  1:5699108/1‑93 (MQ=9)

                 cTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCg                                                                        <  2:2366149/93‑1 (MQ=9)

                 cTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACGCTCGCATGGGGAGCCCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCg                                                                        <  2:5685310/93‑1 (MQ=11)

                  ttttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCgg                                                                       >  1:4858032/1‑93 (MQ=9)

                  ttttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCgg                                                                       >  1:3055950/1‑93 (MQ=9)

                  ttttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCgg                                                                       >  1:1970988/1‑93 (MQ=12)

                   tttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGc                                                                      <  2:2885359/93‑1 (MQ=9)

                   tttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGc                                                                      >  1:1725799/1‑93 (MQ=12)

                    ttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCa                                                                     <  1:6276541/93‑1 (MQ=0)

                    ttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCa                                                                     <  2:5600796/93‑1 (MQ=0)

                    ttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCa                                                                     >  1:1154298/1‑93 (MQ=0)

                    ttGTTTTATTTGATGTCTGCCAGTCCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTCCGGC‑GTTTCACTTCTGAGTTCGGCa                                                                     <  2:3058236/93‑1 (MQ=1)

                                                                   cccACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCttt                      <  2:4964186/93‑1 (MQ=9)

                                                                   cccACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCttt                      <  2:1794561/93‑1 (MQ=9)

                                                                   cccACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCttt                      >  1:2702378/1‑93 (MQ=0)

                                                                    ccACACTACCATCGGCGCTATGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCtttt                     <  1:4203400/93‑1 (MQ=255)

                                                                     cacaCTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTc                    <  2:1059927/93‑1 (MQ=12)

                                                                     cacaCTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTc                    <  1:962756/93‑1 (MQ=12)

                                                                     cacaCTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTc                    <  2:3359044/93‑1 (MQ=12)

                                                                      acacTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCt                   <  2:1909311/93‑1 (MQ=19)

                                                                       cacTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTa                  <  2:119737/93‑1 (MQ=19)

                                                                       cacTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTa                  >  1:2383779/1‑93 (MQ=19)

                                                                        acTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTaa                 <  2:1554006/93‑1 (MQ=255)

                                                                        acTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTaa                 <  2:3160725/93‑1 (MQ=255)

                                                                        acTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTaa                 <  2:2942558/93‑1 (MQ=255)

                                                                         cTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAAt                <  1:5435902/93‑1 (MQ=255)

                                                                         cTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAAt                >  1:4373139/1‑93 (MQ=255)

                                                                          tACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATc               <  1:5944513/93‑1 (MQ=255)

                                                                          tACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATc               <  1:3571448/93‑1 (MQ=255)

                                                                           aCCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCt              <  2:2352940/93‑1 (MQ=255)

                                                                            ccATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTg             <  1:1608026/93‑1 (MQ=255)

                                                                            ccATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTg             >  2:2244692/1‑93 (MQ=255)

                                                                             cATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGc            <  2:1213051/93‑1 (MQ=255)

                                                                             cATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGc            <  2:414512/93‑1 (MQ=255)

                                                                             cATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGc            >  1:5480120/1‑93 (MQ=255)

                                                                              aTCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGcc           >  2:470324/1‑93 (MQ=255)

                                                                              aTCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGcc           <  1:64425/93‑1 (MQ=255)

                                                                                 ggCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGaa        <  1:280420/93‑1 (MQ=255)

                                                                                  gcgcTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAAc       <  1:812163/93‑1 (MQ=255)

                                                                                    gcTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACtt     <  2:6069915/93‑1 (MQ=255)

                                                                                    gcTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACtt     >  1:274350/1‑93 (MQ=255)

                                                                                    gcTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACtt     <  1:2162413/93‑1 (MQ=255)

                                                                                      tACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACtttt   <  2:6049238/93‑2 (MQ=255)

                                                                                      tACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTa   <  1:4104696/93‑1 (MQ=255)

                                                                                      tACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTa   <  2:598632/93‑1 (MQ=255)

                                                                                      tACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTa   <  2:911544/93‑1 (MQ=255)

                                                                                       aCGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTaa  >  1:82623/1‑93 (MQ=255)

                                                                                       aCGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTaa  >  1:492508/1‑93 (MQ=255)

CAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAA  >  NC_000913/3421027‑3421207

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