BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,357,095	A→C	16.5%	V350V (GTT→GTG)	aspA ←	aspartate ammonia‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,357,095	0	A	C	16.5%	582.2 / 145.8	268	V350V (GTT→GTG)	aspA	aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand): ref base A (109/115); new base C (23/21); total (132/136)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

Only 100 of 272 total aligned reads displayed.

ACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAGGCCGGCACGTGGGCCTGAAGAGAG  >  CP009273/4356946‑4357244

acTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTa                                                                                                                                                       >  1:869105/1‑150 (MQ=255)

  ttcttTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAAc                                                                                                                                                     >  1:1264507/1‑150 (MQ=255)

  ttcttTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAAc                                                                                                                                                     <  2:1095199/150‑1 (MQ=255)

     tttGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAgt                                                                                                                                                  <  2:3851018/150‑1 (MQ=255)

      ttGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAgtg                                                                                                                                                 >  1:4261002/1‑150 (MQ=255)

         ttAGCTGTGATGCCGTTAATGCTTTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACTTGGCCTGGCCATTGACCGGCTCCATAACGTTCAACTGCAGCTGTCCTGCTTCTGCTGCCATGGTAACTgtggtg                                                                                                                                              <  1:1270771/150‑1 (MQ=255)

          tAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGgtgt                                                                                                                                             <  2:3011462/150‑1 (MQ=255)

             cAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGt                                                                                                                                          >  2:483231/1‑150 (MQ=255)

              aGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGtt                                                                                                                                         <  2:692306/150‑1 (MQ=255)

              aGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGtt                                                                                                                                         <  2:388351/150‑1 (MQ=255)

                tGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTAc                                                                                                                                       >  2:1644870/1‑150 (MQ=255)

                     ccGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATg                                                                                                                                  <  2:2318763/150‑1 (MQ=255)

                       gTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGAc                                                                                                                                >  1:2710037/1‑150 (MQ=255)

                       gTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGAc                                                                                                                                >  1:4482889/1‑150 (MQ=255)

                        ttAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACt                                                                                                                               <  2:4309010/150‑1 (MQ=255)

                          aaTGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACttt                                                                                                                             <  2:2054253/150‑1 (MQ=255)

                           aTGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTg                                                                                                                            <  1:2741005/150‑1 (MQ=255)

                            tGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGa                                                                                                                           >  2:2462865/1‑150 (MQ=255)

                            tGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGa                                                                                                                           >  2:4476123/1‑150 (MQ=255)

                             gCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGaa                                                                                                                          >  1:1212461/1‑150 (MQ=255)

                             gCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGaa                                                                                                                          >  1:1685637/1‑150 (MQ=255)

                              cATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAg                                                                                                                         <  2:3204329/150‑1 (MQ=255)

                              cATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAg                                                                                                                         <  1:116066/150‑1 (MQ=255)

                               aTTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGc                                                                                                                        <  2:2487984/150‑1 (MQ=255)

                                   ttCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATAc                                                                                                                    >  1:189029/1‑150 (MQ=255)

                                   ttCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATAc                                                                                                                    <  1:1100271/150‑1 (MQ=255)

                                    tCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATAcc                                                                                                                   <  1:2127384/150‑1 (MQ=255)

                                    tCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATAcc                                                                                                                   <  1:2996731/150‑1 (MQ=255)

                                     ccAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCt                                                                                                                  <  1:4051454/150‑1 (MQ=255)

                                     ccAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATANCt                                                                                                                  <  1:3532063/150‑1 (MQ=255)

                                     ccAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCt                                                                                                                  <  1:3172741/150‑1 (MQ=255)

                                       agcagGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTgg                                                                                                                >  1:4198800/1‑150 (MQ=255)

                                       agcagGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTgg                                                                                                                >  1:811439/1‑150 (MQ=255)

                                       agcagGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTgg                                                                                                                <  1:373654/150‑1 (MQ=255)

                                        gcagGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGt                                                                                                               >  2:2391293/1‑150 (MQ=255)

                                        gcagGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGt                                                                                                               <  2:431779/150‑1 (MQ=255)

                                         cagGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAAGGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGtt                                                                                                              <  2:233483/150‑1 (MQ=255)

                                           ggTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTaa                                                                                                            >  2:46839/1‑150 (MQ=255)

                                            gTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAAc                                                                                                           <  1:3354294/150‑1 (MQ=255)

                                              tGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCa                                                                                                         <  2:4100446/150‑1 (MQ=255)

                                               gTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCAc                                                                                                        >  2:3997307/1‑150 (MQ=255)

                                                tAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACt                                                                                                       >  1:362273/1‑150 (MQ=255)

                                                  gCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTc                                                                                                     <  2:1539200/150‑1 (MQ=255)

                                                  gCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTc                                                                                                     >  1:625243/1‑150 (MQ=255)

                                                   cAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTcc                                                                                                    <  2:995579/150‑1 (MQ=255)

                                                       cGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGaa                                                                                                <  2:1272529/150‑1 (MQ=255)

                                                        gTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGaac                                                                                               <  2:1414387/150‑1 (MQ=255)

                                                         ttGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGaaca                                                                                              <  1:2924226/150‑1 (MQ=255)

                                                         ttGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGaaca                                                                                              <  2:4261002/150‑1 (MQ=255)

                                                          tGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGaacaa                                                                                             <  1:2046905/150‑1 (MQ=255)

                                                          tGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGAAACAGTGGTGTCGTTACCGATGACTGTGAAGCATACCTGGTTAACCACTTCCGGaacaa                                                                                             >  1:2401085/1‑150 (MQ=255)

                                                              cAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCgg                                                                                         <  1:1766168/150‑1 (MQ=255)

                                                                 aaTGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGTAGCATACCTGGTTAACCACTTCCGGAACAACCGGGtt                                                                                      <  2:3567987/150‑1 (MQ=255)

                                                                  aTGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGttt                                                                                     >  1:1086855/1‑150 (MQ=255)

                                                                  aTGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGttt                                                                                     <  1:2811226/150‑1 (MQ=255)

                                                                   tgtgAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTa                                                                                    >  1:912602/1‑150 (MQ=255)

                                                                   tgtgAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTa                                                                                    <  2:3222749/150‑1 (MQ=255)

                                                                    gtgAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTAc                                                                                   <  1:1183261/150‑1 (MQ=255)

                                                                      gAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACtt                                                                                 <  1:4124093/150‑1 (MQ=255)

                                                                       aaCGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACttt                                                                                <  2:659345/150‑1 (MQ=255)

                                                                           gATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCt                                                                            <  1:2157729/150‑1 (MQ=255)

                                                                              tCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGc                                                                         <  1:1588702/150‑1 (MQ=255)

                                                                                gAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCat                                                                       >  2:3975823/1‑150 (MQ=255)

                                                                                gAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCat                                                                       <  1:1554579/150‑1 (MQ=255)

                                                                                gAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCat                                                                       <  2:3962630/150‑1 (MQ=255)

                                                                                   cATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCatgat                                                                    <  1:249089/150‑1 (MQ=255)

                                                                                     tggcctggccAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATgg                                                                  >  2:2607220/1‑150 (MQ=255)

                                                                                      ggcctggccAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGa                                                                 <  2:1952876/150‑1 (MQ=255)

                                                                                      ggcctggccAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTGTAGCTGGCATGATGGa                                                                 >  2:2857930/1‑150 (MQ=255)

                                                                                        cctggccAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAag                                                               <  2:283416/150‑1 (MQ=255)

                                                                                         ctggccAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAaga                                                              <  1:2097548/150‑1 (MQ=255)

                                                                                         ctggccAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAaga                                                              >  1:3018862/1‑150 (MQ=255)

                                                                                          tggccAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAagag                                                             <  2:1750983/150‑1 (MQ=255)

                                                                                             ccAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGccc                                                          <  1:2219460/150‑1 (MQ=255)

                                                                                              cAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCg                                                         >  1:1342225/1‑150 (MQ=255)

                                                                                               aaTGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGc                                                        >  2:3503900/1‑150 (MQ=255)

                                                                                                aTGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGcc                                                       <  1:4223819/150‑1 (MQ=255)

                                                                                                aTGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGcc                                                       <  1:537500/150‑1 (MQ=255)

                                                                                                aTGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGcc                                                       >  2:838903/1‑150 (MQ=255)

                                                                                                 tGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCt                                                      >  2:1611316/1‑150 (MQ=255)

                                                                                                   aCCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGc                                                    >  2:1864003/1‑150 (MQ=255)

                                                                                                     cGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAg                                                  >  1:691708/1‑150 (MQ=255)

                                                                                                      ggCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGt                                                 >  2:1953778/1‑150 (MQ=255)

                                                                                                        cTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTc                                               >  1:1161167/1‑150 (MQ=255)

                                                                                                        cTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTc                                               >  1:1568400/1‑150 (MQ=255)

                                                                                                        cTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTc                                               >  1:847161/1‑150 (MQ=255)

                                                                                                          ccATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCg                                             <  1:1461978/150‑1 (MQ=255)

                                                                                                          ccATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCg                                             >  1:2375925/1‑150 (MQ=255)

                                                                                                          ccATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCg                                             >  2:2748195/1‑150 (MQ=255)

                                                                                                                 gTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGtt                                      >  1:4032842/1‑150 (MQ=255)

                                                                                                                    cAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGAt                                   >  1:523140/1‑150 (MQ=255)

                                                                                                                    cAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGAt                                   <  2:3928780/150‑1 (MQ=255)

                                                                                                                         gCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGt                              <  2:1369122/150‑1 (MQ=255)

                                                                                                                         gCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGt                              >  2:728493/1‑150 (MQ=255)

                                                                                                                           aGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTc                            <  1:3872033/150‑1 (MQ=255)

                                                                                                                             cTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAg                          >  1:569033/1‑150 (MQ=255)

                                                                                                                                  cTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAGGCCgg                     >  1:1336972/1‑150 (MQ=255)

                                                                                                                                           ctgcCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAGGCCGGCACGTGGGc            >  1:3331356/1‑150 (MQ=255)

                                                                                                                                                    tAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAGGCCGGCACGTGGGCCTGAagaga   >  2:1936234/1‑150 (MQ=255)

                                                                                                                                                    tAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAGGCCGGCACGTGGGCCTGAagaga   <  1:4528889/150‑1 (MQ=255)

ACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAGGCCGGCACGTGGGCCTGAAGAGAG  >  CP009273/4356946‑4357244

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