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.4%	552.0 / NA	172	A417P (GCA→CCA)	trkH	TrkH
Reads supporting (aligned to +/- strand): ref base . (0/0); major base C (91/80); minor base T (0/1); total (91/81)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 190 total aligned reads displayed.

ATCGCGCACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCG‑CACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGATA  >  NC_000913/4031891‑4032069

aTCGCGCACTGCCGGAACGTATCGTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCg                                                                                         <  2:3389409/93‑1 (MQ=255)

aTCGCGCACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCg                                                                                         >  2:4327969/1‑93 (MQ=255)

aTCGCGCACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCg                                                                                         >  1:3280491/1‑93 (MQ=255)

aTCGCGCACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCg                                                                                         >  2:410171/1‑93 (MQ=255)

aTCGCGCACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCg                                                                                         >  1:2199614/1‑93 (MQ=255)

   gcgcgcTGCCGGAACGTTTCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCG‑cca                                                                                     <  1:1359765/93‑3 (MQ=255)

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

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

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

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

       aCTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACggg                                                                                  <  2:472750/93‑1 (MQ=255)

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

       aCTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACggg                                                                                  <  1:2062145/93‑1 (MQ=255)

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

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

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

              gAACGTATCCTCGAAGCCGTTTGGGTATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGAt                                                                           <  2:1021419/93‑1 (MQ=255)

               aaCGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATg                                                                          <  2:534951/93‑1 (MQ=255)

               aaCGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATg                                                                          <  1:2831117/93‑1 (MQ=255)

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

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

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

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

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

                   tATCCTCGAAGCCATTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACtt                                                                      <  2:328787/93‑1 (MQ=255)

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

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

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

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

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

                         cGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGc                                                                >  2:1932402/1‑93 (MQ=255)

                           aaGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGGGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGACt                                                              >  2:3786432/1‑93 (MQ=255)

                            aGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCtt                                                             <  1:1048684/93‑1 (MQ=255)

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

                             gCCGTTTGGGGATTTTTCTCAGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCttt                                                            <  1:4024485/93‑1 (MQ=255)

                             gCCGTGTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCttt                                                            <  2:4467184/93‑1 (MQ=255)

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

                              ccGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTg                                                           >  1:2034915/1‑93 (MQ=255)

                              ccGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTg                                                           >  2:407421/1‑93 (MQ=255)

                               cGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGTTGACTTTTCTGCCTTTGc                                                          <  2:2379254/93‑1 (MQ=255)

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

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

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

                                       gATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGgttg                                                  >  2:1886592/1‑93 (MQ=255)

                                       gATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGCCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGgttg                                                  <  1:2578388/93‑1 (MQ=255)

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

                                         tttttCTCCGCCTATGCATTCGTGTTTATTGTCAGTATGCTGGAGATTATCGCCACGGTCGTGGATGACTTTTCTGCCTTTGAGTCAgttgtt                                                >  2:2478026/1‑93 (MQ=255)

                                              ctcCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGAc                                           >  2:1854155/1‑93 (MQ=255)

                                               tcCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGAct                                          >  1:4673721/1‑92 (MQ=255)

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

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

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

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

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

                                                    cTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGaa                                     >  2:449554/1‑93 (MQ=255)

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

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

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

                                                      aTGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACCTTGAATa                                   >  1:1913816/1‑93 (MQ=255)

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

                                                       tGCATTGTTGTTTATTGTCAGTATGCTGGCGATTAGCGCCAAGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATaa                                  <  1:1474239/93‑1 (MQ=255)

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

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

                                                          aTTGGTGTTTATTGTCAGTATGCTTGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCt                               >  1:2919374/1‑93 (MQ=255)

                                                          aTTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATACCCt                               >  1:3833895/1‑93 (MQ=255)

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

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

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

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

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

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

                                                             ggTGTTTATTGTCAGTATGCTGGCGATTATCGTCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTggg                            <  2:2105077/93‑1 (MQ=255)

                                                             ggTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACCTTGAATAACCTggg                            >  1:2796945/1‑93 (MQ=255)

                                                              gtgtTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTgggg                           >  2:2274463/1‑93 (MQ=255)

                                                              gtgtTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTgggg                           >  1:4719606/1‑93 (MQ=255)

                                                              gtgtTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTgggg                           >  2:1150793/1‑93 (MQ=255)

                                                                gtTTAGTGTCAGTATTCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGcc                         <  1:1529082/93‑1 (MQ=255)

                                                                 tttATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCa                        >  2:2875382/1‑93 (MQ=255)

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

                                                                  ttATTGTCAGTATGCTGGCGATTATCGCCACGGGCGGGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAg                       >  1:3479578/1‑93 (MQ=255)

                                                                    aTTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAggg                     >  1:2505339/1‑93 (MQ=255)

                                                                    aTTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGATTGTTGCGACATTGAATAACCTGGGGCCgggc                     >  1:2505338/1‑92 (MQ=255)

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

                                                                      tGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCt                   <  1:4558618/93‑1 (MQ=255)

                                                                      tGTCAGAATGCTGGCGTTTGTGG‑CCACGGCGTGGATGCCTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCtt                  <  1:639397/93‑1 (MQ=255)

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

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

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

                                                                          aGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTTAATAACCTGGGGCCAGGGCTTGGc               >  1:2136454/1‑93 (MQ=255)

                                                                            tATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCgg             >  1:3775765/1‑92 (MQ=255)

                                                                             aTGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTg            <  1:1782983/93‑1 (MQ=255)

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

                                                                               gCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTTGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGt          <  1:405526/93‑1 (MQ=255)

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

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

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

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

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

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

                                                                                      aTTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGAt   >  1:903208/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` `< 20 ≤` `ATCG/ATCG` `< 29 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`