BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,890,225	C→T	R345R (CGG→CGA)	b1808 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,890,225	0	C	T	99.5%	731.2 / ‑5.5	202	R345R (CGG→CGA)	b1808	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (113/88); total (114/88)
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.75e-01

Only 100 of 205 total aligned reads displayed.

GAGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGCCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTC  >  NC_000913/1890135‑1890314

gAGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCg                                                                                         >  1:5776347/1‑93 (MQ=255)

 aGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCgt                                                                                        <  2:7587318/93‑1 (MQ=255)

 aGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCgt                                                                                        >  2:570980/1‑93 (MQ=255)

 aGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCgt                                                                                        <  1:7552995/93‑1 (MQ=255)

 aGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGCCgt                                                                                        >  1:378042/1‑93 (MQ=255)

  gcggcgTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCgtg                                                                                       >  2:404390/1‑93 (MQ=255)

  gcggcgTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCgtg                                                                                       >  2:8724111/1‑93 (MQ=255)

   cggcgTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCgtgt                                                                                      <  2:7340268/93‑1 (MQ=255)

    ggcgTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGtgtg                                                                                     <  1:2267239/93‑1 (MQ=255)

       gTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGa                                                                                  >  2:4755469/1‑93 (MQ=255)

        tGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGAt                                                                                 >  1:9116621/1‑93 (MQ=255)

        tGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGAt                                                                                 >  1:6891262/1‑93 (MQ=255)

         gagaGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATa                                                                                <  1:7288458/93‑1 (MQ=255)

          agagCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATAc                                                                               >  1:4460296/1‑93 (MQ=255)

           gagCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACa                                                                              <  1:806768/93‑1 (MQ=255)

            agCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACaa                                                                             <  1:8452748/93‑1 (MQ=255)

            agCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACaa                                                                             <  2:3937215/93‑1 (MQ=255)

             gCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAAc                                                                            <  2:1492702/93‑1 (MQ=255)

              ccAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACg                                                                           >  2:1661392/1‑93 (MQ=255)

               cAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTCCAGTCGTGTGCGATACAACGt                                                                          <  2:6448126/93‑1 (MQ=255)

                agagTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTg                                                                         <  2:6158506/93‑1 (MQ=255)

                agagTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTg                                                                         <  2:852753/93‑1 (MQ=255)

                  agTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGgc                                                                       >  2:8260482/1‑93 (MQ=255)

                    tAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGgcgc                                                                     >  2:2085864/1‑93 (MQ=255)

                        aTGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCta                                                                 >  1:2656873/1‑92 (MQ=255)

                        aTGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTc                                                                 <  1:9578406/93‑1 (MQ=255)

                          ggCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCaa                                                               >  1:4014159/1‑93 (MQ=255)

                           gcgcGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCaaa                                                              <  1:3058919/93‑1 (MQ=255)

                           gcgcGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCaaa                                                              <  1:717241/93‑1 (MQ=255)

                           gcgcGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCaaa                                                              <  2:6647653/93‑1 (MQ=255)

                           gcgcGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCaaa                                                              <  2:4961810/93‑1 (MQ=255)

                           gcgcGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCaaa                                                              <  2:7379966/93‑1 (MQ=255)

                            cgcgAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAAt                                                             <  2:3820219/93‑1 (MQ=255)

                               gAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGcc                                                          >  1:9346212/1‑93 (MQ=255)

                                  gTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCa                                                       <  1:5966714/93‑1 (MQ=255)

                                   tGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTTTGCGATACAACGTGGCGCGCTCAAATGCCGCAt                                                      <  2:3650967/93‑1 (MQ=255)

                                    gCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATc                                                     >  1:9631697/1‑93 (MQ=255)

                                    gCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATc                                                     <  2:3005242/93‑1 (MQ=255)

                                    gCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATc                                                     <  2:9342144/93‑1 (MQ=255)

                                    gCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATc                                                     >  1:5558101/1‑93 (MQ=255)

                                       ttCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAg                                                  >  2:7009157/1‑93 (MQ=255)

                                         cGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCa                                                <  2:4358212/93‑1 (MQ=255)

                                         cGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCa                                                <  2:9899023/93‑1 (MQ=255)

                                         cGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCa                                                >  2:5589815/1‑93 (MQ=255)

                                            aCCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAgg                                             >  1:9776467/1‑93 (MQ=255)

                                            aCCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAgg                                             >  1:5450348/1‑93 (MQ=255)

                                             ccAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAgcc                                            >  1:9242203/1‑91 (MQ=255)

                                             ccAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGc                                            <  1:3583627/93‑1 (MQ=255)

                                              cAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCg                                           >  1:7439613/1‑93 (MQ=255)

                                              cAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCg                                           >  1:4255063/1‑93 (MQ=255)

                                               aGTAGCTGTAGCCCGGCTGATTGATCTCTTTTCGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGGTCAGATGCCGCATCCTGCATGGccg                                          >  2:7038345/1‑91 (MQ=255)

                                               aGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCgg                                          >  1:1683057/1‑93 (MQ=255)

                                               aGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCgg                                          >  2:7038343/1‑93 (MQ=255)

                                                 tAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGaa                                        <  2:1115207/93‑1 (MQ=255)

                                                 tAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGaa                                        <  1:1204015/93‑1 (MQ=255)

                                                 tAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGaa                                        >  1:7944362/1‑93 (MQ=255)

                                                 tAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGCGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGaa                                        <  1:9873459/93‑1 (MQ=255)

                                                 gatcTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGaa                                        <  1:5868897/90‑1 (MQ=255)

                                                      gTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTTCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGccc                                   >  1:5384769/1‑93 (MQ=255)

                                                      gTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGccc                                   >  2:2403205/1‑93 (MQ=255)

                                                      gTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGccc                                   <  2:7231543/93‑1 (MQ=255)

                                                      gTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGccc                                   <  1:445134/93‑1 (MQ=255)

                                                      gTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGccc                                   <  2:6118344/93‑1 (MQ=255)

                                                      gTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGccc                                   >  1:8045834/1‑93 (MQ=255)

                                                       tAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGCTACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGcccc                                  >  1:193884/1‑93 (MQ=255)

                                                       tAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGcccc                                  >  2:3976936/1‑93 (MQ=255)

                                                       tAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGcccc                                  >  1:2506909/1‑93 (MQ=255)

                                                        aGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCa                                 <  1:471348/93‑1 (MQ=255)

                                                         gCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAg                                >  1:891976/1‑93 (MQ=255)

                                                         gCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAg                                >  1:6542592/1‑93 (MQ=255)

                                                           ccGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTg                              <  1:8617629/93‑1 (MQ=255)

                                                           ccGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTg                              >  1:9474726/1‑93 (MQ=255)

                                                            cGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGa                             >  2:3048346/1‑93 (MQ=255)

                                                            cGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGa                             <  2:7731875/93‑1 (MQ=255)

                                                            cGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGa                             >  2:4834483/1‑93 (MQ=255)

                                                              gCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACa                           >  2:9714774/1‑93 (MQ=255)

                                                              gCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACa                           <  1:79089/93‑1 (MQ=255)

                                                              gCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACa                           >  1:6161279/1‑93 (MQ=255)

                                                                tgattgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGt                         >  2:4222403/1‑93 (MQ=255)

                                                                 gattgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGtt                        >  2:3344290/1‑93 (MQ=255)

                                                                  attgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGttt                       >  2:2627852/1‑93 (MQ=255)

                                                                  attgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGttt                       <  1:3752859/93‑1 (MQ=255)

                                                                     gatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGc                    >  2:4549465/1‑93 (MQ=255)

                                                                      atCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGcc                   >  2:5010265/1‑93 (MQ=255)

                                                                      atCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGcc                   >  1:3933488/1‑93 (MQ=255)

                                                                       tctcTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCa                  <  1:1099454/93‑1 (MQ=255)

                                                                        ctctTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCAc                 <  2:5724801/93‑1 (MQ=255)

                                                                        ctctTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCAc                 >  2:2810538/1‑93 (MQ=255)

                                                                        ctctTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCAc                 <  1:3334560/93‑1 (MQ=255)

                                                                         tctTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACg                <  2:512795/93‑1 (MQ=255)

                                                                             ttAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGGTTCGCCACGTCAt            >  2:160421/1‑93 (MQ=255)

                                                                               aGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATaa          >  2:8849336/1‑93 (MQ=255)

                                                                               aGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATaa          <  1:7656669/93‑1 (MQ=255)

                                                                                gCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAAc         >  1:3471655/1‑93 (MQ=255)

                                                                                gCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAAc         >  1:2435324/1‑93 (MQ=255)

                                                                                gCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAAc         <  2:5767356/93‑1 (MQ=255)

                                                                                 ccGCTCCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGACCGCCCCAGTGACAGTTTCGCCACGTCATAACa        <  2:8830015/93‑1 (MQ=255)

                                                                                   gCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACaaa      >  2:3415597/1‑93 (MQ=255)

                                                                                    cTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAAGGCCCCAGTGACAGTTTCGCCACGTTATAACAAAg     <  2:3564714/93‑1 (MQ=255)

                                                                                     ttCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGt    <  1:1584348/93‑1 (MQ=255)

GAGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGCCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTC  >  NC_000913/1890135‑1890314

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: `‑`