BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,552,182	G→C	6.3%	G183R (GGC→CGC)	uxuA →	D‑mannonate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,552,182	0	G	C	6.3%	597.1 / 21.9	175	G183R (GGC→CGC)	uxuA	D‑mannonate dehydratase
Reads supporting (aligned to +/- strand): ref base G (76/88); new base C (3/8); total (79/96)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.50e-01

Only 100 of 183 total aligned reads displayed.

TCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTT  >  NC_000913/4552033‑4552331

tCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCgg                                                                                                                                                       <  1:1312024/150‑1 (MQ=255)

tCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCgg                                                                                                                                                       <  2:2063146/150‑1 (MQ=255)

 cGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCggg                                                                                                                                                      >  1:2203693/1‑150 (MQ=255)

 cGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCggg                                                                                                                                                      >  2:1797926/1‑150 (MQ=255)

  gCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTATACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGgc                                                                                                                                                     <  2:693234/150‑1 (MQ=255)

         tCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGgaag                                                                                                                                              >  2:571609/1‑150 (MQ=255)

          cGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGgaaga                                                                                                                                             >  1:435082/1‑150 (MQ=255)

           gAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGGCAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGgaagaa                                                                                                                                            >  2:300477/1‑150 (MQ=255)

           gAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGCGCCCGgaagac                                                                                                                                            <  2:1002868/150‑2 (MQ=255)

           gAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGCGCCCGgaagac                                                                                                                                            <  2:2183228/150‑2 (MQ=255)

            aaaTGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGCGCCCGgaagacc                                                                                                                                           <  1:2357235/150‑3 (MQ=255)

            aaaTGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCAATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGCGCCCGgaagacc                                                                                                                                           <  2:1015706/150‑3 (MQ=255)

             aaTGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAgg                                                                                                                                          >  1:32597/1‑150 (MQ=255)

             aaTGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAgg                                                                                                                                          <  1:229699/150‑1 (MQ=255)

             aaTGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAgg                                                                                                                                          <  2:27403/150‑1 (MQ=255)

              aTGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGc                                                                                                                                         <  2:2178958/150‑1 (MQ=255)

              aTGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGCGCCCGgaagaccag                                                                                                                                         <  1:272388/150‑5 (MQ=255)

              aTGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGCGCCCGgaagaccag                                                                                                                                         <  2:1311070/150‑5 (MQ=255)

               tGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCt                                                                                                                                        >  1:932675/1‑150 (MQ=255)

                    atCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTACCATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTacaac                                                                                                                                   >  1:2466468/1‑148 (MQ=255)

                      ccTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCTGAAGAAGGCTACACCCt                                                                                                                                 >  2:1546010/1‑150 (MQ=255)

                      ccTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCt                                                                                                                                 >  1:2066524/1‑150 (MQ=255)

                       cTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTc                                                                                                                                <  2:409214/150‑1 (MQ=255)

                          aaaCGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGAc                                                                                                                             <  1:2283168/150‑1 (MQ=255)

                              gCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGt                                                                                                                         <  2:2293651/150‑1 (MQ=255)

                                  aGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCg                                                                                                                     <  2:2283196/150‑1 (MQ=255)

                                  aGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCg                                                                                                                     <  2:715795/150‑1 (MQ=255)

                                   ggCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGt                                                                                                                    <  2:435082/150‑1 (MQ=255)

                                    gcgcGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTa                                                                                                                   <  1:1668274/150‑1 (MQ=255)

                                       cgGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAac                                                                                                                >  2:2342626/1‑150 (MQ=255)

                                       cgGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAac                                                                                                                >  1:1004611/1‑150 (MQ=255)

                                       cgGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAac                                                                                                                <  2:1737699/150‑1 (MQ=255)

                                       cgGAAGCGGATTACACCGAACAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAac                                                                                                                >  2:597987/1‑150 (MQ=255)

                                        ggAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAaca                                                                                                               <  2:2440634/150‑1 (MQ=255)

                                        ggAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAaca                                                                                                               >  1:646331/1‑150 (MQ=255)

                                        ggAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAaca                                                                                                               <  1:713180/150‑1 (MQ=255)

                                        ggAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAaca                                                                                                               >  2:857941/1‑150 (MQ=255)

                                               gATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAg                                                                                                        >  1:2258483/1‑150 (MQ=255)

                                                aTTAGACCGAAGAAGAAATTGCTCAGGCAGCGGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGc                                                                                                       <  2:2266411/150‑1 (MQ=255)

                                                aTTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGC‑GAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCt                                                                                                      <  2:1019911/150‑1 (MQ=255)

                                                   acacCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGt                                                                                                    >  2:639953/1‑150 (MQ=255)

                                                    cacCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTa                                                                                                   <  1:1641877/150‑1 (MQ=255)

                                                       cGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACaa                                                                                                >  2:1418708/1‑150 (MQ=255)

                                                         aagaagaaATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAg                                                                                              >  1:1294953/1‑150 (MQ=255)

                                                         aagaagaaATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAg                                                                                              <  2:2228885/150‑1 (MQ=255)

                                                          agaagaaATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGa                                                                                             <  1:2426589/150‑1 (MQ=255)

                                                          agaagaaATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGa                                                                                             >  2:1199316/1‑150 (MQ=255)

                                                          agaagaaATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGa                                                                                             <  1:778992/150‑1 (MQ=255)

                                                           gaagaaATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGat                                                                                            >  2:1312791/1‑150 (MQ=255)

                                                             agaaATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGatat                                                                                          >  2:1294728/1‑150 (MQ=255)

                                                                       cAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCg                                                                                >  2:952019/1‑150 (MQ=255)

                                                                       cAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGAGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCg                                                                                >  1:1916656/1‑150 (MQ=255)

                                                                              cTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTgc                                                                         >  1:1058516/1‑150 (MQ=255)

                                                                                 aaCGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGcgcg                                                                      <  1:2350804/150‑1 (MQ=255)

                                                                                   cGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCTGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGaa                                                                    >  1:235926/1‑150 (MQ=255)

                                                                                    gCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGaaa                                                                   <  1:466858/150‑1 (MQ=255)

                                                                                     cTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGaaaa                                                                  >  1:2230428/1‑150 (MQ=255)

                                                                                     cTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGaaaa                                                                  <  1:1979691/150‑1 (MQ=255)

                                                                                      ttCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAAc                                                                 <  2:2228567/150‑1 (MQ=255)

                                                                                         gCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACttt                                                              >  2:535692/1‑150 (MQ=255)

                                                                                         gCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACttt                                                              >  1:140723/1‑150 (MQ=255)

                                                                                           cACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGc                                                            <  1:856357/150‑1 (MQ=255)

                                                                                           cACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCATAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGc                                                            <  2:767622/150‑1 (MQ=255)

                                                                                                tGAGCGATGAAGACAAAGCGCGTCTGACCCATAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCt                                                       >  2:334021/1‑150 (MQ=255)

                                                                                                   gCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTcc                                                    <  1:920915/150‑1 (MQ=255)

                                                                                                     gATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTg                                                  <  1:1460480/150‑1 (MQ=255)

                                                                                                      aTGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCCGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGa                                                 >  2:2099504/1‑150 (MQ=255)

                                                                                                       tGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGaa                                                >  1:694408/1‑150 (MQ=255)

                                                                                                           gACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCg                                            >  2:81957/1‑150 (MQ=255)

                                                                                                            aCAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGa                                           >  2:2076521/1‑150 (MQ=255)

                                                                                                              aaaGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCTAAAACTTTGCCGTCTTCCTGAAAGCGatt                                         <  2:224032/150‑1 (MQ=255)

                                                                                                              aaaGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGatt                                         >  1:159256/1‑150 (MQ=255)

                                                                                                              aaaGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGatt                                         <  1:2148717/150‑1 (MQ=255)

                                                                                                               aaGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGatta                                        <  2:804373/150‑1 (MQ=255)

                                                                                                                  cgcgTCTTACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTc                                     <  2:1596163/150‑1 (MQ=255)

                                                                                                                  cgcgTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACACTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTc                                     >  1:593048/1‑150 (MQ=255)

                                                                                                                   gcgTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTcc                                    <  1:1797926/150‑1 (MQ=255)

                                                                                                                   gcgTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTcc                                    <  1:639953/150‑1 (MQ=255)

                                                                                                                    cgTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCa                                   >  1:1560447/1‑150 (MQ=255)

                                                                                                                         gACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGc                              >  2:1769328/1‑150 (MQ=255)

                                                                                                                          aCCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCt                             >  2:1995568/1‑150 (MQ=255)

                                                                                                                          aCCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCt                             <  2:657679/150‑1 (MQ=255)

                                                                                                                           cccGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTg                            <  1:1900664/150‑1 (MQ=255)

                                                                                                                           cccGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTg                            >  2:539793/1‑150 (MQ=255)

                                                                                                                           cccGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACAAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTg                            >  1:25212/1‑150 (MQ=255)

                                                                                                                             cGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTgaa                          >  1:1582585/1‑150 (MQ=255)

                                                                                                                              gTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTgaag                         >  1:37944/1‑150 (MQ=255)

                                                                                                                               tAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTgaagt                        <  1:101833/150‑2 (MQ=255)

                                                                                                                                  catcatTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGaagaagg                     >  2:1734995/1‑149 (MQ=255)

                                                                                                                                  catcatTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGCAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGt                     >  1:1348907/1‑150 (MQ=255)

                                                                                                                                     catTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCgg                  <  1:1602746/150‑1 (MQ=255)

                                                                                                                                     catTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATT‑CAGTTGCTGAAGAAGTCGgc                 <  1:1199316/150‑1 (MQ=255)

                                                                                                                                           ccttCTTCCGCGCCCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGgcgtgcg            >  1:978116/4‑150 (MQ=255)

                                                                                                                                            cttCTTCCGCGCCCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGgcgtgcgt           >  2:109177/3‑150 (MQ=255)

                                                                                                                                               cTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATTTCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATg        >  2:1487361/1‑150 (MQ=255)

                                                                                                                                               cTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATg        >  1:912895/1‑150 (MQ=255)

                                                                                                                                                 tCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGc      >  2:1004890/1‑150 (MQ=255)

                                                                                                                                                 tCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGc      <  1:1261407/150‑1 (MQ=255)

                                                                                                                                                   cGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTg    >  1:459551/1‑150 (MQ=255)

                                                                                                                                                     ggCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGtt  <  1:2099504/150‑1 (MQ=255)

TCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGTT  >  NC_000913/4552033‑4552331

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`