BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,471	A→G	32.4%	S27S (AGT→AGC)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,471	0	A	G	32.4%	187.6 / 116.8	165	S27S (AGT→AGC)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base A (57/51); new base G (24/29); total (83/82)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 170 total aligned reads displayed.

CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGT  >  CP000730/36331‑36609

cGCCATTTGTAATAAGCTTTTTTTTTCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCCGATATATCCCGATaa                                                                                                                                             >  2:654033/1‑141 (MQ=11)

cGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATaa                                                                                                                                             <  1:854560/141‑1 (MQ=18)

 gCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAAc                                                                                                                                            <  2:278061/141‑1 (MQ=39)

  ccctttGTAATAAGCTTTTTTATGCTTTTTCTTCAAAATTTGATATAAAATTGGGGCATATTCTTGAACCTAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACAATACCt                                                                                                                                           >  2:23981/4‑141 (MQ=11)

  ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAACTTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGCCCGTTGAATGATGACCGTTTACACCACGTACCCTTAATATTTCTGATATATCACGATaaac                                                                                                                                           >  2:251514/1‑139 (MQ=11)

  ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATGTTCTTGAACCCAACGGTAGACCGTTGAATGGTGAACGTTTACACCCCGTTCCCTTACTATTTCAGATATATGACGATAACt                                                                                                                                           >  2:582954/1‑141 (MQ=11)

  ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTATTATTTCAGATATATCACGATAGCt                                                                                                                                           >  1:654148/1‑141 (MQ=11)

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

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

  ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACt                                                                                                                                           >  1:44383/1‑141 (MQ=25)

  ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATAACGATAACt                                                                                                                                           >  1:468074/1‑141 (MQ=21)

   cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTc                                                                                                                                          >  1:62229/1‑141 (MQ=11)

   cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTc                                                                                                                                          >  2:563539/1‑141 (MQ=255)

   cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTc                                                                                                                                          >  1:562453/1‑141 (MQ=25)

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

   cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTc                                                                                                                                          <  1:54976/141‑1 (MQ=25)

   cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATCGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATGTCAGATGTATCACGATAACTc                                                                                                                                          <  2:727606/141‑1 (MQ=255)

    aTTTGTAATAATCTTTTTTAGGCTTTTTCTTCCAAATTTGGTATTAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCa                                                                                                                                         <  1:193663/141‑1 (MQ=11)

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

    aTTTGTAATAAGCTTGTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCa                                                                                                                                         <  2:457171/141‑1 (MQ=11)

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

     tttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCACTTAATATTTCAGATATATCACGATAACTCaa                                                                                                                                        >  1:433065/1‑141 (MQ=255)

      ttttaataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAAt                                                                                                                                       <  1:247912/138‑1 (MQ=11)

      ttttaataaGCTTTTTTATGCTTTTTCTTCCAAATTTGAGATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAAt                                                                                                                                       <  2:622074/138‑1 (MQ=11)

      ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCACCGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTCATATTTCAGATATATCACAGTAGCTCact                                                                                                                                       >  2:640572/1‑139 (MQ=11)

      ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGGTATATCACGATAGCTCAAt                                                                                                                                       >  2:206740/1‑141 (MQ=11)

      ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAAt                                                                                                                                       <  1:413385/141‑1 (MQ=255)

      ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAAt                                                                                                                                       <  1:97857/141‑1 (MQ=11)

      ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAAt                                                                                                                                       <  1:347356/141‑1 (MQ=11)

      ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAAt                                                                                                                                       >  2:544905/1‑141 (MQ=255)

      ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAAt                                                                                                                                       <  2:451898/141‑1 (MQ=11)

      ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGACCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAAt                                                                                                                                       <  1:703768/141‑1 (MQ=18)

      ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAAt                                                                                                                                       >  1:985869/1‑141 (MQ=40)

       tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATCTTTCCGATATATCACGATAGCTCAATg                                                                                                                                      >  2:345732/1‑141 (MQ=11)

       tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATg                                                                                                                                      <  1:953806/141‑1 (MQ=11)

       tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATg                                                                                                                                      <  1:1008278/141‑1 (MQ=255)

       tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATg                                                                                                                                      >  1:333149/1‑141 (MQ=25)

       tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATg                                                                                                                                      <  2:366991/141‑1 (MQ=21)

       tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATg                                                                                                                                      <  2:325234/141‑1 (MQ=39)

        ttaataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCACTgc                                                                                                                                     <  2:515413/140‑2 (MQ=11)

        gTAATAAGCTTTTTTATTCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgc                                                                                                                                     <  2:786551/141‑2 (MQ=11)

        gTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgc                                                                                                                                     <  1:1001540/141‑2 (MQ=11)

        gTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATTACGATAACTCAATgc                                                                                                                                     <  1:789018/141‑2 (MQ=12)

         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGTCATATTCTTGATCCCAACGGTAGATCGTTTAATTATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgca                                                                                                                                    <  2:48614/141‑3 (MQ=11)

         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGATCGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCtttttgct                                                                                                                                    >  2:409799/1‑134 (MQ=11)

         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCATTgca                                                                                                                                    >  1:303607/1‑139 (MQ=11)

         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgca                                                                                                                                    <  1:120540/141‑3 (MQ=11)

         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgca                                                                                                                                    <  1:189738/141‑3 (MQ=11)

         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgca                                                                                                                                    >  2:106745/1‑139 (MQ=11)

         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCCCAATgca                                                                                                                                    >  2:107569/1‑139 (MQ=11)

         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAAAGTTTACACAACGATCTCTTAATATTTCAGATATATCACGATAGCTCaaaaca                                                                                                                                    >  2:897913/1‑137 (MQ=11)

         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATgca                                                                                                                                    <  2:504421/141‑3 (MQ=17)

         taataaGCTTTTTTATGCTTTTTCCTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTATTATTTTACATATATCACGATAACTCAATgca                                                                                                                                    <  1:514735/141‑3 (MQ=11)

         taataaGCTTTTTTATGCGTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgca                                                                                                                                    <  2:415527/141‑3 (MQ=11)

          aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTGTTAATATTACGGATATATCCAGATAGCTCGATgcaa                                                                                                                                   >  2:1048947/1‑138 (MQ=11)

          aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTgaatg‑aat                                                                                                                                  >  1:693710/1‑141 (MQ=11)

          aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATGCCTCAATgcaa                                                                                                                                   >  1:20341/1‑138 (MQ=255)

            taaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATgcctat                                                                                                                                 >  1:128654/1‑136 (MQ=0)

                       aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGAtagta                                                                                                                     >  2:860584/1‑141 (MQ=35)

                             tttCCTCCAAATTTGATACAAAATTGGTGCTTATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGccgtt                                                                                                               <  1:35865/141‑4 (MQ=255)

                                  tcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTa                                                                                                          >  1:1025172/1‑141 (MQ=35)

                                       aTTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCAAGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTg                                                                                                     <  1:580186/141‑1 (MQ=25)

                                         ttGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGAt                                                                                                   <  2:506800/141‑1 (MQ=35)

                                          tGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGCTGCCCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGAGAACTCAATG‑TATAGCTTAGATAATAGCCAACGGCTACCGTGATa                                                                                                  >  2:692055/1‑141 (MQ=11)

                                          tGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATa                                                                                                  >  1:228455/1‑141 (MQ=35)

                                            aTACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAAc                                                                                                >  1:86430/1‑141 (MQ=35)

                                               cAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACTTc                                                                                             >  1:656641/1‑141 (MQ=17)

                                                   aTTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTAAAGG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTg                                                                                         <  1:728073/141‑1 (MQ=17)

                                                    ttcgggaagagtCTTTAACCCAACTGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGt                                                                                        <  1:87006/130‑1 (MQ=11)

                                                           atatTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTCAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTCGtacaaacg                                                                                 >  2:630935/1‑134 (MQ=255)

                                                             atTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGtt                                                                               <  1:709239/141‑1 (MQ=32)

                                                             atTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGtt                                                                               <  1:773464/141‑1 (MQ=32)

                                                                      cccAACTGTAGACCGTTGAATGATTAACGTTTACACCACGTTCCCTTAATATTTCAGATATAGCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGa                                                                      <  2:923171/141‑1 (MQ=12)

                                                                          aCGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATa                                                                  >  1:414258/1‑141 (MQ=32)

                                                                             ggtgACCGTTAAATGATTAACGGTTACACCACGTTCCCTTCATATCTCAGATATATCACGACACCTCAATG‑TATATCCTAGATAGTAGCCAACGGCTAAAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGtt                                                               <  1:489085/138‑1 (MQ=11)

                                                                                gACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCAt                                                            >  1:1031138/1‑141 (MQ=32)

                                                                                  ccGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGTTAAGGCAATC‑TCCATTTTATATAGTAGCCAACGGCTACAGTGATAACATCCGTGTTAAATTGTTTATATCTGAAATAGTTCATAc                                                          >  1:508680/1‑141 (MQ=11)

                                                                                  ccGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATAc                                                          <  1:66892/141‑1 (MQ=32)

                                                                                   cGTTGAATTATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACa                                                         <  1:387279/141‑1 (MQ=17)

                                                                                        aaTGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaaga                                                    <  1:312597/141‑1 (MQ=32)

                                                                                          tgatgaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACt                                                  >  2:742127/1‑141 (MQ=32)

                                                                                              gaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCtt                                              <  2:686530/141‑1 (MQ=32)

                                                                                                aCGTTCACACCACGTTCGCTTAATATTACTGATATATCA‑GCCCAGTTAACGTTAGAGGTTAGTTAATATCGATCTGCTACACTGAGTACCTCCTTGTTAAATTGTTTATATCTGAAATATTTCATACAGAAAACTCCtttt                                            >  2:382093/1‑141 (MQ=2)

                                                                                                 cGTTTACACCACGTACCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGGAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGATGAATCCttttt                                           >  2:461282/1‑141 (MQ=255)

                                                                                                   tttACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGt                                         <  2:847029/141‑1 (MQ=32)

                                                                                                    ttACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAACCGTTCAGACGGAAGACTCCTTTTTGtt                                        <  1:1053063/141‑1 (MQ=255)

                                                                                                    tcacatcgcGTTCCCTTTATATTTCAGATATATCACGATAACTCTATG‑GATATCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                        <  2:901211/133‑1 (MQ=1)

                                                                                                     tACAGCACGTTCCCTTAATATTTCAGATATGTCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGttc                                       >  2:646080/1‑140 (MQ=17)

                                                                                                         ccacGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaaaag                                   >  2:645991/1‑140 (MQ=17)

                                                                                                          cacGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAAtt                                  <  1:1003147/141‑1 (MQ=32)

                                                                                                              ttgccTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATAc                              <  2:296580/138‑1 (MQ=17)

                                                                                                                   ttAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATaa                         >  2:955833/1‑141 (MQ=33)

                                                                                                                            cAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACtt                >  1:344295/1‑141 (MQ=35)

                                                                                                                               atatatCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGc             >  1:275111/1‑141 (MQ=35)

                                                                                                                                tatatCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCa            >  2:880522/1‑141 (MQ=35)

                                                                                                                                aatatTACGATAACTCAATG‑TATATCTTAGATAGTAGACAACGGCTAGAGTAATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCa            <  1:248578/140‑1 (MQ=255)

                                                                                                                                 atatCACGATAGCTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCaa           >  1:156417/1‑141 (MQ=25)

                                                                                                                                     cACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGa       >  2:289140/1‑141 (MQ=34)

                                                                                                                                        gATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAcc    <  1:418981/141‑1 (MQ=34)

                                                                                                                                          tAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGt  <  2:441751/141‑1 (MQ=37)

CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGT  >  CP000730/36331‑36609

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