BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,471	A→G	28.3%	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	28.3%	167.4 / 87.7	128	S27S (AGT→AGC)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base A (39/52); new base G (10/26); total (49/79)
Fisher's exact test for biased strand distribution p-value = 1.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-02

Only 100 of 132 total aligned reads displayed.

CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTA  >  CP000730/36331‑36610

cGCCATTCGTAATAAGCTTTATTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATaa                                                                                                                                             >  1:488685/1‑141 (MQ=255)

 gCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAAc                                                                                                                                            >  2:405368/1‑141 (MQ=40)

 gCCATTTGTAATAAGCTTTTTTATGCTTTTTCTGCCAAATTTGATATAAAATTCGTGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAAc                                                                                                                                            <  2:425486/141‑1 (MQ=255)

  ccATTTGTAATAAGCTTTTATATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACGCCACGTTCTCTTAATATTTCAGATATATCACGATAGCt                                                                                                                                           >  1:503297/1‑141 (MQ=11)

  caaTTTGTAATAAGCTTTTTTATGCTTTTTCTTCGAAATTTGATATAAAAGTGGGCTACATTCTTCAACCCAACGCTGGAGCTTCGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCt                                                                                                                                           <  2:192704/139‑1 (MQ=11)

   cATTTGTAATATGCTTTTTTATGCCTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTc                                                                                                                                          <  1:589543/141‑1 (MQ=255)

   cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAACTTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATTTATCACGATAGCTc                                                                                                                                          <  1:428305/141‑1 (MQ=11)

   cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTc                                                                                                                                          <  1:413870/141‑1 (MQ=255)

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

   cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATTAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTc                                                                                                                                          <  2:85226/141‑1 (MQ=25)

   cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTc                                                                                                                                          <  2:141957/141‑1 (MQ=21)

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

   cATTTGTAATAAGCTATTTTATGCTTTTTCTTCCAAAATTGATCCAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAATGTTTACACCACGTTACCTTAATATTTCAGATATCTCACGATGACTc                                                                                                                                          <  1:545325/141‑1 (MQ=11)

     tttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGAGCGCTGAACAACTTAGGTTTACACCACGTTCCCTTACTATTTCAGATATATCACGATAACTCaa                                                                                                                                        >  1:145314/1‑141 (MQ=11)

     tttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGAtcaatcaa                                                                                                                                        >  2:471590/1‑141 (MQ=255)

      ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGTTATATCACGATAGCTCAAt                                                                                                                                       <  1:239669/141‑1 (MQ=11)

      ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGCTATATCACGATAGCTCAAt                                                                                                                                       <  1:411916/141‑1 (MQ=11)

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

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

      ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGACAGCCCAAt                                                                                                                                       <  2:769064/141‑1 (MQ=11)

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

      ttGTAATAAGCTTTTTTATGCTTTTTCTGCCAAATTTGATATAAAATTGGGGCACATTCTTGAACCCAACGGTAGATCGATGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAAt                                                                                                                                       >  1:500010/1‑141 (MQ=11)

      ttGTAATAAGCTTTTTTATGCTTTTTATTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGTATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAAt                                                                                                                                       <  1:472976/141‑1 (MQ=11)

      tagTAATAAGCTTTTTTATGCTTTTTCTTCCGAATTTGATACAAAATTGGGGCATCTTCTTGAACCCAACGGTAGACCGTTAAATGTTGAACGTTTACACCACGTTCCCTTAATATTTCAGATGTATCAAGATAGCTCAAt                                                                                                                                       >  2:168011/3‑141 (MQ=14)

       tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGTACTTTGACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATg                                                                                                                                      <  1:448306/141‑1 (MQ=11)

       tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTTATATTTCAGATATATCACGATAGCTCAATg                                                                                                                                      <  1:202443/141‑1 (MQ=11)

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

       tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATTAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATg                                                                                                                                      <  1:32052/141‑1 (MQ=21)

       tGTAATAAGCGTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATg                                                                                                                                      <  2:112263/141‑1 (MQ=11)

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

        gTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTACTATTTCAGATATATCACTATAACTCAATga                                                                                                                                     >  2:234135/1‑140 (MQ=255)

        gTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGt                                                                                                                                     >  2:241952/1‑141 (MQ=255)

        gTAATAAGCTTTGTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATgc                                                                                                                                     >  2:241184/1‑140 (MQ=17)

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

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

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

         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTTCACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATgca                                                                                                                                    >  1:440904/1‑139 (MQ=12)

         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATgca                                                                                                                                    >  2:246801/1‑139 (MQ=17)

         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACCCAATgca                                                                                                                                    <  1:219323/141‑3 (MQ=12)

         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGGTGACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATgca                                                                                                                                    <  2:566561/141‑3 (MQ=255)

         taataaGCTTTTTTATGCTTTTTCCTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGta                                                                                                                                    <  1:31446/141‑1 (MQ=255)

           ataaGCTTTTTTATGCTTTTTTTTTTTTATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGGCCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGtttt                                                                                                                                  >  2:350404/1‑138 (MQ=14)

           ataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGtata                                                                                                                                  <  2:268758/141‑1 (MQ=40)

                         ccTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCCTAATATTTCAGATCTATCACGATAACGCAATGTATAATTTAAAtagtaac                                                                                                                    >  2:530750/2‑139 (MQ=11)

                            ttttCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCaa                                                                                                                 <  1:133969/141‑1 (MQ=35)

                               tcttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACgg                                                                                                              <  1:694047/141‑1 (MQ=35)

                                cttcCAAACTTGATACAAAATTGGGGCATATTCTTGACCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGc                                                                                                             <  1:390984/141‑1 (MQ=25)

                                 ttcCAAATTTGATACAAAATTGGGGCATATTCTGGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGcc                                                                                                            >  1:263856/1‑140 (MQ=25)

                                   ccAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTAc                                                                                                          >  1:706936/1‑141 (MQ=35)

                                   ccAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTAc                                                                                                          >  1:629450/1‑141 (MQ=35)

                                         ttGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGAt                                                                                                    <  1:282203/141‑1 (MQ=35)

                                            ctACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAAc                                                                                                 <  1:361690/140‑1 (MQ=25)

                                            aTACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATTAACGTTTACACACCGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGGCAACGGCTACGGTGATaaa                                                                                                 >  1:24037/1‑140 (MQ=11)

                                            aTACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAAc                                                                                                 >  2:205029/1‑141 (MQ=35)

                                             tACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTGTATCTTAGATAGTAGCCAACGGCTACAGTGATAACa                                                                                                <  2:708576/141‑1 (MQ=25)

                                                    ttGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTAAATGTATATCGTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGt                                                                                         <  1:729662/141‑1 (MQ=11)

                                                     tGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGtt                                                                                        >  2:701472/1‑141 (MQ=32)

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

                                                         gCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTACCCTTAATATTTCAGATATATCACGAAAACTCAATGTATATCTAAAATAGTAGCCAACGGCTACAGTGATAACATCCTGGTTAAAt                                                                                    >  2:149335/1‑141 (MQ=11)

                                                             atTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACACCCTTGTTAAATTGtt                                                                                >  2:538583/1‑141 (MQ=17)

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

                                                               tCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAAGGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTta                                                                              <  1:308156/141‑1 (MQ=17)

                                                                  tGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTatat                                                                           >  1:224831/1‑141 (MQ=32)

                                                                  tGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTatat                                                                           >  2:146619/1‑141 (MQ=32)

                                                                    aaCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCt                                                                         <  1:179499/141‑1 (MQ=32)

                                                                       ccAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATCTCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGaa                                                                      >  2:403392/1‑141 (MQ=17)

                                                                        cAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGaaa                                                                     <  2:675258/141‑1 (MQ=32)

                                                                        cAACGGTAGACCGTTGAATGATGAACGTTTAAACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGaaa                                                                     <  1:140935/141‑1 (MQ=17)

                                                                         aaCGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAAt                                                                    >  2:488759/1‑141 (MQ=32)

                                                                            ggTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATCACATCCTTGTTAAATTGTTTATATCTGAAATAGt                                                                 >  1:66696/1‑141 (MQ=17)

                                                                            ggTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGt                                                                 >  1:64938/1‑141 (MQ=32)

                                                                               aGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACCGTGATAACATCCTTGTTGAATTGTTTATCGCTGAAATAGTTCa                                                              >  1:281306/1‑141 (MQ=255)

                                                                                       gAATGATGAACCTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaag                                                      <  1:720099/141‑1 (MQ=17)

                                                                                          tgatgaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACt                                                   <  2:263856/141‑1 (MQ=32)

                                                                                               aaCGTTTACACCACGTTCCCTTAATATTTCAGGTATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCATTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttt                                              <  1:530750/141‑1 (MQ=17)

                                                                                               aaCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttt                                              <  1:684453/141‑1 (MQ=32)

                                                                                               aaCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttt                                              <  1:687007/141‑1 (MQ=32)

                                                                                                    ttACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                         >  2:39922/1‑141 (MQ=32)

                                                                                                    ttACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                         <  1:157174/141‑1 (MQ=32)

                                                                                                        accacGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaaaa                                     <  2:461998/141‑1 (MQ=32)

                                                                                                         ccacGTTCCCTTAATATTTCCACTATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACACCCTTGTTAACTTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAAt                                    <  2:316125/141‑1 (MQ=11)

                                                                                                          cacGTTCCCTTAATATTTCAGATATATCACGATACTTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAAtt                                   <  1:596466/141‑1 (MQ=17)

                                                                                                            cGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATtat                                 >  1:214808/1‑141 (MQ=32)

                                                                                                              ttCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATAc                               <  1:620080/141‑1 (MQ=32)

                                                                                                                 ccTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACtat                            >  1:645136/1‑141 (MQ=32)

                                                                                                                     aaTATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAGATTATACTATAAAt                        <  1:403392/141‑1 (MQ=21)

                                                                                                                      atatTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAAtt                       <  1:701472/141‑1 (MQ=34)

                                                                                                                        ttttCCGTTATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCa                     <  1:471590/140‑1 (MQ=25)

                                                                                                                        tgttCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCa                     <  2:267679/139‑1 (MQ=25)

                                                                                                                             aGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACttt                <  2:659495/141‑1 (MQ=35)

                                                                                                                             aGATATATCACGAAAACCCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACttt                <  2:582762/141‑1 (MQ=25)

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

                                                                                                                              aatatatCACCATCAACCAACGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTg               <  2:339164/140‑1 (MQ=11)

                                                                                                                               atatatCACGATAGCTCAATGCATATCTTAGATAGTAGCCAATGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGc              <  1:37123/141‑1 (MQ=255)

                                                                                                                                tatatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACACCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCa             <  1:538583/141‑1 (MQ=25)

                                                                                                                                  tatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTTATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACACAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAAc           >  2:670845/1‑141 (MQ=25)

                                                                                                                                    tCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAg         >  1:263648/1‑141 (MQ=35)

                                                                                                                                       cGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAc      <  1:241952/141‑1 (MQ=35)

                                                                                                                                       cGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAc      <  1:146619/141‑1 (MQ=34)

                                                                                                                                         aTAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCg    <  2:613882/141‑1 (MQ=9)

CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTA  >  CP000730/36331‑36610

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