BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,471	A→G	27.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	27.5%	157.5 / 79.4	129	S27S (AGT→AGC)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base A (46/44); new base G (8/27); total (56/73)
Fisher's exact test for biased strand distribution p-value = 4.84e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.23e-01

Only 100 of 133 total aligned reads displayed.

CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTA  >  CP000730/36331‑36610

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

  ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAATTTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACt                                                                                                                                           <  2:295109/141‑1 (MQ=21)

  ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAATATGGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCt                                                                                                                                           >  1:520323/1‑141 (MQ=11)

  ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTAAGATATATCACGATAGCt                                                                                                                                           >  2:480658/1‑141 (MQ=11)

   cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTc                                                                                                                                          <  1:436454/141‑1 (MQ=11)

    aTTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCa                                                                                                                                         >  1:492423/1‑141 (MQ=255)

    aTTTGTAATAAGCTTTTTTATGCTTCTACTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCa                                                                                                                                         <  2:343036/141‑1 (MQ=11)

     tttGTTATAAGCTTTTTTATGATTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCaa                                                                                                                                        <  1:8243/141‑1 (MQ=11)

     tttGTAATAAGCTTTTTTATGTTTTTTTTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCaa                                                                                                                                        <  1:292667/141‑1 (MQ=11)

     tttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTTATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCaa                                                                                                                                        <  1:87793/141‑1 (MQ=11)

     tttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAATGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCaa                                                                                                                                        >  2:424601/1‑141 (MQ=25)

     tttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCaa                                                                                                                                        >  1:258663/1‑141 (MQ=21)

     tttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCaa                                                                                                                                        <  2:650129/141‑1 (MQ=18)

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

     tttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCaa                                                                                                                                        >  1:285297/1‑141 (MQ=21)

     tttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTAGACCCAACGGTAGACCCCTGAATGATGAACGTTTACAACACGTTCCCTTCATATTTCAGACATATCACGATAACTCaa                                                                                                                                        >  2:36820/1‑141 (MQ=11)

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

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

      ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCATCGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAAt                                                                                                                                       <  2:238182/141‑1 (MQ=255)

       ttgaataaGCTTTTTTATGCGTTGTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAGCCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATg                                                                                                                                      <  1:120111/138‑1 (MQ=11)

       tGTAATAAGCTTTTTTTTGTGTTTTCTTCCAAATTTGATATAAAATGTGGGCATATTCTTGAACCCAACGGTAGAGCGTTGAATGATGAACGTTTACACCACGTTTTCTTAATATTTCAGATATATCACGATAGCTCAATg                                                                                                                                      <  2:51494/141‑1 (MQ=11)

       tGTAATAAGCTTTTTTATGGTGTTTCTTTCTACTTTGATATAAATTTGGGGTATATTTTTGAACCCAACGGTAGACCGTTGAATGATGGACGTTTTCACTCAGTTTCTTTAATTTTTCAGATATATTATGATAACTCTATg                                                                                                                                      >  2:504816/1‑141 (MQ=2)

       tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATACAATTGGGGCATATTCTGGAACCCAACGGTAGATCCTGGAATTATCAACTTTTACACCACGTTCTCTTAATATTGCAGATATATCACGATAGCTCAATg                                                                                                                                      <  2:184949/141‑1 (MQ=11)

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

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

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

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

       tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATg                                                                                                                                      >  2:158804/1‑141 (MQ=21)

       tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTCCCCTTAATATTTCAGATATATCACGATAACTCAATg                                                                                                                                      >  2:585587/1‑141 (MQ=25)

       tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATg                                                                                                                                      >  2:71810/1‑141 (MQ=40)

        gTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGACCGTTTACACCACGTTCTCTTAATATTTCGGATATCTCACGTTACCTCGATgc                                                                                                                                     >  1:583586/1‑140 (MQ=11)

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

        gTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAACTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTGCTCTTAATATTTCAGATATATCACGATAGCTCAATgc                                                                                                                                     <  2:342761/141‑2 (MQ=11)

         taataaGCTTTTTTTTGCTTGTTCTTTCAAATTTGATATGAAATTGGTGCATATTCTTGAACCCAACGGTAGACCTTTGAATGTTGAACCTTTACATCACCTTCCCTTAATATTTCAGATATATCACGATAACTCAATgca                                                                                                                                    <  2:276835/141‑3 (MQ=11)

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

         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgca                                                                                                                                    <  2:261319/141‑3 (MQ=11)

         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCGTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTGTTAATATTTTAGCTATATCACGATTGCTCAAttca                                                                                                                                    >  2:558031/1‑138 (MQ=11)

         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGta                                                                                                                                    >  2:64662/1‑141 (MQ=255)

         taataaGCTTTTTTATGCATTTTCTTCCAAATTTGATATAAAAGTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATgca                                                                                                                                    <  2:567901/141‑3 (MQ=255)

          aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgcct                                                                                                                                   >  2:121832/1‑138 (MQ=11)

          aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACGCCACGTTCCCTTAATATTTCAGCTATATCACGACAACTCAATGtat                                                                                                                                   >  1:51926/1‑141 (MQ=255)

          aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCCGGTATATCACGATAACTCAATGtat                                                                                                                                   >  1:608323/1‑141 (MQ=255)

                 ttttttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAg                                                                                                                            >  1:663701/1‑141 (MQ=35)

                  tttttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGa                                                                                                                           >  1:397241/1‑141 (MQ=35)

                     ttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACCCAATGTATATCTTACAtaa                                                                                                                        >  2:436815/1‑140 (MQ=21)

                     ttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTACCTTAATATTTCATGTATATCACGATAACTCAATGTATATCTTAGAtca                                                                                                                        >  2:464318/1‑139 (MQ=255)

                       aTGCTTTTTCTTCCAAATTTGATAGAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagta                                                                                                                      >  1:424761/1‑141 (MQ=17)

                       aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATCACCCACTGTATATCTTATAttgta                                                                                                                      >  2:48701/1‑141 (MQ=11)

                       aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagta                                                                                                                      <  2:663701/141‑1 (MQ=35)

                          cTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGcc                                                                                                                   <  1:44918/141‑1 (MQ=35)

                           tttttCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCa                                                                                                                  <  2:450838/141‑1 (MQ=35)

                                   ccAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTAc                                                                                                          <  1:417144/141‑1 (MQ=35)

                                      aacttGATACAAAACTGGGGCATATTTTTGACCCCAACGGTAGACCGTTGAATGACCAACGTTTACACCACGTTCCATTAATACTTCATATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAATGGCTACAGt                                                                                                       <  2:114589/138‑1 (MQ=11)

                                      aaTTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGt                                                                                                       <  1:202080/141‑1 (MQ=35)

                                            aTAGAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATACCTCAATGTATATCTTAGACAGTAGCCAACTGCTACAGTGATAAc                                                                                                 <  1:367385/141‑1 (MQ=255)

                                                aaaaTTGGGGCCTATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATcc                                                                                             <  1:439997/141‑1 (MQ=17)

                                                aaaaTTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATcc                                                                                             >  2:111524/1‑141 (MQ=32)

                                                    ttGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTGTATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGt                                                                                         <  1:141151/141‑1 (MQ=17)

                                                     tagggCATATTCTTGAACCCAGCGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGtt                                                                                        <  1:23369/139‑1 (MQ=11)

                                                     tGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGtt                                                                                        <  2:373085/141‑1 (MQ=32)

                                                      ggggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACCCGGCCCGCCCGTAATATTACAGTTAGATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTa                                                                                       >  1:321015/1‑141 (MQ=11)

                                                       gggCATATTCTTGTACCCAAGGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTATTATTTCTGATATATCTCGATACCTCAATGTATATTTTTGATAGTAGCCAACGGCTACAGTGATAACATCCTTGttta                                                                                      >  2:312935/1‑139 (MQ=11)

                                                       gggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGATTACACCACGTTCCCTTAATATTTCAGATATATCACCATAACTCCATCTATATCTTACATAGTAGCCAACGGCTACAGTGATAACACGCTTGTTaa                                                                                      >  2:365537/1‑141 (MQ=11)

                                                        ggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTaaa                                                                                     >  2:255511/1‑141 (MQ=32)

                                                            tatTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGt                                                                                 >  1:660108/1‑141 (MQ=32)

                                                             atTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGtt                                                                                <  1:255970/141‑1 (MQ=32)

                                                               tCTTGCACCTAAAGGTAGACCGTTGAATCATGAACGTTTACACCACGTTCCCTTACTATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAACTGTTta                                                                              <  2:213094/141‑1 (MQ=11)

                                                               tCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTta                                                                              <  1:111524/141‑1 (MQ=32)

                                                                cTTGAACCCAACGGTAGACCGTGGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTtat                                                                             >  2:488425/1‑141 (MQ=17)

                                                                      cccAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGa                                                                       <  2:484848/141‑1 (MQ=32)

                                                                       ccAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGaa                                                                      <  2:430892/141‑1 (MQ=32)

                                                                       ccAACCGTAGACCGTTGCATGACAATCGTTTGCCCCACGTTGTCTTTTTTTATCATATATATCACCTTTACTCAAGGTTTTTTTTACTTAGTATCCGATTTCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGaa                                                                      <  2:57934/141‑1 (MQ=2)

                                                                         caCGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCGTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAAt                                                                    <  1:548927/140‑1 (MQ=17)

                                                                             gTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGtt                                                                <  2:268986/141‑1 (MQ=32)

                                                                               gcaCCGTGGAATGATGAACGTTTACACCACGTTCCCTTAAGATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCa                                                              <  2:100760/139‑1 (MQ=12)

                                                                               aGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCCACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCa                                                              >  1:18724/1‑141 (MQ=17)

                                                                               aGACCGTTAAATGATGAACGTTTACACCACGTTCCCTTAATATATCAGATATACCCCGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGCTAACATCCTTGTTAAATTGTTAATATCTGAAATAGTTcc                                                              >  2:126430/1‑140 (MQ=11)

                                                                                gACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATCTCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCGGAAATGGTTCAt                                                             >  1:19393/1‑141 (MQ=12)

                                                                                gACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCAt                                                             >  2:404312/1‑141 (MQ=32)

                                                                                       gAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTTTATCAGAACTACTTCATCCggaaa                                                      >  1:77724/1‑140 (MQ=11)

                                                                                          tgatgaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACt                                                   >  1:185170/1‑141 (MQ=32)

                                                                                                  gTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTg                                           <  2:321015/141‑1 (MQ=32)

                                                                                                    ttACACCACGTTCCCTTAATATTTCAGATATATCCCGATAACTCGATGTATATCTTAGATAGTAGCCCAAGATGACAGTGATACCACGCGGGTGAAATTTTTTATATCTGAAATAGCTCATACAGAAGACATCTTTTTGtt                                         <  1:488425/141‑1 (MQ=11)

                                                                                                     tACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTa                                        >  1:110000/1‑141 (MQ=32)

                                                                                                           acGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATta                                  <  2:293047/141‑1 (MQ=32)

                                                                                                           acGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATta                                  <  2:519267/141‑1 (MQ=32)

                                                                                                             gTTCCCTTAATTCTTCAGATATACCACGATACCTTAATGTATATCTTAGCTAGAAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATtata                                <  1:349663/141‑1 (MQ=11)

                                                                                                                   ttAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGCAGACTCCTTTTTGTTAAAATTATACTATaa                          <  1:140535/141‑1 (MQ=17)

                                                                                                                           tCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGCTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACt                  <  2:596093/141‑1 (MQ=255)

                                                                                                                              gATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTg               >  1:105887/1‑141 (MQ=35)

                                                                                                                               atatatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGc              >  1:151888/1‑141 (MQ=35)

                                                                                                                               atatatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGc              <  2:151888/141‑1 (MQ=35)

                                                                                                                                  tatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAAc           >  1:199032/1‑141 (MQ=35)

                                                                                                                                    tCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGGAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAg         <  2:429423/141‑1 (MQ=25)

                                                                                                                                       cgctAACTCAATGCATATCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCGTGTGAAAGCGTTTATATCTGAAATAGTTCATACAGAAGCTTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAc      <  2:175517/138‑1 (MQ=1)

                                                                                                                                         aTAGCTGAATGCATAGCTTATATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAATACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCg    <  2:13521/141‑1 (MQ=11)

                                                                                                                                         aTAACTCAATGACTAGCTTAGATAGTAGCCAACAGCTACAGTGATAACATTCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGGCTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAcct    <  2:436827/141‑2 (MQ=11)

                                                                                                                                           aaCTCAATGTATATCTTAGATAGTAGCCAACGGCTCCAGTGCTAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCACCAGAACCGta  >  1:629123/1‑141 (MQ=255)

                                                                                                                                           aaCTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGta  >  2:631284/1‑141 (MQ=38)

                                                                                                                                           aaCTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGta  <  2:185170/141‑1 (MQ=37)

CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTA  >  CP000730/36331‑36610

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