BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	35,997	T→C	31.3%	T185T (ACA→ACG)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	35,997	0	T	C	31.3%	1153.8 / inf	1009	T185T (ACA→ACG)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base T (347/342); new base C (160/154); total (513/496)
Fisher's exact test for biased strand distribution p-value = 8.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 1012 total aligned reads displayed.

CAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTT  >  CP000730/35867‑36137

cAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAGGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATAc                                                                                                                                     <  2:346327/141‑1 (MQ=2)

cAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATAc                                                                                                                                     <  1:357173/141‑1 (MQ=2)

 aGTGTTCGTTTAACTTGTTAGGATGATGCTAATTTCGTGGCATGGTGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACt                                                                                                                                    <  2:217599/141‑1 (MQ=11)

 aGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACt                                                                                                                                    <  1:336484/141‑1 (MQ=11)

 aGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACt                                                                                                                                    <  1:521340/141‑1 (MQ=255)

  gtgtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACtt                                                                                                                                   >  1:430628/1‑141 (MQ=17)

       cgtttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATa                                                                                                                              >  1:499051/4‑141 (MQ=255)

       cgtttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATa                                                                                                                              >  1:139743/4‑141 (MQ=255)

        gtttAGCTGGCTAGCACGATGCTAATGTCGTGGCATGCGGAAAATCCGTAGATCTTAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATat                                                                                                                             <  1:507247/139‑2 (MQ=11)

         tttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATAcc                                                                                                                            <  1:571000/140‑1 (MQ=17)

            aaCTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTg                                                                                                                         <  2:539906/141‑1 (MQ=17)

               tggCTAGCATGATGCTAATTTCCTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGTTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtct                                                                                                                      <  2:422744/139‑1 (MQ=12)

                tGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGtctt                                                                                                                     <  1:148359/141‑1 (MQ=11)

                tGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctt                                                                                                                     <  1:101021/141‑1 (MQ=17)

                tGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctt                                                                                                                     <  1:365024/141‑1 (MQ=17)

                tGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctt                                                                                                                     >  1:418161/1‑141 (MQ=17)

                 gCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGtcttt                                                                                                                    >  2:646660/1‑141 (MQ=255)

                 gCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGtcttt                                                                                                                    >  2:126053/1‑141 (MQ=255)

                 gCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtcttt                                                                                                                    >  1:240026/1‑141 (MQ=17)

                  cTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGCCCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttt                                                                                                                   >  2:392307/1‑140 (MQ=17)

                  cTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttc                                                                                                                   <  2:471382/141‑1 (MQ=17)

                   tAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttct                                                                                                                  <  2:28014/141‑1 (MQ=32)

                   tAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGGATTCTTTGCTGTATTGATACTTTGATACCTTGtctttct                                                                                                                  >  1:527407/1‑141 (MQ=17)

                       atgatgATAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTTCCGTATTGATACTTTGATATCTTGTCTTTCTTACt                                                                                                              <  2:608502/141‑1 (MQ=11)

                        tgatgCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACtt                                                                                                             >  1:439158/1‑141 (MQ=11)

                         gatgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACttt                                                                                                            <  2:291389/141‑1 (MQ=25)

                          atgCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTa                                                                                                           <  1:111775/141‑1 (MQ=11)

                             cTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTCTATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATAGCGTGTCTTTCTTGCTTTAata                                                                                                        >  1:448042/1‑141 (MQ=11)

                              tAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAatat                                                                                                       >  2:511774/1‑141 (MQ=25)

                                aTTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGa                                                                                                     >  1:250744/1‑141 (MQ=35)

                                     ggggCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTACATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTg                                                                                                <  1:218834/139‑1 (MQ=14)

                                     gTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTg                                                                                                <  1:78893/141‑1 (MQ=255)

                                     gTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTg                                                                                                <  2:635871/141‑1 (MQ=255)

                                        gCATGTCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑Tc                                                                                            <  2:595035/141‑1 (MQ=255)

                                          aTGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCtt                                                                                          <  1:564041/141‑1 (MQ=255)

                                           tGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTTCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTg                                                                                         <  2:58702/141‑1 (MQ=255)

                                           tGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTg                                                                                         <  2:496685/141‑1 (MQ=38)

                                           tGGCGAAAATCCGTAGATATGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TATTg                                                                                         <  2:174069/141‑1 (MQ=255)

                                              cGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTc                                                                                      >  1:275268/1‑141 (MQ=38)

                                                aaaaTCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTCAGTCCGCTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCaa                                                                                    <  1:288898/141‑1 (MQ=14)

                                                 aaaTCCGTAGATCTAAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGACTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAAt                                                                                   >  1:369572/1‑141 (MQ=255)

                                                 aaaTCCGTAGAACTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGCTACTTTGCTATCTTGTCTTTCTTCCTTTAATATGACGGTGA‑TTTTGCTCAAt                                                                                   >  2:376348/1‑141 (MQ=11)

                                                  aaTCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATg                                                                                  <  2:3164/141‑1 (MQ=255)

                                                   aTCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGa                                                                                 >  1:577370/1‑141 (MQ=37)

                                                   aTCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGa                                                                                 >  2:148379/1‑141 (MQ=255)

                                                      cGTAGATCTGAAGAGACCTGCGGTTCATTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTCATACTTTGATACCTTGTCTTTCTTACTTTAATATGACTGTGA‑TCTTGCTCAATGAGGt                                                                              <  1:403134/141‑1 (MQ=255)

                                                      cGTAGAACTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGt                                                                              >  2:368680/1‑141 (MQ=255)

                                                       gTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGtt                                                                             <  1:659088/141‑1 (MQ=255)

                                                        tAGATCTGAAGAGACCTACGGTTCTTGTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTa                                                                            <  2:158951/141‑1 (MQ=255)

                                                            tCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTAttaa                                                                        >  2:246451/1‑139 (MQ=37)

                                                            tCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTAATACTTTGATCCCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCa                                                                        >  1:192628/1‑141 (MQ=255)

                                                              tGAAAAGACCTGCGGTTCTTTTTATATTGTGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGTTACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGa                                                                      <  2:567056/141‑1 (MQ=255)

                                                               gAAGAGACCTACGGTTCTTTTTATATTTAGTTTAAATACATTCAATACCTTTTAAAGTATTGTTTGCCGTCTTGATACTTTGTTATCTTTTCTTTCTTCCTTTAATATGAC‑GTGAGACTTGCTCAATATTGTTGTTCAGag                                                                     >  2:409436/1‑140 (MQ=2)

                                                                aaGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCatagt                                                                    <  2:79771/141‑5 (MQ=255)

                                                                aaGAGACCTACGATACTGTTTATATAGAGCGTAAATACATTCAATGCCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATGTTACGGTGG‑TCTTGCTCAATGAGAGTATTCAGata                                                                    >  2:226368/1‑141 (MQ=11)

                                                                   agaCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATAttt                                                                 >  2:485117/1‑141 (MQ=255)

                                                                    gaCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGT‑TCTTTCTCAATGAGGTTATTCATATAtttt                                                                >  1:317812/1‑140 (MQ=255)

                                                                      ccTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGa                                                              <  2:62316/141‑1 (MQ=38)

                                                                      ccTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTGGCTCAATGAGGTTATTCAGATATTTCGa                                                              >  1:123105/1‑141 (MQ=255)

                                                                        tacGGTTCTTTTTATATAGAGCGTAAATACAATCAACACCTTCTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATg                                                            <  1:16470/139‑1 (MQ=11)

                                                                            gTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACa                                                        <  2:269120/141‑1 (MQ=38)

                                                                              tCTTTTTATATAGAGCGTAAATAACATCAATACCTTTTAAATTATTCTTTTCTTTATTGATACTTTGATACCTTGTCTTTCTTACTTGAATATGACGGTTA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAAt                                                      <  2:628657/141‑1 (MQ=11)

                                                                                 ttttATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGAc                                                   >  1:626129/1‑141 (MQ=255)

                                                                                  tttATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGGTATTTCGATGTACCATGGCa                                                  <  2:522653/141‑1 (MQ=255)

                                                                                    tatataGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGt                                                <  1:380637/141‑1 (MQ=255)

                                                                                     atataGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTc                                               >  1:90334/1‑141 (MQ=255)

                                                                                        aagagCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAgg                                            <  1:603517/140‑1 (MQ=255)

                                                                                         agagCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTTATACTTTGCTACCTTGTCTTTCTTACTTTAATATGACGGTGA‑GCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgt                                           >  2:434277/1‑141 (MQ=255)

                                                                                           ggcttAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgttt                                         <  2:657031/137‑1 (MQ=255)

                                                                                           agCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgttt                                         >  1:390023/1‑141 (MQ=255)

                                                                                            gCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgttta                                        <  2:360809/141‑1 (MQ=255)

                                                                                             cGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCCTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaa                                       <  1:459740/141‑1 (MQ=255)

                                                                                               tAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCCGgtttaagt                                     >  2:101851/1‑141 (MQ=255)

                                                                                               tAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagt                                     >  2:446674/1‑141 (MQ=255)

                                                                                                 aatacacgccagaCCTTTTAAAGTATGCTGGGCTGTATTGGTACTTTGATCCCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTAGAATGGTAGTCAGgtttaagttt                                   >  2:576185/13‑141 (MQ=14)

                                                                                                 aaTACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTGTGATATGACGGTGA‑GCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgtttaagttt                                   <  2:554976/141‑1 (MQ=255)

                                                                                                 aaTACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTGTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagttt                                   >  2:189656/1‑141 (MQ=255)

                                                                                                 aaTACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagttt                                   >  1:228990/1‑141 (MQ=255)

                                                                                                 aaTACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagttt                                   >  1:614731/1‑141 (MQ=255)

                                                                                                  aTACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagttta                                  >  1:64149/1‑141 (MQ=255)

                                                                                                   tACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATATTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCTGgtttaggtttaa                                 >  2:332770/1‑141 (MQ=11)

                                                                                                   tACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgtttaagtttaa                                 >  1:175630/1‑141 (MQ=255)

                                                                                                   tACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagtttaa                                 <  1:593294/141‑1 (MQ=255)

                                                                                                     cactCAATACCTTTTAAAGTATTCTTTGCTGTATTGATCCTTTGCTAAATTGCCTTTATTACCTTAATATGATGGGGA‑TCTGGCGCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCGGGTTTAAGTTTaaaa                               >  1:132733/4‑141 (MQ=14)

                                                                                                      actCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATTCCGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCTATGTACAATGGCAGTCAGgttgtagtttacatg                              >  1:207421/3‑137 (MQ=14)

                                                                                                      actCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGCTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCCGGTTTACGTTTAAAAg                              >  1:295053/3‑141 (MQ=255)

                                                                                                          aaTACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCttt                          >  1:490856/1‑141 (MQ=255)

                                                                                                            ttccTTTTAAATTATTCTTTGCTGAATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTaa                        <  2:667934/139‑1 (MQ=255)

                                                                                                             aCCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAAt                       >  1:536084/1‑141 (MQ=255)

                                                                                                             aCCTTTTAAAGTATTCTTTGCAGTTTTAAGGCTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAAt                       >  2:305697/1‑141 (MQ=14)

                                                                                                                ttttAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTAc                    <  2:205563/141‑1 (MQ=255)

                                                                                                                  ttAAAGAATTCATTGCCGTATTGATCCTTTGATATCTTATCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTAttt                  >  2:255487/1‑141 (MQ=14)

                                                                                                                    aaaGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTa                >  2:45613/1‑141 (MQ=255)

                                                                                                                       gTATTCTTTGCTGTATTGCTGCTTTGATACCTTGTCTTTCTTAGTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTGAGATTTTTCGATGTACAATGGGAGTCAGGTTTAAGTTTGAACGCTTTAATTACTTTaatc             >  2:91227/1‑138 (MQ=11)

                                                                                                                        tATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAGAGCTTTAATTACTTTAAcct            >  1:547455/1‑140 (MQ=255)

                                                                                                                         aTTCTTTGCCGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCAt           <  1:262590/141‑1 (MQ=255)

                                                                                                                          ttCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCAtt          >  2:112129/1‑141 (MQ=255)

                                                                                                                            cTTTGCTGTATTGATACTTTGATACCTTGGCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTTAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGc        <  1:155151/141‑1 (MQ=255)

                                                                                                                               tGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTCGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTAc     >  2:382368/1‑141 (MQ=255)

                                                                                                                                 cTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCGTGCTCAATGAGGTTATTCAGATATTGCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCt   <  1:477974/141‑1 (MQ=255)

CAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGA‑TCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTT  >  CP000730/35867‑36137

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