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%	696.2 / ‑5.7	220	A417P (GCA→CCA)	trkH	TrkH
Reads supporting (aligned to +/- strand): ref base . (1/0); major base C (117/101); minor base T (1/0); total (119/101)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 238 total aligned reads displayed.

ATCGCGCACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCG‑CACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGATA  >  NC_000913/4031891‑4032069

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

 tCGCGCACTGCGGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCG‑c                                                                                       <  2:444333/93‑1 (MQ=255)

 tCGCGCACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCG‑c                                                                                       >  2:3620580/1‑93 (MQ=255)

 tCGCGCACTGCCGAAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCG‑c                                                                                       <  2:5710051/93‑1 (MQ=255)

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

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

       aCTGCCGCAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACggg                                                                                  <  1:5097037/93‑1 (MQ=255)

        cTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGc                                                                                 >  1:2710761/1‑93 (MQ=255)

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

          gCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCgg                                                                               >  2:3904696/1‑92 (MQ=255)

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

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

                aCGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGa                                                                         >  1:4627569/1‑93 (MQ=255)

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

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

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

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

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

                     tCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACtttt                                                                    >  1:1760255/1‑93 (MQ=255)

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

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

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

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

                           aaGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCt                                                              <  2:529943/93‑1 (MQ=255)

                           aaGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCt                                                              <  2:2683033/93‑1 (MQ=255)

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

                             gCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCttt                                                            >  2:5614093/1‑93 (MQ=255)

                             gCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGGTTATTGTCAGTATGCTGGCGATTATCGTCACGGGCGTGGATGACTTTTCTGCCttt                                                            >  2:4770055/1‑93 (MQ=255)

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

                               cGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGc                                                          >  2:646008/1‑93 (MQ=255)

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

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

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

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

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

                                 tttGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGt                                                        <  1:1134752/93‑1 (MQ=255)

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

                                  ttGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTc                                                       >  1:3260750/1‑93 (MQ=255)

                                  ttGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTc                                                       >  1:4911805/1‑93 (MQ=255)

                                  ttGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTc                                                       >  2:1034735/1‑93 (MQ=255)

                                   tGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCg                                                      >  1:14005/1‑93 (MQ=255)

                                   tGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCg                                                      <  2:5270725/93‑1 (MQ=255)

                                   tGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCg                                                      <  2:3303764/93‑1 (MQ=255)

                                    tgggATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCgg                                                     <  2:4712691/92‑1 (MQ=255)

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

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

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

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

                                          ttttCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATAGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGttgttg                                               >  1:4018700/1‑93 (MQ=255)

                                             tctcCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGa                                            <  1:2220125/93‑1 (MQ=255)

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

                                              ctcCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGAc                                           <  1:243756/93‑1 (MQ=255)

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

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

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

                                                ccGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGTGTCCGTTGTTGCGACAt                                         >  1:96088/1‑93 (MQ=255)

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

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

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

                                                  gCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTg                                       <  1:859011/93‑1 (MQ=255)

                                                   ccTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGa                                      >  1:5738615/1‑93 (MQ=255)

                                                    cTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGaa                                     <  2:2903534/93‑1 (MQ=255)

                                                    cTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGaa                                     <  2:1527917/93‑1 (MQ=255)

                                                     tATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAAt                                    <  1:1822130/93‑1 (MQ=255)

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

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

                                                       tGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATaa                                  >  1:3485595/1‑93 (MQ=255)

                                                       tGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATaa                                  >  1:1612422/1‑93 (MQ=255)

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

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

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

                                                         gaTTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAAcc                                <  1:283985/92‑1 (MQ=255)

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

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

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

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

                                                          aTTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGGTGTTGCGGCCTTGAATACCCt                               >  2:4546657/1‑93 (MQ=255)

                                                            tGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTAGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTgg                             <  2:1343150/93‑1 (MQ=255)

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

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

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

                                                               tgtTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGGCTTTTCTGCCTTTGCGTCGGTTGTTTCCACCTTGAATAACCTAGGGc                          >  2:3504383/1‑93 (MQ=255)

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

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

                                                                gtTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGcc                         <  2:4852013/93‑1 (MQ=255)

                                                                       gTCAGTATGCGGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCGTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCtt                  <  2:2977619/93‑1 (MQ=255)

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

                                                                          aGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGc               >  1:1882124/1‑93 (MQ=255)

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

                                                                           gTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCg              <  1:251018/93‑1 (MQ=255)

                                                                             aTGCTGGCGATTATCGCCACGGGCGTGGATGACTATTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGCGCTTGGCGTg            >  2:4830623/1‑93 (MQ=255)

                                                                              tGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTgg           >  2:2481699/1‑93 (MQ=255)

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

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

                                                                                cTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGtt         <  1:4751260/93‑1 (MQ=255)

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

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

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

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

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

ATCGCGCACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCG‑CACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGATA  >  NC_000913/4031891‑4032069

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