BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,471	A→G	34.5%	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	34.5%	125.3 / 104.3	130	S27S (AGT→AGC)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base A (48/37); new base G (20/25); total (68/62)
Fisher's exact test for biased strand distribution p-value = 2.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

Only 100 of 135 total aligned reads displayed.

CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACC  >  CP000730/36331‑36607

cGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAACTTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACTCCACGTTTCTTTAATTTTTTAGATATATCACGGTaa                                                                                                                                           >  1:80076/1‑141 (MQ=11)

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

cGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATaa                                                                                                                                           <  2:531112/141‑1 (MQ=21)

 cccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAAc                                                                                                                                          <  2:765457/140‑1 (MQ=40)

  ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGGTATAAAAGTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCt                                                                                                                                         >  1:416904/1‑141 (MQ=11)

  ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATTACGCTAGCt                                                                                                                                         >  2:573776/1‑141 (MQ=11)

  ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCt                                                                                                                                         >  1:565915/1‑141 (MQ=11)

  ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCCCGTTCCCTTAATATTTCGGATATTTTACGGTaatt                                                                                                                                         >  1:264366/1‑139 (MQ=11)

  ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAAcc                                                                                                                                         >  2:202201/1‑140 (MQ=21)

  ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCCTGAACCCAACGGTAGATCGTTGCATGATGAACGTGTAGAGCACGTTCTCTTAATATTTCAGATATATCACGATAGCt                                                                                                                                         >  2:716410/1‑141 (MQ=11)

   cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGCTATAATATTGGGGCTTATTCGTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGCTAACta                                                                                                                                        >  2:479688/1‑140 (MQ=11)

   cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGGTAGCta                                                                                                                                        >  1:496066/1‑140 (MQ=11)

    aTTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCa                                                                                                                                       <  1:283539/141‑1 (MQ=11)

    aTTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCCCTTTATATTTCAGATATATCACGATAGCTCa                                                                                                                                       <  2:437676/141‑1 (MQ=255)

    aTTTGTAATAAGCATATTTATGCTTATACTTCCAAAATTGATATAAAATTGGGGCATATTCTTGAACCAAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCa                                                                                                                                       <  1:638641/141‑1 (MQ=11)

     tttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGCTATATCACGATAGCGCaa                                                                                                                                      >  2:269438/1‑141 (MQ=11)

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

     tttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCaa                                                                                                                                      <  1:257465/141‑1 (MQ=21)

     tttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCaa                                                                                                                                      <  2:633660/141‑1 (MQ=21)

     tttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCaa                                                                                                                                      <  2:540448/141‑1 (MQ=21)

     tttGTAATAAGCTTTCTTATGCTTTTTCTCCCAAATTTGATATATACCTGCGCCTCAGACTCGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCaa                                                                                                                                      <  1:219437/141‑1 (MQ=11)

     tttCTAATAAGCTTTTTCTTGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTCGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCaa                                                                                                                                      <  2:528253/141‑1 (MQ=11)

     gttGTAATAAGCTGTATTATGCTTGTTCTACCAAATTTGATACAAAATTGGGGCATATTCAGGAACCCAACGATAGACCGATGAATGATGAAGGTTTACACACCGTTCCCTTAATATTTCAGATATATCACGATAACTCaa                                                                                                                                      >  2:353346/2‑141 (MQ=11)

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

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

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

      ttGTAATAAGCTTTTTTATGCTTTTTCTCCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAAt                                                                                                                                     <  2:55673/141‑1 (MQ=11)

       tGTAATAAGCTTTTTTATGCTTTTTCTTCCACATTTGATATAAAATTGGGGCATATTTTTTAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATg                                                                                                                                    >  2:167028/1‑141 (MQ=255)

       tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGTACGTTTACACCACGTTCTTTTAATATTTCAGATATATCACGATAGCTCAATg                                                                                                                                    >  1:595859/1‑141 (MQ=11)

       tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTATACCACGTTCTCTTAATATTTCAGATATATTACGATAGCTCAATg                                                                                                                                    <  1:431441/141‑1 (MQ=11)

       tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATg                                                                                                                                    >  1:786455/1‑141 (MQ=11)

       tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTGTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATg                                                                                                                                    >  1:105735/1‑141 (MQ=11)

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

       tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATg                                                                                                                                    <  2:690831/141‑1 (MQ=18)

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

        ttaataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgc                                                                                                                                   <  1:233503/140‑2 (MQ=11)

        ttaataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgc                                                                                                                                   <  2:225225/140‑2 (MQ=11)

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

        gTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATgc                                                                                                                                   >  2:642231/1‑140 (MQ=18)

         ttataaGCTTTTTTATGCTTTTTCTGTCAAATTTGATATAAAATGGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGCTAGCTCAATgca                                                                                                                                  <  2:214525/139‑3 (MQ=11)

         taataaGCTTTTTTATGCTTTTTCTTCCAACTTTGATATAAAATTGGGGCATAATCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATgca                                                                                                                                  >  1:389432/1‑139 (MQ=255)

         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCACTgca                                                                                                                                  >  1:635242/1‑139 (MQ=11)

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

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

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

         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgca                                                                                                                                  >  2:103200/1‑139 (MQ=255)

         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACtgcctgac                                                                                                                                  >  1:369917/1‑134 (MQ=255)

         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGtt                                                                                                                                  >  2:535678/1‑140 (MQ=255)

          aataaGCTTTTTTTTTTTTTTTTTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTTAGATATATCACGATAACTCAAttcat                                                                                                                                 >  2:549609/1‑137 (MQ=11)

          aataaGCTTTTTTATTTTTTTTCTTCCAAATTTGATATAAAATTGTGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAAtccat                                                                                                                                 >  2:635205/1‑137 (MQ=11)

          aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACTCCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgcat                                                                                                                                 <  2:129491/141‑4 (MQ=11)

          aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgcat                                                                                                                                 <  2:47209/141‑4 (MQ=11)

          aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgcat                                                                                                                                 <  2:147434/141‑4 (MQ=11)

          aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGATCTCTTAATATTTCAGATATATCACGATAGCTCAATgcaa                                                                                                                                 >  2:685558/1‑138 (MQ=11)

          aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGGACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATgcat                                                                                                                                 >  1:495088/1‑138 (MQ=2)

          aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCCCGATAACTCATTGtaa                                                                                                                                 >  2:268896/1‑140 (MQ=255)

           ataaGCTTTTTTATGCTTTATATTGCACATTTGATATAAACTTGGGGGAGATTAGAGAACCCGAAGATGGATCGGTGAGTGATGAACGTTTAAGAGACGTTATCTTAATATTTCAGATACATCACGATAGCGCaagatatt                                                                                                                                >  2:204848/1‑135 (MQ=255)

                     ttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACAGCACGTTCACTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtag                                                                                                                      >  1:582588/1‑141 (MQ=25)

                     ttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtag                                                                                                                      >  2:337639/1‑141 (MQ=35)

                      tATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagt                                                                                                                     >  1:96053/1‑141 (MQ=35)

                      tATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagt                                                                                                                     >  1:467483/1‑141 (MQ=35)

                      tATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACAGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagt                                                                                                                     >  1:512139/1‑141 (MQ=25)

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

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

                         gCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGc                                                                                                                  >  1:564426/1‑141 (MQ=35)

                           tcttgcttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCa                                                                                                                <  1:241332/136‑1 (MQ=25)

                                   ccAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTTTATCTTAGATAGTAGCCAACGGCTAc                                                                                                        >  1:599114/1‑141 (MQ=21)

                                    ttttttttATACAAAATTGGGGCATAATCTTGTACCCAACGGTAGACCGTTGAATGGTGAACGTTTACACCACTTTTTTTTATTATTTCAGATATATGATGATAACTCAATTTATATTTTAGATAGTAGCCAACgggtgaa                                                                                                       >  2:495916/5‑136 (MQ=11)

                                           gATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGAAGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAAGGGCTACAGTgctca                                                                                                >  1:521911/1‑137 (MQ=255)

                                                  aaTTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATGCtt                                                                                        >  2:313297/1‑141 (MQ=17)

                                                    ttGGGGCATATTCTTGAACCCAACGGTATACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACAGCTACAGTG‑ATAACATCCTTGt                                                                                      <  1:447623/141‑1 (MQ=17)

                                                    tgggggCATATTCTTGAACCCACCCGTAGACCGTTGAATGATGAACCTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGt                                                                                      <  1:162038/139‑1 (MQ=12)

                                                       gggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTaa                                                                                   >  1:320265/1‑141 (MQ=32)

                                                       gggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTGCCCTTAATATTTCAGATATATCACGATAACTCAATTTAGATCTTAGATAGTAGCCAACGGCTACAGTG‑CTACGATCCTTGTTaa                                                                                   >  1:22184/1‑141 (MQ=11)

                                                        ggCCTATTCTTGAACCCAACGGTAGCCCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATACCACGATAACTCAATGTATGTCTTAGATAGTAGCCACCGGCCACAGTG‑ATAACCTCCCTGTTaca                                                                                  >  2:92202/1‑139 (MQ=11)

                                                        agCATATTCTTGAACCCAACGGTAGACAGTTGAATGATGAACGTTTACACCACGATCCCTTCATATTTCAGAAATATCACGATAACTCACTGTAAATCTTAGATCATGGACAACGCCT‑CAATGAATAACATCCTTGTTcaa                                                                                  >  2:92214/2‑141 (MQ=11)

                                                                cTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTCTATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGGTtat                                                                          >  2:2533/1‑141 (MQ=17)

                                                                    aaCCCATCCGTAGATTGTTTAATGATGAACGTTTACACGACGTCCTCCTAATATTTCAGATATATCACGATAGCTCAATGTATATCTTAGATAGTAGCCAACGTCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCt                                                                      <  2:537633/141‑1 (MQ=11)

                                                                      cccAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGa                                                                    >  2:394832/1‑141 (MQ=32)

                                                                                gACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCAt                                                          >  1:303520/1‑141 (MQ=32)

                                                                                     ttGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACaga                                                     >  2:773191/1‑141 (MQ=32)

                                                                                                aCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGATTCCtttt                                          >  2:581495/1‑141 (MQ=17)

                                                                                                    ttACACCACGTTCCTTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                      <  2:610387/141‑1 (MQ=17)

                                                                                                    ttACACCACCTTCCCTTAACATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                      <  1:28769/141‑1 (MQ=17)

                                                                                                    tgacacCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                      <  1:14448/139‑1 (MQ=17)

                                                                                                      acacCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaa                                    <  1:124833/141‑1 (MQ=32)

                                                                                                      acacCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaa                                    <  1:129676/141‑1 (MQ=32)

                                                                                                         ccacGTTCCCTTAATATTGCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAAt                                 <  1:23867/141‑1 (MQ=17)

                                                                                                            cGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATtat                              >  2:18857/1‑141 (MQ=32)

                                                                                                               tCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACt                           >  2:497075/1‑141 (MQ=32)

                                                                                                               tCCCTTAACATATCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACt                           <  2:507431/141‑1 (MQ=17)

                                                                                                                     aaTATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAAt                     <  1:632673/141‑1 (MQ=35)

                                                                                                                           tCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACt               >  1:440501/1‑141 (MQ=35)

                                                                                                                            cAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACtt              >  2:773679/1‑141 (MQ=35)

                                                                                                                              aatatatCACGATAACTCAATGGATATCTTAGATAGTAGCCAACAGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTg            <  2:144553/140‑1 (MQ=17)

                                                                                                                                tatatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAATTTATAATATAAATTTAACTTTGCa          >  2:471689/1‑141 (MQ=255)

                                                                                                                                tatatCACGATAACTCAACGGATATCGTAGATAGAAGCCAACAGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATGCAGAAGACTCCTTTTTGTTAAAATTATACTCTAAATTCAACTTTGCa          <  1:202110/141‑1 (MQ=11)

                                                                                                                                taaatCACGATAACTAAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCa          <  1:353346/138‑1 (MQ=21)

                                                                                                                                  tatCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATAAAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAAc        <  1:416965/141‑1 (MQ=255)

                                                                                                                                        gATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAccc  >  2:193791/1‑141 (MQ=25)

CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACC  >  CP000730/36331‑36607

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