BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,052,436	G→A	5.6%	I324I (ATC→ATT)	sspB ←	Staphopain B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,052,436	0	G	A	5.6%	365.1 / 6.8	119	I324I (ATC→ATT)	sspB	Staphopain B
Reads supporting (aligned to +/- strand): ref base G (39/73); new base A (3/4); total (42/77)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

Only 100 of 121 total aligned reads displayed.

GCATCTTGGATTGATAATTCTGTATCCCAAGGATTCCAGTAAATAAGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGTATTCAATCATTTGATTAGGGAATGTTGCACAATTAGGAAGGT  >  USA300TCH1516_ALE/1052296‑1052575

gCATCTTGGATTGATAATTCTGTATCCCAAGGATTCCAGTAAATAAGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAgg                                                                                                                                             <  2:727739/141‑1 (MQ=255)

 cATCTTGGATTGATAATTCTGTATCCCAAGGATTCCAGTAAATAAGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGa                                                                                                                                            <  3:394447/141‑1 (MQ=255)

  aTCTTGGATTGATAATTCTGTATCCCAAGGATTCCAGTAAATAAGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGAt                                                                                                                                           <  3:15229/141‑1 (MQ=255)

    cTTGGATTGATAATTCTGTATCCCAAGGATTCCAGTAAATAAGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCa                                                                                                                                         <  2:879020/141‑1 (MQ=255)

    cTTGGATTGATAATTCTGTATCCCAAGGATTCCAGTAAATAAGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCa                                                                                                                                         <  1:627604/141‑1 (MQ=255)

     ttGGATTGATAATTCTGTATCCCAAGGATTCCAGTAAATAAGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCAt                                                                                                                                        <  3:232514/141‑1 (MQ=255)

       ggATTGATAATTCTGTATCCCAAGGATTCCAGTAAATAAGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATaa                                                                                                                                      <  4:172129/141‑1 (MQ=255)

       ggATTGATAATTCTGTATCCCAAGGATTCCAGTAAATAAGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATaa                                                                                                                                      <  3:309281/141‑1 (MQ=255)

           tGATAATTCTGTAGCCCAAGGATTCCAGTAAATAAGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTcc                                                                                                                                  >  2:154156/1‑141 (MQ=255)

            gATAATTCTGTATCCCAAGGATTCCAGTAAATAAGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCt                                                                                                                                 <  2:121356/141‑1 (MQ=255)

              tAATTCTGTATCCCAAGGATTCCAGTAAATAAGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTAc                                                                                                                               <  3:102303/141‑1 (MQ=255)

              tAATTCTGTATCCCAAGGATTCCAGTAAATAAGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTAc                                                                                                                               <  1:163718/141‑1 (MQ=255)

              tAATACTGTATCCCAAGGATTCCAGTAAATAAGTGTTTCTTGGTCATTAATTTTAGCATTACCAACAGCTGCTAGCGCATGTCATAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTAc                                                                                                                               <  2:218798/141‑1 (MQ=255)

               aaTTCTGTATCCCAAGGATTCCAGTAAATAAGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACa                                                                                                                              <  4:343919/141‑1 (MQ=255)

                aTTCTGTATCCCAAGGATTCCAGTAAATAAGTTTTTTTTTTTTTTTTATTTTATCATTACCAACAACTGTTAGTGCATGTCCTAATTGTTGATCATTTGTGTTTTGCGATACAATTTGTGAAAGAATTGTAATTCCtttgt                                                                                                                             >  4:367578/1‑137 (MQ=255)

                aTTCTGTATCCCAAGGATTCCAGTAAATAAGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACAt                                                                                                                             <  2:908269/141‑1 (MQ=255)

                 ttCTGTATCCCAAGGATTCCAGTAAATAAGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACAtt                                                                                                                            <  2:452355/141‑1 (MQ=255)

                   cTGTATCCCAAGGATTCCAGTAAATAAGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTAt                                                                                                                          >  2:650552/1‑141 (MQ=255)

                     gTATCCCAAGGATTCCAGTAAATAAGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCt                                                                                                                        >  2:207098/1‑141 (MQ=255)

                      tATCCCAAGGATTCCAGTAAATAAGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCtt                                                                                                                       <  1:730089/141‑1 (MQ=255)

                      tATCCCAAGGATTCCAGTAAATAAGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCtt                                                                                                                       <  2:1023058/141‑1 (MQ=255)

                       aTCCCAAGGATTCCAGTAAATAAGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGGGCATGTCCTAAATGTGGATCATTAGGGTTTTGAGATACACTTTGTGCAAGAATCATAATTCCTACATTATCttt                                                                                                                      <  2:788744/141‑1 (MQ=255)

                        tCCCAAGGATTCCAGTAAATAAGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCtttt                                                                                                                     <  1:584351/141‑1 (MQ=255)

                         cccAAGGATTCCAGTAAATAAGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTg                                                                                                                    <  1:891469/141‑1 (MQ=255)

                             aGGATTCCAGTAAATAAGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAg                                                                                                                <  2:557270/141‑1 (MQ=255)

                                 ttCCAGTAAATAAGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGa                                                                                                            >  3:103249/1‑141 (MQ=255)

                                   ccAGTAAATAAGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATc                                                                                                          >  1:552258/1‑141 (MQ=255)

                                   ccAGTAAATAAGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATc                                                                                                          <  1:356108/141‑1 (MQ=255)

                                   ccAGTAAATAAGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATc                                                                                                          <  4:21102/141‑1 (MQ=255)

                                   ccAGTAAATAAGTTTTTCTTGGTCATTAATGTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGCTTTGAGATACACTTTGTGCAAGGATCATAATTGCTACATTATCTTTTGTAAGTGGATc                                                                                                          >  4:294103/1‑141 (MQ=255)

                                    cAGTAAATAAGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCa                                                                                                         <  2:307225/141‑1 (MQ=255)

                                    cAGTAAATAAGTTTTTCTAGGCCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTAGCTTTTGGAAGTTGATCa                                                                                                         >  2:237403/1‑141 (MQ=255)

                                      gTAAATAAGTTTTTCTTGGTCATTAATTTTAGCAGTACCAACAACTGCGAGCGCTTGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAAc                                                                                                       <  2:1059499/141‑1 (MQ=255)

                                           tAAGTTTTTCTTGGTCATTAATTTTAGCATTACCTACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGtt                                                                                                  <  2:495905/141‑1 (MQ=255)

                                             aGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCa                                                                                                >  2:1059262/1‑141 (MQ=255)

                                             aGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCa                                                                                                <  4:149791/141‑1 (MQ=255)

                                             aGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGGTCa                                                                                                >  2:137268/1‑141 (MQ=255)

                                              gTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAACGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCat                                                                                               >  1:934865/1‑141 (MQ=255)

                                                 tttCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGGGCATGTCCTAAATGTGGATCATTAGGGTTTTGAGATACACTTTGTGCAAGAATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATatg                                                                                            <  3:341970/141‑1 (MQ=255)

                                                     ttGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATgg                                                                                        >  2:67917/1‑141 (MQ=255)

                                                      tGGTCATTAATTTTAGCATTACCAACAACTGCTAGGGCATGTCCTAAATGTGGATCATTAGGGTTTTGAGATACACTTTGTGCAAGAATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGt                                                                                       >  1:669741/1‑141 (MQ=255)

                                                      tGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGt                                                                                       >  2:603826/1‑141 (MQ=255)

                                                       ggTCATTAATTTTAGCATTACCAACAACTGCTAGGGCATGTCCTAAATGTGGATCATTAGGGTTTTGAGATACACTTTGTGCAAGAATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTa                                                                                      <  1:534245/141‑1 (MQ=255)

                                                          cATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGc                                                                                   <  3:211722/141‑1 (MQ=255)

                                                             tAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCctt                                                                                <  3:11475/141‑1 (MQ=255)

                                                              aaTTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCcttc                                                                               <  4:159952/141‑1 (MQ=255)

                                                              aaTTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCcttc                                                                               <  4:179230/141‑1 (MQ=255)

                                                              aaTTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCcttc                                                                               >  2:825008/1‑141 (MQ=255)

                                                                ttttAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTATTCATATGATGGTACGCattctt                                                                             >  2:267781/1‑141 (MQ=255)

                                                                     gCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATaa                                                                        >  1:948371/1‑141 (MQ=255)

                                                                     gCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATaa                                                                        <  3:393935/141‑1 (MQ=255)

                                                                         tACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGaa                                                                    <  1:49851/141‑1 (MQ=255)

                                                                            caacaaCTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAatat                                                                 <  2:975540/141‑1 (MQ=255)

                                                                             aacaacTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATAtc                                                                >  3:177555/1‑141 (MQ=255)

                                                                             aacaacTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATAtc                                                                >  2:347905/1‑141 (MQ=255)

                                                                               caacTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATAtctc                                                              <  4:303272/141‑1 (MQ=255)

                                                                               caacTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATAtctc                                                              >  1:436523/1‑141 (MQ=255)

                                                                               caacTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATAtctc                                                              <  2:704942/141‑1 (MQ=255)

                                                                                 acTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTg                                                            <  1:895962/141‑1 (MQ=255)

                                                                                   tGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGcc                                                          >  2:877265/1‑141 (MQ=255)

                                                                                   tGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTATTGATAATGAAGCTCTCTGcc                                                          <  4:204810/141‑1 (MQ=255)

                                                                                    gCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCt                                                         <  3:378272/141‑1 (MQ=255)

                                                                                     cTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCtt                                                        >  2:953084/1‑141 (MQ=255)

                                                                                         cgcATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGa                                                    <  1:218619/141‑1 (MQ=255)

                                                                                           cATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGAtt                                                  <  3:205253/141‑1 (MQ=255)

                                                                                               tCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGAACATAATTCCTACATTCTCTTTTGTAAGTTGATCAACTTGTTCATATTATGGTACGCCTTCTTGATAATGAGTATCTCTGCCTTGTGATTTAcc                                              >  1:291867/1‑141 (MQ=255)

                                                                                                 cTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGt                                            <  1:913635/141‑1 (MQ=255)

                                                                                                    aaTGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGTAtt                                         <  3:394696/141‑1 (MQ=255)

                                                                                                     aTGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGTATTc                                        <  2:166962/141‑1 (MQ=255)

                                                                                                      tgtgGGTCATCAGGTTGTTGAGACACACTTTGTGCAAGGAGCATAATGCCTACATTATCTTTTGTAAGTTGATTGACTTGTTCATCTGAGGGTACGCCGTCGTGATAATGAATGTCGCTGGCTTGGGATTTACCGTATTCa                                       <  4:255810/141‑1 (MQ=255)

                                                                                                        tgGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGTATTCAAt                                     >  1:406383/1‑141 (MQ=255)

                                                                                                         gggTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGAGCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGTATTCAATc                                    <  2:1048166/141‑1 (MQ=255)

                                                                                                          ggTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGTATTCAATCa                                   <  4:378756/141‑1 (MQ=255)

                                                                                                          ggTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGTATTCAATCa                                   >  2:588528/1‑141 (MQ=255)

                                                                                                          ggTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGTATTCAATCa                                   >  2:74086/1‑141 (MQ=255)

                                                                                                           gTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGTATTCAATCAt                                  <  2:1053927/141‑1 (MQ=255)

                                                                                                           gTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATACCTCTGCCTTGTGATTTACCGTATTCAATCAt                                  <  1:630067/141‑1 (MQ=255)

                                                                                                              aTTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGTATTCAATCATTTg                               >  2:304453/1‑141 (MQ=255)

                                                                                                                 aGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGGTTTACCGTATTCAATCATTTGAtt                            <  3:58055/141‑1 (MQ=255)

                                                                                                                 aGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGTATTCAATCATTTGAtt                            <  1:1013341/141‑1 (MQ=255)

                                                                                                                   ggTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGTATTTAATCATTTGATTAg                          >  1:1060426/1‑141 (MQ=255)

                                                                                                                   ggTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGTATTCAATCATTTGATTAg                          >  1:776454/1‑141 (MQ=255)

                                                                                                                   ggTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGTATTCAATCATTTGATTAg                          <  1:853334/141‑1 (MQ=255)

                                                                                                                     ttttGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGTATTCAATCATTTGATTAggg                        >  2:827108/1‑141 (MQ=255)

                                                                                                                     ttttGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGTATTCAATCATTTGATTAggg                        <  3:372948/141‑1 (MQ=255)

                                                                                                                        tGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGTATTCAATCATTTTCTTACGGAAt                     >  2:129747/1‑141 (MQ=255)

                                                                                                                           gaTACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGTATTCAATCATTTGATTAGGGAATGtt                  >  1:349786/1‑141 (MQ=255)

                                                                                                                           gaTACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGTATTCAATCATTTGATTAGGGAATGtt                  >  1:334448/1‑141 (MQ=255)

                                                                                                                              acacTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGTATTCAATCATTTGATTAGGGAATGTTGca               <  4:230686/141‑1 (MQ=255)

                                                                                                                               cacTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGTATTCAATCATTTGATTAGGGAATGTTGcac              <  4:134642/141‑1 (MQ=255)

                                                                                                                                   ttGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGTATTCAATCATTTGATTAGGGAATGTTGCACAAtt          <  3:279948/141‑1 (MQ=255)

                                                                                                                                   ttGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGTATTCAATCATTTGATTAGGGAATGTTGCACAAtt          >  1:742040/1‑141 (MQ=255)

                                                                                                                                    tgtgCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGTATTCAATCATTTGATTAGGGAATGTTGCACAATTa         >  1:230083/1‑141 (MQ=255)

                                                                                                                                    tgtgCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGTATTCAATCATTTGATTAGGGAATGTTGCACAATTa         <  2:606831/141‑1 (MQ=255)

                                                                                                                                     gtgCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGTATTCAATCATTTGATTAGGGAATGTTGCACAATTAg        >  2:684374/1‑141 (MQ=255)

                                                                                                                                       gCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGTATTCAATCATTTGATTAGGGAATGTTGCACAATTAGGa      <  2:635351/141‑1 (MQ=255)

                                                                                                                                       gCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGTATTCAATCATTTGATTAGGGAATGTTGCACAATTAGGa      <  1:436664/141‑1 (MQ=255)

                                                                                                                                         aaGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGTATTCAATCATTTGATTAGGGAATGTTGCACCATTAGGAAg    <  2:427688/141‑1 (MQ=255)

                                                                                                                                         aaGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGTATTCAATCATTTGATTAGGGAATGTTGCACAATTAGGAAg    >  2:680871/1‑141 (MQ=255)

                                                                                                                                           ggATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGTATTCAATCATTTGATTAGGGAATGTTGCACAATTAGGAAGGt  <  1:650552/141‑1 (MQ=255)

GCATCTTGGATTGATAATTCTGTATCCCAAGGATTCCAGTAAATAAGTTTTTCTTGGTCATTAATTTTAGCATTACCAACAACTGCTAGCGCATGTCCTAAATGTGGGTCATTAGGGTTTTGAGATACACTTTGTGCAAGGATCATAATTCCTACATTATCTTTTGTAAGTTGATCAACTTGTTCATATGATGGTACGCCTTCTTGATAATGAATATCTCTGCCTTGTGATTTACCGTATTCAATCATTTGATTAGGGAATGTTGCACAATTAGGAAGGT  >  USA300TCH1516_ALE/1052296‑1052575

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