BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	35,997	T→C	26.7%	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	26.7%	843.8 / inf	594	T185T (ACA→ACG)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base T (230/203); new base C (92/66); total (324/270)
Fisher's exact test for biased strand distribution p-value = 3.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 594 total aligned reads displayed.

CAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCT  >  CP000730/35867‑36136

 aGTGTTCGTTCAACCTGCTAGCATGATGCTAATTCCCTAGCCAAGCGAAAACAACTCGTTCTGTAGAGATCTGCGGTCCTTTTTATATAGAGCGTAACTACACTCAATACCGTTTAAAGTATTCTTTGCTGTAATGATACt                                                                                                                                   >  2:110901/1‑141 (MQ=2)

  gtgtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGGAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACtt                                                                                                                                  >  1:211052/1‑141 (MQ=17)

    gctCGCTTACCTTGGTAGCAGTATGGTAATTTCGTCGCATGGCGAAAACCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCCGAAACACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTg                                                                                                                                <  2:21708/139‑1 (MQ=11)

    gatccgttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTg                                                                                                                                <  2:135715/135‑1 (MQ=12)

      tCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATGCTTTGAt                                                                                                                              >  1:29749/1‑141 (MQ=17)

         tttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATatc                                                                                                                           <  1:17934/140‑3 (MQ=255)

           tAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCtt                                                                                                                         >  2:619131/1‑141 (MQ=17)

           tAACATGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCtt                                                                                                                         >  1:471522/1‑141 (MQ=17)

             aCTCGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGt                                                                                                                       <  2:68554/141‑1 (MQ=11)

              cTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGtc                                                                                                                      <  1:130829/141‑1 (MQ=11)

               tggCTAGCTTCATGCTAATTTCGTGGCATGTCGAAACTCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtct                                                                                                                     <  1:317054/139‑1 (MQ=11)

                  cTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGtctttt                                                                                                                  >  2:23780/1‑140 (MQ=255)

                  cTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCGATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttc                                                                                                                  <  1:427960/141‑1 (MQ=17)

                  cTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttc                                                                                                                  >  2:655252/1‑141 (MQ=17)

                     gCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTa                                                                                                               >  1:278941/1‑141 (MQ=17)

                     gCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACATGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGGTACCTTGTCTTTCTTa                                                                                                               >  1:246290/1‑141 (MQ=12)

                         gatgCTAATTTCGTGGCATGGGGGAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACttt                                                                                                           <  2:344277/141‑1 (MQ=11)

                         gatgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACttt                                                                                                           >  1:638843/1‑141 (MQ=21)

                         gatgCTAACTTCGGCGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACttt                                                                                                           <  1:38257/141‑1 (MQ=255)

                            gCTAATTTCGTGGCATGGGGAAAATCCGTAGATATGCACAAACTACATACTACTTGCGCAAGGAGCGTAAATCAATTCAATACCTTTTAAAGTATTTTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAat                                                                                                        <  2:304882/141‑1 (MQ=2)

                             cTAATTGCGTGGCATGGCGAAAAGCCGTAGATCTGAAGAGGCCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTGTGATACCTTGTCTTTCTTACTTTAata                                                                                                       <  2:504750/141‑1 (MQ=255)

                              tAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGCTATCTTGTCTTTCTTAATTTAatat                                                                                                      >  1:361332/1‑141 (MQ=11)

                                    cGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACAGt                                                                                                >  2:593156/1‑141 (MQ=11)

                                     gTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTg                                                                                               <  1:255109/141‑1 (MQ=39)

                                      tGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATAACGACTCTCTTACTTTAATATGACGGTGa                                                                                              <  1:324777/141‑1 (MQ=14)

                                         cATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTCACTTTAATCTGACGGTGATCt                                                                                           <  1:348161/141‑1 (MQ=255)

                                          aTGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGAGAGCtt                                                                                          >  2:543198/1‑141 (MQ=255)

                                           tGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTg                                                                                         <  2:373322/141‑1 (MQ=38)

                                           tGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTg                                                                                         >  2:568524/1‑141 (MQ=37)

                                              ggAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTCATATGACGGTGATCTTG‑CTc                                                                                     >  1:152169/2‑141 (MQ=255)

                                               gAAAATCCGTAGATCCCAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCa                                                                                    <  1:50338/141‑1 (MQ=14)

                                                 aaaTCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTC‑CTCAAt                                                                                  >  2:53616/1‑141 (MQ=255)

                                                  aaTCCGTAGAGCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATg                                                                                 <  2:481024/141‑1 (MQ=255)

                                                      cGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGAACTTG‑CTCAATGCGGt                                                                             >  2:380495/1‑141 (MQ=255)

                                                          gATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTAtt                                                                         <  1:412785/141‑1 (MQ=255)

                                                           aTCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTc                                                                        <  2:56202/141‑1 (MQ=38)

                                                            tCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCa                                                                       <  1:529422/141‑1 (MQ=39)

                                                            tCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCa                                                                       <  2:348242/141‑1 (MQ=255)

                                                             cTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAg                                                                      <  2:259566/141‑1 (MQ=38)

                                                             cTGAACAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAg                                                                      <  2:426997/141‑1 (MQ=255)

                                                              tGAAGTGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGa                                                                     <  2:101323/141‑1 (MQ=255)

                                                              tGAAGAGCCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGa                                                                     <  1:265197/141‑1 (MQ=255)

                                                              tGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGa                                                                     <  1:116087/141‑1 (MQ=255)

                                                              cgAAGAGACCTGCGGTTCTTTTTATATAGCGCGTAACTACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGa                                                                     <  1:16063/140‑1 (MQ=255)

                                                              agAAGAGACCTACTGTTCTTTTTATATAGAGCGTAAATACATTCAATACATTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGCCGGTGATCTTG‑CTCAATGAGGTTATCCAGa                                                                     <  2:406750/140‑1 (MQ=14)

                                                                  gagaCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGCTCTTG‑CTCAATGAGGTTATTCAGATAtt                                                                 <  2:207203/141‑1 (MQ=255)

                                                                  gagaCCTACGGTTCTTTTTATATTGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATAtt                                                                 <  1:49907/141‑1 (MQ=255)

                                                                   agaCCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATAttt                                                                >  2:246740/1‑141 (MQ=255)

                                                                       cTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGAt                                                            <  1:356837/141‑1 (MQ=255)

                                                                        tGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATg                                                           >  2:634511/1‑141 (MQ=255)

                                                                         tctgTTTTTTTTATATATAGCGTAAACACCCTTAATTCCTTTTAAATTATTCTTTGCTGTATTCATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGACCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGt                                                          <  2:363872/138‑1 (MQ=11)

                                                                         gCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGt                                                          >  1:240526/1‑141 (MQ=255)

                                                                         acGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTCAAGTATTCTTTGCCGTATTGATACTTTGATAGCTTTTCTGTCTTACATTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGt                                                          >  1:89011/2‑141 (MQ=11)

                                                                           ggTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTAc                                                        >  1:269062/1‑141 (MQ=255)

                                                                             ttCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CCCAATGAGGTTATTCAGATATTTCGATGTACaa                                                      <  1:594438/141‑1 (MQ=255)

                                                                              tCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAAt                                                     <  1:567393/141‑1 (MQ=255)

                                                                              tCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACATTGTCTTTCTTACTTTAATATGACGGTGATCTTA‑CTCAATGAGGTTATTCAGATATTTCGATGTACaag                                                     >  2:629279/1‑140 (MQ=255)

                                                                                tttttATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTTATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATgg                                                   <  2:94197/141‑2 (MQ=255)

                                                                                 ttttATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGTTACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATggc                                                  >  1:58133/1‑139 (MQ=255)

                                                                                  tttATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATggca                                                 <  1:368798/141‑4 (MQ=255)

                                                                                  tttATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCa                                                 >  1:193318/1‑141 (MQ=255)

                                                                                    tatttaTAGTGTAAATACTTTCAATACCTTTTTAAGTATTCTTTTCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGt                                               <  2:291867/141‑1 (MQ=14)

                                                                                    tatataGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATAACTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGt                                               <  1:73568/141‑1 (MQ=255)

                                                                                     atataAAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTc                                              >  2:554874/1‑141 (MQ=255)

                                                                                           agCGTAAATACATTCAATTCCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgttt                                        <  2:85571/141‑1 (MQ=255)

                                                                                            gCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgttta                                       >  1:14153/1‑141 (MQ=255)

                                                                                                aaaTACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgtttaagtt                                   <  1:101004/141‑1 (MQ=255)

                                                                                                aaaTACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGggttaagtt                                   >  1:340563/1‑141 (MQ=255)

                                                                                                 aaTACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATCCCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATCTTTCGCTTTACCATCGCCGTCAGgcttacggtt                                  >  2:124711/1‑136 (MQ=14)

                                                                                                  aTACATTCAATACCTTTTAACGTATTCTTTGCCGTATTGATACTTTGATATCTTATCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATTAGGTTTTTCAGCTAGTTCGATGTACAATGACAGACAGgtttaagttta                                 >  2:192070/1‑141 (MQ=14)

                                                                                                   tACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATAACGGAGCTCTTG‑CTCAATGAGATTAATGAGATATTTCGATGTACAATAACAGTCAGgtttaagtttaa                                >  1:452122/1‑141 (MQ=11)

                                                                                                   tACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACACTGGCAGTCAGgtgtaagtttaa                                >  2:388575/1‑141 (MQ=255)

                                                                                                   tACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagtttaa                                >  1:226612/1‑141 (MQ=255)

                                                                                                    acactCAATACCTTTTAAAGTACTCTTTGCTCAATTGATACTTTGATACCATGTCTTACTTACTTTAATATGACGGTGATCTTG‑CTCAATGATGTTATTCAGATATGTCGCTCTACAACGACAGTAAGgtttcaatttaca                               >  2:145038/5‑139 (MQ=255)

                                                                                                    aCATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTaaa                               >  2:575787/1‑141 (MQ=255)

                                                                                                     cactCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTaaaa                              >  2:601273/4‑141 (MQ=255)

                                                                                                      actCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGACGTAAGGTTTAAGTTTACAAg                             >  1:511290/3‑141 (MQ=14)

                                                                                                      actCAATACCTTTTAAAGTATTCTTTGCTGTATTGAGACTTTCCTACATAGTCTTTATTACATTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAg                             >  1:461766/3‑141 (MQ=14)

                                                                                                        tCAATACCTTTTAAAGTATTCTTTGCTGTAGTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGGTATTCCGATATTTCGATGCACCACCGCAGTCAGGTTTAAGTTTAAAAGcg                           <  2:491088/141‑2 (MQ=14)

                                                                                                         cAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCtt                          <  1:399339/141‑1 (MQ=255)

                                                                                                           aTACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTa                        >  1:656483/1‑141 (MQ=255)

                                                                                                             aCCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAAt                      <  1:547318/141‑1 (MQ=255)

                                                                                                             aCCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAAt                      >  1:617988/1‑141 (MQ=255)

                                                                                                              ccTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAAtt                     >  2:442877/1‑141 (MQ=255)

                                                                                                              ccTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAAtt                     >  1:603487/1‑141 (MQ=40)

                                                                                                              ccTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACACTGACAGTCAGGTTTAAGTTTAAAAGCTTTAAtt                     >  2:322767/1‑141 (MQ=255)

                                                                                                               cTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTa                    >  2:179078/1‑141 (MQ=255)

                                                                                                               cTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCCTCAGCTA‑AAGAAGCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTa                    >  1:485854/1‑141 (MQ=14)

                                                                                                                  ttAAAGTTTTCCTTGCTGTATTGATACGTTGATACCTTGTCTTTCTTACTTTAATATGCCGGTGATCTTG‑CTGAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACtt                 <  2:401615/141‑1 (MQ=14)

                                                                                                                  ttAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGCTCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTaattaatt                 >  1:289520/1‑141 (MQ=255)

                                                                                                                  ttAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACtt                 <  1:437817/141‑1 (MQ=255)

                                                                                                                   tAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACttt                >  1:444929/1‑141 (MQ=255)

                                                                                                                      aGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGc             >  1:202535/1‑141 (MQ=255)

                                                                                                                        tATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAGGTTTAAAAGCTTTAATTACTTTAACCa           >  1:55593/1‑141 (MQ=255)

                                                                                                                         aTTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCAt          >  2:513772/1‑141 (MQ=255)

                                                                                                                         aTTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCCGGTTTCAGCTTAAAAGCTTTAATTACTTTAGCCAt          >  2:468064/1‑141 (MQ=255)

                                                                                                                            cTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAACAGCTTTAATTACTTTATCCATTGc       <  2:647479/141‑1 (MQ=255)

                                                                                                                             tttGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCt      >  1:538265/1‑141 (MQ=255)

                                                                                                                                gCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTAcc   <  2:460323/141‑1 (MQ=255)

                                                                                                                                gCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTAcc   >  2:282688/1‑141 (MQ=255)

CAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCT  >  CP000730/35867‑36136

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