BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	791,797	T→C	G405G (GGA→GGG)	modF ←	ModF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	791,797	0	T	C	99.6%	857.9 / NA	250	G405G (GGA→GGG)	modF	ModF
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (138/111); total (139/111)
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 269 total aligned reads displayed.

CAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCTGGC  >  NC_000913/791707‑791885

  gCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCcccc                                                                                      <  1:4830929/93‑1 (MQ=255)

  gCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCcccc                                                                                      <  1:5368727/93‑1 (MQ=255)

  gCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCcccc                                                                                      >  2:8457539/1‑93 (MQ=255)

    gtagcAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAg                                                                                    <  2:8086323/89‑1 (MQ=255)

    ggATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGCCAATCAGCGCCAGACGCTGCTGCCCCCAg                                                                                    >  2:7938550/1‑93 (MQ=255)

    ggATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAg                                                                                    >  1:4615403/1‑93 (MQ=255)

    ggATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAg                                                                                    <  1:6609728/93‑1 (MQ=255)

      aTCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGa                                                                                  <  2:6247428/93‑1 (MQ=255)

       tCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGaa                                                                                 <  2:5492058/93‑1 (MQ=255)

         aaGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAg                                                                               >  1:6714703/1‑93 (MQ=255)

           gCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGAc                                                                             >  2:5860391/1‑93 (MQ=255)

           gCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGAc                                                                             >  1:8187219/1‑93 (MQ=255)

              ccTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTAt                                                                          >  2:8666106/1‑93 (MQ=255)

              ccTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTAt                                                                          >  2:6445734/1‑93 (MQ=255)

               cTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATg                                                                         >  1:4949996/1‑93 (MQ=255)

               cTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATg                                                                         <  2:1657864/93‑1 (MQ=255)

                tGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATgg                                                                        >  2:3484364/1‑93 (MQ=255)

                tGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATgg                                                                        >  1:8659209/1‑93 (MQ=255)

                 gTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCCGACGCTGCTGTCCCCCGGGAAGAATATGGa                                                                       >  2:5129591/1‑93 (MQ=255)

                 gTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGa                                                                       >  2:153847/1‑93 (MQ=255)

                  tAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGaa                                                                      <  2:3941236/93‑1 (MQ=255)

                  tAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGaa                                                                      <  2:5040055/93‑1 (MQ=255)

                  tAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGaa                                                                      >  1:3682366/1‑93 (MQ=255)

                    gTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGCCTATGGAACg                                                                    >  2:7721145/1‑93 (MQ=255)

                    gTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACg                                                                    >  2:1256576/1‑93 (MQ=255)

                    gTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACg                                                                    >  1:4026607/1‑93 (MQ=255)

                     tGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCGGCCCCCAGGAAAGACTATGGAACgg                                                                   <  1:3038465/93‑1 (MQ=255)

                     aggTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACgg                                                                   <  1:6705153/92‑1 (MQ=255)

                      ggTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGa                                                                  >  2:566055/1‑93 (MQ=255)

                      ggTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGa                                                                  <  1:2621497/93‑1 (MQ=255)

                       gcccATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAg                                                                 >  1:666104/4‑93 (MQ=255)

                       gTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAg                                                                 >  2:8822338/1‑93 (MQ=255)

                        ttCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGc                                                                >  2:8707875/1‑93 (MQ=255)

                              gagaATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGc                                                          >  2:25973/1‑93 (MQ=255)

                                gaATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCg                                                        >  2:1800424/1‑93 (MQ=255)

                                 aaTAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGt                                                       <  1:5127985/93‑1 (MQ=255)

                                 aaTAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGt                                                       >  2:526090/1‑93 (MQ=255)

                                 aaTAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGt                                                       >  1:2104860/1‑93 (MQ=255)

                                    aaGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCg                                                    <  1:4325141/93‑1 (MQ=255)

                                    aaGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCg                                                    >  1:6173871/1‑93 (MQ=255)

                                    aaGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCg                                                    <  2:8146903/93‑1 (MQ=255)

                                     aGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGt                                                   >  1:7364166/1‑93 (MQ=255)

                                     aGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGt                                                   >  2:439595/1‑93 (MQ=255)

                                     aGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGt                                                   >  2:3239987/1‑93 (MQ=255)

                                      gCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCCGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGtt                                                  >  2:2363855/1‑93 (MQ=255)

                                        agcGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTa                                                <  1:3873800/91‑1 (MQ=255)

                                        aaCGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTa                                                <  2:1216734/93‑1 (MQ=255)

                                         aCGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTAt                                               >  1:2699643/1‑93 (MQ=255)

                                         aCGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTAt                                               >  2:3119608/1‑93 (MQ=255)

                                         aCGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTAt                                               >  1:291965/1‑93 (MQ=255)

                                           gTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCa                                             >  1:1595311/1‑93 (MQ=255)

                                            tCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCaa                                            <  1:1418329/93‑1 (MQ=255)

                                             cGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAAt                                           >  2:661674/1‑93 (MQ=255)

                                               gATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGc                                         <  2:3279584/93‑1 (MQ=255)

                                                    ttCACCAGTGCGCGGACAATCAGCGACAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGGCgaga                                    <  1:2226235/93‑1 (MQ=255)

                                                     tCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGaa                                   <  2:4792958/93‑1 (MQ=255)

                                                      cACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAat                                  >  2:4847051/1‑93 (MQ=255)

                                                      cACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAat                                  >  2:6151079/1‑93 (MQ=255)

                                                       aCCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAata                                 >  1:5622231/1‑93 (MQ=255)

                                                       aCCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAata                                 <  1:4753678/93‑1 (MQ=255)

                                                        ccAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAatat                                >  1:8035660/1‑93 (MQ=255)

                                                        ccAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAatat                                >  1:858256/1‑93 (MQ=255)

                                                         cAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATc                               <  2:5609958/93‑1 (MQ=255)

                                                          aGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATcc                              <  2:8345584/93‑1 (MQ=255)

                                                           gTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCa                             >  2:4021620/1‑93 (MQ=255)

                                                           gTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCa                             >  2:7388625/1‑93 (MQ=255)

                                                              cgcgGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGcc                          <  2:8636092/93‑1 (MQ=255)

                                                              cgaagACATTCAGCGCCAGACGCTGCGGCCCCTAGGAAAGACTATGGAACGGAACGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGcc                          <  2:3381647/89‑1 (MQ=255)

                                                               gcgGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCa                         >  2:3995239/1‑93 (MQ=255)

                                                                cgGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCAc                        <  1:7960527/93‑1 (MQ=255)

                                                                cgGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCAc                        >  1:8443418/1‑93 (MQ=255)

                                                                 ggACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCAct                       <  1:5542111/93‑1 (MQ=255)

                                                                  gACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCActg                      >  2:7008113/1‑93 (MQ=255)

                                                                  gACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCActg                      >  1:4087037/1‑93 (MQ=255)

                                                                  gACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCActg                      >  2:1338603/1‑93 (MQ=255)

                                                                   gcAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCActgc                     <  1:2814042/92‑1 (MQ=255)

                                                                     aaTCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCActgctg                   >  1:1124369/1‑93 (MQ=255)

                                                                     aaTCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCActgctg                   >  2:5903092/1‑93 (MQ=255)

                                                                     aaTCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCActgctg                   >  2:3591055/1‑93 (MQ=255)

                                                                      aTCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGGGTTTATCAATGCCGAGAATATCCAGCCACtgctgc                  >  1:2950078/1‑93 (MQ=255)

                                                                      aTCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACtgctgc                  >  1:6748815/1‑93 (MQ=255)

                                                                       tCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCa                 <  2:7642246/93‑1 (MQ=255)

                                                                       tCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGGGCGGTTATCAATGCCGGGAATATCCAGCCACTGCTGCa                 >  2:6598307/1‑93 (MQ=255)

                                                                         aGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCAcc               <  2:537955/93‑1 (MQ=255)

                                                                             ccAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGtt           <  2:6194875/93‑1 (MQ=255)

                                                                             ccAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGtt           >  1:966831/1‑93 (MQ=255)

                                                                               aGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGtttt         >  1:1929638/1‑93 (MQ=255)

                                                                               aGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGtttt         >  1:3230477/1‑93 (MQ=255)

                                                                               aGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGtttt         >  1:7741543/1‑93 (MQ=255)

                                                                                gACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGttttt        <  2:8237244/93‑1 (MQ=255)

                                                                                gACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGttttt        >  2:3024074/1‑93 (MQ=255)

                                                                                gACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGttttt        >  1:5412465/1‑93 (MQ=255)

                                                                                 aCGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTg       <  2:797290/93‑1 (MQ=255)

                                                                                  cGCTGCTGCCCCCAGGAAAGACTATGGAGCGGAGCGTCAGCCTTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGc      <  1:5818818/93‑1 (MQ=255)

                                                                                  cGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGc      >  1:8548885/1‑93 (MQ=255)

                                                                                   gctgctGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCt     <  1:4504772/93‑1 (MQ=255)

                                                                                   gctgctGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCt     >  1:6027481/1‑93 (MQ=255)

                                                                                    ctgctgCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCCCCAGTTTTTGCTg    <  2:6483478/93‑1 (MQ=255)

                                                                                    ctgctgCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCTg    <  1:6049613/93‑1 (MQ=255)

                                                                                      gctgCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCTGGc  <  1:6978675/93‑1 (MQ=255)

CAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCTGGC  >  NC_000913/791707‑791885

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