BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,758,859	A→G	R6R (CGT→CGC)	selA ←	selenocysteine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,758,859	0	A	G	99.4%	1084.6 / ‑2.7	313	R6R (CGT→CGC)	selA	selenocysteine synthase
Reads supporting (aligned to +/- strand): ref base A (0/2); new base G (172/139); total (172/141)
Fisher's exact test for biased strand distribution p-value = 2.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

Only 100 of 321 total aligned reads displayed.

CACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAG  >  NC_000913/3758769‑3758948

 acGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAgcgc                                                                                        <  2:5272629/93‑1 (MQ=255)

 acGCGGGTGTGACAATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAgcgc                                                                                        <  1:285650/93‑1 (MQ=255)

  cgcgGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGcgcg                                                                                       <  1:6109955/93‑1 (MQ=255)

  cgcgGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGcgcg                                                                                       >  2:133350/1‑93 (MQ=255)

   gcgGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCCCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGt                                                                                      <  2:4561467/93‑1 (MQ=255)

    cgGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGtt                                                                                     <  2:2500731/93‑1 (MQ=255)

    cgGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGtt                                                                                     <  1:636337/93‑1 (MQ=255)

    cgGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGtt                                                                                     <  1:3003980/93‑1 (MQ=255)

     gggTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGttt                                                                                    >  2:3647564/1‑93 (MQ=255)

       gtgtGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCg                                                                                  >  1:2296197/1‑93 (MQ=255)

       gtgtGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCg                                                                                  <  2:2182564/93‑1 (MQ=255)

         gtgACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGt                                                                                <  1:4766551/93‑1 (MQ=255)

         gtgACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGt                                                                                <  1:4486848/93‑1 (MQ=255)

          tgACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGtt                                                                               <  2:4457566/93‑1 (MQ=255)

          tgACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGtt                                                                               <  2:1498948/93‑1 (MQ=255)

            aCCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGt                                                                             >  2:3390222/1‑93 (MQ=255)

             ccATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTc                                                                            >  1:5826613/1‑93 (MQ=255)

             ccATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTc                                                                            >  2:1460980/1‑93 (MQ=255)

              cATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCa                                                                           <  2:5755343/93‑1 (MQ=255)

              cATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCa                                                                           >  1:6042379/1‑93 (MQ=255)

               aTAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCAt                                                                          >  2:2094893/1‑93 (MQ=255)

               aTAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCAt                                                                          <  1:5075885/93‑1 (MQ=255)

               aTAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCAt                                                                          <  2:4717540/93‑1 (MQ=255)

                  aGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAgg                                                                       <  2:4988287/93‑1 (MQ=255)

                   gTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGa                                                                      <  1:3089883/93‑1 (MQ=255)

                      tCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAAc                                                                   <  1:3204916/93‑1 (MQ=255)

                      tCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAAc                                                                   <  2:3340954/93‑1 (MQ=255)

                        aCGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCt                                                                 <  2:5162041/93‑1 (MQ=255)

                           cAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGc                                                              <  2:1051716/93‑1 (MQ=255)

                            aaagaaagGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCt                                                             <  1:4405123/93‑1 (MQ=255)

                            aaagaaagGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCt                                                             <  2:3886902/93‑1 (MQ=255)

                             aagaaagGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTg                                                            >  2:4658940/1‑93 (MQ=255)

                             aagaaagGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTg                                                            >  1:3265770/1‑93 (MQ=255)

                             aagaaagGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTg                                                            >  2:2600638/1‑93 (MQ=255)

                             aagaaagGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTg                                                            >  1:2471020/1‑93 (MQ=255)

                              agaaagGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGa                                                           >  1:4679401/1‑93 (MQ=255)

                              agaaagGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGa                                                           >  1:5988130/1‑93 (MQ=255)

                              agaaagGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGa                                                           <  1:353139/93‑1 (MQ=255)

                               gaaagGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGAt                                                          >  2:295649/1‑93 (MQ=255)

                                aaagGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATc                                                         >  2:1571801/1‑93 (MQ=255)

                                  aggaaggaGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCaa                                                       >  1:3847435/1‑93 (MQ=255)

                                  aggaaggaGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCaa                                                       >  1:1584414/1‑93 (MQ=255)

                                   ggaaggaGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAg                                                      <  2:485812/93‑1 (MQ=255)

                                   ggaaggaGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAg                                                      <  2:5420415/93‑1 (MQ=255)

                                   ggaaggaGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAg                                                      <  2:1840620/93‑1 (MQ=255)

                                    gaaggaGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAgg                                                     <  2:5086762/93‑1 (MQ=255)

                                      aggaGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGcc                                                   >  1:4333557/1‑93 (MQ=255)

                                       ggaGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGccc                                                  <  2:3799962/93‑1 (MQ=255)

                                       ggaGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGccc                                                  >  1:4916746/1‑93 (MQ=255)

                                       ggaGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGccc                                                  <  2:6134235/93‑1 (MQ=255)

                                       ggaGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGGTGTGATAGGAAACCTGGCTGATCAAgtacc                                                  >  1:136257/1‑89 (MQ=255)

                                         aGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCctc                                                >  1:5195013/1‑93 (MQ=255)

                                         aGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCctc                                                <  2:6417124/93‑1 (MQ=255)

                                          gCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCctct                                               <  2:6266624/93‑1 (MQ=255)

                                           cTATCACGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCtctc                                              <  1:6127046/93‑1 (MQ=255)

                                            tATCTCTCAATAAGCGATCAATAGCCGGAATTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTca                                             <  2:31388/93‑1 (MQ=255)

                                             aTCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTcac                                            >  1:4171157/1‑93 (MQ=255)

                                                gcgcAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACg                                         <  1:2372064/93‑1 (MQ=255)

                                                  gcAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGga                                       >  2:813309/1‑93 (MQ=255)

                                                    aaTAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGgaga                                     <  1:3715536/93‑1 (MQ=255)

                                                    aaTAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGgaga                                     >  2:5899413/1‑93 (MQ=255)

                                                    aaTAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGgaga                                     <  1:1400423/93‑1 (MQ=255)

                                                     aTAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGaa                                    >  1:5312777/1‑93 (MQ=255)

                                                     aTAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGaa                                    >  1:2800184/1‑93 (MQ=255)

                                                      tAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGCAg                                   >  1:1217314/1‑93 (MQ=255)

                                                      tAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAg                                   <  2:3695960/93‑1 (MQ=255)

                                                      tAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAg                                   <  1:5593702/93‑1 (MQ=255)

                                                       aaGCGATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGTGAAgg                                  <  1:2029916/93‑1 (MQ=255)

                                                           gATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGt                              >  1:5264648/1‑93 (MQ=255)

                                                           gATCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGt                              >  2:5544475/1‑93 (MQ=255)

                                                            aTCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGtt                             >  2:3452661/1‑93 (MQ=255)

                                                            aTCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGtt                             <  1:3146214/93‑1 (MQ=255)

                                                            aTCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGtt                             >  2:6091632/1‑93 (MQ=255)

                                                             tCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGttt                            >  2:243054/1‑93 (MQ=255)

                                                             tCAATAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCACTCACACGGAGAAGGGCGttt                            <  1:1502569/93‑1 (MQ=255)

                                                               aaTAGCCGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGACCTCTCACACGGAGAAGGGCGTTTaa                          >  1:6419839/1‑93 (MQ=255)

                                                                    ccGGAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCCCACGGAGAAGGGCGTTTAACATaa                     >  2:4166126/1‑93 (MQ=255)

                                                                      ggAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAAcc                   >  1:4798122/1‑93 (MQ=255)

                                                                      ggAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAAcc                   <  1:529957/93‑1 (MQ=255)

                                                                      ggAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGCCCCTCTCACACGGAGAAGGGCGTTTAACATAAcc                   <  1:3333469/93‑1 (MQ=255)

                                                                       gAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCa                  >  2:4070188/1‑93 (MQ=255)

                                                                       gAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCa                  >  2:3563312/1‑93 (MQ=255)

                                                                       gAAGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCa                  <  2:1227580/93‑1 (MQ=255)

                                                                        aaGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCAc                 <  2:2354978/93‑1 (MQ=255)

                                                                        aaGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCAc                 <  2:264374/93‑1 (MQ=255)

                                                                         aGTTGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACg                >  2:4342636/1‑93 (MQ=255)

                                                                           ttGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGa              <  2:4603145/93‑1 (MQ=255)

                                                                           ttGACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGa              <  1:82746/93‑1 (MQ=255)

                                                                             gACTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGAtt            >  2:1605490/1‑93 (MQ=255)

                                                                              aCTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTg           >  1:5371009/1‑93 (MQ=255)

                                                                              aCTATAGAGGGAGCGCGCTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTg           >  2:5186353/1‑93 (MQ=255)

                                                                               cTATAGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGt          <  1:4985635/93‑1 (MQ=255)

                                                                                tataGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTa         >  1:5442141/1‑93 (MQ=255)

                                                                                 ataGAGGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTaa        >  1:3642586/1‑93 (MQ=255)

                                                                                     agGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTg    <  1:2599318/93‑1 (MQ=255)

                                                                                     agGGAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTg    >  1:174373/1‑93 (MQ=255)

                                                                                      gggAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTga   >  1:6222740/1‑93 (MQ=255)

                                                                                      gggAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTga   >  1:2709375/1‑93 (MQ=255)

                                                                                      gggAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTga   <  1:1041704/93‑1 (MQ=255)

                                                                                       ggAGCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTgag  >  1:925588/1‑93 (MQ=255)

CACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAG  >  NC_000913/3758769‑3758948

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