BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,031,984	+C	coding (1250/1299 nt)	trkH →	TrkH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,031,983	1	.	C	99.6%	1664.0 / ‑3.2	490	coding (1249/1299 nt)	trkH	TrkH
Reads supporting (aligned to +/- strand): ref base . (0/2); new base C (267/221); total (267/223)
Fisher's exact test for biased strand distribution p-value = 2.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 538 total aligned reads displayed.

ATCGCGCACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCG‑CACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGATAA  >  NC_000913/4031891‑4032070

   gcgcACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCG‑cca                                                                                      >  1:1000869/1‑91 (MQ=255)

   gcgcACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCG‑cca                                                                                      >  1:4290688/1‑91 (MQ=255)

   gcgcACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCG‑cca                                                                                      >  1:6877078/1‑91 (MQ=255)

   gcgcACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCG‑cca                                                                                      <  1:5699009/93‑3 (MQ=255)

    cgcACTGCCGGAACGTATCCTCGAAGTCGTTTGGGGATTTTTCTCCGCCTATGCATTGGCGTTTATTGTCAGTATGCTGGCGATTATCG‑ccac                                                                                     <  1:6433162/93‑4 (MQ=255)

    cgcACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCG‑ccac                                                                                     <  1:1252194/93‑4 (MQ=255)

    cgcACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCG‑ccac                                                                                     <  2:3754156/93‑4 (MQ=255)

     gcACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCG‑ccacg                                                                                    >  2:7149951/1‑89 (MQ=255)

     gcACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCG‑ccacg                                                                                    <  1:4135932/93‑5 (MQ=255)

       aCTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACggg                                                                                   >  2:1636509/1‑93 (MQ=255)

         tGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCg                                                                                 <  2:4501671/93‑1 (MQ=255)

         tGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCg                                                                                 <  1:3228799/93‑1 (MQ=255)

         tGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCg                                                                                 >  2:1329383/1‑93 (MQ=255)

            cGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTgg                                                                              <  2:7336576/93‑1 (MQ=255)

             ggAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGa                                                                             >  1:3597993/1‑93 (MQ=255)

             ggAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGa                                                                             <  1:3702694/93‑1 (MQ=255)

             ggAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGa                                                                             >  2:2434178/1‑93 (MQ=255)

             ggAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGa                                                                             <  1:2605438/93‑1 (MQ=255)

             ggAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGa                                                                             >  2:860487/1‑93 (MQ=255)

                aCGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGa                                                                          <  1:1985443/93‑1 (MQ=255)

                aCGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGa                                                                          <  2:6032916/93‑1 (MQ=255)

                  gTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACt                                                                        >  2:2819445/1‑93 (MQ=255)

                   tATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACtt                                                                       <  2:2312061/93‑1 (MQ=255)

                    aTCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACttt                                                                      <  2:2179165/93‑1 (MQ=255)

                    aTCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACttt                                                                      <  1:3803875/93‑1 (MQ=255)

                      ccTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTc                                                                    <  1:2431414/93‑1 (MQ=255)

                      ccTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTc                                                                    <  2:5693827/93‑1 (MQ=255)

                      ccTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTc                                                                    <  1:536594/93‑1 (MQ=255)

                       cTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCt                                                                   >  2:3467541/1‑93 (MQ=255)

                          gAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGcc                                                                <  2:5826989/93‑1 (MQ=255)

                           aaGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCt                                                               >  2:4140900/1‑93 (MQ=255)

                            aGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCtt                                                              >  1:1166036/1‑93 (MQ=255)

                             gCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCttt                                                             <  1:5444689/93‑1 (MQ=255)

                             gCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCttt                                                             <  2:3327065/93‑1 (MQ=255)

                              ccGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTg                                                            <  1:2968389/93‑1 (MQ=255)

                              ccGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTg                                                            <  1:3130937/93‑1 (MQ=255)

                              ccGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTg                                                            <  2:21057/93‑1 (MQ=255)

                               cGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGc                                                           <  1:1894744/93‑1 (MQ=255)

                                 tttGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGt                                                         >  2:4336211/1‑93 (MQ=255)

                                    ggggATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCgg                                                      >  1:4589440/1‑93 (MQ=255)

                                     gggATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGgt                                                     >  1:5788806/1‑93 (MQ=255)

                                       gATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGgttg                                                   <  1:2673258/93‑1 (MQ=255)

                                        aTTTTTCTCCGCCTATGCTTTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGgttgt                                                  <  1:2236398/93‑1 (MQ=255)

                                        aTTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGgttgt                                                  <  1:2140996/93‑1 (MQ=255)

                                        aTTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGgttgt                                                  <  1:6602882/93‑1 (MQ=255)

                                        aTTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGgttgt                                                  >  2:7471732/1‑93 (MQ=255)

                                         tttttCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGgttgtt                                                 >  1:1428929/1‑93 (MQ=255)

                                         tttttCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGgttgtt                                                 >  2:2586926/1‑93 (MQ=255)

                                          ttttCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGttgttg                                                <  1:1547592/93‑1 (MQ=255)

                                               tcCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACa                                           >  1:1805471/1‑93 (MQ=255)

                                               tcCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACa                                           <  2:6721996/93‑1 (MQ=255)

                                               tcCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACa                                           >  1:6763751/1‑93 (MQ=255)

                                               tcCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACa                                           >  1:3350885/1‑93 (MQ=255)

                                               tcCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGGCGGTTGTTGCGACa                                           <  2:6450629/93‑1 (MQ=255)

                                                ccGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACAt                                          <  2:6165985/93‑1 (MQ=255)

                                                 cGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACAtt                                         <  2:3545212/93‑1 (MQ=255)

                                                   ccTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGa                                       <  1:4399779/93‑1 (MQ=255)

                                                     tATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAAt                                     >  2:1693365/1‑93 (MQ=255)

                                                      aTGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATa                                    >  1:5247754/1‑93 (MQ=255)

                                                      aTGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATa                                    <  1:7482638/93‑1 (MQ=255)

                                                      aTGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATa                                    >  2:3104238/1‑93 (MQ=255)

                                                        gCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAAc                                  <  1:7234863/93‑1 (MQ=255)

                                                         cATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAAcc                                 >  1:3549351/1‑93 (MQ=255)

                                                         cATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAAcc                                 <  2:5221866/93‑1 (MQ=255)

                                                         cATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAAcc                                 <  1:3419862/93‑1 (MQ=255)

                                                          aTTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCt                                >  2:5511223/1‑93 (MQ=255)

                                                           ttGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTg                               >  1:6537506/1‑93 (MQ=255)

                                                           ttGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTg                               <  1:6033278/93‑1 (MQ=255)

                                                           ttGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTg                               <  1:5983349/93‑1 (MQ=255)

                                                             ggTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTggg                             <  2:2007567/93‑1 (MQ=255)

                                                             ggTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTggg                             <  2:754588/93‑1 (MQ=255)

                                                             ggTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTggg                             >  1:4841744/1‑93 (MQ=255)

                                                              gtgtTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTgggg                            <  2:3846719/93‑1 (MQ=255)

                                                               tgtTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGc                           <  2:4703380/93‑1 (MQ=255)

                                                                gtTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGcc                          >  1:4773370/1‑93 (MQ=255)

                                                                gtTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGcc                          >  1:2902300/1‑93 (MQ=255)

                                                                 tttATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCa                         <  2:14904/93‑1 (MQ=255)

                                                                   tATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAgg                       >  1:5904698/1‑93 (MQ=255)

                                                                     ttGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGc                     >  2:1147613/1‑93 (MQ=255)

                                                                     ttGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGc                     <  2:6642644/93‑1 (MQ=255)

                                                                      tGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCt                    >  2:5455368/1‑93 (MQ=255)

                                                                       gTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCtt                   <  1:5523146/93‑1 (MQ=255)

                                                                       gTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCtt                   <  2:5058411/93‑1 (MQ=255)

                                                                         cAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTgg                 <  2:5939331/93‑1 (MQ=255)

                                                                          aGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGc                <  2:1035849/93‑1 (MQ=255)

                                                                          aGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGc                <  2:4233446/93‑1 (MQ=255)

                                                                              tGCTGGCGATTATCGCCCCGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTgg            >  2:4902709/1‑93 (MQ=255)

                                                                              tGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTgg            <  2:2743677/93‑1 (MQ=255)

                                                                               gCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGt           <  1:1620110/93‑1 (MQ=255)

                                                                                 tGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTg         >  1:4190145/1‑93 (MQ=255)

                                                                                  ggCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGc        >  2:6322702/1‑93 (MQ=255)

                                                                                  ggCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGc        >  1:519220/1‑93 (MQ=255)

                                                                                   gCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCt       <  1:848864/93‑1 (MQ=255)

                                                                                     gATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGa     <  2:2538550/93‑1 (MQ=255)

                                                                                     gATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGa     >  1:3913987/1‑93 (MQ=255)

                                                                                     gATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGa     <  2:6931473/93‑1 (MQ=255)

                                                                                       ttATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGATa   >  2:5402816/1‑93 (MQ=255)

                                                                                       ttATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGATa   <  1:4973369/93‑1 (MQ=255)

                                                                                       ttATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGATa   >  1:6966144/1‑93 (MQ=255)

                                                                                        tATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGATaa  >  2:1784465/1‑93 (MQ=255)

ATCGCGCACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCG‑CACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGATAA  >  NC_000913/4031891‑4032070

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