BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,049,164	C→A	8.8%	E311D (GAG→GAT)	ligA ←	DNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,049,164	0	C	A	8.8%	347.6 / 19.4	102	E311D (GAG→GAT)	ligA	DNA ligase
Reads supporting (aligned to +/- strand): ref base C (55/38); new base A (5/4); total (60/42)
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.92e-01

Only 100 of 106 total aligned reads displayed.

TACTTTTACTGGTTCTAAAATAGCAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATCCCATCTCATCCTGTTGATCTAAATCATTAACCTTAAT  >  USA300TCH1516_ALE/2049074‑2049253

tACTTTTACTGGTTCTAAAATAGCAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTcc                                                                                           >  1:2256844/1‑91 (MQ=255)

  cTTTTACTGGTTCTAAAATAGCAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTc                                                                                         <  1:3281769/91‑1 (MQ=255)

  cTTTTACTGGTTCTAAAATAGCAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTc                                                                                         >  1:3529950/1‑91 (MQ=255)

   ttttACTGGTTCTAAAATAGCAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCa                                                                                        >  1:776699/1‑91 (MQ=255)

    tttACTGGTTCTAAAATAGCAGTAGGTGTGACTACACCTGTTCGTCCCATACTTAATTCAATATCTAATAATTTAGTCACTACTTCCTCAg                                                                                       >  1:4148760/1‑91 (MQ=255)

     ttACTGGTTCTAAAATAGCAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGc                                                                                      >  1:1911789/1‑91 (MQ=255)

     ttACTGGTTCTAAAATAGCAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGc                                                                                      <  1:1168332/91‑1 (MQ=255)

     ttACTGGTTCTAAAATAGCAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGc                                                                                      >  1:2862753/1‑91 (MQ=255)

      tACTGGTTCTAAAATAGCAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCt                                                                                     >  1:3213395/1‑91 (MQ=255)

        cTGGTTCTAAAATAGCAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTgg                                                                                   >  1:1680417/1‑91 (MQ=255)

         tGGTTCTAAAATAGCAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGa                                                                                  <  1:1764639/91‑1 (MQ=255)

         tGGTTCTAAAATAGCAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGa                                                                                  <  1:2275511/91‑1 (MQ=255)

           gTTCTAAAATAGCAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCATCAGCTGGaaa                                                                                >  1:1658724/1‑91 (MQ=255)

            ttCTAAAATAGCAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAAt                                                                               >  1:3825680/1‑91 (MQ=255)

            ttCTAAAATAGCAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCATCAGCTGGAAAt                                                                               >  1:855482/1‑91 (MQ=255)

             tCTAAAATAGCAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAAtt                                                                              >  1:1540051/1‑91 (MQ=255)

                 aaaTAGCAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTtata                                                                          <  1:1528806/91‑1 (MQ=255)

                  aaTAGCAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATaa                                                                         >  1:2354808/1‑91 (MQ=255)

                  aaTAGCAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATaa                                                                         >  1:2970288/1‑91 (MQ=255)

                    tAGCAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGTAAATTTATAAGc                                                                       <  1:3189140/91‑1 (MQ=255)

                    tAGCAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAga                                                                       >  1:2974492/1‑90 (MQ=255)

                     aGCAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCa                                                                      >  1:921935/1‑91 (MQ=255)

                      gCAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCaa                                                                     <  1:1391808/91‑1 (MQ=255)

                      gCAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCaa                                                                     >  1:536920/1‑91 (MQ=255)

                       cAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAAt                                                                    >  1:2957446/1‑91 (MQ=255)

                       cAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAAt                                                                    >  1:3894432/1‑91 (MQ=255)

                       cAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAAt                                                                    <  1:3930781/91‑1 (MQ=255)

                       cAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAAt                                                                    >  1:3895082/1‑91 (MQ=255)

                       cAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAAt                                                                    >  1:3894929/1‑91 (MQ=255)

                       cAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAAt                                                                    >  1:3894367/1‑91 (MQ=255)

                       cAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAAt                                                                    >  1:16241/1‑91 (MQ=255)

                       cAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGAAAATTTATAAGCAAt                                                                    >  1:159687/1‑91 (MQ=255)

                        aGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATg                                                                   >  1:2713634/1‑91 (MQ=255)

                         gTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATgg                                                                  >  1:760152/1‑91 (MQ=255)

                          tAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCATCAGCTGGAAATTTATAAGCAATGGc                                                                 >  1:2743264/1‑91 (MQ=255)

                           aGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCATCAGCTGGAAATTTATAAGCAATGGcc                                                                <  1:3368444/91‑1 (MQ=255)

                            ggTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGccc                                                               >  1:3111590/1‑91 (MQ=255)

                            ggTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGccc                                                               >  1:3110013/1‑91 (MQ=255)

                            ggTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGccc                                                               >  1:3664551/1‑91 (MQ=255)

                             gtgtGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCa                                                              <  1:897910/91‑1 (MQ=255)

                              tgtgACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCAt                                                             <  1:3183647/91‑1 (MQ=255)

                                tgACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCt                                                           <  1:3608726/91‑1 (MQ=255)

                                tgACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCt                                                           <  1:3308771/91‑1 (MQ=255)

                                  aCTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAg                                                         >  1:1828891/1‑91 (MQ=255)

                                   cTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAgg                                                        <  1:1182600/91‑1 (MQ=255)

                                    tACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGga                                                       <  1:56574/91‑1 (MQ=255)

                                    tACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGga                                                       <  1:1968884/91‑1 (MQ=255)

                                      cacCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGgaga                                                     <  1:2564779/91‑1 (MQ=255)

                                      cacCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGgaga                                                     >  1:199035/1‑91 (MQ=255)

                                      cacCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGgaga                                                     <  1:3096516/91‑1 (MQ=255)

                                       acCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGAt                                                    >  1:2493464/1‑91 (MQ=255)

                                        ccTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGAtt                                                   >  1:3748578/1‑91 (MQ=255)

                                         cTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGAttt                                                  <  1:3155824/91‑1 (MQ=255)

                                         cTGTTCGTCCAATACTTAATTCAATATCTAAGAATTTAGTTACTGCTTCATCAGCTGCAAATTTATAAGCAATGGCCCATCTAGGAGAttt                                                  <  1:3310419/91‑1 (MQ=255)

                                         cTGTTCGTCCAATACTTAATTCAATATCTAAGAATTTAGTTACTGCTTCATCACCTGCAAATTTATAAGCAATGGCCCATCTAGGAGAttt                                                  <  1:3310449/91‑1 (MQ=255)

                                           gTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCATCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGAttttt                                                >  1:2924415/1‑91 (MQ=255)

                                             tCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTtgt                                              <  1:1820719/91‑1 (MQ=255)

                                             tCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTtgt                                              <  1:1821476/91‑1 (MQ=255)

                                                tCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTtgtgtg                                           >  1:815240/1‑91 (MQ=255)

                                                tCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTtgtgtg                                           >  1:2249831/1‑91 (MQ=255)

                                                tCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTtgtgtg                                           >  1:1691250/1‑91 (MQ=255)

                                                tCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTtgtgtg                                           >  1:1691305/1‑91 (MQ=255)

                                                 ccAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGa                                          <  1:394154/91‑1 (MQ=255)

                                                 ccAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGa                                          <  1:394155/91‑1 (MQ=255)

                                                  cAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGaa                                         >  1:1624383/1‑91 (MQ=255)

                                                  cAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGaa                                         >  1:3448551/1‑91 (MQ=255)

                                                    aTACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATc                                       >  1:2903674/1‑91 (MQ=255)

                                                    aTACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATc                                       >  1:1368348/1‑91 (MQ=255)

                                                     tACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATcc                                      >  1:4104700/1‑91 (MQ=255)

                                                      aCTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATccc                                     >  1:938843/1‑91 (MQ=255)

                                                         tAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATCCCAtc                                  >  1:2401409/1‑91 (MQ=255)

                                                          aaTTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATCCCAtct                                 <  1:709548/91‑1 (MQ=255)

                                                             tCAATATCTAATAATTTAGTTACTAATTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATCCCATCTCAt                              <  1:3846283/91‑1 (MQ=255)

                                                                atatCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATCCCATCTCATCCt                           >  1:3131674/1‑91 (MQ=255)

                                                                atatCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGAATTTTGTGTGAATCCCATCTCATCCt                           >  1:3853620/1‑91 (MQ=255)

                                                                 tatCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATCCCATCTCATCCTg                          >  1:537693/1‑91 (MQ=255)

                                                                 tatCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATCCCATCTCATCCTg                          >  1:825366/1‑91 (MQ=255)

                                                                 tatCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATCCCATCTCATCCTg                          <  1:1095897/91‑1 (MQ=255)

                                                                 tatCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATCCCATCTCATCCTg                          >  1:3891859/1‑91 (MQ=255)

                                                                 tatCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATCCCATCTCATCCTg                          <  1:3029227/91‑1 (MQ=255)

                                                                 tatCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATCCCATCTCATCCTg                          >  1:2164777/1‑91 (MQ=255)

                                                                   tCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATCCCATCTCATCCTGtt                        <  1:2909270/91‑1 (MQ=255)

                                                                    cTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATCCCATCTCATCCTGTTg                       >  1:1331375/1‑91 (MQ=255)

                                                                    cTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATCCCATCTCATCCTGTTg                       <  1:2452913/91‑1 (MQ=255)

                                                                    cTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATCCCATCTCATCCTGTTg                       <  1:1541972/91‑1 (MQ=255)

                                                                     taataaTTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATCCCATCTCATCCTGTTGa                      <  1:1033782/91‑1 (MQ=255)

                                                                     taataaTTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATCCCATCTCATCCTGTTGa                      >  1:3011185/1‑91 (MQ=255)

                                                                              aGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATCCCATCTCATCCTGTTGATCTAAATCa             >  1:2466501/1‑91 (MQ=255)

                                                                               gTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATCCCATCTCATCCTGTTGATCTAAATCAt            >  1:1170876/1‑91 (MQ=255)

                                                                               gTTACTAATTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATCCCATCTCATCCTGTTGATCTAAATCAt            <  1:2073626/91‑1 (MQ=255)

                                                                                ttACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATCCCATCTCATCCTGTTGATCTAAATCAtt           <  1:638941/91‑1 (MQ=255)

                                                                                 tactacTTCATCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATCCCATCTCATCCTGTTGATCTAAATCATTa          <  1:3595610/91‑1 (MQ=255)

                                                                                  actactTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATCCCATCTCATCCTGTTGATCTAAATCATTaa         <  1:373474/91‑1 (MQ=255)

                                                                                  actactTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATCCCATCTCATCCTGTTGATCTAAATCATTaa         <  1:1464429/91‑1 (MQ=255)

                                                                                     actTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATCCCATCTCATCCTGTTGATCTAAATCATTAACCt      >  1:530463/1‑91 (MQ=255)

                                                                                      ctTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATCCCATCTCATCCTGTTGATCTAAATCATTAACCtt     <  1:252625/91‑1 (MQ=255)

                                                                                       ttCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATCCCATCTCATCCTGTTGATCTAAATCATTAACCTta    <  1:1821990/91‑1 (MQ=255)

                                                                                         ccTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATCCCATCTCATCCTGTTGATCTAAATCATTAACCTtaat  >  1:2179635/1‑91 (MQ=255)

                                                                                         ccTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATCCCATCTCATCCTGTTGATCTAAATCATTAACCTtaat  <  1:1350252/91‑1 (MQ=255)

                                                                                         ccTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATCCCATCTCATCCTGTTGATCTAAATCATTAACCTtaat  >  1:3916933/1‑91 (MQ=255)

TACTTTTACTGGTTCTAAAATAGCAGTAGGTGTGACTACACCTGTTCGTCCAATACTTAATTCAATATCTAATAATTTAGTTACTACTTCCTCAGCTGGAAATTTATAAGCAATGGCCCATCTAGGAGATTTTTGTGTGAATCCCATCTCATCCTGTTGATCTAAATCATTAACCTTAAT  >  USA300TCH1516_ALE/2049074‑2049253

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