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.5%	686.3 / ‑3.8	194	coding (1249/1299 nt)	trkH	TrkH
Reads supporting (aligned to +/- strand): ref base . (1/0); new base C (102/91); total (103/91)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01

Only 100 of 216 total aligned reads displayed.

ATCGCGCACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCG‑CACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGATAACTT  >  NC_000913/4031891‑4032073

aTCGCGCACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCg                                                                                             <  1:2987298/93‑1 (MQ=255)

aTCGCGCACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCg                                                                                             <  2:1031521/93‑1 (MQ=255)

  cgcgCACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCG‑cc                                                                                          <  1:1142694/93‑2 (MQ=255)

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

   gcgcACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCG‑cca                                                                                         <  2:4349889/93‑3 (MQ=255)

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

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

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

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

          gCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGt                                                                                   <  2:2124662/93‑1 (MQ=255)

          gCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGATGGCGATTATCGCCACGGGCGt                                                                                   >  1:2252364/1‑93 (MQ=255)

           ccGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTg                                                                                  <  2:5417809/93‑1 (MQ=255)

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

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

               aaCGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATg                                                                              >  2:1474159/1‑93 (MQ=255)

                 cGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGAc                                                                            <  2:2047340/93‑1 (MQ=255)

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

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

                   tATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACtt                                                                          >  1:2889370/1‑93 (MQ=255)

                    ttCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACttt                                                                         <  1:4303074/92‑1 (MQ=255)

                    aTCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACttt                                                                         >  2:2664472/1‑93 (MQ=255)

                     tCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGCGGATTACtttt                                                                        >  1:5394846/1‑93 (MQ=255)

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

                        tCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTg                                                                     >  1:5421767/1‑93 (MQ=255)

                        tCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTg                                                                     >  2:1616166/1‑93 (MQ=255)

                        tCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTg                                                                     >  2:4956982/1‑93 (MQ=255)

                          gAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGcc                                                                   >  1:2873653/1‑93 (MQ=255)

                          gAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGcc                                                                   >  1:435465/1‑93 (MQ=255)

                            agacGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCTTGGATGACTTTTCTGCAtt                                                                 <  1:1937209/90‑1 (MQ=255)

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

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

                                gTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCg                                                             >  2:4563578/1‑93 (MQ=255)

                                gTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCg                                                             <  1:5050093/93‑1 (MQ=255)

                                gTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCg                                                             >  2:3891282/1‑93 (MQ=255)

                                 gttGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGt                                                            <  1:2290024/92‑1 (MQ=255)

                                    ggtgATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCgg                                                         <  1:2124879/90‑1 (MQ=255)

                                    ggggATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCgg                                                         <  2:4014346/93‑1 (MQ=255)

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

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

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

                                     gggATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGgt                                                        <  1:2968943/93‑1 (MQ=255)

                                      ggATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGgtt                                                       <  2:3597209/93‑1 (MQ=255)

                                      ggATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGgtt                                                       <  1:3591400/93‑1 (MQ=255)

                                          ttttCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGttgttg                                                   >  1:3936813/1‑93 (MQ=255)

                                           tttCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGc                                                  <  2:5148316/93‑1 (MQ=255)

                                           tttCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGc                                                  <  2:3918546/93‑1 (MQ=255)

                                             tctcCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGa                                                >  2:2808058/1‑93 (MQ=255)

                                             tctcCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGa                                                >  2:1940126/1‑93 (MQ=255)

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

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

                                                  gCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTg                                           >  2:5240435/1‑93 (MQ=255)

                                                  gCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTg                                           >  1:4942762/1‑93 (MQ=255)

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

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

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

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

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

                                                       tGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATaa                                      <  2:3678724/93‑1 (MQ=255)

                                                        gCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAAc                                     >  1:968090/1‑93 (MQ=255)

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

                                                          aTTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCt                                   <  2:3758770/93‑1 (MQ=255)

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

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

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

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

                                                            tGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTgg                                 <  2:3573952/93‑1 (MQ=255)

                                                            tGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTgg                                 >  1:1224549/1‑93 (MQ=255)

                                                            tGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTgg                                 <  2:1107115/93‑1 (MQ=255)

                                                               tgtTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGc                              >  2:255443/1‑93 (MQ=255)

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

                                                                 tttATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGGCa                            >  1:1400918/1‑93 (MQ=255)

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

                                                                  ttATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAg                           >  1:4201551/1‑93 (MQ=255)

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

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

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

                                                                        tCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTg                     <  1:848476/93‑1 (MQ=255)

                                                                        tCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTg                     <  2:2509142/93‑1 (MQ=255)

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

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

                                                                           gTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCg                  >  1:5871779/1‑93 (MQ=255)

                                                                           gTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCg                  >  1:128844/1‑93 (MQ=255)

                                                                               gCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGt              >  1:2010595/1‑93 (MQ=255)

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

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

                                                                                cTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTgggt             >  1:1632596/1‑91 (MQ=255)

                                                                                  ggCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGc           <  1:6083899/93‑1 (MQ=255)

                                                                                   gCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCt          >  1:2041989/1‑93 (MQ=255)

                                                                                   gCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCt          >  1:801510/1‑93 (MQ=255)

                                                                                   gCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCt          >  2:5073278/1‑93 (MQ=255)

                                                                                    cGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTg         <  1:722119/93‑1 (MQ=255)

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

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

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

                                                                                      aTTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGAt       >  2:3590641/1‑93 (MQ=255)

                                                                                      aTTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGAt       >  1:485340/1‑93 (MQ=255)

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

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

                                                                                        tATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGATaa     <  2:2873653/93‑1 (MQ=255)

                                                                                          tCG‑CGCGGGGGTGGATGACTTTTCTGCCTTTGCGTCGGTTGGTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGATAACtt  <  1:4717837/93‑1 (MQ=255)

ATCGCGCACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCG‑CACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGATAACTT  >  NC_000913/4031891‑4032073

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