BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	269,671	T→G	9.6%	T82P (ACG→CCG)	fadA ←	3‑ketoacyl‑CoA thiolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	269,671	0	T	G	9.6%	553.0 / 48.4	182	T82P (ACG→CCG)	fadA	3‑ketoacyl‑CoA thiolase
Reads supporting (aligned to +/- strand): ref base T (69/95); new base G (16/2); total (85/97)
Fisher's exact test for biased strand distribution p-value = 2.34e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.38e-01

Only 100 of 195 total aligned reads displayed.

CTCCTTGACCAGCCATAATTTGATTGGCTGCAATCGCGATGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCGTACCAACAATGACATCTTC  >  USA300TCH1516_ALE/269581‑269761

ctcCTTGACCAGCCATAATTTGATTGGCTGCAATCGCGATGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCgt                                                                                            <  2:718175/91‑1 (MQ=255)

ctcCTTGACCAGCCATAATTTGATTGGCTGCAATCGCGATGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCgt                                                                                            >  1:880432/1‑91 (MQ=255)

 tcCTTGACCAGCCATAATTTGATTGGCTGCAATCGCGATGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCgtg                                                                                           >  1:3291114/1‑91 (MQ=255)

  ccTTGACCAGCCATAATTTGATTGGCTGCAATCGCGATGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCgtgt                                                                                          <  2:1112954/91‑1 (MQ=255)

   cTTGACCAGCCATAATTTGATTGGCTGCAATCGCGATGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTc                                                                                         <  2:1542709/91‑1 (MQ=255)

    ttGACCAGCCATAATTTGATTGGCTGCAATCGCGATGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACACGGTACCGGGTCa                                                                                        <  2:3696095/91‑1 (MQ=255)

    tggACCGGCCATAATTTGATTGGCTGCAATCGCGATGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCa                                                                                        <  2:2247875/89‑1 (MQ=255)

        ccAGCCATAATTTGATTGGCTGCAATCGCGATGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATa                                                                                    <  2:257272/91‑1 (MQ=255)

         cAGCCATAATTTGATTGGCTGCAATCGCGATGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATaa                                                                                   <  2:59733/91‑1 (MQ=255)

         cAGCCATAATTTGATTGGCTGCAATCGCGAAGGATCGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATaa                                                                                   >  2:235974/1‑91 (MQ=255)

          aGCCATAATTTGATTGGCTGCAATCGCGATGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAAt                                                                                  >  2:937836/1‑91 (MQ=255)

           gCCATAATTTGATTGGCTGCAATCGCGATGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATc                                                                                 <  2:1811332/91‑1 (MQ=255)

            ccATAATTTGATTTGCTGCAATCGCGATGGTTTGTAATACTGATGAGCAGTAGCGATTCACTGTTTGACCAGGTACCGTGTCAGATAATcc                                                                                >  2:529653/1‑91 (MQ=255)

             cATAATTTGATTGGCTGCAATCGCGATGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATccc                                                                               <  2:407587/91‑1 (MQ=255)

               tAATTTGATTGGCTGCAATCGCGATGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCcgc                                                                             <  1:3748368/91‑1 (MQ=255)

                 aTTTGATTGGCTGCAATCGCGATGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCcgcac                                                                           >  2:2952604/1‑91 (MQ=255)

                   ttGATTGGCTGCAATCGCGATGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCcgcacgc                                                                         <  1:139952/91‑1 (MQ=255)

                   ttGAGTGGCTGCAATCGCGATGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCcgcacgc                                                                         <  2:3370308/91‑1 (MQ=255)

                       ttGGCTGCAATCGCGATGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGGCCCGGTACCGTGTCAGATAATCCCGCACgcaatg                                                                     >  2:483837/1‑91 (MQ=255)

                       ttGGCTGCAATCGCGATGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACgcaatg                                                                     <  1:3196600/91‑1 (MQ=255)

                        tGGCTGCAATCGCGATGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACgcaatgc                                                                    >  1:1342346/1‑91 (MQ=255)

                          gCTGCAATCGCGATGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACgcaatgcaa                                                                  <  1:1450790/91‑1 (MQ=255)

                           cTGCAATCGCGATGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACgcaatgcaat                                                                 <  1:3821127/91‑1 (MQ=255)

                              cAATCGCGATGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAAtcgt                                                              >  1:2929864/1‑91 (MQ=255)

                                aTCGCGATGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAAtcgttc                                                            >  2:2756782/1‑91 (MQ=255)

                                 tCGCGATGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAAtcgttcg                                                           <  1:2859922/91‑1 (MQ=255)

                                 tCGCGATGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAAtcgttcg                                                           >  2:2183154/1‑91 (MQ=255)

                                    cgATGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGTAATCGTTCGTGc                                                        <  1:2969270/91‑1 (MQ=255)

                                    cgATGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGc                                                        <  1:3590593/91‑1 (MQ=255)

                                      aTGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCaa                                                      <  1:2693467/91‑1 (MQ=255)

                                      aTGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCaa                                                      >  2:1794559/1‑91 (MQ=255)

                                       tGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAAt                                                     <  1:3793356/91‑1 (MQ=255)

                                       tGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAAt                                                     >  2:713770/1‑91 (MQ=255)

                                        ggTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATg                                                    <  2:3486381/91‑1 (MQ=255)

                                          tttGTAATCCTGATGAGCAGTAGCGATACACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGtt                                                  >  2:1273568/1‑91 (MQ=255)

                                           ttGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGttt                                                 <  1:2488495/91‑1 (MQ=255)

                                             gTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTg                                               >  1:3665142/1‑91 (MQ=255)

                                               aaTCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGc                                             <  1:194554/91‑1 (MQ=255)

                                                  ccTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCtt                                          <  2:3254008/91‑1 (MQ=255)

                                                   cTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTg                                         <  2:2903565/91‑1 (MQ=255)

                                                   cTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTg                                         >  1:2723596/1‑91 (MQ=255)

                                                    tgatgaGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGt                                        <  1:3398325/91‑1 (MQ=255)

                                                    tgatgaGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGt                                        <  2:3134584/91‑1 (MQ=255)

                                                     gatgaGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTa                                       >  1:918657/1‑91 (MQ=255)

                                                        gaGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATc                                    >  2:3883783/1‑91 (MQ=255)

                                                           cAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCtt                                 <  1:1923744/91‑1 (MQ=255)

                                                           cAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCtt                                 <  2:2298066/91‑1 (MQ=255)

                                                            aGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTc                                <  1:1215394/91‑1 (MQ=255)

                                                             gTAGCGATTCACTGTTTGGCCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCt                               <  2:3818795/91‑1 (MQ=255)

                                                             gTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCt                               >  2:1973952/1‑91 (MQ=255)

                                                               aGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTgg                             <  2:2317858/91‑1 (MQ=255)

                                                               aGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTgg                             <  1:1004967/91‑1 (MQ=255)

                                                               aGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTgg                             <  1:3119723/91‑1 (MQ=255)

                                                                gCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGa                            >  2:3665382/1‑91 (MQ=255)

                                                                gCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGa                            <  1:1380649/91‑1 (MQ=255)

                                                                 cGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGaa                           >  1:1568731/1‑91 (MQ=255)

                                                                 cGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGaa                           <  2:3222069/91‑1 (MQ=255)

                                                                 cGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGaa                           <  1:235974/91‑1 (MQ=255)

                                                                 cGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGaa                           <  1:1198502/91‑1 (MQ=255)

                                                                 cGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGaa                           <  1:1064139/91‑1 (MQ=255)

                                                                  gATTCACTTTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGAGCAATGTTTTGGCCTTGTAATCCTTCTGGaaa                          <  2:3801730/91‑1 (MQ=255)

                                                                  gATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGaaa                          >  2:1396877/1‑91 (MQ=255)

                                                                  gATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGaaa                          >  2:1742767/1‑91 (MQ=255)

                                                                  gATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGaaa                          <  2:849271/91‑1 (MQ=255)

                                                                  gATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGaaa                          <  2:849288/91‑1 (MQ=255)

                                                                   aTTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGaaaa                         >  1:3431028/1‑91 (MQ=255)

                                                                    ttCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAg                        <  2:1600991/91‑1 (MQ=255)

                                                                       aCTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCg                     >  1:1897566/1‑91 (MQ=255)

                                                                       aCTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCg                     >  2:3121631/1‑91 (MQ=255)

                                                                        cTGTTTGACCCGGTACCGTGTCAGATGGTCCCGCACGCAATGCAATCGTTCGTGCAATTTTTTGGCCTTGTAATCCTTCTGGAAATGCCGt                    <  1:1403226/91‑1 (MQ=255)

                                                                        cTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCGt                    <  1:1633595/91‑1 (MQ=255)

                                                                        cTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCGt                    <  1:426865/91‑1 (MQ=255)

                                                                         tGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCGTa                   >  1:2353064/1‑91 (MQ=255)

                                                                         tGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCCTGTAATCCTTCTGGAAAAGCCGTa                   >  2:3602272/1‑91 (MQ=255)

                                                                            ttGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCGTACca                <  1:3555830/91‑1 (MQ=255)

                                                                              gACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCGTACcaac              >  1:2266240/1‑91 (MQ=255)

                                                                              gACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCGTACcaac              >  1:3704986/1‑91 (MQ=255)

                                                                               acacGGTACCGGGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCGTACcaaca             >  2:650069/4‑91 (MQ=255)

                                                                                cccGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCGTACcaacaa            <  1:1994921/91‑1 (MQ=255)

                                                                                cccGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCGTACcaacaa            >  2:24635/1‑91 (MQ=255)

                                                                                 ccGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCGTACCAACAAt           <  1:3760117/91‑1 (MQ=255)

                                                                                 acGGTACCGGGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCGTACCAACAAt           >  2:445617/2‑91 (MQ=255)

                                                                                 acGGTACCGGGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCGTACCAACAAt           >  1:2008538/2‑91 (MQ=255)

                                                                                   ggTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCGTACCAACAATGa         >  1:3736493/1‑91 (MQ=255)

                                                                                   ggTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCGTACCAACAATGa         <  1:1740380/91‑1 (MQ=255)

                                                                                   ggTACCGGGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCGTACCAACAATGa         >  1:1797481/1‑91 (MQ=255)

                                                                                   ggTACCGGGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCGTACCAACAATGa         >  2:2069702/1‑91 (MQ=255)

                                                                                    gTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCGTACCAACAATGAc        >  1:2886543/1‑91 (MQ=255)

                                                                                    gTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCGTACCAACAATGAc        <  1:2092530/91‑1 (MQ=255)

                                                                                     tACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCGTACCAACAATGACa       <  2:766176/91‑1 (MQ=255)

                                                                                     tACCGGGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCGTACCAACAATGACa       >  2:3306184/1‑91 (MQ=255)

                                                                                     tACCGGGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCGTACCAACAATGACa       >  1:785062/1‑91 (MQ=255)

                                                                                      aCCGGGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCGTACCAACAATGACAt      >  1:2750359/1‑91 (MQ=255)

                                                                                      aCCGGGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCGTACCAACAATGACAt      >  2:3046756/1‑91 (MQ=255)

                                                                                        catgtCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCGTACCAACAATGACATCt    >  2:3561661/3‑91 (MQ=255)

                                                                                        cGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCGTACCAACAATGACATCt    >  2:131025/1‑91 (MQ=255)

                                                                                          tgtCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCGTACCAACAATGACATCTTc  >  1:3570211/1‑91 (MQ=255)

                                                                                          tgtCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCGTACCAACAATGACATCTTc  >  1:532107/1‑91 (MQ=255)

                                                                                          tgtCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCGTACCAACAATGACATCTTc  >  2:3109637/1‑91 (MQ=255)

                                                                                          tgtCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCGTACCAACAATGACATCTTc  >  2:2155363/1‑91 (MQ=255)

CTCCTTGACCAGCCATAATTTGATTGGCTGCAATCGCGATGGTTTGTAATCCTGATGAGCAGTAGCGATTCACTGTTTGACCCGGTACCGTGTCAGATAATCCCGCACGCAATGCAATCGTTCGTGCAATGTTTTGGCCTTGTAATCCTTCTGGAAAAGCCGTACCAACAATGACATCTTC  >  USA300TCH1516_ALE/269581‑269761

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