BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,378	C→T	88.2%	L58L (TTG→TTA)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,378	0	C	T	88.2%	1076.7 / 135.8	589	L58L (TTG→TTA)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base C (37/33); new base T (257/258); total (296/293)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 590 total aligned reads displayed.

CGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACA  >  CP000730/36238‑36516

cGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATAc                                                                                                                                            <  1:409666/141‑1 (MQ=40)

  cAACCAAATATCTAATGTATGTACCACGGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCT‑CATCAATACGCC‑GTTTGTACTAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATaa                                                                                                                                          >  1:372882/1‑141 (MQ=11)

     ccAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCT‑CATCAATATGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAAt                                                                                                                                       >  1:11288/1‑141 (MQ=255)

      cAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAAtt                                                                                                                                      >  2:100306/1‑141 (MQ=38)

        aaTATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATGTACAATTgg                                                                                                                                    >  2:748882/1‑141 (MQ=11)

         atatCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTggg                                                                                                                                   <  1:306414/141‑1 (MQ=37)

          tcatcaaTGCATGTCCCTCTGTATGAATGGCATGCTAGAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTgggg                                                                                                                                  <  2:141063/136‑1 (MQ=11)

           atCTTATTTATGTCCCNCAGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGc                                                                                                                                 <  1:655441/141‑1 (MQ=11)

            tCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGCTATAAAATTGGGGCa                                                                                                                                >  1:282192/1‑141 (MQ=255)

             ttagtgtatgtCCCTCTGCATCAATCTGACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCat                                                                                                                               <  2:36956/137‑1 (MQ=11)

             ttAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCat                                                                                                                               <  1:447290/140‑1 (MQ=255)

             gtAATGTAGGGCCCTCTGTATCAATGGCAGGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCat                                                                                                                               <  2:728442/140‑1 (MQ=11)

             cTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCat                                                                                                                               <  2:184415/141‑1 (MQ=255)

             cTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCat                                                                                                                               <  2:209807/141‑1 (MQ=25)

              tAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCata                                                                                                                              <  2:288630/141‑1 (MQ=37)

                atgtatgtCCCTCTGCATCAATGACACGACATAAATAGCTCCATTATCCTTTTATTTTGATGTACGTCT‑CATCAATACGCC‑ATATGTAATAAGCTATTTTATGCTTTTTCTTCCAAAATTGATATAAAATAGGGACAAAtt                                                                                                                            >  2:761739/1‑141 (MQ=11)

                  gtatgtCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTTTTTTGATAGACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCt                                                                                                                          <  1:145028/141‑1 (MQ=11)

                  gtatgtCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCt                                                                                                                          <  1:169535/141‑1 (MQ=255)

                   tatgtCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCATTTTTTTTTGATGTACGTCT‑CATCATTATGAC‑ATTTGTAATATGATTTTTTATGCTTTTTATTGCACATTTGATATAAAATTAGGGattattatt                                                                                                                         >  1:543934/1‑141 (MQ=11)

                    atgtCCCTCTGCATCAGTGGCACGATATAAATATCTCCATTTTCCTTTTATTTTGATGTACGTCT‑CTTCATTACGCC‑ATTTATAATATGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATGGTGGCAGATTCTTg                                                                                                                        >  1:383313/1‑141 (MQ=11)

                     tgtCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCACATTTGATCTAAAATTTGGGCATATTCTTGa                                                                                                                       >  2:352180/1‑141 (MQ=11)

                     tgtCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGa                                                                                                                       <  1:395110/141‑1 (MQ=255)

                     tgtCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGa                                                                                                                       >  1:738968/1‑141 (MQ=255)

                      gtCCCTTTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGACATATGATTGTGGCATATTCTTGaa                                                                                                                      <  2:544892/141‑1 (MQ=11)

                           tctGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCT‑CATTAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTTCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCaa                                                                                                                 >  1:476005/1‑141 (MQ=255)

                           cctGCATCACTGGCACGATATAAATAGCTCCCTTTTCCTTTTATTTTGCTGTACGTCT‑CATCAATACGCC‑ATTTGTAATAATCTTTTTTATGCTTTTTCTTCAAAATTTGATATAAAATTGGGGCATATTCTTGAACCCaa                                                                                                                 <  2:117472/140‑1 (MQ=11)

                             tGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCT‑CATCAATACGCC‑ATTTGTAATAAACTTTTTTATGTTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACg                                                                                                               >  1:266748/1‑141 (MQ=255)

                              gtagcAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACgg                                                                                                              <  1:220136/137‑1 (MQ=255)

                              gtaacAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACgg                                                                                                              <  2:481548/137‑1 (MQ=255)

                              gtaTCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACgg                                                                                                              >  2:511280/3‑141 (MQ=255)

                                aTCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTa                                                                                                            <  2:631215/141‑1 (MQ=16)

                                 ttaaTGTCACGATATCAATAGCTCCATTTTCCTTTTATTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAg                                                                                                           <  2:152694/139‑1 (MQ=255)

                                 tCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAg                                                                                                           <  1:588632/141‑1 (MQ=39)

                                  cAATGGCACGATATAAATAGCTCCATTTTCTTTTTATTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGa                                                                                                          <  1:331870/141‑1 (MQ=255)

                                   aaTGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGAc                                                                                                         >  1:350903/1‑141 (MQ=40)

                                     tGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCg                                                                                                       <  1:100306/141‑1 (MQ=21)

                                     tGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGatcg                                                                                                       <  2:488246/141‑4 (MQ=255)

                                        cACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTg                                                                                                    <  1:156675/141‑1 (MQ=18)

                                        cACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCT‑CATCAATACTCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGTGGCATATTCTTGAACCCAACGGTAGATCGTTg                                                                                                    >  1:512468/1‑141 (MQ=11)

                                         aCGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGa                                                                                                   >  1:506705/1‑141 (MQ=21)

                                         aCGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGAATATTCTTGAACCCAACGGTAGACCGTTGa                                                                                                   >  1:193486/1‑141 (MQ=21)

                                           gATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTTATAAGCTTTTTTATGCTTTTTTTTCCAAATTTGATATAAAATTGGGGAATATTCTTGAACACAACAGTAGATCGTTTAat                                                                                                 >  2:218008/1‑141 (MQ=11)

                                             tataAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGTGCATATTCTTGAACCCAACGGTAGACCGTTGAatga                                                                                               >  2:446952/1‑141 (MQ=21)

                                             attaAAAAGCTCCATTTTCCTTTTATTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAatga                                                                                               <  2:314305/139‑1 (MQ=11)

                                              ataAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAatgat                                                                                              >  2:431976/1‑141 (MQ=40)

                                              ataAATAGCTCCATTTTCCTTTTATTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAatgat                                                                                              >  1:565593/1‑141 (MQ=255)

                                               taAATAGCTCCATTTTCCTTTTATTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACTCAACGGTAGATAGTTGTatgttg                                                                                             >  2:660025/1‑141 (MQ=11)

                                                  aTAGCTCCATTTTCCTTTTATTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAAc                                                                                          >  1:403164/1‑141 (MQ=255)

                                                    ggcttctttttCCTTTTATTTTGATGTACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGt                                                                                        <  1:8606/134‑1 (MQ=255)

                                                    aGCTCCATTTTCCTTTTATTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATCAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGt                                                                                        >  2:693031/1‑141 (MQ=255)

                                                    aGCTCCATTTTCCTCTTATTTTGATGTACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGt                                                                                        >  2:471335/1‑141 (MQ=21)

                                                     gCTCCATTTTTCTTTTATTTTGACATCCGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGtt                                                                                       <  1:456135/141‑1 (MQ=11)

                                                     gCTCCATTTTCCTTTTATTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGtt                                                                                       <  2:420670/141‑1 (MQ=255)

                                                     gCTCCATGTTCCTTTTATTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGtt                                                                                       <  2:271496/141‑1 (MQ=255)

                                                         cATTTTCCTTTTATTTTGATGTACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTaca                                                                                   >  1:748316/1‑141 (MQ=18)

                                                          aTTTTCCTTTTATTTTGATGTACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTacac                                                                                  <  1:701054/141‑1 (MQ=18)

                                                          aTTTTCCTTTTATTTTCATGTACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTacac                                                                                  <  1:241539/141‑1 (MQ=18)

                                                           ttttCCTTTTATTTTGATG‑GCGTTTCCATCAATACGCC‑ATTTGTAATAAGCTTGGTTATGCGTTTTCTTCCAAATGCGATATAAAATTGGGTCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTAcacc                                                                                 <  2:594307/141‑1 (MQ=11)

                                                           ttttCCTTTTATTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTGTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTAcacc                                                                                 <  2:498982/141‑1 (MQ=11)

                                                           tttcCCTTTTATTCTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTAcacc                                                                                 <  2:406504/141‑1 (MQ=11)

                                                            tttCCTTTTGTTCTCATGTACGTCT‑CATCAATACGCC‑ATTTGTTATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTAcacca                                                                                <  2:607230/141‑1 (MQ=11)

                                                              tCCTTTTATTTTGATGTACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACg                                                                              <  2:707438/141‑1 (MQ=21)

                                                              tCCTTTTATTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACg                                                                              <  2:331030/141‑1 (MQ=255)

                                                                cTTTTATTTTGATGTACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTACACCGTTCAATGATGACCGTTTACACCACGtc                                                                            >  2:56175/1‑140 (MQ=255)

                                                                cTTTTATTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTGCACCACGtt                                                                            >  1:122485/1‑141 (MQ=255)

                                                                cTTTTATTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGtt                                                                            <  2:716110/141‑1 (MQ=255)

                                                                cTTTTATTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATAGTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGtt                                                                            >  1:572877/1‑141 (MQ=255)

                                                                 ttttATTTTGATGTACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTc                                                                           <  1:136002/141‑1 (MQ=18)

                                                                    tATTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATGTTCGTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTctct                                                                        >  1:61216/1‑141 (MQ=11)

                                                                     tttttGATCTACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCtt                                                                       <  1:622304/140‑1 (MQ=11)

                                                                     aTTTTGTTGTACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCtt                                                                       <  2:198214/141‑1 (MQ=255)

                                                                     aTTTTGATGTACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCtt                                                                       <  2:584272/141‑1 (MQ=21)

                                                                     aTTTTGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCtt                                                                       <  2:171367/141‑1 (MQ=11)

                                                                      ttttGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTctttta                                                                      >  2:383671/1‑136 (MQ=11)

                                                                        ttGATATACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAat                                                                    <  2:620471/141‑1 (MQ=11)

                                                                         tGATATACGTCT‑CATCAATACCCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAata                                                                   <  2:556183/141‑1 (MQ=11)

                                                                          gatatACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAatat                                                                  >  1:208481/5‑141 (MQ=255)

                                                                           aTGTACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGATTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACTGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATAtt                                                                 >  1:313224/1‑141 (MQ=255)

                                                                            tatACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATAttt                                                                >  2:241130/3‑141 (MQ=255)

                                                                            tGTACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATAttt                                                                >  1:676196/1‑141 (MQ=40)

                                                                              tgcGTCT‑CATCAATACGCCGAAGTGT‑ATAAGCGTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGGACGTTTACACCACGTTCTCTTAATATTTCa                                                              <  2:535325/139‑1 (MQ=11)

                                                                              tACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCa                                                              >  1:692127/1‑141 (MQ=255)

                                                                                cGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGa                                                            >  1:131612/1‑141 (MQ=255)

                                                                                 gTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGat                                                           <  2:305079/141‑1 (MQ=40)

                                                                                  tct‑cATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGata                                                          >  1:527599/1‑141 (MQ=40)

                                                                                   ct‑cATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGatat                                                         <  1:728404/141‑1 (MQ=18)

                                                                                        tcaATACGCC‑ATTTGTAATAAGCTTTTTTTTTCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTCAATATTTCAGATATATCa                                                     >  2:12830/1‑141 (MQ=255)

                                                                                        tcaATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCa                                                     >  2:375495/1‑141 (MQ=255)

                                                                                           aTACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGa                                                  <  2:91709/141‑1 (MQ=255)

                                                                                            tACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGaa                                                 >  1:64208/1‑140 (MQ=255)

                                                                                              cGCC‑ATTTGTAATGAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATaa                                               <  1:1837/141‑1 (MQ=17)

                                                                                                cc‑ATTTGTAATAAGCTTTTTTATGCTTTTTTTTTTATATTTGATATAAAATTGGGACATATTGTATAACCCAACGGTAGGTGGTTGAATGATTAACGTTTACACCACGTCCTCTTAATATTTCAGATATATCACgctgtct                                             >  2:234020/1‑135 (MQ=11)

                                                                                                cc‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACt                                             >  1:423351/1‑141 (MQ=25)

                                                                                                 c‑ATTTATAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATAGAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTc                                            <  1:367767/141‑1 (MQ=11)

                                                                                                      tGTAATATGCTTTTTTATGCTTTTTCTTTCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATg                                        <  1:251246/141‑1 (MQ=255)

                                                                                                      tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACTGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAAGATTTCAGATATATCACGATAACTCAATg                                        <  1:481338/141‑1 (MQ=255)

                                                                                                      tGTAATAAGCGTTTTTATGCCTTTTCGTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATg                                        <  2:678825/141‑1 (MQ=11)

                                                                                                       gTAATAAGCTTTTTTATGCTTTTTCTTCCACACTTGATATCAATTTGGGTCATATTCTTGAACCCAATGAGGGATGGTTCACAGTTTTGAAGTTACTCCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgc                                       <  1:754291/141‑2 (MQ=2)

                                                                                                        taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTAAATgca                                      >  1:435072/1‑139 (MQ=11)

                                                                                                                      aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagta                        >  1:492918/1‑141 (MQ=35)

CGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCT‑CATCAATACGCC‑ATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACA  >  CP000730/36238‑36516

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