BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,413	C→T	55.6%	V47I (GTC→ATC)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,413	0	C	T	55.6%	140.6 / inf	735	V47I (GTC→ATC)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base C (158/165); new base T (214/192); total (375/360)
Fisher's exact test for biased strand distribution p-value = 3.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.61e-04

Only 100 of 748 total aligned reads displayed.

AATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATAC  >  CP000730/36273‑36553

 aTGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGAcc                                                                                                                                             >  1:333404/1‑141 (MQ=21)

 aTGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGAcc                                                                                                                                             <  1:190249/141‑1 (MQ=18)

 aTGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATCCAATTGGGCCATCTTCTTCAACCCAACGGTAGAct                                                                                                                                             <  2:552432/141‑2 (MQ=11)

  tGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCg                                                                                                                                            >  2:273274/1‑141 (MQ=25)

  tGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCg                                                                                                                                            >  1:355605/1‑141 (MQ=255)

   ggCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGt                                                                                                                                           <  1:473560/141‑1 (MQ=40)

    gCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGtt                                                                                                                                          >  2:490359/1‑141 (MQ=255)

     cccgctataAATGGCTCTATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGTGCATATTCTTGAACCCAACGGTAGACCGTTg                                                                                                                                         <  1:785556/136‑1 (MQ=11)

      aCGATATTAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGa                                                                                                                                        <  2:723382/141‑1 (MQ=11)

      aCGATATAAATAGCTCCATTTTCCGTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGa                                                                                                                                        <  2:279600/141‑1 (MQ=255)

       cGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTTTTCTAATTTTGAGATAAAAATGGGGCATATTCTTGAACCCAACGGTTGTCCGTGGaa                                                                                                                                       >  2:506816/1‑141 (MQ=11)

       cGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTTGACCGTTGaa                                                                                                                                       >  2:701511/1‑141 (MQ=21)

       cGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGaa                                                                                                                                       >  1:104890/1‑141 (MQ=255)

       cGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGaa                                                                                                                                       >  2:458377/1‑141 (MQ=255)

        gATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAat                                                                                                                                      <  2:20117/141‑1 (MQ=18)

        gATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCCACGGTAGATCGTTGAat                                                                                                                                      >  2:824895/1‑141 (MQ=255)

         atatAAATAGCTTCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAAGTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAatg                                                                                                                                     <  2:666155/141‑1 (MQ=255)

         atatAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAatg                                                                                                                                     >  2:372072/1‑141 (MQ=25)

         atatAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAatg                                                                                                                                     >  2:604860/1‑141 (MQ=255)

         atatAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAatg                                                                                                                                     >  1:837972/1‑141 (MQ=255)

         ataagaaTGGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAatg                                                                                                                                     <  2:751579/136‑1 (MQ=255)

          tataAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTTTTTTTTTTTTTCCAAATTTGATATAAAATTGGGGCATATTTTTGAACCCAACGGTAGATCGTTGAatga                                                                                                                                    >  2:34861/1‑141 (MQ=11)

          aataAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAatga                                                                                                                                    <  2:218890/140‑1 (MQ=18)

           ataAATAGCTCCATTTTCCTTTTATTTTTATGTACTTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAatgat                                                                                                                                   <  1:775638/141‑1 (MQ=255)

           ataAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAatgat                                                                                                                                   <  1:659594/141‑1 (MQ=255)

            taAATAGCTCCCTGTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAatgatg                                                                                                                                  <  2:792580/141‑1 (MQ=255)

             aaaTAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAAtgatga                                                                                                                                 >  1:946972/1‑141 (MQ=255)

              aacaGCTCCTGTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGaa                                                                                                                                <  2:172345/138‑1 (MQ=11)

               aTAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAAc                                                                                                                               <  1:229680/141‑1 (MQ=25)

                 aGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGt                                                                                                                             <  1:597354/141‑1 (MQ=255)

                  gCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGtt                                                                                                                            >  2:856249/1‑141 (MQ=255)

                      cATTTTCCTTTTATTTTGTTATACGTCTCATCAATACGCCATTTGTAATACGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTaca                                                                                                                        <  2:912970/141‑1 (MQ=11)

                      cATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGGACGTTTacc                                                                                                                        >  1:498296/1‑140 (MQ=11)

                      cATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTaca                                                                                                                        >  1:460377/1‑141 (MQ=255)

                      cATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGATTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTGGAATGATGAAAGTTTaca                                                                                                                        >  1:131072/1‑141 (MQ=11)

                      cATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCAATTTGTAATAAGCTTTTTTATGCTTTTTCTTACAAATTTGATATAAAATTGGAGCATATTCTTGAACCCAAAGGTAGATCGGTGGataatacatgtctgac                                                                                                                        >  2:563499/1‑130 (MQ=11)

                      cATTTTCCTTTTATTTTGATAGACGTCTCAGCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTaca                                                                                                                        <  2:308938/141‑1 (MQ=11)

                       aTTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTacac                                                                                                                       <  1:709670/141‑1 (MQ=255)

                       aTTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTacac                                                                                                                       <  1:543288/141‑1 (MQ=255)

                       aTTTTCCTTTTATTTTGATATACGGCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTacac                                                                                                                       >  1:936828/1‑141 (MQ=11)

                        ttttCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTAcacc                                                                                                                      <  1:729012/141‑1 (MQ=18)

                         tttCCTTTTATTTTGATGTACGCCTCATCACTACGCCATTTGTAATAAGCTTTTAAATTCTTGTGCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTAcacca                                                                                                                     <  1:447249/141‑1 (MQ=11)

                           tctttttACTTCGATATACGTCTCATTAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATCAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACg                                                                                                                   <  2:413384/138‑1 (MQ=11)

                           tCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACg                                                                                                                   >  1:366002/1‑141 (MQ=255)

                            ccTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGt                                                                                                                  >  1:383021/1‑141 (MQ=25)

                            ccTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGt                                                                                                                  >  1:208945/1‑141 (MQ=255)

                             cTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCCCGtt                                                                                                                 >  1:398641/1‑141 (MQ=21)

                             cTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGtt                                                                                                                 <  1:919728/141‑1 (MQ=255)

                             cTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGtt                                                                                                                 >  2:802102/1‑141 (MQ=255)

                             cTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGtt                                                                                                                 >  2:480079/1‑141 (MQ=255)

                               tttATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATACAATTGGGGCATATCCTTGAACCCACAGGTACATAGTTGAATGATTAATGTTTACAACACGTTct                                                                                                               >  2:780051/1‑140 (MQ=11)

                                tttttttGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTccc                                                                                                              <  2:25322/138‑1 (MQ=18)

                                ttATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTctc                                                                                                              <  2:882408/141‑3 (MQ=255)

                                 tATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTctct                                                                                                             <  1:563895/141‑1 (MQ=255)

                                  tttttGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCtt                                                                                                            <  2:520826/140‑1 (MQ=21)

                                  tttttGATCTACGTCTAATCAATACGCCATTTGTAATCAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCtt                                                                                                            <  1:757351/140‑1 (MQ=255)

                                  tttttGATATACGTCTTATCAATACGCCATTTCTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATACTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCtt                                                                                                            <  1:630642/140‑1 (MQ=11)

                                  aTTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCtt                                                                                                            >  2:896341/1‑141 (MQ=21)

                                  aTTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCtt                                                                                                            <  1:808264/141‑1 (MQ=255)

                                    tttGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTaa                                                                                                          <  1:529540/141‑1 (MQ=18)

                                       gatatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAatat                                                                                                       >  1:299343/5‑141 (MQ=255)

                                       gatatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAatat                                                                                                       >  1:863613/5‑141 (MQ=255)

                                       gatatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCATGAACCCTACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAatat                                                                                                       >  1:696382/5‑141 (MQ=11)

                                       gATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATGTAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCATTGAATGATGAACGTTTACACCACGTTCCCTTAatat                                                                                                       >  1:618617/1‑141 (MQ=255)

                                       gATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCGTAatat                                                                                                       >  2:112482/1‑141 (MQ=25)

                                           tACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCa                                                                                                   <  1:721874/141‑1 (MQ=255)

                                           tACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCa                                                                                                   >  1:510911/1‑141 (MQ=21)

                                           tACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCa                                                                                                   <  2:323496/141‑1 (MQ=18)

                                            aCGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAg                                                                                                  >  2:234751/1‑141 (MQ=255)

                                            aCGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCaa                                                                                                  >  2:439776/1‑140 (MQ=21)

                                             cGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGa                                                                                                 >  2:148372/1‑141 (MQ=255)

                                             cGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGa                                                                                                 >  2:874259/1‑141 (MQ=255)

                                              gTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGat                                                                                                <  1:425013/141‑1 (MQ=255)

                                               tctcATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGata                                                                                               >  2:938378/1‑141 (MQ=255)

                                                 tcatcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCACATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGatata                                                                                             <  2:682925/141‑1 (MQ=255)

                                                   atcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGatatatt                                                                                           >  2:625123/1‑140 (MQ=255)

                                                   atcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATc                                                                                           >  1:176909/1‑141 (MQ=21)

                                                    tcaATACGCCATTTGTAATAAGCTTTTTTTTGTTTTTTCTTCCAAAATTCATCTAACATTGGAGAGTATTCTTAAACCCAACGGTAGGTCGTTGAATGATTACTGTTTGCACCACGTTCTCTTAATATTTCAGATATATCa                                                                                          >  2:292067/1‑141 (MQ=2)

                                                      aaTACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACg                                                                                        >  2:836248/1‑141 (MQ=255)

                                                        tACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGAt                                                                                      >  2:398591/1‑141 (MQ=255)

                                                        tACGACATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTATAGCACGTTCCCTTAATATTTCAGATATATCACGAt                                                                                      >  1:641276/1‑141 (MQ=255)

                                                          cGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCATTACTATTTCAGATATATCACGATaa                                                                                    >  2:763279/1‑141 (MQ=255)

                                                            ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCt                                                                                  >  1:582/1‑141 (MQ=255)

                                                            ccATTTGTAAGAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACTTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACt                                                                                  >  1:682164/1‑141 (MQ=255)

                                                             cttttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTc                                                                                 <  1:256056/139‑1 (MQ=18)

                                                             cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTc                                                                                 >  1:66601/1‑141 (MQ=21)

                                                              aTTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCa                                                                                <  1:828474/141‑1 (MQ=21)

                                                               tttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCaa                                                                               <  1:311516/141‑1 (MQ=11)

                                                                ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCCCAAt                                                                              >  1:900912/1‑141 (MQ=11)

                                                                 tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCaagg                                                                             >  1:846773/1‑139 (MQ=11)

                                                                  gTAATAAGGTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATgc                                                                            <  2:643960/141‑2 (MQ=17)

                                                                  gTAATAAGCTTTTTTATGCTTTTGCTTCCAAATTTGATATAAAATGAGGGCATGTTCTTGAACCCAACGGTTGCTCGATCAATGATGAAAGCTTACACCACGATCTCTTACTCTTTAAGATATATCACGATAGCTCactcc                                                                            >  2:185216/1‑137 (MQ=2)

                                                                   taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAGTGATGAGCGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgca                                                                           <  1:409986/141‑3 (MQ=11)

                                                                   taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCCCGTTCCCTTAATATTTCAGATATATCAGGATAACTCAATgca                                                                           >  2:385080/1‑139 (MQ=255)

                                                                                          cttcCAAAGTTGATACAAAAGTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGc                                                    <  1:277769/141‑1 (MQ=25)

                                                                                                    tGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATa                                          >  1:795137/1‑141 (MQ=35)

                                                                                                            aaTTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCtt                                  <  1:883063/141‑1 (MQ=32)

                                                                                                                           ttGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTtata                   >  2:571846/1‑141 (MQ=32)

                                                                                                                            tGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATCTTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTatat                  >  1:153761/1‑141 (MQ=17)

                                                                                                                               aCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTg               >  1:819449/1‑141 (MQ=32)

AATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATAC  >  CP000730/36273‑36553

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