BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	473,816	A→G	D190D (GAT→GAC)	tesB ←	acyl‑CoA thioesterase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	473,816	0	A	G	98.6%	894.6 / 0.6	288	D190D (GAT→GAC)	tesB	acyl‑CoA thioesterase II
Reads supporting (aligned to +/- strand): ref base A (0/4); new base G (156/128); total (156/132)
Fisher's exact test for biased strand distribution p-value = 4.30e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.69e-01

Only 100 of 311 total aligned reads displayed.

GAGAAAACCGATGCCGTGCGGCTGTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACG  >  NC_000913/473726‑473905

gagaAAACCGATGCCGTGCGGCTGTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGgtcgtc                                                                                         >  1:9220408/1‑93 (MQ=255)

gagaAAACCGATGCCGGGCGGCTGTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGgtcgtc                                                                                         <  1:4486312/93‑1 (MQ=255)

 agaAAACCGATGCCGTGCGGCTGTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTcc                                                                                        >  1:6423343/1‑93 (MQ=255)

 agaAAACCGATGCCGTGCGGCTGTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTcc                                                                                        >  1:7794865/1‑93 (MQ=255)

 agaAAACCGATGCCGTGCGGCTGTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTcc                                                                                        >  1:2363428/1‑93 (MQ=255)

   aaaaCCGATGCCGTGCGGCTGTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACGGATGAACGCACAGGTCGTCCgg                                                                                      >  1:6847710/1‑93 (MQ=255)

    aaaCCGATGCCGTGCGGCTGTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGc                                                                                     <  2:9406912/93‑1 (MQ=255)

    aaaCCGATGCCGTGCGGCTGTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGc                                                                                     <  1:9387624/93‑1 (MQ=255)

     aaCCGATGCCGTGCGGCTGTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCa                                                                                    <  1:1092907/93‑1 (MQ=255)

     aaCCGATGCCGTGCGGCTGTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCa                                                                                    >  1:3095014/1‑93 (MQ=255)

     aaCCGATGCCGTGCGGCTGTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCa                                                                                    <  2:4532858/93‑1 (MQ=255)

       ccGATGCCGTGCGGCTGTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACg                                                                                  >  1:8771717/1‑93 (MQ=255)

        cGATGCCGTGCGGCTGTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGc                                                                                 >  2:917137/1‑93 (MQ=255)

         gATGCCGTGCGGCTGTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCt                                                                                >  1:8108913/1‑93 (MQ=255)

         gATGCCGTGCGGCTGTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCt                                                                                >  1:7938581/1‑93 (MQ=255)

           tGCCGTGCGGCTGTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTAc                                                                              <  1:8919712/93‑1 (MQ=255)

              cGTGCGGCTGTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCAt                                                                           <  2:2889324/93‑1 (MQ=255)

               gTGCGGCTGTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCAtt                                                                          <  2:6821017/93‑1 (MQ=255)

               gTGCGGCTGTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCAtt                                                                          <  2:3004847/93‑1 (MQ=255)

                tGCGGCTGTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCAttt                                                                         >  1:1271008/1‑93 (MQ=255)

                    gCTGTAGAGCTACCGGCAGGCAGTTAAGATCGGAAACGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTccgg                                                                     >  1:8187892/1‑92 (MQ=255)

                    gCTGTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTgcgc                                                                     >  2:912624/1‑93 (MQ=255)

                       gTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGa                                                                  >  1:351299/1‑93 (MQ=255)

                       gTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGa                                                                  >  1:9967444/1‑93 (MQ=255)

                       gTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGCGTAAACGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCggc                                                                  >  1:8280814/1‑92 (MQ=255)

                        tAGAGCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGAt                                                                 <  1:1198112/93‑1 (MQ=255)

                        tAGAGCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGAt                                                                 <  2:5655625/93‑1 (MQ=255)

                          gagCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATcc                                                               <  1:3686965/93‑1 (MQ=255)

                           agCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGGAGGGCGTCCGGCACGCTACCATTTGCGCGGATccc                                                              >  1:2131795/1‑92 (MQ=255)

                           agCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCca                                                              <  1:1996820/93‑1 (MQ=255)

                            gCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCcac                                                             <  1:7277478/93‑1 (MQ=255)

                              tACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCacac                                                           >  1:4162786/1‑93 (MQ=255)

                              tACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCacac                                                           >  1:3162145/1‑93 (MQ=255)

                               aCCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACAcc                                                          <  2:8612970/93‑1 (MQ=255)

                               aCCGGCAGGAAGGTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGAACGCTACCATTTGCGCGGATCCCCAcc                                                          >  2:5911174/1‑93 (MQ=255)

                                ccGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCt                                                         <  2:6606197/93‑1 (MQ=255)

                                ccGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCt                                                         <  1:5669225/93‑1 (MQ=255)

                                 cGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGACGTCCGGCACGCTACCATTTGCGCGGATCCACCCCTg                                                        >  1:8415237/1‑93 (MQ=255)

                                  ggCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGa                                                       <  2:2210365/93‑1 (MQ=255)

                                  ggCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGa                                                       <  2:8752401/93‑1 (MQ=255)

                                  ggCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGa                                                       <  1:1777873/93‑1 (MQ=255)

                                    cAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACg                                                     >  2:2563413/1‑93 (MQ=255)

                                     aGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGa                                                    >  1:5701181/1‑93 (MQ=255)

                                       gaggTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTCCCATTTGCGCGGATCCACCCCTGACGAtg                                                  <  1:348388/90‑1 (MQ=255)

                                       gAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGAtg                                                  <  2:4433477/93‑1 (MQ=255)

                                          gTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTgg                                               <  2:3894320/93‑1 (MQ=255)

                                           ttAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGt                                              >  2:2164043/1‑93 (MQ=255)

                                             aaGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTc                                            <  1:2621711/93‑1 (MQ=255)

                                             aaGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTc                                            >  2:1623929/1‑93 (MQ=255)

                                              aGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCt                                           >  1:7049522/1‑93 (MQ=255)

                                              aGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCt                                           <  2:4054381/93‑1 (MQ=255)

                                                aTCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGc                                         <  2:1127268/93‑1 (MQ=255)

                                                aTCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGc                                         >  1:6041414/1‑93 (MQ=255)

                                                aTCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGc                                         >  1:6566990/1‑93 (MQ=255)

                                                    gAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACg                                     <  2:8492066/93‑1 (MQ=255)

                                                     aaGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGt                                    >  1:9608061/1‑93 (MQ=255)

                                                     aaGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGt                                    <  2:7430060/93‑1 (MQ=255)

                                                      aGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTg                                   <  2:4732844/93‑1 (MQ=255)

                                                        cGTAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGAc                                 >  2:1357386/1‑93 (MQ=255)

                                                          tAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCt                               >  2:5284108/1‑93 (MQ=255)

                                                          tAACCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCt                               >  2:731217/1‑93 (MQ=255)

                                                           aaCCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCtt                              <  2:5155497/93‑1 (MQ=255)

                                                            aCCGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCttt                             <  1:2978121/93‑1 (MQ=255)

                                                              cGAGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCa                           <  2:8342706/93‑1 (MQ=255)

                                                                aGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAgg                         >  2:8945151/1‑92 (MQ=255)

                                                                aGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGt                         >  1:5428116/1‑93 (MQ=255)

                                                                aGCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGt                         <  2:5028671/93‑1 (MQ=255)

                                                                 gCAGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTg                        <  1:8498997/93‑1 (MQ=255)

                                                                   aGATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTggg                      <  1:7459561/93‑1 (MQ=255)

                                                                    gATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGt                     <  1:471334/93‑1 (MQ=255)

                                                                    gATACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGt                     >  2:5063432/1‑93 (MQ=255)

                                                                     aTACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGtt                    >  2:6245828/1‑93 (MQ=255)

                                                                      tACTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTAAGTGGGTTa                   <  1:3350794/93‑1 (MQ=255)

                                                                       aCTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCCCACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTAt                  <  1:4454953/93‑1 (MQ=255)

                                                                       aCTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTAt                  >  2:1489164/1‑93 (MQ=255)

                                                                       aCTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTAt                  >  2:354437/1‑93 (MQ=255)

                                                                        cTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATg                 >  1:5902368/1‑93 (MQ=255)

                                                                        cTGATGAACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATg                 >  2:7965058/1‑93 (MQ=255)

                                                                         tgatgaACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGa                >  2:4178388/1‑93 (MQ=255)

                                                                         tgatgaACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGa                >  1:3947039/1‑93 (MQ=255)

                                                                          gatgaACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGaa               >  1:6142460/1‑93 (MQ=255)

                                                                           atgaACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGaaa              <  2:1429604/93‑1 (MQ=255)

                                                                           atgaACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGaaa              >  2:2400968/1‑93 (MQ=255)

                                                                            tgaACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGGCCTTTCAGTGGGTTATGAAAc             <  2:8765256/93‑1 (MQ=255)

                                                                            tgaACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGaact             >  1:148727/1‑91 (MQ=255)

                                                                            tgaACGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAAc             <  1:9777236/93‑1 (MQ=255)

                                                                            taaaCGCGCAGGTCGTCCGGCACGCTACCATTTACGCGGATCCACACCTGACGATGAGGATCTGCTCCGTGACCTTTCAGTGGGGTATGaacc             >  2:3698989/3‑91 (MQ=255)

                                                                               aCGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTcc          >  2:3249363/1‑93 (MQ=255)

                                                                               aCGCGCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTcc          <  2:8086566/93‑1 (MQ=255)

                                                                                cgcgCAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCa         <  2:463498/93‑1 (MQ=255)

                                                                                  cgcAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCAcc       <  1:2693875/93‑1 (MQ=255)

                                                                                  cgcAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCAcc       >  2:1830569/1‑93 (MQ=255)

                                                                                  cgcAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCAcc       <  2:9609525/93‑1 (MQ=255)

                                                                                  cgcAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCAcc       >  2:2678592/1‑93 (MQ=255)

                                                                                   gcAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACcg      >  2:643111/1‑93 (MQ=255)

                                                                                   gcAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACcg      >  1:5317915/1‑93 (MQ=255)

                                                                                    cAGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACcgg     <  2:7251219/93‑1 (MQ=255)

                                                                                     aGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACcgga    <  1:4999117/93‑1 (MQ=255)

                                                                                     aGGTCGTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACcgga    <  2:5385914/93‑1 (MQ=255)

                                                                                       gtcgtcCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACcggacg  <  2:2108485/93‑1 (MQ=255)

GAGAAAACCGATGCCGTGCGGCTGTAGAGCTACCGGCAGGAAGTTAAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACG  >  NC_000913/473726‑473905

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