BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,444	C→T	47.0%	R36R (AGG→AGA)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,444	0	C	T	47.0%	89.8 / inf	341	R36R (AGG→AGA)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base C (97/83); new base T (78/82); total (176/165)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.01e-01

Only 100 of 343 total aligned reads displayed.

TTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑CAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAA  >  CP000730/36304‑36575

tttATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACATTTACACGCCGTTcc                                                                                                                                     >  2:387680/1‑141 (MQ=11)

  tATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGATTTTTCTTCAAAATTTGATA‑TAAAATTGGGGAATATTCTTGAACCCAAAGGTAGACCGTTGAATGATGAACGTTTACAACTCGTTctct                                                                                                                                   >  2:651188/1‑141 (MQ=11)

   aTTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCtt                                                                                                                                  <  1:374041/141‑1 (MQ=11)

    ttttGATGTACGTCTCATCAATACTCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTa                                                                                                                                 >  1:425888/1‑141 (MQ=21)

    ttttGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTa                                                                                                                                 >  2:562608/1‑141 (MQ=25)

    ttttGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTa                                                                                                                                 <  1:397903/141‑1 (MQ=25)

     tttGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTaa                                                                                                                                <  2:505249/141‑1 (MQ=255)

     tttGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTATTTCTTCCAAATATGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACCCCACGTTCTCTTaa                                                                                                                                <  2:455169/141‑1 (MQ=11)

       tGATGTACGTCTCATCAATACTCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑CAAAATTGGGGCATATTCTTGAACCAAACGTTTGGCCGTTGGATGATGAACGTTTACACCACGTTCCCTTAatt                                                                                                                              >  1:345998/1‑140 (MQ=14)

        tatatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTCTGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAatat                                                                                                                             <  2:583887/137‑1 (MQ=11)

        gatatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAAAAAGACTCACATT‑TTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCCCGTTCTCTTAatat                                                                                                                             >  1:484322/5‑141 (MQ=11)

        gatatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGCTGAACGTTTACACCACGTTCTCTTActat                                                                                                                             >  2:532171/5‑141 (MQ=11)

        gatatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGACCGTTTACACCACGTTCTCTTAatat                                                                                                                             >  1:120386/5‑141 (MQ=11)

        gatatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCtt                                                                                                                                  >  1:220700/5‑136 (MQ=11)

        gatatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCtt                                                                                                                                  <  2:220700/132‑1 (MQ=11)

        gatatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAattt                                                                                                                             >  2:129083/5‑139 (MQ=11)

        gatatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAatat                                                                                                                             >  1:390802/5‑141 (MQ=255)

        gatatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAatat                                                                                                                             >  1:537424/5‑141 (MQ=255)

        gatatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAatat                                                                                                                             >  1:282564/5‑141 (MQ=255)

        gatatACGTCTCATCAATAAGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATGGTTGAATGATAAACGTTTAAACAACGTTCTCTTAatat                                                                                                                             >  2:361192/5‑141 (MQ=11)

        gATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAatat                                                                                                                             >  1:63959/1‑141 (MQ=21)

        gATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAatat                                                                                                                             >  1:112187/1‑141 (MQ=25)

        gATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAatat                                                                                                                             >  1:580205/1‑141 (MQ=25)

          tGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATAttt                                                                                                                           <  1:362340/141‑1 (MQ=18)

            tACGTCTCATCAATACGCCAGTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑CAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACTCCACTTTCCCTTAATATTTCa                                                                                                                         <  1:340930/141‑1 (MQ=255)

              cGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGCAGCCGACTGGATTCCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGa                                                                                                                       >  1:19385/1‑141 (MQ=11)

              cGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGa                                                                                                                       >  1:262620/1‑141 (MQ=255)

              cGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGa                                                                                                                       >  2:424581/1‑141 (MQ=25)

              cGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGCTAGATCGTTGAATGACGAACGTTTACACCACGTTTTCTTATATTTTCAGa                                                                                                                       >  2:456605/1‑141 (MQ=11)

              cGTCTCATCAATACCCCATTTGTAATAAGGTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGa                                                                                                                       <  2:154611/141‑1 (MQ=255)

               ggctcAGCAATACGCCATTTGTAATAAGCTGTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTGGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGat                                                                                                                      <  1:370730/139‑1 (MQ=255)

               gTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TACAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATTTTTCAGat                                                                                                                      >  1:166976/1‑141 (MQ=11)

               gTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGat                                                                                                                      <  1:635427/141‑1 (MQ=255)

               gTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGat                                                                                                                      <  2:638204/141‑1 (MQ=21)

               gTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑CAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCCCGTTCCCTTAATATTTCAGat                                                                                                                      >  2:610854/1‑141 (MQ=255)

               gTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTATTTTTATTTTGATA‑CAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATTTTTCAGat                                                                                                                      >  2:356719/1‑141 (MQ=14)

               ctctcATCAATACGCCATTTGTAATAAGCTTTTTTATTCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGat                                                                                                                      <  1:171406/140‑1 (MQ=11)

                 ctcATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAGCGTTTACACCACGTTCACTTAATATTTCAGatat                                                                                                                    >  2:14215/1‑141 (MQ=255)

                 ctcATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGatat                                                                                                                    <  1:54999/141‑1 (MQ=255)

                 ctcATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGatat                                                                                                                    <  1:310625/141‑1 (MQ=255)

                 ctcATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGatat                                                                                                                    >  2:614662/1‑141 (MQ=21)

                 ctcATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGatat                                                                                                                    >  2:575922/1‑141 (MQ=21)

                  tcatcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGatata                                                                                                                   <  1:603169/141‑1 (MQ=255)

                   catcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTTTTTCTTTTTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGatattt                                                                                                                  >  2:251408/1‑139 (MQ=11)

                   catcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGatatat                                                                                                                  >  1:32051/1‑141 (MQ=255)

                   catcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGatatat                                                                                                                  <  2:620900/141‑1 (MQ=255)

                   catcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGatatat                                                                                                                  >  1:555505/1‑141 (MQ=255)

                    atcaATACGCGATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCTACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATc                                                                                                                 <  1:246839/141‑1 (MQ=11)

                    atcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATc                                                                                                                 >  1:213816/1‑141 (MQ=21)

                    atcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATc                                                                                                                 >  1:460211/1‑141 (MQ=18)

                     tcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCa                                                                                                                <  1:32539/141‑1 (MQ=255)

                     tcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCa                                                                                                                >  1:46622/1‑141 (MQ=255)

                      cagtACGCCATTTGGAATGAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAGCGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCAc                                                                                                               <  1:589635/138‑1 (MQ=11)

                      caATACGCCATTTTTAATAAGCATTTTTAATCTTATTCTTCCAAAATTGATA‑TAAAATTGGGACATATTCTTAAACCCAAAAATAGAGCGTAAAATAATAAACGTTTACACCACGTTCTCTTAATATTTCAGATATATccc                                                                                                               >  2:238465/1‑139 (MQ=2)

                       aaTACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACg                                                                                                              <  2:166287/141‑1 (MQ=21)

                       aaTACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTATTGCACCCACCGGTAGATCGTAGGAAGATGAACGTTTACACCACGCTATCTTCATATTCCAGatatcataca                                                                                                              >  2:19142/1‑135 (MQ=11)

                        aTACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGa                                                                                                             <  1:563745/141‑1 (MQ=255)

                         tatgtctatcgTAATAAGCTTTTTTATGCTTTTTCCTCCAGATGTGATA‑TTAGAGTTGGGCCTATCCTTGAACACAAAGGTTGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGAt                                                                                                            <  2:329426/131‑1 (MQ=2)

                         tACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACAGTAGGTCGTTGTATGATTAACGATTTCACCACGTTCTCTTAATATTTTAGATATATCAAGaa                                                                                                            >  2:638953/1‑140 (MQ=11)

                         gaCGCCATGTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGAt                                                                                                            <  1:148813/140‑1 (MQ=11)

                          aCGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑CAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATa                                                                                                           <  2:106734/141‑1 (MQ=39)

                           cGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATag                                                                                                          >  1:143213/1‑140 (MQ=255)

                            gCCGTTTGTAATAAGCTTTTTTATTCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATagc                                                                                                         <  1:421008/141‑3 (MQ=11)

                            gCCATTTGTAATAATCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGACCCCAACGGTAGACCGTTGAATAATGAACGTTTACACCACGTTCCCTTAGTATTCCAGATATATCACGATAAc                                                                                                         <  1:128837/141‑1 (MQ=11)

                            gCCATTTGTAATAAGCTTTTTTAGGCGTTTTCTTCCAAATTTGATA‑TAAGATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATagc                                                                                                         <  1:299116/141‑3 (MQ=11)

                            cccATTTCTAATAAGCTTTTTTATGCTTCTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAATGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAAc                                                                                                         <  1:179016/140‑1 (MQ=255)

                             ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATagct                                                                                                        >  2:179311/1‑138 (MQ=255)

                             ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCt                                                                                                        >  1:268128/1‑141 (MQ=255)

                             ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACt                                                                                                        >  1:352206/1‑141 (MQ=25)

                              cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTc                                                                                                       <  2:66244/141‑1 (MQ=11)

                              cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTc                                                                                                       <  2:31454/141‑1 (MQ=255)

                              cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTc                                                                                                       <  1:228831/141‑1 (MQ=18)

                               aTTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAACTTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCa                                                                                                      >  2:206630/1‑141 (MQ=25)

                               aTTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTGGATA‑TAAAATTGGGGCATTTTCTTGAACCCAACGGTAGACCGTTGCATGATGAACGTTTACCCCACGGTCCTTTACTATTTCCGATATATCACGATAACTCa                                                                                                      >  1:94804/1‑141 (MQ=11)

                                tttGTAATCAGCTTTTTTATGCTTTTTCTTCCCAATTTGATC‑TACTATTATGGCATACTCCTGAGCCCAACGGTTGACCGTTGAATGATGAACGGTTACACCACGCTCCCTTAAAATTTCAGATATATCACGATAACTcca                                                                                                     >  2:304558/1‑139 (MQ=11)

                                 ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAAt                                                                                                    <  1:210572/141‑1 (MQ=21)

                                 ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAAt                                                                                                    <  1:445465/141‑1 (MQ=21)

                                 ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAAt                                                                                                    >  1:134731/1‑141 (MQ=21)

                                   gTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCCCGTTCCCTTAATATTTCAGATATATCACGATAACTCAATgc                                                                                                  >  2:121550/1‑140 (MQ=16)

                                   gTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTACTATTTCAGATATATCACGATAACTCAATgc                                                                                                  >  2:86735/1‑140 (MQ=12)

                                    taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgca                                                                                                 <  2:299070/141‑3 (MQ=11)

                                    taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATgca                                                                                                 >  2:541173/1‑139 (MQ=17)

                                      ataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑CAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGtcta                                                                                               >  2:378104/1‑141 (MQ=255)

                                             tttttATGCTTTTTCTTCCAAATTTGATA‑CAAAATTGGGGCATATTCTTGAACAAAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGa                                                                                        >  1:625013/1‑141 (MQ=25)

                                                ttATGCTTTTTCTTCCAAATTTGATA‑CAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtag                                                                                     <  1:100280/141‑1 (MQ=35)

                                                        tttATTGCAAATTTGATA‑CAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTACTATGTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAAc                                                                             <  1:157664/141‑1 (MQ=255)

                                                          tcttcCAAATTTGATA‑CAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACgg                                                                           <  2:77742/141‑1 (MQ=35)

                                                          tcttcCAAATTTGATA‑CAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACgg                                                                           >  2:547771/1‑141 (MQ=35)

                                                              ccAAATTTGATA‑CAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAAATAGTAGCCAACggatgc                                                                       >  2:478196/1‑137 (MQ=25)

                                                              ccAAATTTGATA‑CAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTAAAATATATCACGATAACTCAATGTATATCTTAGATGGTAGCCAACGGCTAc                                                                       >  1:286692/1‑141 (MQ=255)

                                                               cAAATTTGATA‑CAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACa                                                                      >  2:519093/1‑141 (MQ=35)

                                                                    ttGATA‑CAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTATAGTGAt                                                                 <  2:368021/141‑1 (MQ=25)

                                                                        tA‑CAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATAATTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCCACGGCTACAGTGCTAACa                                                             >  2:482998/1‑141 (MQ=255)

                                                                              aaTTGGTGCATATTCTTGAACCCAACGGTAGACCGTTGAATGAAGACCGTTTTCACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCtt                                                        <  2:270466/141‑1 (MQ=255)

                                                                                               gagcccAGCGGTAGCTCGTGGAGTGTGTAACGTTTACACTACGTTTTCTTAGTATGCTCCATATACCACGATAACTCAAGTTTTATCTTATATAGGAGCCAACAGCTACAGTGATAACATCCTTGTGAAATTGttggttgt                                       <  2:276829/138‑7 (MQ=2)

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

                                                                                                        ggTAGAGCGTTGAATGATGAACGTTTACAGCCCGTGCCCTAAATATTTCAGATATATCACGATAACACAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTCGTTAAATTGTTTATATCTGAAATAGt                              >  1:489431/1‑141 (MQ=11)

                                                                                                                gTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACag                      >  1:667046/1‑141 (MQ=32)

                                                                                                                gTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACag                      >  2:591781/1‑141 (MQ=32)

                                                                                                                                  acacCACGTTCCCTTAATATTTCAGATATATGCCGATAACTCAGTGTATGTCTTAGATAGTAGCCAACGGCTACAGTGATGACAGCCTTGTGAAATTGTTTATATCTGCAATAGTTCATACAGAATCAGCATTTTTGTTaa    >  2:282174/1‑141 (MQ=11)

TTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑CAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAA  >  CP000730/36304‑36575

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