BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	418,825	A→T	10.9%	N53I (AAT→ATT)	rpsF →	30S ribosomal protein S6

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	418,825	0	A	T	10.9%	353.5 / 29.9	110	N53I (AAT→ATT)	rpsF	30S ribosomal protein S6
Reads supporting (aligned to +/- strand): ref base A (49/49); new base T (6/6); total (55/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

Only 100 of 111 total aligned reads displayed.

TGAACGTTTCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTAAAATCA  >  USA300TCH1516_ALE/418736‑418914

tGAACGTTTCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAAt                                                                                          <  1:3491174/91‑1 (MQ=255)

tGAACGTTTCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAAt                                                                                          <  1:3496848/91‑1 (MQ=255)

  aaCGTTTCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCATTGa                                                                                        <  1:2566606/91‑1 (MQ=255)

  aaCGTTTCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCATTGa                                                                                        <  1:3151301/91‑1 (MQ=255)

  aaCGTTTCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGa                                                                                        <  1:706863/91‑1 (MQ=255)

  aaCGTTTCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAAAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCATTGa                                                                                        <  1:3316180/91‑1 (MQ=255)

   aCGTTTCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCATTGAt                                                                                       <  1:4088065/91‑1 (MQ=255)

   aCGTTTCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAAtaga                                                                                       >  1:1651062/1‑88 (MQ=255)

   aCGTTTCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGAt                                                                                       >  1:3422913/1‑91 (MQ=255)

    cGTTTCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCATTGAtt                                                                                      <  1:488265/91‑1 (MQ=255)

     cgttCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGAttt                                                                                     >  1:521953/3‑91 (MQ=255)

      tttCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTTTGAAATCATTGATTTc                                                                                    <  1:2316469/91‑1 (MQ=255)

      tttCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTc                                                                                    >  1:2060161/1‑91 (MQ=255)

      tttCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTc                                                                                    <  1:2738836/91‑1 (MQ=255)

       ttCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCa                                                                                   >  1:1364429/1‑91 (MQ=255)

       ttCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCa                                                                                   >  1:1364702/1‑91 (MQ=255)

       ttCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCa                                                                                   >  1:269842/1‑91 (MQ=255)

       ttCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCa                                                                                   >  1:879840/1‑91 (MQ=255)

       ttCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATCTCa                                                                                   >  1:3343597/1‑91 (MQ=255)

        tCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCaa                                                                                  <  1:1767720/91‑1 (MQ=255)

        tCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTAGGGTAAACGTCGCGTAGCTTATGAAATCAATGACTTCaa                                                                                  >  1:3860912/1‑91 (MQ=255)

         cAACGGTATCTTAGCTACTGAAGGTGCAGAACTTTTAGAAGCAAAAGACTGGGGTAATCGTCGCCTAGCTTATGAAATCAATGATTTCaaa                                                                                 >  1:3681098/1‑91 (MQ=255)

           aCGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGa                                                                               <  1:3403559/91‑1 (MQ=255)

           aCGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGa                                                                               >  1:1076120/1‑91 (MQ=255)

           aCGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGa                                                                               <  1:1245653/91‑1 (MQ=255)

           aCGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGa                                                                               <  1:920033/91‑1 (MQ=255)

           aCGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGa                                                                               >  1:1076195/1‑91 (MQ=255)

            cGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGAt                                                                              >  1:1969704/1‑91 (MQ=255)

            cGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGAt                                                                              <  1:2729965/91‑1 (MQ=255)

                 tCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCtt                                                                         >  1:2865921/1‑91 (MQ=255)

                  cTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTc                                                                        >  1:687933/1‑91 (MQ=255)

                  cTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTc                                                                        >  1:687978/1‑91 (MQ=255)

                   ttAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCt                                                                       >  1:2805127/1‑91 (MQ=255)

                    tAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTa                                                                      <  1:2985906/91‑1 (MQ=255)

                       cTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTacaa                                                                   <  1:3962671/91‑1 (MQ=255)

                       cTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTacaa                                                                   >  1:3955732/1‑91 (MQ=255)

                        tACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTacaac                                                                  <  1:2250942/91‑1 (MQ=255)

                         aCTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTacaaca                                                                 <  1:405798/91‑1 (MQ=255)

                              aGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTa                                                            <  1:2350694/91‑1 (MQ=255)

                                gTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACg                                                          <  1:1148297/91‑1 (MQ=255)

                                 tGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACgt                                                         <  1:1822650/91‑1 (MQ=255)

                                 tGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACgt                                                         <  1:2006673/91‑1 (MQ=255)

                                  gCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACgtg                                                        <  1:2838434/91‑1 (MQ=255)

                                   cAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATTGCTTCTACAACATCGTACgtgt                                                       <  1:20016/91‑1 (MQ=255)

                                   cAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACgtgt                                                       <  1:1289611/91‑1 (MQ=255)

                                   cAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACgtgt                                                       <  1:1289644/91‑1 (MQ=255)

                                   cAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACgtgt                                                       <  1:1281681/91‑1 (MQ=255)

                                    aGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGtt                                                      <  1:2296367/91‑1 (MQ=255)

                                     gAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTa                                                     >  1:3933938/1‑91 (MQ=255)

                                      aaGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTaa                                                    <  1:2895829/91‑1 (MQ=255)

                                       aGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTaaa                                                   <  1:2552659/91‑1 (MQ=255)

                                       aGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTaaa                                                   <  1:2552643/91‑1 (MQ=255)

                                       aGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTaaa                                                   <  1:2552626/91‑1 (MQ=255)

                                         ttttAGAAGAAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTATTAAATc                                                 >  1:2080411/1‑91 (MQ=255)

                                           ttAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTg                                               <  1:3400782/91‑1 (MQ=255)

                                            tAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGa                                              <  1:1529575/91‑1 (MQ=255)

                                             aGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGAt                                             >  1:37913/1‑91 (MQ=255)

                                              gAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATa                                            <  1:803864/91‑1 (MQ=255)

                                              gAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATa                                            <  1:2908729/91‑1 (MQ=255)

                                                  cAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATaacaa                                        >  1:1486347/1‑91 (MQ=255)

                                                     aaGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAAcaacaa                                     >  1:1795551/1‑91 (MQ=255)

                                                      aGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACaaa                                    >  1:3511407/1‑91 (MQ=255)

                                                        aCTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGc                                  >  1:2031268/1‑91 (MQ=255)

                                                         cTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCt                                 <  1:1343267/91‑1 (MQ=255)

                                                          tGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTa                                <  1:618216/91‑1 (MQ=255)

                                                          tGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTa                                >  1:3694876/1‑91 (MQ=255)

                                                           ggggTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTAc                               >  1:2465437/1‑91 (MQ=255)

                                                            gggTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACt                              >  1:890956/1‑91 (MQ=255)

                                                            gggTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACt                              >  1:3750355/1‑91 (MQ=255)

                                                             ggTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTg                             <  1:1789773/91‑1 (MQ=255)

                                                             ggTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTg                             <  1:774062/91‑1 (MQ=255)

                                                             ggTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTg                             >  1:2440591/1‑91 (MQ=255)

                                                              gTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGa                            <  1:118204/91‑1 (MQ=255)

                                                              gTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGa                            <  1:116745/91‑1 (MQ=255)

                                                                   cGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAAt                       >  1:697496/1‑91 (MQ=255)

                                                                   cGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAAt                       <  1:2568243/91‑1 (MQ=255)

                                                                   cGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAAt                       >  1:1872654/1‑91 (MQ=255)

                                                                     tCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTc                     >  1:2473945/1‑91 (MQ=255)

                                                                     tCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTc                     >  1:855871/1‑91 (MQ=255)

                                                                     tCGCCTAGCTTATGAAATAAATGATTTCAAAGATGGCTTCTACTACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTc                     >  1:684210/1‑91 (MQ=255)

                                                                      cGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTcc                    >  1:963809/1‑91 (MQ=255)

                                                                      cGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTcc                    <  1:2513680/91‑1 (MQ=255)

                                                                      cGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTcc                    >  1:2618740/1‑91 (MQ=255)

                                                                      cGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTcc                    <  1:2628093/91‑1 (MQ=255)

                                                                         cTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAAc                 >  1:3628641/1‑91 (MQ=255)

                                                                          tAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACg                >  1:327565/1‑91 (MQ=255)

                                                                          tAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACg                >  1:330540/1‑91 (MQ=255)

                                                                            gCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTc              >  1:2252042/1‑91 (MQ=255)

                                                                             cTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCt             >  1:542227/1‑91 (MQ=255)

                                                                                   aaaTCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTaa       <  1:3073650/91‑1 (MQ=255)

                                                                                    aaTCATTGATTTCATAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTaaa      >  1:3071241/1‑91 (MQ=255)

                                                                                    aaTCATTGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTaaa      >  1:2188101/1‑91 (MQ=255)

                                                                                    aaTCATTGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTaaa      >  1:851151/1‑91 (MQ=255)

                                                                                    aaTCATTGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTaaa      >  1:1159506/1‑91 (MQ=255)

                                                                                    aaTCATTGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTaaa      >  1:393377/1‑91 (MQ=255)

                                                                                    aaTCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTaaa      <  1:887822/91‑1 (MQ=255)

                                                                                     aTCATTGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTaaaa     >  1:2183218/1‑91 (MQ=255)

                                                                                     aTCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTaaaa     <  1:2008271/91‑1 (MQ=255)

                                                                                      tCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTAAAAt    >  1:2697819/1‑91 (MQ=255)

                                                                                        aaTGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTAAAATCa  >  1:2992525/1‑91 (MQ=255)

TGAACGTTTCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTAAAATCA  >  USA300TCH1516_ALE/418736‑418914

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