BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	111,560	T→C	12.8%	R210R (CGT→CGC)	USA300HOU_RS00520 →	putative lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	111,560	0	T	C	12.8%	616.1 / 61.4	225	R210R (CGT→CGC)	USA300HOU_RS00520	putative lipoprotein
Reads supporting (aligned to +/- strand): ref base T (105/90); new base C (12/17); total (117/108)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 227 total aligned reads displayed.

TTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATAT‑‑AAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATGAAAATATTTACTTTACTGA  >  USA300TCH1516_ALE/111420‑111700

    aGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGCaaaa                                                                                                                                            >  2:1406087/1‑141 (MQ=255)

    aGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTaaaa                                                                                                                                            >  2:1465530/1‑141 (MQ=255)

    aGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCAGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTaaaa                                                                                                                                            >  4:461140/1‑141 (MQ=255)

     gTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAg                                                                                                                                           <  1:295785/141‑1 (MQ=255)

      ttttttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGCAAAAGa                                                                                                                                          <  1:918950/141‑1 (MQ=255)

        ttttGCGAAATACACCAATTGTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGCAAAAGata                                                                                                                                        <  1:1334967/141‑1 (MQ=14)

         tttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGCAAAAGatat                                                                                                                                       >  2:230655/1‑141 (MQ=255)

          ttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGCAAAAGATATg                                                                                                                                      >  2:655024/1‑141 (MQ=255)

          ttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATg                                                                                                                                      <  1:20494/141‑1 (MQ=255)

           tGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGa                                                                                                                                     >  3:595766/1‑141 (MQ=255)

           tGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGa                                                                                                                                     >  3:92312/1‑141 (MQ=255)

             cgcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGCAAAAGATATGata                                                                                                                                   <  2:806678/141‑3 (MQ=255)

             cgcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGCAAAAGATATGata                                                                                                                                   <  2:723142/141‑3 (MQ=255)

              gcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGCAAAAgatatgatat                                                                                                                                  >  4:164281/1‑141 (MQ=255)

              gcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGCAAAAGATATGatat                                                                                                                                  <  4:544542/141‑4 (MQ=14)

              gcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACAt                                                                                                                                  <  2:1143904/141‑1 (MQ=255)

               cAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATa                                                                                                                                 <  1:565492/141‑1 (MQ=255)

               cAATATAGCAATTTTAAAGATTTAATGAATAATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCGCAAGATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATa                                                                                                                                 >  2:570338/1‑141 (MQ=255)

                  tataGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGCAAAAGATATGATATAccg                                                                                                                              <  2:387292/141‑2 (MQ=14)

                  tataGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACca                                                                                                                              <  2:1626902/141‑1 (MQ=255)

                  tataGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACca                                                                                                                              <  2:1188643/141‑1 (MQ=255)

                   ataGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGCAAAAGATATGATATAccga                                                                                                                             <  3:508147/141‑3 (MQ=14)

                    taGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTCACTAATGATGATTATAATGT‑AAAACAATTACGCAAAAGATATGATATACcgac                                                                                                                            <  1:845055/141‑1 (MQ=11)

                    taGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGCAAAAGATATGATATAccgac                                                                                                                            <  1:1317223/141‑4 (MQ=14)

                     aGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATGATATAATCCAGAGGTGCAGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGCAAAAGATATGATATAccgacg                                                                                                                           <  4:453738/141‑5 (MQ=11)

                     aGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGCAAAAGATATGATATAccgacg                                                                                                                           <  2:46789/141‑5 (MQ=11)

                      gCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGCAAAAGATATGATATACcgacga                                                                                                                          >  2:417367/1‑139 (MQ=11)

                      gCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACcaacaa                                                                                                                          >  3:699214/1‑141 (MQ=255)

                       cAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACaaa                                                                                                                         >  2:1067589/1‑141 (MQ=255)

                       aaaTTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGCAAAAGATATGATATAccgacgag                                                                                                                         <  1:626315/140‑7 (MQ=11)

                        aaTTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGCAAAAGATATGATATAccgacgagt                                                                                                                        <  2:662059/141‑8 (MQ=11)

                           tttAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATaaa                                                                                                                     <  2:962928/141‑1 (MQ=255)

                               aaGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAgcgc                                                                                                                 >  3:215956/1‑141 (MQ=255)

                               aaGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATAACAATTAACTAATGATGATTATAATGT‑ACAACAATTAGGTAACAGATATGACAGACCAACAAATAAAgcgc                                                                                                                 >  1:59733/1‑141 (MQ=18)

                                aGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGcc                                                                                                                <  3:440127/141‑1 (MQ=255)

                                 gATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCg                                                                                                               <  1:1373747/141‑1 (MQ=255)

                                  aTTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGATTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGa                                                                                                              <  4:414608/141‑1 (MQ=255)

                                        tGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtgt                                                                                                        <  3:263745/141‑1 (MQ=255)

                                             tataAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGttgttg                                                                                                   <  2:1217690/141‑1 (MQ=255)

                                             tataAAGATGGAGATATATCATATAACCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGttgttg                                                                                                   <  2:1313883/141‑1 (MQ=255)

                                              ataAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGa                                                                                                  <  4:358477/141‑1 (MQ=255)

                                               taAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGa                                                                                                  >  3:458472/1‑141 (MQ=255)

                                               taAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGaa                                                                                                 <  1:1174006/141‑1 (MQ=255)

                                                  aGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAgg                                                                                              <  4:320763/141‑1 (MQ=255)

                                                     tGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTAc                                                                                           >  1:1532639/1‑141 (MQ=255)

                                                     tGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTAc                                                                                           >  1:172168/1‑141 (MQ=255)

                                                         gaTATATCATATAATCCAGAGTTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAggg                                                                                       >  1:791296/1‑141 (MQ=255)

                                                         gaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGa                                                                                                                                          >  4:653197/1‑90 (MQ=255)

                                                         gaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAggg                                                                                       <  4:402328/141‑1 (MQ=255)

                                                         gaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAggg                                                                                       <  2:308767/141‑1 (MQ=255)

                                                                cATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTaa                                                                                <  1:1127121/141‑1 (MQ=255)

                                                                 atatAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTaaa                                                                               >  1:1049823/1‑141 (MQ=255)

                                                                     aaTCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGt                                                                           <  1:1067589/141‑1 (MQ=255)

                                                                     aaTCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATGAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGt                                                                           >  1:1062679/1‑141 (MQ=255)

                                                                      tttatgagGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGtt                                                                          <  2:1531214/136‑1 (MQ=255)

                                                                      aTCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGtt                                                                          >  1:1126133/1‑141 (MQ=255)

                                                                        ccAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtca                                                                        <  1:525849/141‑1 (MQ=255)

                                                                         cAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtcat                                                                       >  1:756479/1‑141 (MQ=255)

                                                                          agagGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGcc                                                                                                                >  3:148937/1‑99 (MQ=255)

                                                                          agagGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtcatc                                                                      >  1:1559669/1‑141 (MQ=255)

                                                                           gagGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtcatca                                                                     >  2:922078/1‑141 (MQ=255)

                                                                             ggTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTAAGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGt                                                                   <  4:245191/141‑1 (MQ=255)

                                                                                gCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTgg                                                                >  2:1344811/1‑141 (MQ=255)

                                                                                gCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTgg                                                                >  3:296712/1‑141 (MQ=255)

                                                                                gCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAAATTAAAAGGTTCATCAGTTgg                                                                >  4:440830/1‑141 (MQ=255)

                                                                                  cGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGat                                                              >  4:74488/1‑141 (MQ=255)

                                                                                   gAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGata                                                             <  4:616663/141‑1 (MQ=255)

                                                                                     gTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGAtat‑‑a                                                         >  3:469563/1‑141 (MQ=255)

                                                                                       tATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATaaa                                                                                                                     <  2:45388/81‑1 (MQ=255)

                                                                                       tATTCAGCTCAATATCAATTAACTAATGATGATTATAATGG‑AAAACAATTACGTAAAAGATATGACATACCAACAAATaaa                                                                                                                     >  1:45388/1‑81 (MQ=255)

                                                                                                 aaTATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATAT‑‑AAAAAAATTGAAt                                             >  2:1295098/1‑141 (MQ=255)

                                                                                                       aaTTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATTCCAACAAATAAAGCGCCGAATCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATAT‑‑AAAAAAATTGAATTTACtt                                       <  3:136057/141‑1 (MQ=255)

                                                                                                           aaCTAATGATTATGATAATTT‑ATAAGAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAATCTGTTGTTGTAAGGTACAGGTTATTTAATAGTTTTATGAGTTTGATAT‑‑AAAAAAATTGAATGTACtttttt                                   <  2:984836/141‑1 (MQ=17)

                                                                                                            aCTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATAT‑‑AAAAAAATTGAATTTACTTTTTTa                                  >  2:1329233/1‑141 (MQ=255)

                                                                                                            aCTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATAT‑‑AAAAAAATTGAATTTACTTTTTTa                                  >  2:87231/1‑141 (MQ=255)

                                                                                                             cTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATAT‑‑AAAAAAATTg                                                <  3:341355/126‑1 (MQ=255)

                                                                                                              aaaagatgatTATAATGT‑AAAACAATTACGAAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTT‑TTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATAT‑‑AAAAAAATTGAATTTACTTTTTTagag                               <  1:363772/137‑1 (MQ=255)

                                                                                                                atgatgATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATAT‑‑AAAAAAATTGAATTTACTTTTTTAgaga                              <  1:1513821/141‑1 (MQ=255)

                                                                                                                 tgatgaTTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATAT‑‑AAAAAAATTGAATTTACTTTTTTAGAGaa                             >  4:57246/1‑141 (MQ=255)

                                                                                                                 tgatgaTTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATAT‑‑AAAAAAATTGAATTTACTTTTTTAGAGaa                             >  1:397108/1‑141 (MQ=255)

                                                                                                                  gatgatTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATAT‑‑AAAAAAATTGAATTTACTTTTTTAGAGAAt                            <  1:263673/141‑1 (MQ=255)

                                                                                                                    tgatTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGAGAAGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATAT‑‑AAAAAAATTGAATTTAATTTTTTAGAGAATaa                          >  2:908437/1‑141 (MQ=255)

                                                                                                                     gatTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCGGTTGGAGAT‑‑AAAAAAATTGAATTTACTTTTTTAGAGAATaaa                         <  1:1423025/141‑1 (MQ=255)

                                                                                                                     gatTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATAT‑‑AAAAAAATTGAATTTACTTTTTTAGAGAATaaa                         >  1:1152889/1‑141 (MQ=255)

                                                                                                                       ttATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATAT‑‑AAAAAAATTGAATTTACTTTTGTAGAGAATaaaaa                       >  2:747388/1‑141 (MQ=255)

                                                                                                                        tataATGT‑AAAACAAGTACGTAAACGATATGACATACCAACAAATAAAGCGCCGAAGCTGCTGTTGAAAGGTACAGGGAATATAGAAGGTTCATCAGTTGGATAT‑‑AAAAAAATTGAATTTACTTTTTTAGAGAATAAAAAt                      >  2:1444027/1‑141 (MQ=18)

                                                                                                                         ataATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATAT‑‑AAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATg                     >  1:1038937/1‑141 (MQ=255)

                                                                                                                           aaTGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATg                     >  1:382319/1‑141 (MQ=255)

                                                                                                                             tGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATAT‑‑AAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATGaaaa                 <  2:887189/141‑1 (MQ=255)

                                                                                                                                    aaaaTTACGTAAAAGATATGACATAACAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATAT‑‑AAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATGAAAATAttta           <  2:1004885/139‑1 (MQ=255)

                                                                                                                                     cAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAggg                                                                                       >  2:1608582/1‑66 (MQ=255)

                                                                                                                                     cAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATAT‑‑AAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATGAAAATAtttac          <  3:358080/141‑1 (MQ=255)

                                                                                                                                      aaTTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATAT‑‑AAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATGAAAATAtttact         <  2:1054584/141‑1 (MQ=255)

                                                                                                                                       aTTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATAT‑‑AAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATGAAAATAtttactt        >  2:1606738/1‑141 (MQ=255)

                                                                                                                                        ttACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATAT‑‑AAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATGAAAATAtttacttt       >  3:210601/1‑141 (MQ=255)

                                                                                                                                        ttACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATAT‑‑AAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATGAAAATAtttacttt       >  4:514951/1‑141 (MQ=255)

                                                                                                                                         tACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATAT‑‑AAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATGAAAATAtttacttta      >  3:662134/1‑141 (MQ=255)

                                                                                                                                          aCGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATAT‑‑AAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATGAAAATAtttactttac     >  2:536276/1‑141 (MQ=255)

                                                                                                                                          aCGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATAT‑‑AAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATGAAAATAtttactttac     <  2:181131/141‑1 (MQ=255)

                                                                                                                                           cGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATAT‑‑AAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATGAAAATATttactttact    >  2:547572/1‑141 (MQ=255)

TTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGT‑AAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATAT‑‑AAAAAAATTGAATTTACTTTTTTAGAGAATAAAAATGAAAATATTTACTTTACTGA  >  USA300TCH1516_ALE/111420‑111700

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