BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,479	T→C	21.0%	T25A (ACA→GCA)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,479	0	T	C	21.0%	165.7 / 41.4	104	T25A (ACA→GCA)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base T (40/41); new base C (11/11); total (51/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

TAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAAT  >  CP000730/36340‑36619

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

 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)

                     ttcttcCAAATTTGATACAAAATTGGGGCATATTCTTGCACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACg                                                                                                                        >  1:429423/1‑141 (MQ=25)

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

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

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

                              aTTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTg                                                                                                               <  1:395460/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)

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

                                               ggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACAACCTTGGTaaa                                                                                              >  2:484384/1‑141 (MQ=17)

                                                   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)

                                                              ccAACGGCAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGaa                                                                               <  2:430878/141‑1 (MQ=17)

                                                              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)

                                                                      aGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATAGCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTATTTATATGTGCAATAGTTCa                                                                       >  2:170783/1‑141 (MQ=255)

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

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

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

                                                                           tttGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACag                                                                  <  2:495732/140‑1 (MQ=17)

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

                                                                                atgatgAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGAc                                                             >  2:641631/1‑141 (MQ=32)

                                                                                 tgatgaACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACt                                                            <  2:147724/141‑1 (MQ=17)

                                                                                 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)

                                                                                                       cccTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACATAAGACTCCTTTTTGTTAAAATTATACta                                      <  1:269732/141‑1 (MQ=17)

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

                                                                                                           tAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATaaa                                  <  1:636726/141‑1 (MQ=33)

                                                                                                                 ttCAGTTATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAAc                            <  2:53165/141‑1 (MQ=255)

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

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

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

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

                                                                                                                       aatatCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTTATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCa                      <  2:442934/140‑1 (MQ=255)

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

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

                                                                                                                         aatcacgttagtccAATGTATATCTTAGATAGTACCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATATTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAAc                    <  2:137670/128‑1 (MQ=11)

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

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

                                                                                                                              cGATAACTTAATGTATATCTTATATATGAGCCAACGGTTACAGTGATAACATCATTGTTAAATTGTTTATATCTGAAATAGTTCAAACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAc               <  2:31324/141‑1 (MQ=11)

                                                                                                                              cGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCGTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTCTACTATAAATTCAACTTTGCAACAGacc               >  2:11810/1‑139 (MQ=19)

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

                                                                                                                                agagctaaaTGCAGATCTTAGAAAGGAGCCCACGGCAACCGCTATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTCCATACAGAAGACGCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCg             <  1:123033/134‑1 (MQ=11)

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

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

                                                                                                                                 taggtCAATGTATATATTAGACAGTAGCCAACGGCTCCAGTGCTAAAATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctc            <  1:302909/137‑3 (MQ=11)

                                                                                                                                 tagcgcAATGCATATCTTAGATAGTAGCCAACGGCTACAGTTATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt            <  2:316804/136‑1 (MQ=255)

                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt            <  1:72600/138‑1 (MQ=12)

                                                                                                                                 cagcttaaTGCATATATTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTTCTTTTTGTTCAAATTATATTATAAATTCAACTTTGCAACAGAACCGt            <  1:24872/135‑1 (MQ=11)

                                                                                                                                 aagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAATCGt            <  1:480658/138‑1 (MQ=255)

                                                                                                                                 aagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTATTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt            <  1:399491/138‑1 (MQ=255)

                                                                                                                                  agctaaaTGCATATATTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTTAAATTATATTATAAATTCAACTTTGCAACATAACCGta           >  2:637716/6‑141 (MQ=11)

                                                                                                                                  agcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTATATTTATATTATAAATACAACTTTTCAACATAACCGta           >  2:170807/3‑141 (MQ=11)

                                                                                                                                  agcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGta           >  2:279844/3‑141 (MQ=9)

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

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

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

                                                                                                                                   gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACGTCCTTGTTAAATTGTTTATATCTGAAATAGGTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAAAACCGtat          >  1:155309/2‑141 (MQ=11)

                                                                                                                                   gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat          >  1:306638/2‑141 (MQ=12)

                                                                                                                                   gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat          >  2:220571/2‑141 (MQ=17)

                                                                                                                                   gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat          >  1:222228/2‑141 (MQ=12)

                                                                                                                                   gcTCAATCCATATCTTAGATGGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTAAACTTAGCAACAGAACCGtat          >  1:417661/2‑141 (MQ=11)

                                                                                                                                    cTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAATTTATACTATAAATTCAACTTTGCAACAGAACCGtatt         >  2:567772/1‑141 (MQ=255)

                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCACCAGAACCGtatt         >  1:125607/1‑141 (MQ=12)

                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt         >  1:183136/1‑141 (MQ=17)

                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATAGTTTATATCTGCAATAGTTCATACAGACGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGATCCGtatt         >  2:74901/1‑141 (MQ=11)

                                                                                                                                       aaTGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTTTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATg      >  2:245581/1‑141 (MQ=255)

                                                                                                                                        aTGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTCATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAAGTATGCTATAAGTTCAACTTTTCACCAGAACCGTATTATgg     >  2:23050/1‑141 (MQ=14)

                                                                                                                                          gTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCACCTTTGCAACAGAACCGCATTATGGaa   >  1:220429/1‑141 (MQ=255)

                                                                                                                                           tataTCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAAt  <  1:37897/141‑1 (MQ=255)

TAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAAT  >  CP000730/36340‑36619

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