BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,407,307	A→T	13.2%	H180L (CAC→CTC)	yjfP →	acyl CoA esterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,407,307	0	A	T	13.2%	1007.3 / 123.5	334	H180L (CAC→CTC)	yjfP	acyl CoA esterase
Reads supporting (aligned to +/- strand): ref base A (147/141); new base T (18/26); total (166/168)
Fisher's exact test for biased strand distribution p-value = 2.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

Only 100 of 342 total aligned reads displayed.

GGTGAGATGTACCGCCAGCATGATGGGATCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTCATGGCAGCCAGGCGTGC  >  CP009273/4407158‑4407456

   gagaTGTACCGCCAGCATGATGGGATCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCt                                                                                                                                                    >  1:2460640/1‑150 (MQ=255)

     gaTGTACCGCCAGCATGATGGGATCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTgg                                                                                                                                                  <  1:254846/150‑1 (MQ=255)

      ttGTACCGCCAGCATGATGGGATCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGa                                                                                                                                                 <  1:2294682/149‑1 (MQ=255)

       tGTACCGGCATCATGTTTGGTTCGGGCCTTTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGaa                                                                                                                                                <  2:4812693/150‑1 (MQ=255)

       tGTACCGCCAGCATGATGGGATCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCTCCTGGaa                                                                                                                                                <  2:1825230/150‑1 (MQ=255)

       tGTACCGCCAGCATGATGGGATCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGaa                                                                                                                                                <  1:2712396/150‑1 (MQ=255)

        gTACCGCCAGCATGATGGGATCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGaac                                                                                                                                               >  2:2885541/1‑150 (MQ=255)

          aCCGCCAGCATGATGGGATCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGaacna                                                                                                                                             <  1:2968371/150‑1 (MQ=255)

           ccGCCAGCATGATGGGATCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGaacaac                                                                                                                                            <  2:2207378/150‑1 (MQ=255)

            cGCCAGCATGATGGGATCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCTCCTGGAACAACt                                                                                                                                           >  2:869102/1‑150 (MQ=255)

            cGCCAGCATGATGGGATCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACt                                                                                                                                           >  2:22581/1‑150 (MQ=255)

             gCCAGCATGATGGGATCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACtt                                                                                                                                          >  1:1972122/1‑150 (MQ=255)

             gCCAGCATGATGGGATCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACtt                                                                                                                                          >  1:1567429/1‑150 (MQ=255)

              ccAGCATGATGGGATCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTa                                                                                                                                         >  1:360326/1‑150 (MQ=255)

                aGCATGATGGGATCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGt                                                                                                                                       >  2:3052764/1‑150 (MQ=255)

                    tgatgGGATCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGaaa                                                                                                                                   <  1:1419081/150‑3 (MQ=255)

                    tgatgGGATCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACa                                                                                                                                   <  2:3543865/150‑1 (MQ=255)

                     gatgGGATCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCTCCTGGAACAACTTAGTGACAg                                                                                                                                  <  1:4804140/150‑1 (MQ=255)

                       tgGGATCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCTCCTGGAACAACTTAGTGACAGAc                                                                                                                                >  1:4886817/1‑150 (MQ=255)

                       tgGGATCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGAc                                                                                                                                <  1:2610033/150‑1 (MQ=255)

                        gggATCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGAcc                                                                                                                               >  2:522887/1‑150 (MQ=255)

                           aTCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACctct                                                                                                                            <  1:4167637/150‑1 (MQ=255)

                              gggCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTAct                                                                                                                         >  1:4391099/1‑150 (MQ=255)

                                gCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTActgc                                                                                                                       >  2:4008959/1‑150 (MQ=255)

                                        aCATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCAt                                                                                                               <  1:2931662/150‑1 (MQ=255)

                                            cACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGcc                                                                                                           <  2:3838716/150‑1 (MQ=255)

                                              cTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTc                                                                                                         >  2:2828450/1‑150 (MQ=255)

                                               tCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTcg                                                                                                        <  1:3074618/150‑1 (MQ=255)

                                                 gCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTcgac                                                                                                      <  1:1099797/150‑1 (MQ=255)

                                                 gCCCGGTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGACGGAATTCCCTAACATTGTCGCGCCCCGGGCAGAGTGGGACGCGACACACCCCCTGGACCCACTGCGTGACCGACCTCTACTGCTGTGGCCTGGCCTcgac                                                                                                      >  2:4533512/1‑150 (MQ=255)

                                                     gTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATg                                                                                                  <  1:4406720/150‑1 (MQ=255)

                                                       tCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGAc                                                                                                <  2:4652496/150‑1 (MQ=255)

                                                       tCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGAc                                                                                                <  2:3395768/150‑1 (MQ=255)

                                                       tCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGAc                                                                                                <  2:4860076/150‑1 (MQ=255)

                                                        cACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCTCCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACg                                                                                               <  2:4734566/150‑1 (MQ=255)

                                                        cACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACg                                                                                               <  2:1502562/150‑1 (MQ=255)

                                                         aCTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACGTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGt                                                                                              >  2:3611000/1‑150 (MQ=255)

                                                           tGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTtg                                                                                            >  2:4377195/1‑150 (MQ=255)

                                                            gTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTtgt                                                                                           <  1:1692906/150‑1 (MQ=255)

                                                                ccACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCTCCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGtgcct                                                                                       <  1:3443030/150‑1 (MQ=255)

                                                                ccACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGATGtgcct                                                                                       >  1:618263/1‑150 (MQ=255)

                                                                   ccGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGtgcctgcc                                                                                    <  1:785103/150‑1 (MQ=255)

                                                                      cTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCcgac                                                                                 >  1:5036491/1‑150 (MQ=255)

                                                                      cTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCcgac                                                                                 >  1:3283196/1‑150 (MQ=255)

                                                                            ccTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCa                                                                           <  1:2816561/150‑1 (MQ=255)

                                                                               gAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCTCCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTa                                                                        >  2:1004709/1‑150 (MQ=255)

                                                                                 aaCGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACg                                                                      <  1:4104087/150‑1 (MQ=255)

                                                                                   cGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGtt                                                                    <  2:3126865/150‑1 (MQ=255)

                                                                                   cGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGtt                                                                    >  1:3755722/1‑150 (MQ=255)

                                                                                      cagcacagcaGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTgc                                                                 <  1:3218886/150‑1 (MQ=255)

                                                                                       agcacagcaGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTgca                                                                >  2:3253817/1‑150 (MQ=255)

                                                                                       agcacagcaGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTgca                                                                >  1:2075799/1‑150 (MQ=255)

                                                                                              cagAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGcc                                                         >  2:3373204/1‑150 (MQ=255)

                                                                                              cagAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGcc                                                         >  2:2791918/1‑150 (MQ=255)

                                                                                              cagAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGcc                                                         >  2:546496/1‑150 (MQ=255)

                                                                                                gaatggatTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCtt                                                       <  1:166803/150‑1 (MQ=255)

                                                                                                gaatgaatTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCtt                                                       <  1:3781391/150‑1 (MQ=255)

                                                                                                 aatgaatTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTa                                                      >  2:1083882/1‑150 (MQ=255)

                                                                                                 aatgaatTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTa                                                      >  2:1420398/1‑150 (MQ=255)

                                                                                                  atgaatTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTaa                                                     <  2:1567429/150‑1 (MQ=255)

                                                                                                   tgaatTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAg                                                    <  2:3934855/150‑1 (MQ=255)

                                                                                                    gaatTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGc                                                   >  1:388730/1‑150 (MQ=255)

                                                                                                     aatTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCg                                                  <  1:2157029/150‑1 (MQ=255)

                                                                                                       ttCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCgag                                                <  2:704433/150‑1 (MQ=255)

                                                                                                       ttCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCgag                                                >  2:4918902/1‑150 (MQ=255)

                                                                                                        tCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCgaga                                               <  2:1387976/150‑1 (MQ=255)

                                                                                                         cAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAgac                                              >  2:3294780/1‑150 (MQ=255)

                                                                                                          aaTAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAgacg                                             >  2:2481137/1‑150 (MQ=255)

                                                                                                            tAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAgacggg                                           >  1:3019181/1‑150 (MQ=255)

                                                                                                              ccagtGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTCCTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAgacgggac                                         <  2:2351567/146‑1 (MQ=255)

                                                                                                               cagtGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAgacgggacg                                        <  2:3711392/147‑1 (MQ=255)

                                                                                                                 ttGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCTCCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGacgggacggg                                      >  1:1006160/1‑150 (MQ=255)

                                                                                                                  tGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTCCTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGAcgggacggga                                     <  2:2527541/150‑1 (MQ=255)

                                                                                                                      gcgcCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATaaa                                 >  1:4527275/1‑150 (MQ=255)

                                                                                                                       cgcCACTGGCAGAGTGGGAAGCGACAAACCACCTGGACCAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGACGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACTGGATAAAc                                <  2:2617796/150‑1 (MQ=255)

                                                                                                                       cgcCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAAc                                >  1:1131625/1‑150 (MQ=255)

                                                                                                                        gcCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACt                               >  1:4729042/1‑150 (MQ=255)

                                                                                                                        gcCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACt                               >  1:511214/1‑150 (MQ=255)

                                                                                                                          cACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGc                             >  1:2334290/1‑150 (MQ=255)

                                                                                                                          cACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGc                             >  1:1581045/1‑150 (MQ=255)

                                                                                                                           aCTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCt                            >  1:4997322/1‑150 (MQ=255)

                                                                                                                           aCTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCt                            >  1:850077/1‑150 (MQ=255)

                                                                                                                            cTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTa                           >  1:10782/1‑150 (MQ=255)

                                                                                                                             tGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTaa                          >  2:2130715/1‑150 (MQ=255)

                                                                                                                             tGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTaa                          >  2:4518769/1‑150 (MQ=255)

                                                                                                                                 agagTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTg                      <  1:2379771/150‑1 (MQ=255)

                                                                                                                                 agagTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTg                      <  1:2135005/150‑1 (MQ=255)

                                                                                                                                    gTGGGAAGCGACAAACCTCCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTc                   <  1:1859587/150‑1 (MQ=255)

                                                                                                                                    gTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTc                   <  2:2031652/150‑1 (MQ=255)

                                                                                                                                     tGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTCa                  <  1:2120181/150‑1 (MQ=255)

                                                                                                                                      gggAAGCGACAAACCTCCTGGAACAACTTAGGGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGCCGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTCAt                 <  2:3798940/150‑1 (MQ=255)

                                                                                                                                          aGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTCATGGCa             >  1:2677964/1‑150 (MQ=255)

                                                                                                                                          aGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTCATGGCa             <  2:1170310/150‑1 (MQ=255)

                                                                                                                                             gACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTCATGGCAGcc          >  2:1890521/1‑150 (MQ=255)

                                                                                                                                              aCAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTCATGGCAGCCa         >  2:4347690/1‑150 (MQ=255)

                                                                                                                                               cAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTCATGGCAGCCAg        <  2:4156274/150‑1 (MQ=255)

                                                                                                                                                aaaCCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTCATGGCAGCCAgg       <  1:4691376/150‑1 (MQ=255)

                                                                                                                                                 aaCCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTCATGGCAGCCAGgc      >  2:1733963/1‑150 (MQ=255)

                                                                                                                                                 aaCCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTCATGGCAGCCAGgc      <  1:1061807/150‑1 (MQ=255)

                                                                                                                                                     accTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTCATGGCAGCCAGgcgtgc  <  1:4008959/150‑1 (MQ=255)

GGTGAGATGTACCGCCAGCATGATGGGATCGGGCTATTTTACATCACTCGCCCGTTCACTGTTTCCACCGCTGATACCTGAAACGGCAGCACAGCAGAATGAATTCAATAACATTGTCGCGCCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTCATGGCAGCCAGGCGTGC  >  CP009273/4407158‑4407456

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