BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	418,825	A→T	7.2%	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	7.2%	796.5 / 36.5	235	N53I (AAT→ATT)	rpsF	30S ribosomal protein S6
Reads supporting (aligned to +/- strand): ref base A (127/91); new base T (7/10); total (134/101)
Fisher's exact test for biased strand distribution p-value = 2.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 242 total aligned reads displayed.

TGAACGTTTCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTAAAATCAG  >  USA300TCH1516_ALE/418736‑418915

   aCGTTTCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGAt                                                                                         <  1:817372/91‑1 (MQ=255)

   aCGTTTCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGAt                                                                                         <  1:2599480/91‑1 (MQ=255)

   aCGTTTCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTCAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGAt                                                                                         <  2:215818/91‑1 (MQ=255)

    cGTTTCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGAtt                                                                                        >  1:3623783/1‑91 (MQ=255)

     cgttCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGAttt                                                                                       >  2:3196974/3‑91 (MQ=255)

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

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

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

        tCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCa‑a                                                                                   >  2:2648446/1‑91 (MQ=255)

         cAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCa‑aa                                                                                  <  2:572168/91‑1 (MQ=255)

         cAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCa‑aa                                                                                  >  1:3122497/1‑91 (MQ=255)

         cAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCa‑aa                                                                                  <  1:3556618/91‑1 (MQ=255)

          aaCGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAg                                                                                 >  1:1643973/1‑91 (MQ=255)

          aaCGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAg                                                                                 >  1:1622138/1‑91 (MQ=255)

          aaCGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAg                                                                                 <  1:3680536/91‑1 (MQ=255)

           aCGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGa                                                                                >  1:3055347/1‑91 (MQ=255)

           aCGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGa                                                                                >  2:3189019/1‑91 (MQ=255)

           aCGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGa                                                                                >  1:1063701/1‑91 (MQ=255)

           aCGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGa                                                                                >  1:1788259/1‑91 (MQ=255)

            cGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGAt                                                                               <  1:1363445/91‑1 (MQ=255)

            cGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGAt                                                                               <  1:1517178/91‑1 (MQ=255)

              gTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATgg                                                                             <  2:1072637/91‑1 (MQ=255)

               tATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGc                                                                            >  1:2814102/1‑91 (MQ=255)

                   ttAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCt                                                                        >  1:2460925/1‑91 (MQ=255)

                   ttAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCt                                                                        >  1:799621/1‑91 (MQ=255)

                     aGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTac                                                                      >  2:3563379/1‑91 (MQ=255)

                     aGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTac                                                                      <  2:1460326/91‑1 (MQ=255)

                      gCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTaca                                                                     >  1:2286722/1‑91 (MQ=255)

                        tACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTacaac                                                                   <  2:3132284/91‑1 (MQ=255)

                        tACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTacaac                                                                   <  2:2819680/91‑1 (MQ=255)

                          cTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACAt                                                                 >  1:662425/1‑91 (MQ=255)

                          cTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTCCTACAACAt                                                                 <  1:1948823/91‑1 (MQ=255)

                           tGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATc                                                                <  2:1781499/91‑1 (MQ=255)

                           tGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATc                                                                <  2:140690/91‑1 (MQ=255)

                           tGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATc                                                                <  2:1335315/91‑1 (MQ=255)

                           tGAAGGCGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATc                                                                <  2:3545479/91‑1 (MQ=255)

                                gTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACg                                                           >  2:2729844/1‑91 (MQ=255)

                                gTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACg                                                           <  2:1675222/91‑1 (MQ=255)

                                gTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACg                                                           <  2:2236109/91‑1 (MQ=255)

                                 tGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACgt                                                          <  2:2139330/91‑1 (MQ=255)

                                 tGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACgt                                                          >  1:2576530/1‑91 (MQ=255)

                                   cAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACgtgt                                                        >  2:1766313/1‑91 (MQ=255)

                                   cAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACgtgt                                                        <  2:243970/91‑1 (MQ=255)

                                   cAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACgtgt                                                        <  1:1349198/91‑1 (MQ=255)

                                   cAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACgtgt                                                        <  1:411895/91‑1 (MQ=255)

                                    aGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGtt                                                       <  1:3541189/91‑1 (MQ=255)

                                    aGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGtt                                                       >  1:2789456/1‑91 (MQ=255)

                                      aaGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTaa                                                     >  1:3469852/1‑91 (MQ=255)

                                       aGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTaaa                                                    >  1:212217/1‑91 (MQ=255)

                                            tAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGa                                               >  1:3428552/1‑91 (MQ=255)

                                            tAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGa                                               >  2:520176/1‑91 (MQ=255)

                                             aGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGAt                                              >  1:2723774/1‑91 (MQ=255)

                                             aGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGAt                                              >  1:1830844/1‑91 (MQ=255)

                                              gAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATa                                             <  2:2336931/91‑1 (MQ=255)

                                                aGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATaac                                           <  1:2028751/91‑1 (MQ=255)

                                                 gCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATaaca                                          >  2:3276917/1‑91 (MQ=255)

                                                   aaaaGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATaacaac                                        <  2:799621/91‑1 (MQ=255)

                                                    aaaGACTGGGGTGAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTa                                                                       <  1:3080598/59‑1 (MQ=255)

                                                    aaaGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAacaaca                                       >  2:3505983/1‑91 (MQ=255)

                                                    aaaGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAacaaca                                       >  2:2511845/1‑91 (MQ=255)

                                                      aGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACaaa                                     <  1:628450/91‑1 (MQ=255)

                                                      aGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACaaa                                     >  2:406513/1‑91 (MQ=255)

                                                       gACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAGTCTGATAACAACAAAg                                    >  2:2977823/1‑91 (MQ=255)

                                                       gACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAg                                    >  1:1738907/1‑91 (MQ=255)

                                                       gACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAg                                    >  2:3341488/1‑91 (MQ=255)

                                                        aCTGGGGTAAACGTCGCCTAGCTTATGAAATCCATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGc                                   >  2:3815190/1‑91 (MQ=255)

                                                        aCTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGc                                   >  2:932702/1‑91 (MQ=255)

                                                        aCTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGc                                   >  2:2035628/1‑91 (MQ=255)

                                                        aCTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGc                                   >  1:3231970/1‑91 (MQ=255)

                                                        aCTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGc                                   >  1:3286204/1‑91 (MQ=255)

                                                        aCTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGc                                   >  2:2093202/1‑91 (MQ=255)

                                                          tGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTa                                 <  1:2753455/91‑1 (MQ=255)

                                                          tGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTa                                 >  1:3151781/1‑91 (MQ=255)

                                                           ggggTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTAc                                <  2:2443191/91‑1 (MQ=255)

                                                            gggTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACt                               >  1:232189/1‑91 (MQ=255)

                                                            gggTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACt                               <  1:3766822/91‑1 (MQ=255)

                                                             ggTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTg                              >  2:1177948/1‑91 (MQ=255)

                                                              gTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGa                             <  2:2003615/91‑1 (MQ=255)

                                                                aaaCGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACg                           <  1:222106/91‑1 (MQ=255)

                                                                aaaCGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACg                           >  2:1961494/1‑91 (MQ=255)

                                                                     tCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTc                      >  1:3601310/1‑91 (MQ=255)

                                                                     tCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTc                      >  2:3744706/1‑91 (MQ=255)

                                                                     tCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTc                      >  1:869545/1‑91 (MQ=255)

                                                                      cGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTcc                     >  2:2869114/1‑91 (MQ=255)

                                                                       gCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCa                    >  1:2410296/1‑91 (MQ=255)

                                                                        ccTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCaa                   >  1:1869189/1‑91 (MQ=255)

                                                                          tAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACg                 >  2:2929165/1‑91 (MQ=255)

                                                                            gCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTc               <  1:1118114/91‑1 (MQ=255)

                                                                            gCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTc               >  1:3354669/1‑91 (MQ=255)

                                                                             cTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCGACTGACGAATTCCAACGTCt              <  2:2765011/91‑1 (MQ=255)

                                                                                 tGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCt          <  1:166705/91‑1 (MQ=255)

                                                                                   aaaTCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTaa        >  2:2359346/1‑91 (MQ=255)

                                                                                    aaTCATTGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTaaa       >  1:70571/1‑91 (MQ=255)

                                                                                     aTCATTGATTTCATAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTaaa       >  1:2927926/1‑91 (MQ=255)

                                                                                     aTCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTaaaa      >  1:562963/1‑91 (MQ=255)

                                                                                       cAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTAAAATc    >  1:2563761/1‑91 (MQ=255)

                                                                                        aaTGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTAAAATCa   >  1:1181002/1‑91 (MQ=255)

                                                                                         aTGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTAAAATCAg  >  2:1295376/1‑91 (MQ=255)

                                                                                         aTGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTAAAATCAg  >  2:2233921/1‑91 (MQ=255)

                                                                                         aTGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTAAAATCAg  >  2:3042396/1‑91 (MQ=255)

TGAACGTTTCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCA‑AAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTAAAATCAG  >  USA300TCH1516_ALE/418736‑418915

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