BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,253,781	T→C	E171G (GAG→GGG)	yhaNM ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,253,781	0	T	C	99.4%	475.1 / ‑4.8	153	E171G (GAG→GGG)	yhaNM	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (70/82); total (70/83)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

Only 100 of 164 total aligned reads displayed.

ATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACAC  >  NC_000913/3253690‑3253870

aTCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCccct                                                                                         <  2:4386748/93‑1 (MQ=255)

 tCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTg                                                                                        >  2:707107/1‑93 (MQ=255)

  cgccgcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGc                                                                                       >  2:3729421/1‑93 (MQ=255)

   gccgcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGcc                                                                                      >  1:3329270/1‑93 (MQ=255)

   gccgcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGcc                                                                                      >  2:2936/1‑93 (MQ=255)

    ccgcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCca                                                                                     >  1:3048813/1‑93 (MQ=255)

    ccgcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCca                                                                                     <  1:4156111/93‑1 (MQ=255)

    ccgcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCca                                                                                     <  2:3987867/93‑1 (MQ=255)

    ccgcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCca                                                                                     <  1:2185175/93‑1 (MQ=255)

       cAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGTTCCCTTCcacac                                                                                  >  2:3865396/1‑93 (MQ=255)

       cAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCcacac                                                                                  >  1:4326720/1‑93 (MQ=255)

       cAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCcacac                                                                                  <  1:2243885/93‑1 (MQ=255)

        aaaCGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCcacaca                                                                                 <  2:976136/93‑1 (MQ=255)

         accGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCCCCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCacacac                                                                                <  1:1911347/91‑1 (MQ=255)

         aaCGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCacacac                                                                                <  2:1160772/93‑1 (MQ=255)

         aaCGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCacacac                                                                                >  2:2971096/1‑93 (MQ=255)

         aaCGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCacacac                                                                                >  1:2425122/1‑93 (MQ=255)

            gggACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGcc                                                                             <  1:3357827/93‑1 (MQ=255)

              gACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGGCACACACGCctg                                                                           >  2:2224465/1‑93 (MQ=255)

               aCTGCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCctgc                                                                          <  1:99056/93‑1 (MQ=255)

                cTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCctgct                                                                         >  2:4364004/1‑93 (MQ=255)

                cTTCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCctgct                                                                         >  1:2722896/1‑93 (MQ=255)

                  tCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCCCCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTgg                                                                       <  2:4186229/93‑1 (MQ=255)

                  tCATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTgg                                                                       >  1:1704000/1‑93 (MQ=255)

                   cATTGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTggg                                                                      >  1:4594485/1‑93 (MQ=255)

                     ttGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTa                                                                    >  1:3815020/1‑93 (MQ=255)

                      tGACGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTaa                                                                   <  2:4653281/93‑1 (MQ=255)

                       taCGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTaaa                                                                  <  1:1045294/92‑1 (MQ=255)

                       taCGAACTTCAGGATCTCAACCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTaaa                                                                  <  1:3033180/92‑1 (MQ=255)

                         cGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAaca                                                                >  2:4559086/1‑93 (MQ=255)

                         cGAACTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAaca                                                                >  1:1496086/1‑93 (MQ=255)

                           accTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAATGCAGCGTCAGCGGAGATTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAacacc                                                              <  2:1443562/91‑1 (MQ=255)

                            aCTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAacacca                                                             <  2:4422/93‑1 (MQ=255)

                             cTTCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAacaccac                                                            <  2:3836670/93‑1 (MQ=255)

                               tCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAacaccacac                                                          >  1:1630661/1‑93 (MQ=255)

                               tCAGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAacaccacac                                                          >  2:384129/1‑93 (MQ=255)

                                 aGGATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAAcaccacacca                                                        <  2:1901215/93‑1 (MQ=255)

                                  ggATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCAt                                                       >  1:3939645/1‑93 (MQ=255)

                                  ggATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCAt                                                       >  1:658740/1‑93 (MQ=255)

                                  ggATCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCAt                                                       <  1:2374358/93‑1 (MQ=255)

                                   gATCTCAGCCAGCGTC‑GTTCTGGAAATCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATc                                                      <  1:1408047/93‑1 (MQ=255)

                                    aTCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCg                                                     >  1:3210175/1‑93 (MQ=255)

                                    aTCTCAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCg                                                     <  2:3096878/93‑1 (MQ=255)

                                      ctcAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGACCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTg                                                   >  1:4497050/1‑93 (MQ=255)

                                       tcAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGc                                                  <  2:3573152/93‑1 (MQ=255)

                                       tcAGCCAGCGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGc                                                  <  1:1177722/93‑1 (MQ=255)

                                         agccagcGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGt                                                >  2:608163/1‑93 (MQ=255)

                                         agccagcGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGt                                                >  1:1718483/1‑93 (MQ=255)

                                          gccagcGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGtc                                               <  1:4232207/93‑1 (MQ=255)

                                          gccagcGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGtc                                               >  1:2076115/1‑93 (MQ=255)

                                          gccagcGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGtc                                               <  2:1325881/93‑1 (MQ=255)

                                           ccagcGTC‑GTGCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGtct                                              <  2:307149/93‑1 (MQ=255)

                                            cagcGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGtctc                                             >  1:1396540/1‑93 (MQ=255)

                                              gcGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGa                                           >  2:1475465/1‑93 (MQ=255)

                                              gcGTC‑GTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGa                                           <  2:2830369/93‑1 (MQ=255)

                                              gcGTC‑GTTCGGGAAAGCACCGTCAGCGTAGACTCTGGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGa                                           <  1:4025055/93‑1 (MQ=255)

                                                  tcgTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGatat                                        <  2:2408878/91‑1 (MQ=255)

                                                  tcgTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGatat                                        <  1:1016117/91‑1 (MQ=255)

                                                  tcgTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGatat                                        <  2:3121769/91‑1 (MQ=255)

                                                  tcgTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGatat                                        >  2:3932612/3‑93 (MQ=255)

                                                   cgTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATg                                       >  2:1913077/2‑93 (MQ=255)

                                                    gtgcTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGc                                      <  2:1843933/90‑1 (MQ=255)

                                                    gTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGc                                      >  1:603800/1‑93 (MQ=255)

                                                    gTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGc                                      >  2:3288273/1‑93 (MQ=255)

                                                     ttCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGca                                     >  2:491144/1‑93 (MQ=255)

                                                            aaGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGtt                              >  2:3218912/1‑93 (MQ=255)

                                                            aaGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGtt                              <  1:3513203/93‑1 (MQ=255)

                                                            aaGCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGtt                              <  2:3685974/93‑1 (MQ=255)

                                                              gCACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTgg                            <  1:3253277/93‑1 (MQ=255)

                                                               cACCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATggtggt                           >  1:4369287/1‑93 (MQ=255)

                                                                aCCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTa                          <  1:2787369/93‑1 (MQ=255)

                                                                aCCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTa                          <  2:1273893/93‑1 (MQ=255)

                                                                 ccGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTAt                         >  1:2643919/1‑93 (MQ=255)

                                                                   gTCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGc                       <  2:1440233/93‑1 (MQ=255)

                                                                    tCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGcc                      <  1:4168165/93‑1 (MQ=255)

                                                                    tCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGcc                      >  2:57695/1‑93 (MQ=255)

                                                                    tCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGcc                      >  1:2994875/1‑93 (MQ=255)

                                                                    tCAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGcc                      >  1:790379/1‑93 (MQ=255)

                                                                     cAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGccc                     <  1:1613879/93‑1 (MQ=255)

                                                                     cAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGccc                     <  2:4317108/93‑1 (MQ=255)

                                                                     cAGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGccc                     <  2:4639055/93‑1 (MQ=255)

                                                                     cAGCGAAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATCTGCACAATGTTGGTATGccc                     <  2:2666969/93‑1 (MQ=255)

                                                                      aGCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCccg                    >  2:4046220/1‑93 (MQ=255)

                                                                       gCGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCccgc                   <  2:739207/93‑1 (MQ=255)

                                                                        cGGAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACCCCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCccgcc                  <  2:2542568/93‑1 (MQ=255)

                                                                         ggAGACTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCccgccg                 >  2:2629828/1‑93 (MQ=255)

                                                                           ggaCTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCcgac               <  2:3030841/92‑1 (MQ=255)

                                                                            gaCTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCcgacg              <  2:564555/93‑1 (MQ=255)

                                                                             aCTCTTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCcgacga             >  2:4358672/1‑93 (MQ=255)

                                                                              ctctTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGAt            <  1:1751754/93‑1 (MQ=255)

                                                                              ctctTGCTAGCCCCCTGCCACACACGCCTGCTGAGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGAt            <  2:1751226/93‑1 (MQ=255)

                                                                               tctTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATg           >  1:3097081/1‑93 (MQ=255)

                                                                               tctTGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGACGACGATg           >  1:2027789/1‑93 (MQ=255)

                                                                                 ttGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGt         >  2:4218550/1‑93 (MQ=255)

                                                                                 ttGCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGt         >  1:1553224/1‑93 (MQ=255)

                                                                                   gCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGa       <  1:1202970/93‑1 (MQ=255)

                                                                                   gCTCGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGa       <  1:2581119/93‑1 (MQ=255)

                                                                                      cGCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGacac    >  1:3623277/1‑93 (MQ=255)

                                                                                       gCCCCCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGacaca   <  1:1756936/93‑1 (MQ=255)

                                                                                        ccccCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGacacac  <  2:1725487/93‑1 (MQ=255)

ATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACAC  >  NC_000913/3253690‑3253870

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