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	.	95.7%	297.4 / 11.5	117	noncoding (61/121 nt)	rrfF	5S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base G (3/2); new base . (54/58); total (57/60)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-02

Only 100 of 129 total aligned reads displayed.

CAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAA  >  NC_000913/3421027‑3421207

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

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

  gTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgtt                                                                                        <  1:2663882/93‑3 (MQ=38)

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

   tCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTNCGGCGttt                                                                                       >  1:5298609/1‑93 (MQ=255)

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

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

   tCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGGtt                                                                                       >  1:323530/1‑93 (MQ=255)

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

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

     ttCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGAttca                                                                                     >  1:3445869/1‑91 (MQ=255)

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

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

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

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

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

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

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

         gACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTc                                                                                >  2:3723162/1‑93 (MQ=16)

         gACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTc                                                                                >  2:1262761/1‑93 (MQ=9)

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

         gACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTc                                                                                >  1:4409309/1‑93 (MQ=16)

         gACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTc                                                                                >  2:6454468/1‑93 (MQ=16)

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

          aCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCt                                                                               >  2:5704654/1‑93 (MQ=9)

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

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

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

             gAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAg                                                                            >  1:4788449/1‑93 (MQ=12)

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

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

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

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

               gCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGtt                                                                          <  1:2226496/93‑1 (MQ=9)

                ccTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTc                                                                         <  2:2897623/93‑1 (MQ=9)

                ccTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTc                                                                         <  2:2749007/93‑1 (MQ=9)

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

                ccTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTc                                                                         <  1:1780018/93‑1 (MQ=9)

                ccTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTc                                                                         >  2:2903486/1‑93 (MQ=9)

                 cTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGG‑TTCACTTCTGAGTTCg                                                                        >  1:1534205/1‑93 (MQ=255)

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

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

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

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

                 cTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GATTGACTTCTGAGTTCg                                                                        >  1:1555122/1‑93 (MQ=11)

                  ttttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTTgg                                                                       >  1:1663997/1‑93 (MQ=9)

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

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

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

                  ttttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCgg                                                                       <  1:3062279/93‑1 (MQ=9)

                  tcttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTCAGTTcag                                                                       >  1:5274983/3‑91 (MQ=1)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                       cacTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCAGGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTa                  >  2:4681439/1‑93 (MQ=255)

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

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

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

                                                                          tACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATc               >  2:5291716/1‑93 (MQ=255)

                                                                          tACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATc               >  2:541686/1‑93 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

                                                                               tCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCg          >  2:710093/1‑93 (MQ=255)

                                                                               tCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCg          >  2:4664964/1‑93 (MQ=255)

                                                                               tCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCg          <  1:4537108/93‑1 (MQ=255)

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

                                                                                  gcgcTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCCCCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAAc       <  1:4056680/93‑1 (MQ=255)

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

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

                                                                                  gcgcTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAAc       >  1:136795/1‑93 (MQ=255)

                                                                                   cgcTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGactt      >  2:2106130/1‑90 (MQ=255)

                                                                                   cgcTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACt      <  1:3985051/93‑1 (MQ=255)

                                                                                   cgcTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACt      >  2:1821942/1‑93 (MQ=255)

                                                                                   cgcTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACt      <  1:5448153/93‑1 (MQ=255)

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

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

                                                                                     cTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACttt    >  2:847706/1‑93 (MQ=255)

                                                                                      tACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTa   >  1:1230195/1‑93 (MQ=255)

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

CAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAA  >  NC_000913/3421027‑3421207

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