BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,463,242	T→A	8.2%	V373E (GTG→GAG) ‡	paaK →	phenylacetyl‑CoA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,463,242	0	T	A	8.2%	332.5 / 13.6	127	V373E (GTG→GAG) ‡	paaK	phenylacetyl‑CoA ligase
Reads supporting (aligned to +/- strand): ref base T (56/59); new base A (6/5); total (62/65)
Fisher's exact test for biased strand distribution p-value = 7.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.22e-01

Only 100 of 141 total aligned reads displayed.

CATCAGCGGACGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAG  >  NC_000913/1463093‑1463391

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

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

   cAGCGGACGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGaa                                                                                                                                                    <  2:1136815/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)

         aCGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGa                                                                                                                                              >  2:1415629/1‑150 (MQ=255)

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

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

           gCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCACAGAAAGTgaatc                                                                                                                                            <  1:1268211/150‑5 (MQ=255)

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

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

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

               cGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTgaa                                                                                                                                        <  2:192596/150‑1 (MQ=255)

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

                    aCATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTgaaagaaa                                                                                                                                   <  2:411245/150‑1 (MQ=255)

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

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

                       tGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTa                                                                                                                                >  1:1710274/1‑150 (MQ=255)

                          tGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGcc                                                                                                                             >  1:545847/1‑150 (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)

                                ttCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACAt                                                                                                                       <  1:885617/150‑1 (MQ=255)

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

                                     ggTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACg                                                                                                                  <  1:456789/150‑1 (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)

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

                                                      gTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCa                                                                                                 <  1:46293/150‑1 (MQ=255)

                                                      gTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCa                                                                                                 >  1:1622440/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)

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

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

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

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

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

                                                                   gaagaGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAg                                                                                    <  2:460263/150‑1 (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)

                                                                                     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)

                                                                                            tATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGAATGCCACCAGCTGCGCCATCGGATTAAGTCGAtgg                                                           <  2:579304/150‑1 (MQ=255)

                                                                                               cgccgcATTACCAACTGTAGTTGAACCCCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGAAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGAtggtgg                                                        <  2:1854390/150‑1 (MQ=255)

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

                                                                                                    cATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGAtc                                                   >  2:512014/1‑150 (MQ=255)

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

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

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

                                                                                                           aaCTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCg                                            >  1:1525578/1‑150 (MQ=255)

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

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

                                                                                                                gAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGAtctg                                       >  1:990840/1‑150 (MQ=255)

                                                                                                                  ggTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGCGGAGTTGAAAGAAAGTAGCCTCACATTGACGCATGAGCAGCGTTGACAGGTATGCCACCAGCTGCGCCATCGGATTAAGTAGATGGTGGGGAACTCTACCGATCtgat                                     >  1:1569024/1‑150 (MQ=255)

                                                                                                                     gAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTgatgat                                  >  2:1792790/1‑150 (MQ=255)

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

                                                                                                                       aCCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCg                                >  1:1604725/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)

                                                                                                                                        actccACTTTCACAGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATc               >  2:1674700/5‑150 (MQ=255)

                                                                                                                                         ctccACTTTCACAGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATcc              >  2:598584/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:1415629/150‑1 (MQ=255)

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

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

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

                                                                                                                                                    gtgAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGaa   >  2:500454/1‑150 (MQ=255)

                                                                                                                                                    gtgAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGAGGCAGTATCCCGCGTTCAGaa   >  2:1658295/1‑150 (MQ=255)

                                                                                                                                                     tgAAAGTGGAGTTGAAAGAACGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGaag  <  2:68530/150‑1 (MQ=255)

                                                                                                                                                     tgAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGaag  >  2:1530890/1‑150 (MQ=255)

                                                                                                                                                     tgAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGaag  <  1:628469/150‑1 (MQ=255)

CATCAGCGGACGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAG  >  NC_000913/1463093‑1463391

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