BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,463,240	T→A	7.6%	S372S (TCT→TCA)	paaK →	phenylacetyl‑CoA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,463,240	0	T	A	7.6%	367.9 / 13.1	127	S372S (TCT→TCA)	paaK	phenylacetyl‑CoA ligase
Reads supporting (aligned to +/- strand): ref base T (57/60); new base A (5/5); total (62/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

Only 100 of 141 total aligned reads displayed.

CGCATCAGCGGACGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGA  >  NC_000913/1463091‑1463389

 gCATCAGCGGACGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCtg                                                                                                                                                      >  2:451235/1‑150 (MQ=255)

 gCATCAGCGGACGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATAACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCtg                                                                                                                                                      >  1:1672097/1‑150 (MQ=255)

 gCATCAGCGGACGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTCATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCtg                                                                                                                                                      <  1:1550758/150‑1 (MQ=255)

  cATCAGCGGACGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCtgt                                                                                                                                                     <  1:551164/150‑1 (MQ=255)

   aTCAGCGGACGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCtgtg                                                                                                                                                    <  1:1635732/150‑1 (MQ=255)

   aTCAGCGGAAGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCtgtg                                                                                                                                                    <  2:1394755/150‑1 (MQ=255)

     cAGCGGACGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGaa                                                                                                                                                  >  2:10040/1‑150 (MQ=255)

      aGCGGACGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGaaa                                                                                                                                                 >  2:1817154/1‑150 (MQ=255)

         ggACGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCACAGAAAGTg                                                                                                                                              <  2:1251827/150‑1 (MQ=255)

            cGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGTAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAg                                                                                                                                           <  1:497343/150‑1 (MQ=255)

            cGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAg                                                                                                                                           <  2:485348/150‑1 (MQ=255)

            cGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAg                                                                                                                                           >  2:1791803/1‑150 (MQ=255)

              cAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGtt                                                                                                                                         >  2:1578940/1‑150 (MQ=255)

               cgCGATGACATGTGGATTATTCGCGGTGTAAATGTCTTTCCGGCACAACTGGAAGAAGAGAGTGTCATATTCGAACATTTATCGCCGCATTCCCAACTGGAGGGGACCCGCCGGGGGCATCTTGATTCACTTTCACAGAAAGTgaatcaa                                                                                                                                        <  1:952014/149‑7 (MQ=255)

                gCGATGACATGTTGTTTATTCGCGGTGTAAATGTCTTTCCGTCACAACTTGAAGAAGAGATTGGCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGTGCATCTTGATTCACTTTCACAGAAAGTgaatcaag                                                                                                                                       <  1:906892/150‑8 (MQ=255)

                gCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAACTGGAGTTga                                                                                                                                       >  1:1561829/1‑150 (MQ=255)

                   aTGACATGTTGATTATTCTCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGTGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTgaaag                                                                                                                                    <  2:472891/150‑1 (MQ=255)

                       cATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGaaagaaag                                                                                                                                >  2:552202/1‑150 (MQ=255)

                       cATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGaaagaaag                                                                                                                                <  2:1758565/150‑1 (MQ=255)

                       cATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGaaagaaag                                                                                                                                <  2:874258/150‑1 (MQ=255)

                       aaTGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGaaagaaag                                                                                                                                <  1:867153/149‑1 (MQ=255)

                             gATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCt                                                                                                                          >  1:685298/1‑150 (MQ=255)

                              attattCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTAAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTg                                                                                                                         <  2:787165/150‑1 (MQ=255)

                               ttattCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGa                                                                                                                        <  2:1451535/150‑1 (MQ=255)

                                 attCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACa                                                                                                                      >  1:846412/1‑150 (MQ=255)

                                 attCGCGGTGTAAATGTCTTTCAGTAACAACTTGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGTGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGCAAGTAGCCTGACa                                                                                                                      <  1:1068744/150‑1 (MQ=255)

                                  ttCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACAt                                                                                                                     >  1:1244100/1‑150 (MQ=255)

                                      cgGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGAc                                                                                                                 <  2:1839895/150‑1 (MQ=255)

                                       ggTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACg                                                                                                                <  1:456789/150‑1 (MQ=255)

                                        gtgtAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGc                                                                                                               <  2:1859463/150‑1 (MQ=255)

                                         tgtAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCa                                                                                                              >  1:1583527/1‑150 (MQ=255)

                                           tAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATg                                                                                                            <  2:407183/150‑1 (MQ=255)

                                           tAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATg                                                                                                            <  1:1254781/150‑1 (MQ=255)

                                             aaggTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGag                                                                                                          <  1:798976/147‑1 (MQ=255)

                                              aTGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCTTCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGagc                                                                                                         <  2:874545/150‑1 (MQ=255)

                                              aTGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGagc                                                                                                         <  2:1405213/150‑1 (MQ=255)

                                                  cTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCg                                                                                                     >  2:523672/1‑150 (MQ=255)

                                                       cGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGcc                                                                                                >  1:289310/1‑150 (MQ=255)

                                                        gTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCa                                                                                               >  2:933146/1‑150 (MQ=255)

                                                          cacaACTTGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCTACTGGAGGTGAACCGCCGTGGGCATATTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAgg                                                                                             <  1:1103629/150‑1 (MQ=255)

                                                           acaACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGAGAAAGTGGAGTTGAAAGAAAGTCGCCTGACATTGACGCATGAGCAGCGTTGCCAGGt                                                                                            >  1:1153158/1‑150 (MQ=255)

                                                            caACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTa                                                                                           >  1:723191/1‑150 (MQ=255)

                                                            caACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTa                                                                                           >  2:642061/1‑150 (MQ=255)

                                                              aCTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATg                                                                                         <  2:915098/150‑1 (MQ=255)

                                                              aCTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATg                                                                                         <  2:1001986/150‑1 (MQ=255)

                                                                tGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCCGCGTTGCCAGGTATGcc                                                                                       >  2:1813053/1‑150 (MQ=255)

                                                                  gaagaaGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGctac                                                                                     >  1:1508900/1‑150 (MQ=255)

                                                                    agaagaGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGcctcca                                                                                   >  1:1546482/1‑150 (MQ=255)

                                                                    agaagaGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGccacca                                                                                   >  1:601836/1‑150 (MQ=255)

                                                                    agaagaGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGccacca                                                                                   >  1:1241384/1‑150 (MQ=255)

                                                                     gaagaGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAg                                                                                  >  1:389175/1‑150 (MQ=255)

                                                                       agagATTGTCAAATTCGAACATTTATCGCCGCATTACAAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCt                                                                                >  2:543144/1‑150 (MQ=255)

                                                                            ttGTCAAATTCGAACATTTATCGCCGCATTACCAACTGCAGGTGAACCGCCTGGGGCATATTGATTCACTATCTGTGAAAGTGGAGTTGAAAGCAAGTAGCCTGACATTGACGCATGAGCAGCGTAGCCAGGTATGCCACCAGCAGcccc                                                                           >  1:1011429/1‑150 (MQ=255)

                                                                               tCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATc                                                                        <  1:552202/150‑1 (MQ=255)

                                                                                cAAATTCGAACATTTATCTCCGCATTACCAACTGTAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGGGGAGTTGAAAGAAAGTAGCCTGACATAGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCg                                                                       <  1:98010/150‑1 (MQ=255)

                                                                                cAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCg                                                                       <  1:151532/150‑1 (MQ=255)

                                                                                  aaTTCGAACATTTATCGCCGCATTACCATCTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGa                                                                     <  1:846835/150‑1 (MQ=255)

                                                                                   aTTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGGGTTGCCAGGTATGCCACCAGCTGCGCCATCGGAt                                                                    <  2:762563/150‑1 (MQ=255)

                                                                                     tCGAAAATTTATCGCCGCATTACCAACTCGCGGTGAACCGCCTGGTGAATCTTGATTCACTTTATGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAACAGCGTTGCCAGGTAAGCCACCAGCTGCGCCATCGGATTa                                                                  <  2:314941/150‑1 (MQ=255)

                                                                                      cGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTaa                                                                 >  2:233643/1‑150 (MQ=255)

                                                                                       gAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCAGCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAg                                                                <  2:78480/150‑1 (MQ=255)

                                                                                        aaCATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGt                                                               <  2:1066669/150‑1 (MQ=255)

                                                                                        aaCATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTCCGCCATCGGATTAAGt                                                               >  1:517439/1‑150 (MQ=255)

                                                                                          cATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCg                                                             >  2:1191804/1‑150 (MQ=255)

                                                                                           aTTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGa                                                            <  1:1578940/150‑1 (MQ=255)

                                                                                                 cgccgcATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGACTTAATCCGAtggcgc                                                      <  1:21652/150‑4 (MQ=255)

                                                                                                  gccgcATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTggg                                                     <  2:1404684/150‑1 (MQ=255)

                                                                                                  gccgcATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCAATCGGATTAAGTCGATGGTtgg                                                     >  1:1615341/1‑150 (MQ=255)

                                                                                                   ccgcATTCCCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACATTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTgggg                                                    <  2:1079284/150‑1 (MQ=255)

                                                                                                       aTTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGTGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGAtct                                                >  2:820774/1‑150 (MQ=255)

                                                                                                         tACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCCTGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGAGTAAGTCGATGGTGGGGATctct                                              >  1:1638442/1‑150 (MQ=255)

                                                                                                          aCCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTa                                             >  2:1608897/1‑150 (MQ=255)

                                                                                                           ccAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTAc                                            >  1:1607217/1‑150 (MQ=255)

                                                                                                            cAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTAcc                                           <  2:442318/150‑1 (MQ=255)

                                                                                                             aaCTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGCCGCATGATCAGCGTTGCCAGGTATGCCACCCGCTGCGCCATCGGATTAAGTATATGGTGGGGATCTCTACCg                                          >  1:465181/1‑150 (MQ=255)

                                                                                                               cTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGAt                                        >  1:1005743/1‑150 (MQ=255)

                                                                                                                tGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGAtg                                       <  1:1817154/150‑1 (MQ=255)

                                                                                                                        aaCCGCCGGGGGCAGCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATc                               <  1:984446/150‑1 (MQ=255)

                                                                                                                         aCCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCg                              <  1:1254141/150‑1 (MQ=255)

                                                                                                                         aCCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCg                              >  1:1301952/1‑150 (MQ=255)

                                                                                                                         aCCGCCGGGGGCAGCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCg                              <  2:1861238/150‑1 (MQ=255)

                                                                                                                                ggggCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACtg                       <  1:1000085/150‑1 (MQ=255)

                                                                                                                                 gggCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGAAGCGTTGCCAGGTATGCAACCAGCTGCGCCATAGGATGAAGTCGATGGTGGGGATCTCTACCGAAGTCATGATCGTTAAAtgt                      >  1:231387/1‑150 (MQ=255)

                                                                                                                                    cATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCTTCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGc                   <  1:477061/150‑1 (MQ=255)

                                                                                                                                    cATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGc                   <  2:765760/150‑1 (MQ=255)

                                                                                                                                      tCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAg                 >  2:1040025/1‑150 (MQ=255)

                                                                                                                                         tGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTAt              >  2:1449958/1‑150 (MQ=255)

                                                                                                                                         aactccACTTTCACAGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGCTGTGATGATCGTTAACTGTGGCAGTAt              >  1:1283920/6‑150 (MQ=255)

                                                                                                                                           ctccACTTTCACAGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATcc            >  2:590763/4‑150 (MQ=255)

                                                                                                                                            ttCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTATGTGGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATccc           >  2:1148443/1‑150 (MQ=255)

                                                                                                                                            ttCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATccc           >  2:696286/1‑150 (MQ=255)

                                                                                                                                              cACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGCCGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCcgc         <  2:685298/150‑1 (MQ=255)

                                                                                                                                              cACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCcgc         >  2:59427/1‑150 (MQ=255)

                                                                                                                                               aCTTTTTGTGAAAGTGGAGTTTAAAGAAAGTAGCCTTACTTTGACGAATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCcgcg        <  2:1241384/150‑1 (MQ=255)

                                                                                                                                                atttCTGTGAAAGTGGAGGTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGt       <  1:475783/149‑1 (MQ=255)

                                                                                                                                                    cTGTGAAAGTGGAGTTGACAGAAACTCGCCTGACTTTTAGGCTTGAGCAGCGTTGCCAGGTATGCCACCAGCAGCGCCATCGGATTAATTCGATGGTGGGGATCTCTTCCGATGTGATGCTCGTTAACTGGGGCAGTATCCCGCGTTCAg   <  1:139571/150‑1 (MQ=255)

                                                                                                                                                    cTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGCCATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAg   >  2:1553123/1‑150 (MQ=255)

                                                                                                                                                     tgtgAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGa  >  1:1790262/1‑150 (MQ=255)

                                                                                                                                                     tgtgAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGa  >  1:924050/1‑150 (MQ=255)

                                                                                                                                                     tgtgAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGa  >  1:803055/1‑150 (MQ=255)

CGCATCAGCGGACGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGA  >  NC_000913/1463091‑1463389

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