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%	733.1 / ‑4.3	234	coding (1249/1299 nt)	trkH	TrkH
Reads supporting (aligned to +/- strand): ref base . (0/1); new base C (129/104); total (129/105)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 250 total aligned reads displayed.

ATCGCGCACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCG‑C‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGATA  >  NC_000913/4031891‑4032069

   gcgcACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCTTTGGTGTTTATTGTCAGTATGCTGGCGATTATCG‑cca                                                                                      <  2:3011777/93‑3 (MQ=255)

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

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

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

       aCTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACggg                                                                                  <  1:285314/93‑1 (MQ=255)

       aCTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACggg                                                                                  >  2:1115857/1‑93 (MQ=255)

        cTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGc                                                                                 <  1:1287933/93‑1 (MQ=255)

         tGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCg                                                                                >  2:4298506/1‑93 (MQ=255)

         tGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCg                                                                                >  2:1890775/1‑93 (MQ=255)

          gCCGGAACGTATCCTGGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGt                                                                               <  2:3706133/93‑1 (MQ=255)

           ccGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTg                                                                              >  1:149394/1‑93 (MQ=255)

           ccGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTg                                                                              <  1:3325635/93‑1 (MQ=255)

            cGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTgg                                                                             <  1:533600/93‑1 (MQ=255)

             ggAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGa                                                                            <  1:2047357/93‑1 (MQ=255)

             cgAACGTATCCTCGAAGCCGTTTGGGGAGTTGTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGGGGa                                                                            <  1:3172678/92‑1 (MQ=255)

              gAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGAt                                                                           >  2:3531357/1‑93 (MQ=255)

                aCGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGa                                                                         >  1:2472962/1‑93 (MQ=255)

                  gTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACt                                                                       >  1:4101993/1‑93 (MQ=255)

                  gTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACt                                                                       >  1:913354/1‑93 (MQ=255)

                   tATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACtt                                                                      >  1:3977044/1‑93 (MQ=255)

                   tATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACtt                                                                      <  2:1960776/93‑1 (MQ=255)

                   tATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACtt                                                                      <  1:1193429/93‑1 (MQ=255)

                    aTCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACttt                                                                     >  2:3221992/1‑93 (MQ=255)

                     tCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACtttt                                                                    >  1:3231650/1‑93 (MQ=255)

                      ccTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTc                                                                   >  1:1565264/1‑93 (MQ=255)

                      ccTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTc                                                                   <  1:600385/93‑1 (MQ=255)

                      ccTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTc                                                                   <  1:4475312/93‑1 (MQ=255)

                        tCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTg                                                                 >  2:1396880/1‑93 (MQ=255)

                          gAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGcc                                                               <  2:378260/93‑1 (MQ=255)

                           aaGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCt                                                              <  2:694293/93‑1 (MQ=255)

                            aGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCtt                                                             <  1:3900123/93‑1 (MQ=255)

                            aGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCtt                                                             >  1:834375/1‑93 (MQ=255)

                            aGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCtt                                                             >  2:1225742/1‑93 (MQ=255)

                             gCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCttt                                                            >  2:1801234/1‑93 (MQ=255)

                             gCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCttt                                                            <  2:3067589/93‑1 (MQ=255)

                             gCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCttt                                                            >  2:659529/1‑93 (MQ=255)

                             gCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCttt                                                            <  2:902453/93‑1 (MQ=255)

                             gCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCttt                                                            <  1:1920388/93‑1 (MQ=255)

                               cGTCTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGc                                                          <  1:3848805/93‑1 (MQ=255)

                                gTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCg                                                         >  1:2828133/1‑93 (MQ=255)

                                 tttGGGGTTTTTTCTCCGCCTATGCGTTGGTGTTTCTTGTCAGTATGCTGGCGATTATCG‑CTACGGGCGGGGCTGACTTTTCTGCCTTTGCGt                                                        >  2:3953885/1‑93 (MQ=255)

                                 tttGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGt                                                        >  2:1093073/1‑93 (MQ=255)

                                  ttGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTc                                                       <  2:1259582/93‑1 (MQ=255)

                                  ttGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTc                                                       <  1:688564/93‑1 (MQ=255)

                                      ggATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGgtt                                                   <  1:2423748/93‑1 (MQ=255)

                                         tttttCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGgttgtt                                                <  2:2289613/93‑1 (MQ=255)

                                          ttttCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGttgttg                                               >  2:2704199/1‑93 (MQ=255)

                                          ttttCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGttgttg                                               >  2:827935/1‑93 (MQ=255)

                                            ttCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCg                                             >  1:2119969/1‑93 (MQ=255)

                                            ttCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCg                                             <  2:4219007/93‑1 (MQ=255)

                                               tcCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACa                                          >  1:4571081/1‑93 (MQ=255)

                                               tcCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACa                                          >  2:3496719/1‑93 (MQ=255)

                                                 cGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACAtt                                        <  2:1622363/93‑1 (MQ=255)

                                                  gCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTg                                       >  1:2950824/1‑93 (MQ=255)

                                                   ccTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGa                                      >  1:4445511/1‑93 (MQ=255)

                                                   catATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGa                                      <  2:555858/91‑1 (MQ=255)

                                                    cTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGaa                                     <  2:2853523/93‑1 (MQ=255)

                                                    cTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGaa                                     >  1:1761295/1‑93 (MQ=255)

                                                     tATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAAt                                    >  2:241814/1‑93 (MQ=255)

                                                     tATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAAt                                    >  1:2083433/1‑93 (MQ=255)

                                                     tATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTGGTTGCGACATTGAAt                                    <  2:1513324/93‑1 (MQ=255)

                                                      aTGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATa                                   >  2:4440932/1‑93 (MQ=255)

                                                      aTGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATa                                   <  1:2157825/93‑1 (MQ=255)

                                                       tGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATaa                                  >  1:1028508/1‑93 (MQ=255)

                                                        gCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAAc                                 <  1:4157474/93‑1 (MQ=255)

                                                        gCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAAc                                 <  2:2473358/93‑1 (MQ=255)

                                                        gCATTGATGTTTATTGTCAGTATCCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAAc                                 <  2:4081553/93‑1 (MQ=255)

                                                         cATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAAcc                                <  2:3397837/93‑1 (MQ=255)

                                                         cATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAAcc                                >  2:617487/1‑93 (MQ=255)

                                                         cATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTCCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAAcc                                <  1:2205903/93‑1 (MQ=255)

                                                          aTTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCt                               <  1:1718186/93‑1 (MQ=255)

                                                          aTTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCt                               <  1:2994921/93‑1 (MQ=255)

                                                          aTTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCt                               <  1:3069316/93‑1 (MQ=255)

                                                           ttGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTg                              >  1:2731024/1‑93 (MQ=255)

                                                           ttGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTg                              >  1:4630588/1‑93 (MQ=255)

                                                           ttGGTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTg                              <  1:4599815/93‑1 (MQ=255)

                                                             ggTGTTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTggg                            <  1:3049618/93‑1 (MQ=255)

                                                               tgtTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGc                          >  1:3379436/1‑93 (MQ=255)

                                                               tgtTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGc                          <  2:4631426/93‑1 (MQ=255)

                                                                gtTTATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGcc                         >  2:2934689/1‑93 (MQ=255)

                                                                 tttATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCa                        <  1:4004842/93‑1 (MQ=255)

                                                                  ttATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAg                       <  2:1834740/93‑1 (MQ=255)

                                                                  ttATTGTCACTATGCTGGCGATTATCGCC‑ACGGGCGTAGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCAACATTGAATAACCTGGGGCCAg                       <  1:2352199/93‑1 (MQ=255)

                                                                   tATTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAgg                      >  1:647353/1‑93 (MQ=255)

                                                                    aTTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAggg                     >  2:2635735/1‑93 (MQ=255)

                                                                    aTTGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAggg                     <  2:2046550/93‑1 (MQ=255)

                                                                     ttGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGc                    >  1:232462/1‑93 (MQ=255)

                                                                      tGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCt                   <  2:4553039/93‑1 (MQ=255)

                                                                      tGTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGAGGCCAGGGCt                   >  1:182535/1‑93 (MQ=255)

                                                                       gTCAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCtt                  <  2:4196019/93‑1 (MQ=255)

                                                                         cAGTATGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTgg                >  1:3539016/1‑93 (MQ=255)

                                                                            tATGCTGGCGATTATCGCC‑ACGGGCGAAAATGACTTTTCTGCCTTTGCATCGATTGTTGCGACATTGAATAAACTGGGGCCAGGGCTTGGCGt             <  1:4680229/93‑1 (MQ=255)

                                                                             aTGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTg            <  1:4206291/93‑1 (MQ=255)

                                                                              tGCTGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTgg           >  1:257065/1‑93 (MQ=255)

                                                                                 tGGCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTCGGGCCAGGGCTTGGCGTGGTTg        >  2:2260596/1‑93 (MQ=255)

                                                                                  ggCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGc       >  1:1428557/1‑93 (MQ=255)

                                                                                  ggCGATTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGc       >  1:1066735/1‑93 (MQ=255)

                                                                                      aTTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGAt   >  2:3305544/1‑93 (MQ=255)

                                                                                      aTTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGAt   >  2:3204455/1‑93 (MQ=255)

                                                                                      aTTATCGCC‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGAt   >  1:1406917/1‑93 (MQ=255)

ATCGCGCACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCG‑C‑ACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGATA  >  NC_000913/4031891‑4032069

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