BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	931,267	T→C	T3A (ACG→GCG)	trxB ←	thioredoxin reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	931,267	0	T	C	92.5%	794.9 / 42.7	296	T3A (ACG→GCG)	trxB	thioredoxin reductase
Reads supporting (aligned to +/- strand): ref base T (10/12); new base C (160/113); total (171/125)
Fisher's exact test for biased strand distribution p-value = 2.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01

Only 100 of 311 total aligned reads displayed.

CAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGTGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAAGAGTATGACGATTTTGTTA  >  NC_000913/931177‑931356

  ggttgcaggttgGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGccc                                                                                       >  2:2493835/1‑93 (MQ=255)

    ttgcaggttgGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCAt                                                                                     >  2:337684/1‑93 (MQ=255)

       caggttgGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATgag                                                                                  <  2:274865/93‑1 (MQ=255)

         ggttgGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGAt                                                                                >  2:6160443/1‑93 (MQ=255)

         ggttgGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGAt                                                                                >  2:8784216/1‑93 (MQ=255)

          gttgGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATc                                                                               >  2:6267428/1‑93 (MQ=255)

           ttgGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATcc                                                                              <  2:7429271/93‑1 (MQ=255)

             ggCGCCCGCCGCGTAGACAGCAGCGGTGCATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATccgc                                                                            >  1:5308300/1‑91 (MQ=255)

              gcgcgcgcCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGTGCCCATGAGATCCCCa                                                                           <  1:59681/93‑1 (MQ=255)

                 cgcgcCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAg                                                                        <  1:6753711/93‑1 (MQ=255)

                 cgcgcCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAg                                                                        >  1:9408834/1‑93 (MQ=255)

                 cgcgcCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAg                                                                        <  1:3512930/93‑1 (MQ=255)

                   cgccgcGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAgtt                                                                      >  2:4093826/1‑93 (MQ=255)

                   cgccgcGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAgtt                                                                      >  1:7714954/1‑93 (MQ=255)

                    gccgcGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAgttg                                                                     >  2:9457424/1‑93 (MQ=255)

                     ccgcGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAgttgt                                                                    <  2:7713851/93‑1 (MQ=255)

                       gcgTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCCTAGttgttg                                                                  >  1:2473871/1‑93 (MQ=255)

                        cgTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTgg                                                                 >  1:1998285/1‑93 (MQ=255)

                         gTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGc                                                                >  2:575034/1‑93 (MQ=255)

                          tAGACAGCAGCGGTGTTTCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCa                                                               >  2:8952694/1‑93 (MQ=255)

                          tAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCa                                                               <  2:730255/93‑1 (MQ=255)

                             aCAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGAc                                                            >  2:7656424/1‑93 (MQ=255)

                               agcagcGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACaa                                                          <  2:5453074/93‑1 (MQ=255)

                               agcagcGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACaa                                                          <  2:6894803/93‑1 (MQ=255)

                                gcagcGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAAt                                                         >  1:3269587/1‑93 (MQ=255)

                                gcagcGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAAt                                                         <  2:2988331/93‑1 (MQ=255)

                                 cagcGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATg                                                        <  1:7148231/93‑1 (MQ=255)

                                    cGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGc                                                     >  2:9183966/1‑93 (MQ=255)

                                      gtgtATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAg                                                   >  1:9478403/1‑93 (MQ=255)

                                        gtATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGa                                                 >  1:5526634/1‑93 (MQ=255)

                                        gtATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGa                                                 >  1:9850102/1‑93 (MQ=255)

                                         tATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGAt                                                >  2:4673570/1‑93 (MQ=255)

                                         tATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGAt                                                <  1:6726141/93‑1 (MQ=255)

                                           tCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTg                                              >  1:6951097/1‑93 (MQ=255)

                                             ccgccgGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTa                                            >  1:3938510/1‑93 (MQ=255)

                                              cgccgGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAg                                           <  1:2275473/93‑1 (MQ=255)

                                               gccgTGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAgg                                          <  2:1254046/93‑1 (MQ=255)

                                                 cgGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGTGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGa                                        <  2:6002224/93‑1 (MQ=255)

                                                  gggCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGaa                                       >  1:8738590/1‑93 (MQ=255)

                                                  gggCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGaa                                       >  1:8063056/1‑93 (MQ=255)

                                                   ggCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAAt                                      >  1:4385274/1‑93 (MQ=255)

                                                    gCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAAtt                                     >  1:6471459/1‑93 (MQ=255)

                                                       tGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTAc                                  >  2:9647214/1‑93 (MQ=255)

                                                       tGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTAc                                  >  1:8443992/1‑93 (MQ=255)

                                                        gAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGGATTTACa                                 >  2:3875094/1‑93 (MQ=255)

                                                        gAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACa                                 <  1:3123280/93‑1 (MQ=255)

                                                         aaCCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGGTTGTAGGGAATTTACAg                                >  1:9047950/1‑93 (MQ=255)

                                                         aaCCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAg                                <  1:4312341/93‑1 (MQ=255)

                                                         aaCCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAg                                <  2:581807/93‑1 (MQ=255)

                                                         aaCCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAg                                <  2:6560768/93‑1 (MQ=255)

                                                         aaCCCAGGATAAGCAGTTTACTGTGTTTGATCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAg                                >  1:1042138/1‑93 (MQ=255)

                                                          aCCCAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGa                               <  1:3051248/93‑1 (MQ=255)

                                                           cccAGGATAAGCAGTTTACTGTGTTTGGTCGTGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGAc                              >  1:1380458/1‑93 (MQ=255)

                                                           cccAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGAc                              >  2:9874018/1‑93 (MQ=255)

                                                           cccAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGAc                              >  1:2660719/1‑93 (MQ=255)

                                                           cccAGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGAc                              >  1:4339286/1‑93 (MQ=255)

                                                              aGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTa                           <  2:414761/93‑1 (MQ=255)

                                                              aGGATAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTa                           >  1:7908085/1‑93 (MQ=255)

                                                                  tAAGCAGTTTACTGTGTTTGGTCGTGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTaaaaa                       >  2:8552960/1‑93 (MQ=255)

                                                                  tAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTaaaaa                       >  2:3673043/1‑93 (MQ=255)

                                                                  tAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTaaaaa                       >  1:9596033/1‑93 (MQ=255)

                                                                  tAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTaaaaa                       >  1:1244102/1‑93 (MQ=255)

                                                                  tAAGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTaaaaa                       <  2:2002199/93‑1 (MQ=255)

                                                                   aaGCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTaaaaaa                      <  2:4481511/93‑1 (MQ=255)

                                                                    ttttgTTTATTGTGTGTGGTCGCGCCCATGAGATCCCAATAGTTGTGGGAAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTaaaaaaa                     <  2:9541254/89‑1 (MQ=255)

                                                                     gCAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAag                    <  2:8182658/93‑1 (MQ=255)

                                                                      cAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAaga                   >  1:2791682/1‑93 (MQ=255)

                                                                      cAGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAaga                   <  2:5185781/93‑1 (MQ=255)

                                                                       aGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAagag                  >  1:468062/1‑93 (MQ=255)

                                                                       aGTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAagag                  >  2:8089835/1‑93 (MQ=255)

                                                                        gTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAAGAGt                 >  2:2890572/1‑93 (MQ=255)

                                                                        gTTTACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAAGAGt                 <  1:6392370/93‑1 (MQ=255)

                                                                         tttACTGTGTTTGGTCGTGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAAGAGTa                >  1:1890577/1‑93 (MQ=255)

                                                                         tttACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAAGAGTa                <  1:8256565/93‑1 (MQ=255)

                                                                         tttACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAAGAGTa                >  2:5624080/1‑93 (MQ=255)

                                                                         tttACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAAGAGTa                >  1:1099481/1‑93 (MQ=255)

                                                                         tttACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGAAGTAAAAAAAGAGTa                >  2:2260112/1‑93 (MQ=255)

                                                                          ttACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAAGAGTAt               <  1:2704657/93‑1 (MQ=255)

                                                                          ttACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAAGAGTAt               <  1:7559398/93‑1 (MQ=255)

                                                                          ttACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAAGAGTAt               <  1:619127/93‑1 (MQ=255)

                                                                          ttACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAAGAGTAt               >  1:4906505/1‑93 (MQ=255)

                                                                           tACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAAGAGTATg              <  2:2374238/93‑1 (MQ=255)

                                                                           tACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAAGAGTATg              <  1:478403/93‑1 (MQ=255)

                                                                           tACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAAGAGTATg              >  2:1494790/1‑93 (MQ=255)

                                                                           tACTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAAGAGGATg              >  1:2012926/1‑93 (MQ=255)

                                                                            aCTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAAGAGTATGa             <  1:4135991/93‑1 (MQ=255)

                                                                             cTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAAGAGTATGAc            <  1:160942/93‑1 (MQ=255)

                                                                             cTGTGTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAAGAGTATGAc            >  2:7040390/1‑93 (MQ=255)

                                                                              tgtgTTTGGTCGTGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAAGAGTATGACg           >  1:8403865/1‑93 (MQ=255)

                                                                              tgtgTTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAAGAGTATGACg           >  2:2178383/1‑93 (MQ=255)

                                                                               gtgtTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAAGAGTATGACGa          <  2:7413222/93‑1 (MQ=255)

                                                                               gtgtTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAAGAGTATGACGa          >  1:5850411/1‑93 (MQ=255)

                                                                               gtgtTGGGTCGCGCACATGAGTTCCACATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGTAATTTACAGACGTAAAAAAAGAGTATGACGa          <  2:800144/93‑1 (MQ=255)

                                                                                 gtTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGTGAATTTACAGACGTAAAAAAAGAGTATGACGAtt        >  2:8401736/1‑93 (MQ=255)

                                                                                 gtTTGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACATAAAAAAAGAGTATGACGAtt        >  1:9944034/1‑93 (MQ=255)

                                                                                  tttGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAAGAGTATGACGAttt       <  1:770594/93‑1 (MQ=255)

                                                                                   ttGGTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAAGAGTATGACGAtttt      <  2:3565929/93‑1 (MQ=255)

                                                                                     ggTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAAGAGTATGACGATTTTGt    >  2:50768/1‑93 (MQ=255)

                                                                                     ggTCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAAGAGTATGACGATTTTGt    <  1:7388911/93‑1 (MQ=255)

                                                                                       tCGCGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAAGAGTATGACGATTTTGTTa  >  1:2298356/1‑93 (MQ=255)

CAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGTGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAAGAGTATGACGATTTTGTTA  >  NC_000913/931177‑931356

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