BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	139,326	C→T	G634S (GGC→AGC)	gcd ←	glucose dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	139,326	0	C	T	98.9%	638.7 / NA	188	G634S (GGC→AGC)	gcd	glucose dehydrogenase
Reads supporting (aligned to +/- strand): ref base C (1/1); new base T (87/99); total (88/100)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 191 total aligned reads displayed.

TTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTT  >  NC_000913/139231‑139421

tttACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACCCCATACGGTACACCGTACTGTGGCTAAATGCCGGATTCCGTACCCGTGCtttt                                                                                              <  2:2508455/99‑1 (MQ=255)

tttACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCtttt                                                                                              <  1:2531376/99‑1 (MQ=255)

 ttACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTg                                                                                             >  1:2979304/1‑99 (MQ=255)

 ttACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTg                                                                                             <  2:4501243/99‑1 (MQ=255)

  tACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTgg                                                                                            <  2:4528598/99‑1 (MQ=255)

   aCATGGCAGACCAAATGGTGAGAGGAACGGGTTGGGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGc                                                                                           <  2:1808275/99‑1 (MQ=255)

     aTGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCAt                                                                                         >  2:3657632/1‑99 (MQ=255)

     aTGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCAt                                                                                         >  1:2850499/1‑99 (MQ=255)

      tGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATc                                                                                        >  1:2304760/1‑99 (MQ=255)

       ggCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCt                                                                                       >  1:3640167/1‑99 (MQ=255)

        gccgACCAAAGGGTGAGGGGACCGCGTTGAGCGTGACACCACACGGTACACCGCCCTGTGGCTGAACGTCGGATCCCATACCCGTGCTTTTGGCATCtt                                                                                      >  2:3070792/4‑99 (MQ=255)

         taGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCttt                                                                                     <  1:1805240/98‑1 (MQ=255)

         cAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCttt                                                                                     >  2:4543248/1‑99 (MQ=255)

         cAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCttt                                                                                     <  1:637362/99‑1 (MQ=255)

          aGACCAAATGGTGAGAGGAGCGGGTTGAGCGTGACACCATACGGTTTACCGTACTGTGGCTGAATGCTGGTTTCCGTACTCGTGCTTTTGGCATATTTc                                                                                    <  2:4462106/99‑1 (MQ=255)

          aGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTc                                                                                    <  2:4266994/99‑1 (MQ=255)

          aGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTc                                                                                    <  1:3198180/99‑1 (MQ=255)

          aGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTc                                                                                    >  2:4091918/1‑99 (MQ=255)

           gACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTca                                                                                   >  2:1770830/1‑98 (MQ=255)

           gACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCg                                                                                   <  1:3903013/99‑1 (MQ=255)

              cAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGgct                                                                                >  1:523935/1‑99 (MQ=255)

               aaaTGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGgctg                                                                               >  1:5151351/1‑99 (MQ=255)

                aaTGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGgctgc                                                                              <  2:49538/99‑1 (MQ=255)

                aaTGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGgctgc                                                                              <  1:4371837/99‑1 (MQ=255)

                 aTGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGgctgct                                                                             <  2:781528/99‑1 (MQ=255)

                  tGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTc                                                                            >  2:399299/1‑99 (MQ=255)

                    gTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCa                                                                          >  1:3456100/1‑99 (MQ=255)

                     tGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCAt                                                                         <  1:3753713/99‑1 (MQ=255)

                      gagaGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATc                                                                        <  2:3441679/99‑1 (MQ=255)

                         agGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCggg                                                                     >  1:509865/1‑99 (MQ=255)

                          ggAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGt                                                                    <  1:3699613/99‑1 (MQ=255)

                          ggAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACACGTGCTCTTGGCCTCTTCCGGCCTCTCAACCGGGt                                                                    >  2:3936563/1‑99 (MQ=255)

                           gAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGtt                                                                   <  2:3850001/99‑1 (MQ=255)

                           gAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGtt                                                                   >  2:1602243/1‑99 (MQ=255)

                           gAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGtt                                                                   <  1:2991657/99‑1 (MQ=255)

                            aaCGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTg                                                                  >  1:1648162/1‑99 (MQ=255)

                            aaCGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTg                                                                  >  1:3071333/1‑99 (MQ=255)

                             cgtggTTGAGCGTGGCACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGc                                                                 <  1:411532/96‑1 (MQ=255)

                              cGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGcc                                                                <  1:3574098/99‑1 (MQ=255)

                              cGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGcc                                                                <  2:3609246/99‑1 (MQ=255)

                                ggTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAg                                                              <  1:2063345/99‑1 (MQ=255)

                                ggTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAg                                                              <  1:1498688/99‑1 (MQ=255)

                                 gTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAgg                                                             >  2:849356/1‑99 (MQ=255)

                                 gTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAgg                                                             >  2:3896412/1‑99 (MQ=255)

                                 gTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAgg                                                             >  1:603347/1‑99 (MQ=255)

                                  ttGGGCGTGACACCATACGGTAGACCGTACCGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCGCCATCGGGTTGCCAGGa                                                            <  2:5023774/99‑1 (MQ=255)

                                  ttGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGa                                                            <  2:2947939/99‑1 (MQ=255)

                                  gtGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGa                                                            <  2:2576319/98‑1 (MQ=255)

                                   tGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGAc                                                           <  1:2686041/99‑1 (MQ=255)

                                         tGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCAcgcg                                                     >  2:73537/1‑99 (MQ=255)

                                          gACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCCTCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCgg                                                    >  2:1395508/1‑99 (MQ=255)

                                          gACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCgg                                                    >  1:3769901/1‑99 (MQ=255)

                                             acCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGAt                                                 >  2:1897187/1‑99 (MQ=255)

                                              ccATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTCCCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACAACGCCGGaaa                                                >  2:2726409/1‑97 (MQ=255)

                                              ccATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATc                                                <  2:1995383/99‑1 (MQ=255)

                                              ccATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATc                                                <  2:1856106/99‑1 (MQ=255)

                                               cATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCa                                               >  1:4983540/1‑99 (MQ=255)

                                                  aCGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGATGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGtt                                            >  1:2573421/1‑99 (MQ=255)

                                                   cGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGttt                                           >  1:1315583/1‑99 (MQ=255)

                                                   cGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGttt                                           <  1:717982/99‑1 (MQ=255)

                                                    ggTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTc                                          <  2:1314291/99‑1 (MQ=255)

                                                    ggTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTc                                          <  1:4499845/99‑1 (MQ=255)

                                                    ggTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTc                                          <  1:4776040/99‑1 (MQ=255)

                                                     gTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCg                                         <  1:1994144/99‑1 (MQ=255)

                                                      tACACCGTACTTTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCGAGGACCACGCGGGATCAGTTTCGa                                        <  2:4046692/99‑1 (MQ=255)

                                                       acacCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGaa                                       <  1:1084601/99‑1 (MQ=255)

                                                       acacCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGaa                                       <  2:1564874/99‑1 (MQ=255)

                                                        cacCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGaaa                                      <  1:868596/99‑1 (MQ=255)

                                                           cGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGaaacaa                                   >  1:2741127/1‑99 (MQ=255)

                                                             tACTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGaaacaaac                                 <  2:2794302/99‑1 (MQ=255)

                                                               cTGTGGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACgg                               <  2:3811432/99‑1 (MQ=255)

                                                                  tgGCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAg                            <  1:4614228/99‑1 (MQ=255)

                                                                   ggCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGt                           <  2:210530/99‑1 (MQ=255)

                                                                    gCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTGGGCTGCTTCATCGGGTTGCCAGGACCACGCGGGATCAGTTTTGAAACAAACGGCAGTg                          >  1:4928074/1‑99 (MQ=255)

                                                                    gCTGAATGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTg                          >  2:4637609/1‑99 (MQ=255)

                                                                          tGCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTg                    >  2:3667782/1‑99 (MQ=255)

                                                                           gCCGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCAGCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCAATgg                   <  2:4280359/99‑1 (MQ=255)

                                                                            ccGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTggg                  <  2:1106182/99‑1 (MQ=255)

                                                                            ccGGATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTggg                  >  1:88563/1‑99 (MQ=255)

                                                                              ggATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGACAGTGCCATTGGGtt                <  1:2825725/99‑1 (MQ=255)

                                                                               gATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTg               <  1:3203749/99‑1 (MQ=255)

                                                                               gATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTg               <  2:3303030/99‑1 (MQ=255)

                                                                               gATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTg               <  2:2227279/99‑1 (MQ=255)

                                                                               gATTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTg               <  2:2462112/99‑1 (MQ=255)

                                                                                aTTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTgg              <  2:4820090/99‑1 (MQ=255)

                                                                                aTTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTgg              <  2:898034/99‑1 (MQ=255)

                                                                                aTTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTgg              <  2:3599059/99‑1 (MQ=255)

                                                                                aTTCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTgg              <  1:2196471/99‑1 (MQ=255)

                                                                                 ttCCGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGc             <  1:770028/99‑1 (MQ=255)

                                                                                  gccGCACCCGTGCCTTCGGCATCTTCCGTCTTCTCCATCGGGTTGCCAGTACCACGCGGGATCAGGTTCGAAACAAACGGCAGTGCCATTGGGTTTGCa            <  1:4146241/98‑1 (MQ=255)

                                                                                   ccGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCaa           <  1:4846989/99‑1 (MQ=255)

                                                                                   ccGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCaa           <  1:4443320/99‑1 (MQ=255)

                                                                                    cGTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAt          >  1:3499591/1‑99 (MQ=255)

                                                                                     gTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATc         <  1:334041/99‑1 (MQ=255)

                                                                                     gTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATc         >  1:2079444/1‑99 (MQ=255)

                                                                                     gTACCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATc         <  1:1011640/99‑1 (MQ=255)

                                                                                       aCCCGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGc       <  1:4694090/99‑1 (MQ=255)

                                                                                          cGTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCAc    >  1:1131740/1‑99 (MQ=255)

                                                                                           gTGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACt   >  2:1021888/1‑99 (MQ=255)

                                                                                            tGCTTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACtt  >  1:1721346/1‑99 (MQ=255)

TTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTT  >  NC_000913/139231‑139421

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