BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	977,289	T→A	9.0%	I190N (ATT→AAT)	clpB_1 →	Chaperone protein ClpB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	977,289	0	T	A	9.0%	436.2 / 20.6	129	I190N (ATT→AAT)	clpB_1	Chaperone protein ClpB
Reads supporting (aligned to +/- strand): ref base T (54/63); new base A (2/10); total (56/73)
Fisher's exact test for biased strand distribution p-value = 6.69e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 132 total aligned reads displayed.

AGCTAAATATGGCCGCGACTTAGTAGAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAGGTGTTGGTAAAACTGCAATTGTTGAAGGATTAGCGC  >  USA300TCH1516_ALE/977200‑977378

aGCTAAATATGGCCGCGACTTAGTAGAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGAtt                                                                                          <  1:1905836/91‑1 (MQ=255)

aGCTAAATATGGCCGCGACTTAGTAGAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGAtt                                                                                          >  1:2517859/1‑91 (MQ=255)

 gCTAAATATGGCCGCGACTTATGAGAAGCAGTTGTCGAAGGTAAAATGGATCATGTTATAGTAAGCGATGAAGAAATTCGAAATACGAATc                                                                                         <  1:3743058/91‑1 (MQ=255)

 gCTAAATATGGCCGCGACTTAGTAGAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTc                                                                                         >  1:1272553/1‑91 (MQ=255)

 gCTAAATATGGCCGCGACTTAGTAGAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGAATc                                                                                         <  1:1045976/91‑1 (MQ=255)

  cTAAATATGGCCGCGACTTAGTAGAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCg                                                                                        <  1:1488257/91‑1 (MQ=255)

  cTAAATATGGCCGCGACTTAGTAGAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCg                                                                                        >  1:17276/1‑91 (MQ=255)

  cTAAATATGGCCGCGACTTAGTAGAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGAATCg                                                                                        <  1:3714302/91‑1 (MQ=255)

  cTAAATATGGCCGCGACTTAGTAGAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGAATCg                                                                                        <  1:3714411/91‑1 (MQ=255)

  cTAAATATGGCCGCGACTTAGTAGAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGAATCg                                                                                        <  1:2941278/91‑1 (MQ=255)

  cTAAATATGGCCGCGACTTAGTAGAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGAATCg                                                                                        <  1:2837302/91‑1 (MQ=255)

   tAAATATGGCCGCGACTTAGTAGAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGAATCGt                                                                                       <  1:2806679/91‑1 (MQ=255)

          ggCCGCGACTTAGTAGAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTaa                                                                                >  1:3108710/1‑91 (MQ=255)

          ggCCGCGACTTAGTAGAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTaa                                                                                >  1:260233/1‑91 (MQ=255)

          ggCCGCGACTTAGTAGAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTaa                                                                                >  1:261507/1‑91 (MQ=255)

           gCCGCGACTTAGTAGAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTGAAg                                                                               >  1:3341621/1‑91 (MQ=255)

            ccGCGACTTAGTAGAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGt                                                                              <  1:2427243/91‑1 (MQ=255)

            ccGCGACTTAGTAGAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGt                                                                              <  1:2950604/91‑1 (MQ=255)

            ccGCGACTGAGTAGAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGt                                                                              <  1:2426895/91‑1 (MQ=255)

             cgcgACTTAGTAGAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTc                                                                             <  1:1629578/91‑1 (MQ=255)

             cgcgACTTAGTAGAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTc                                                                             >  1:1434636/1‑91 (MQ=255)

             cgcgACTTAGTAGAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTc                                                                             >  1:359685/1‑91 (MQ=255)

               cgACTTAGTAGAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGt                                                                           >  1:1149724/1‑91 (MQ=255)

                 aCTTAGTAGAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTaa                                                                         <  1:3674875/91‑1 (MQ=255)

                 aCTGAGTAGAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTaa                                                                         <  1:954396/91‑1 (MQ=255)

                  cTTAGTAGAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTaaa                                                                        >  1:2898589/1‑91 (MQ=255)

                    tagtagAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAAc                                                                      <  1:748504/91‑1 (MQ=255)

                     agtagAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACt                                                                     >  1:2157795/1‑91 (MQ=255)

                     agtagAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACt                                                                     >  1:2157120/1‑91 (MQ=255)

                     agtagAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACt                                                                     >  1:2156635/1‑91 (MQ=255)

                      gtagAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTa                                                                    >  1:2688219/1‑91 (MQ=255)

                      gtagAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTa                                                                    >  1:1073111/1‑91 (MQ=255)

                      gtagAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTa                                                                    >  1:3986040/1‑91 (MQ=255)

                      gtagAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGGAAAACTa                                                                    >  1:2080971/1‑91 (MQ=255)

                        agaagaAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTaaa                                                                  <  1:3610198/91‑1 (MQ=255)

                        agaagaAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTaaa                                                                  <  1:3609774/91‑1 (MQ=255)

                            gaagTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAaaca                                                              >  1:1392170/1‑91 (MQ=255)

                             aagTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAaacaa                                                             <  1:2334538/91‑1 (MQ=255)

                               gTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAAcc                                                           >  1:2668227/1‑91 (MQ=255)

                                  aGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCtg                                                        <  1:2132936/91‑1 (MQ=255)

                                  aGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCtg                                                        >  1:1197389/1‑91 (MQ=255)

                                  aGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCtg                                                        <  1:2132767/91‑1 (MQ=255)

                                  aGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCtg                                                        >  1:1523526/1‑91 (MQ=255)

                                  aGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCtg                                                        >  1:69184/1‑91 (MQ=255)

                                  aGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCtg                                                        <  1:2132964/91‑1 (MQ=255)

                                  aGACAAGGTAAAATGGATCCGGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCtg                                                        <  1:3289095/91‑1 (MQ=255)

                                     cAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGc                                                     <  1:288181/91‑1 (MQ=255)

                                        ggTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCa                                                  <  1:1241533/91‑1 (MQ=255)

                                        ggTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCa                                                  <  1:1241534/91‑1 (MQ=255)

                                        ggTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCa                                                  >  1:2496933/1‑91 (MQ=255)

                                        ggTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCa                                                  <  1:1373342/91‑1 (MQ=255)

                                          tAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAATCCTGTGCTCAtt                                                <  1:3201313/91‑1 (MQ=255)

                                          tAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAAGCCTGTGCTCAtt                                                <  1:3201247/91‑1 (MQ=255)

                                          tAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCAtt                                                <  1:2576617/91‑1 (MQ=255)

                                          tAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCAtt                                                <  1:594435/91‑1 (MQ=255)

                                          tAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCAtt                                                <  1:594471/91‑1 (MQ=255)

                                                 gATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATAGGTGAAc                                         >  1:2070445/1‑91 (MQ=255)

                                                  aTCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAAcc                                        >  1:736690/1‑91 (MQ=255)

                                                  aTCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAAcc                                        <  1:347078/91‑1 (MQ=255)

                                                   tCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCa                                       >  1:2026506/1‑91 (MQ=255)

                                                   tCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCa                                       >  1:4073721/1‑91 (MQ=255)

                                                   tCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCa                                       <  1:2859927/91‑1 (MQ=255)

                                                   tCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCa                                       >  1:4075883/1‑91 (MQ=255)

                                                    ccTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAg                                      >  1:1806222/1‑91 (MQ=255)

                                                    ccTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAg                                      <  1:794858/91‑1 (MQ=255)

                                                     cTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATATTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCCgg                                     >  1:779287/1‑91 (MQ=255)

                                                      tGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAGgt                                    >  1:2592211/1‑91 (MQ=255)

                                                         tataGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAGGTGtt                                 >  1:2897625/1‑91 (MQ=255)

                                                           taGGAGGAGATGAAGAAATTCGAAATACGATTCGTATTTGAAGTCGTAAAACTGAAAACAACCCTGTGCTCATTGGTGAACCAGGTGTTgg                               <  1:3861360/91‑1 (MQ=255)

                                                           taGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAGGTGTTgg                               >  1:430336/1‑91 (MQ=255)

                                                            aGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAGGTGTTGGt                              >  1:1047730/1‑91 (MQ=255)

                                                             ggAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAGGTGTTGGTa                             >  1:2637627/1‑91 (MQ=255)

                                                               aaGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAGGTGTTGGTaaa                           >  1:2296270/1‑91 (MQ=255)

                                                                   gaTGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAGGTGTTGGTAAAACTg                       <  1:488160/91‑1 (MQ=255)

                                                                   gaTGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAAATAAAAACAACCCTGTGCTCATTGGTGAACCAGGTGTTGGTAAAACTg                       <  1:4150608/91‑1 (MQ=255)

                                                                    aTGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAGGTGTTGGTAAAACTGc                      <  1:1487048/91‑1 (MQ=255)

                                                                    aTGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAGGTGTTGGTAAAACTGc                      <  1:79069/91‑1 (MQ=255)

                                                                     tGAAGAAATTCGAAATACGATTCGTATTTTATGTCGTAAAACTAGAAACAACCATGTGCTCATTGGTGAACCAGGTGTTGGTAAAACTGCa                     <  1:3312648/91‑1 (MQ=255)

                                                                     tGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAGGTGTTGGTAAAACTGCa                     >  1:3568304/1‑91 (MQ=255)

                                                                     tGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAGGTGTTGGTAAAACTGCa                     >  1:2954649/1‑91 (MQ=255)

                                                                     tGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAGGTGTTGGTAAAACTGCa                     <  1:1645517/91‑1 (MQ=255)

                                                                      gaagaaATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAGGTGTTGGTAAAACTGCaa                    <  1:2756831/91‑1 (MQ=255)

                                                                      gaagaaATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAGGTGTTGGTAAAACTGCaa                    <  1:937625/91‑1 (MQ=255)

                                                                      gaagaaATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAGGTGTTGGTAAAACTGCaa                    >  1:3662768/1‑91 (MQ=255)

                                                                        agaaATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAGGTGTTGGTAAAACTGCAAtt                  >  1:2060951/1‑91 (MQ=255)

                                                                        agaaATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAGGTGTTGGTAAAACTGCAAtt                  >  1:2060988/1‑91 (MQ=255)

                                                                         gaaATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAGGTGTTGGTAAAACTGCAAttg                 <  1:1961059/91‑1 (MQ=255)

                                                                         gaaATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAGGTGTTGGTAAAACTGCAAttg                 <  1:1961702/91‑1 (MQ=255)

                                                                            aTTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAGGTGTTGGTAAAACTGCAAttgttg              >  1:2251327/1‑91 (MQ=255)

                                                                               cGAAATACGATTCGTATTTTAAGTCGTACAACTAAAAACAACCCTGTGCTCATTGGTGAACTAGGTGTTGGTAAAACTGCGATTGTTGAAg           >  1:3950542/1‑91 (MQ=255)

                                                                                  aaTACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAGGTGTTGGTAAAACTGCAATTGTTGAAGGAt        >  1:872296/1‑91 (MQ=255)

                                                                                  aaTACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAGGTGTTGGTAAAACTGCAATTGTTGAAGGAt        <  1:2515686/91‑1 (MQ=255)

                                                                                  aaTACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAGGTGTTGGTAAAACTGCAATTGTTGAAGGAt        <  1:2515427/91‑1 (MQ=255)

                                                                                   aTACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAGGTGTTGGTAAAACTGCAATTGTTGAAGGAtt       <  1:2959865/91‑1 (MQ=255)

                                                                                   aTACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAGGTGTTGGTAAAACTGCAATTGTTGAAGGAtt       <  1:1158678/91‑1 (MQ=255)

                                                                                    tACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAGGTGTTGGTAAAACTGCAATTGTTGAAGGATTa      <  1:1039159/91‑1 (MQ=255)

                                                                                    tACGAATCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAGGTGTTGGTAAAACTGCAATTGTTGAAGGATTa      >  1:1787443/1‑91 (MQ=255)

                                                                                      cGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAGGTGTTGGTAAAACTGCAATTGTTGAAGGATTAgc    <  1:556404/91‑1 (MQ=255)

                                                                                       gATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAGGTGTTGGTAAAACTGCAATTGTTGAAGGATTAgcg   <  1:2606286/91‑1 (MQ=255)

                                                                                        aTTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAGGTGTTGGTAAAACTGCAATTGTTGAAGGATTAgcgc  <  1:3550141/91‑1 (MQ=255)

AGCTAAATATGGCCGCGACTTAGTAGAAGAAGTTAGACAAGGTAAAATGGATCCTGTTATAGGAAGAGATGAAGAAATTCGAAATACGATTCGTATTTTAAGTCGTAAAACTAAAAACAACCCTGTGCTCATTGGTGAACCAGGTGTTGGTAAAACTGCAATTGTTGAAGGATTAGCGC  >  USA300TCH1516_ALE/977200‑977378

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