BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,593,190	T→A	6.2%	D194V (GAC→GTC)	ypfH ←	esterase YpfH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,593,190	0	T	A	6.2%	506.5 / 14.9	160	D194V (GAC→GTC)	ypfH	esterase YpfH
Reads supporting (aligned to +/- strand): ref base T (74/76); new base A (3/7); total (77/83)
Fisher's exact test for biased strand distribution p-value = 3.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 160 total aligned reads displayed.

CGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGC  >  NC_000913/2593042‑2593338

  cgTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGa                                                                                                                                                   >  1:710428/1‑150 (MQ=255)

  cgTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGa                                                                                                                                                   <  1:107982/150‑1 (MQ=255)

       tATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCa                                                                                                                                              >  2:361162/1‑150 (MQ=255)

       tATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCa                                                                                                                                              >  2:2465810/1‑150 (MQ=255)

          cAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATCGACAATCGCagta                                                                                                                                           <  1:1890971/150‑4 (MQ=255)

          cAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATCGACAATCGCagta                                                                                                                                           <  2:687471/150‑4 (MQ=255)

           aGGCCAGGTAATCTGTCTTTTCCATCTTTTCAATCACATCATGATCGCCAGGCTTTCCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGAGCGACAATCGCagtat                                                                                                                                          <  2:2059651/150‑5 (MQ=255)

           aGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATCGACAATCGCagtat                                                                                                                                          <  1:1960845/150‑5 (MQ=255)

           aGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTAACTCCACTTAGCGCTTCTTCGAAATAATGTTTCGGAATGGTATAACTCAAGTGATCCAGGTAAAATTGCATACTGCGATCGACAATCGCagtat                                                                                                                                          <  2:2051135/150‑5 (MQ=255)

            ggCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGcc                                                                                                                                         >  2:2199555/1‑150 (MQ=255)

            ggCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGcc                                                                                                                                         >  2:2236691/1‑150 (MQ=255)

             gCCAGCTAATCTGTCTTATACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTAAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCt                                                                                                                                        >  2:1000149/1‑150 (MQ=255)

                 ggTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAAt                                                                                                                                    >  1:1937405/1‑150 (MQ=255)

                    aaTCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAAtcct                                                                                                                                 >  1:1144074/1‑150 (MQ=255)

                      tCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCAt                                                                                                                                             <  2:869509/136‑1 (MQ=255)

                      tCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAAtcctcc                                                                                                                               <  2:344532/150‑1 (MQ=255)

                       cTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCa                                                                                                                              >  1:346148/1‑150 (MQ=255)

                        tGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCAc                                                                                                                             <  1:1027740/150‑1 (MQ=255)

                           cTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGat                                                                                                                          <  1:1729739/150‑1 (MQ=255)

                              ttACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATc                                                                                                                       >  1:802598/1‑150 (MQ=255)

                                 catcatTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAg                                                                                                                    >  2:1147587/1‑150 (MQ=255)

                                 catcatTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAg                                                                                                                    >  1:931632/1‑150 (MQ=255)

                                  atcatTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGc                                                                                                                   >  1:874277/1‑150 (MQ=255)

                                   tcatTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGcg                                                                                                                  <  1:2167571/150‑1 (MQ=255)

                                   tcatTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGcg                                                                                                                  <  2:2074759/150‑1 (MQ=255)

                                    catTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGcgt                                                                                                                 <  2:1246774/150‑1 (MQ=255)

                                        tCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGcgtcacg                                                                                                             <  2:931632/150‑1 (MQ=255)

                                         cAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGcgtcacgt                                                                                                            >  1:1070193/1‑150 (MQ=255)

                                          aaTCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGcgtcacgtc                                                                                                           <  1:990161/150‑1 (MQ=255)

                                           aTCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGcgtcacgtca                                                                                                          <  1:1694110/150‑1 (MQ=255)

                                            tCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCgtcacgtcac                                                                                                         <  2:858747/150‑1 (MQ=255)

                                              acaTCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTcacca                                                                                                       >  2:223381/1‑150 (MQ=255)

                                                  catcatcGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCgg                                                                                                   <  2:78847/150‑1 (MQ=255)

                                                  catcatcGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCgg                                                                                                   <  1:1564640/150‑1 (MQ=255)

                                                  catcatcGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCgg                                                                                                   <  1:678080/150‑1 (MQ=255)

                                                   atcatcGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGc                                                                                                  >  1:1982805/1‑150 (MQ=255)

                                                    tcatcGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCa                                                                                                 >  2:727972/1‑150 (MQ=255)

                                                      atcGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACt                                                                                               >  1:758594/1‑150 (MQ=255)

                                                       ccGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATACCGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTg                                                                                              <  2:2329399/149‑1 (MQ=255)

                                                        cGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCCCCGGCACTGa                                                                                             >  1:2432731/1‑150 (MQ=255)

                                                        cGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGa                                                                                             >  2:1845219/1‑150 (MQ=255)

                                                         gCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGAt                                                                                            >  2:74943/1‑150 (MQ=255)

                                                         gCCAGGCTTACCGCCACTTAGCGCTTCATCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGAt                                                                                            <  2:1709273/150‑1 (MQ=255)

                                                              gCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACg                                                                                       >  1:1027209/1‑150 (MQ=255)

                                                              gCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACg                                                                                       <  1:848379/150‑1 (MQ=255)

                                                               cTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGc                                                                                      >  2:1286110/1‑150 (MQ=255)

                                                                  aCCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTc                                                                                   <  1:1845219/150‑1 (MQ=255)

                                                                   ccGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTcc                                                                                  <  2:1937405/150‑1 (MQ=255)

                                                                   ccGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTcc                                                                                  <  1:1387590/150‑1 (MQ=255)

                                                                   ccGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTcc                                                                                  <  2:1164398/150‑1 (MQ=255)

                                                                    cGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCt                                                                                 <  1:2366207/150‑1 (MQ=255)

                                                                    cGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCt                                                                                 <  1:2374891/150‑1 (MQ=255)

                                                                      ccACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCGCCACCGGCACTGATTAACGCTTCCTgc                                                                               >  2:962035/1‑150 (MQ=255)

                                                                      ccACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTgc                                                                               <  1:1433873/150‑1 (MQ=255)

                                                                       cACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTgcg                                                                              <  1:2465810/150‑1 (MQ=255)

                                                                        aCTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTgcgc                                                                             >  2:2389671/1‑150 (MQ=255)

                                                                          ttAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGcgccg                                                                           <  2:1980805/150‑1 (MQ=255)

                                                                               gcTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAAcc                                                                      <  2:752493/150‑1 (MQ=255)

                                                                                     tcgAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCCACCGCATgc                                                                >  2:1818400/1‑150 (MQ=255)

                                                                                     tcgAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATgc                                                                <  2:1959218/150‑1 (MQ=255)

                                                                                        aaaTAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCGTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAAGCGCATGCGcc                                                             <  2:111689/150‑1 (MQ=255)

                                                                                        aaaTAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGcc                                                             <  1:639214/150‑1 (MQ=255)

                                                                                         aataatGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCa                                                            >  1:965483/1‑150 (MQ=255)

                                                                                          ataatGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAg                                                           >  2:617156/1‑150 (MQ=255)

                                                                                           taatGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCGACCGCATGCGCCAGa                                                          >  1:2432876/1‑150 (MQ=255)

                                                                                              tGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCg                                                       >  1:2484196/1‑150 (MQ=255)

                                                                                                     gAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCg                                                <  1:1617503/150‑1 (MQ=255)

                                                                                                       aTGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGa                                              >  2:1431957/1‑150 (MQ=255)

                                                                                                         ggTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATc                                            >  2:831954/1‑150 (MQ=255)

                                                                                                          gTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCt                                           >  1:1588549/1‑150 (MQ=255)

                                                                                                          gTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCt                                           <  2:188982/150‑1 (MQ=255)

                                                                                                             taACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCa                                        >  2:2372610/1‑150 (MQ=255)

                                                                                                               aCGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCAcc                                      <  1:769148/150‑1 (MQ=255)

                                                                                                               aCGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACAACCGGCACTGATTAAAGCTTACTGCGCCGCAACCGCATGCGCAAGATCGATAACCGGATCTTCAcc                                      >  1:2056742/1‑150 (MQ=255)

                                                                                                                cGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCg                                     >  2:1466558/1‑150 (MQ=255)

                                                                                                                cGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCg                                     >  1:2257390/1‑150 (MQ=255)

                                                                                                                cGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCg                                     <  1:2005882/150‑1 (MQ=255)

                                                                                                                    aGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCAt                                 >  1:89661/1‑150 (MQ=255)

                                                                                                                    aGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCAt                                 <  2:1017407/150‑1 (MQ=255)

                                                                                                                       taTCCAGGGCAAATTGCATACTGCTATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGACAGATCTATAACCGGATCTTCACCGCCATGaa                              <  2:2121735/149‑1 (MQ=255)

                                                                                                                       gATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGaa                              >  1:268635/1‑150 (MQ=255)

                                                                                                                            aGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCTACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAgg                         <  1:734449/150‑1 (MQ=255)

                                                                                                                              ggCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATTCTCCAGATCTATAACCTGATCTTCACCGCCATGAATCAGGTg                       <  2:802598/150‑1 (MQ=255)

                                                                                                                                cAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGaa                     >  1:721171/1‑150 (MQ=255)

                                                                                                                                cAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGaa                     >  2:1784156/1‑150 (MQ=255)

                                                                                                                                 aaaTTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAAt                    <  2:1529355/150‑1 (MQ=255)

                                                                                                                                   aTTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTg                  >  2:1558586/1‑150 (MQ=255)

                                                                                                                                     cacgttcTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTc                <  2:277393/144‑1 (MQ=255)

                                                                                                                                       cATACTTCGATTGTAGATAGCATGGCATAAATCATCCACGATATCCATCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGt              >  2:2440652/1‑150 (MQ=255)

                                                                                                                                        gccatGCGATCGACAATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTg             >  1:1201913/5‑150 (MQ=255)

                                                                                                                                         tACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTgg            <  1:361162/150‑1 (MQ=255)

                                                                                                                                         ccatGCGATCGACAATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTgg            >  1:228984/4‑150 (MQ=255)

                                                                                                                                          aCTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGc           <  1:2139154/150‑1 (MQ=255)

                                                                                                                                           cTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCt          >  1:22348/1‑150 (MQ=255)

                                                                                                                                           cTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCt          >  2:94003/1‑150 (MQ=255)

                                                                                                                                             gCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGt        >  1:1346232/1‑150 (MQ=255)

                                                                                                                                              cGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTa       >  1:148581/1‑150 (MQ=255)

                                                                                                                                                aTTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGa     <  1:1294565/150‑1 (MQ=255)

                                                                                                                                                   gTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCATCCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGCCGTGGCTGTAGACGc  >  2:1640180/1‑150 (MQ=255)

                                                                                                                                                   gTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGc  >  1:83660/1‑150 (MQ=255)

CGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTTTCGGAATGGTATAACGCAAGTGATCCAGGGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGC  >  NC_000913/2593042‑2593338

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