BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,060,235	T→A	10.1%	K568* (AAA→TAA) ‡	atl_1 ←	Bifunctional autolysin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,060,235	0	T	A	10.1%	678.6 / 55.0	218	K568* (AAA→TAA) ‡	atl_1	Bifunctional autolysin
Reads supporting (aligned to +/- strand): ref base T (94/102); new base A (5/17); total (99/119)
Fisher's exact test for biased strand distribution p-value = 2.55e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 229 total aligned reads displayed.

GTGTAGGCGTAGGTTTAGCAGTATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGACACACTACCAGCAACTTGTTTAGATGTACCCCAAGGTAC  >  USA300TCH1516_ALE/1060145‑1060325

gtgtAGGCGTAGGTTTAGCAGTATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGAttt                                                                                            >  1:3645294/1‑91 (MQ=255)

 tgtAGGCGTAGGTTTAGCAGTATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTa                                                                                           <  2:530784/91‑1 (MQ=255)

   tAGGCGTAGGTTTAGCAGTATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATc                                                                                         <  1:2696699/91‑1 (MQ=255)

     ggCGTAGGTTTAGCAGTATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCaa                                                                                       <  2:2219475/91‑1 (MQ=255)

      gCGTAGGTTTCGCAGTATCCACTAAATCTGCTTTACTTACCCAACCAGATTTCCCATTCACAGAGCCATATAACTAACTTGATTTATCatt                                                                                      >  2:3595461/1‑89 (MQ=21)

      gCGTAGGTTTAGCAGTATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAAt                                                                                      <  2:843206/91‑1 (MQ=255)

      gCGTAGGTTTAGCAGTATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAAt                                                                                      >  1:2274213/1‑91 (MQ=255)

        gTAGGTTTAGCAGTATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAAttt                                                                                    <  2:1871516/91‑1 (MQ=255)

        gTAGGTTTAGCAGTATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATAAATCAAttt                                                                                    <  1:1672628/91‑1 (MQ=255)

        gTAGGTTTAGCAGTATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATAAATCAAttt                                                                                    <  2:2724196/91‑1 (MQ=255)

         tAGGTTTAGCAGTATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATttg                                                                                   >  1:2609837/1‑91 (MQ=255)

         tAGGTTTAGCAGTATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATttg                                                                                   >  2:1337816/1‑91 (MQ=255)

         tAGGTTTAGCAGTATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATttg                                                                                   >  1:3191538/1‑91 (MQ=255)

         tAGGTTTAGCAGTATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATAAATCAAttta                                                                                   <  1:1916402/91‑2 (MQ=255)

         tAGGTTTAGCAGTATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATAAATCAAttta                                                                                   <  1:3246338/91‑2 (MQ=255)

         tAGGTTTAGCAGTATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATAAATCAAttta                                                                                   <  1:3257923/91‑2 (MQ=255)

         tAGGTTTAGCAGTATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATAAATCAAttta                                                                                   <  1:1390400/91‑2 (MQ=255)

         tAGGTTTAGCAGTATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATAAATCAAttta                                                                                   <  1:552266/91‑2 (MQ=255)

          aGGTTTAGCAGTATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATttgt                                                                                  <  1:582038/91‑1 (MQ=255)

          aGGTTTAGCAGTATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATttgt                                                                                  <  2:1641468/91‑1 (MQ=255)

           ggTTTAGCAGTATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATttgtt                                                                                 <  1:3558665/91‑1 (MQ=255)

           ggTTTAGCAGTATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATttgtt                                                                                 >  2:2187081/1‑91 (MQ=255)

             tttAGCAGTATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATttgttgt                                                                               <  1:919480/91‑1 (MQ=255)

             tttAGCAGTATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATttgttgt                                                                               <  2:480696/91‑1 (MQ=255)

              ttAGCAGTATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATttgttgtt                                                                              >  2:1252049/1‑91 (MQ=255)

              ttAGCAGTATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATttgttgtt                                                                              <  1:2039225/91‑1 (MQ=255)

              ttAGCAGTATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATttgttgtt                                                                              >  1:2892040/1‑91 (MQ=255)

               tAGCAGTATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATttgttgttg                                                                             >  1:552947/1‑91 (MQ=255)

                 gCAGTATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCt                                                                           >  2:944462/1‑91 (MQ=255)

                  cAGTATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCtt                                                                          <  2:3376914/91‑1 (MQ=255)

                  cAGTATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCtt                                                                          <  1:1337816/91‑1 (MQ=255)

                     tATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGa                                                                       >  2:499386/1‑91 (MQ=255)

                     tATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGa                                                                       <  1:1186786/91‑1 (MQ=255)

                          aCTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCt                                                                  <  2:3028244/91‑1 (MQ=255)

                          aCTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCt                                                                  >  2:2666819/1‑91 (MQ=255)

                           cTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCtt                                                                 <  2:2274213/91‑1 (MQ=255)

                           cTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCtt                                                                 <  2:5430/91‑1 (MQ=255)

                            tAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTa                                                                >  1:167808/1‑91 (MQ=255)

                             aaaTATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTaa                                                               <  1:1443163/91‑1 (MQ=255)

                              aaTATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTaaa                                                              >  1:2405760/1‑91 (MQ=255)

                                tatGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATg                                                            >  2:802581/1‑91 (MQ=255)

                                tatGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATg                                                            <  1:2488738/91‑1 (MQ=255)

                                  tGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGtt                                                          <  2:1790090/91‑1 (MQ=255)

                                  tGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGtt                                                          >  2:1698440/1‑91 (MQ=255)

                                    cTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTg                                                        <  2:2940969/91‑1 (MQ=255)

                                       tacttacCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGtt                                                     >  1:1782749/1‑91 (MQ=255)

                                       tacttacCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGtt                                                     <  2:151213/91‑1 (MQ=255)

                                         cttacCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTc                                                   <  1:1851684/91‑1 (MQ=255)

                                          ttacCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTcc                                                  >  1:333967/1‑91 (MQ=255)

                                           tacCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCa                                                 <  2:3631027/91‑1 (MQ=255)

                                             cccAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCaga                                               >  1:519700/1‑91 (MQ=255)

                                              ccAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCagag                                              >  2:3157461/1‑91 (MQ=255)

                                              ccAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCagag                                              >  2:3811024/1‑91 (MQ=255)

                                              ccAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCagag                                              >  1:1205998/1‑91 (MQ=255)

                                               cAACCGGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGc                                             <  2:3142132/91‑1 (MQ=255)

                                               cAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGc                                             <  2:2058608/91‑1 (MQ=255)

                                               cAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGc                                             <  2:2321468/91‑1 (MQ=255)

                                               cAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGc                                             <  2:2055192/91‑1 (MQ=255)

                                                aaCCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGttt                                                    >  1:2104928/1‑83 (MQ=255)

                                                 aCCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCa                                           >  2:699455/1‑91 (MQ=255)

                                                 aCCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCa                                           >  2:1978771/1‑91 (MQ=255)

                                                  ccAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAg                                          <  1:3771223/91‑1 (MQ=255)

                                                   cAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGa                                         <  1:360426/91‑1 (MQ=255)

                                                     gATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGaca                                       >  1:2535414/1‑91 (MQ=255)

                                                       tttACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGacaca                                     >  2:1599547/1‑91 (MQ=255)

                                                       tttACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGacaca                                     >  2:1455466/1‑91 (MQ=255)

                                                       tttACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGacaca                                     >  1:410273/1‑91 (MQ=255)

                                                        ttACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGacacac                                    >  2:694927/1‑91 (MQ=255)

                                                         tACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGACACACt                                   <  2:1051913/91‑1 (MQ=255)

                                                         tACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGACACACt                                   <  2:2632010/91‑1 (MQ=255)

                                                         tACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGACACACt                                   >  1:3737269/1‑91 (MQ=255)

                                                           ccATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGACACACTAc                                 <  1:2434528/91‑1 (MQ=255)

                                                              ttCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGACACACTACCAg                              <  2:3703579/91‑1 (MQ=255)

                                                               tCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGACACACTACCAGc                             >  2:1936040/1‑91 (MQ=255)

                                                               tCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGACACACTACCAGc                             >  1:1228606/1‑91 (MQ=255)

                                                                cacaGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGACACACTACCAGCa                            <  2:1721920/91‑1 (MQ=255)

                                                                   agagCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGACACACTACCAGCAACt                         >  2:3535952/1‑91 (MQ=255)

                                                                    gagCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCag                                                <  2:3417450/67‑1 (MQ=255)

                                                                    gagCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCag                                                >  1:3417450/1‑67 (MQ=255)

                                                                       ccATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGACACACTACCAGCAACTTGtt                     >  1:1000669/1‑91 (MQ=255)

                                                                         atatAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGACACACTACCAGCAACTTGTTTa                   <  1:997065/91‑1 (MQ=255)

                                                                         atatAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGACACACTACCAGCAACTTGTTTa                   >  2:3667514/1‑91 (MQ=255)

                                                                             aaataaatTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGACACACTACCAGCAACTTGTTTAGATg               <  2:1754003/91‑1 (MQ=255)

                                                                              aataaatTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGACACACTACCAGCAACTTGTTTAGATGt              >  1:27766/1‑91 (MQ=255)

                                                                              aataaatTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGACACACTACCAGCAACTTGTTTAGATGt              <  1:967343/91‑1 (MQ=255)

                                                                                taaatTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGACACACTACCAGCAACTTGTTTAGATGTAc            <  1:219832/91‑1 (MQ=255)

                                                                                caaatTGATAAATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGACACACTACCAGCAACTTGTTTAGATGTAc            >  1:3808858/2‑91 (MQ=255)

                                                                                  aatTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGTGCCTGACACACTACCAGCAACTTGTTTAGATGTAccc          <  1:3095848/91‑1 (MQ=255)

                                                                                  aatTGATAAATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGACACACTACCAGCAACTTGTTTAGATGTAccc          >  1:602871/1‑91 (MQ=255)

                                                                                   atTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGACACACTACCAGCAACTTGTTTAGATGTAcccc         <  1:3868866/91‑1 (MQ=255)

                                                                                      gATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGACACACTACCAGCAACTTGTTTAGATGTACCCCAAg      <  1:1593545/91‑1 (MQ=255)

                                                                                      gATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGACACACTACCAGCAACTTGTTTAGATGTACCCCAAg      >  2:2912737/1‑91 (MQ=255)

                                                                                       aTTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGACACACTACCAGCAACTTGTTTAGATGTACCCCAAgg     <  2:2378824/91‑1 (MQ=255)

                                                                                        tttATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGACACACTACCAGCAACTTGTTTAGATGTACCCCAAGGt    <  2:2956146/91‑1 (MQ=255)

                                                                                        tttATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGACACACTACCAGCAACTTGTTTAGATGTACCCCAAGGt    >  2:1608419/1‑91 (MQ=255)

                                                                                        tttATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGACACACTACCAGCAACTTGTTTAGATGTACCCCAAGGt    <  1:717992/91‑1 (MQ=255)

                                                                                         ttATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGACACACTACCAGCAACTTGTTTAGATGTACCCCAAGGTa   <  2:823140/91‑1 (MQ=255)

                                                                                          tATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGACACACTACCAGCAACTTGTTTAGATGTACCCCAAGGTAc  <  2:1166121/91‑1 (MQ=255)

                                                                                          tATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGACACACTACCAGCAACTTGTTTAGATGTACCCCAAGGTAc  >  1:1830186/1‑91 (MQ=255)

                                                                                          tATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGACACACTACCAGCAACTTGTTTAGATGTACCCCAAGGTAc  >  2:1078010/1‑91 (MQ=255)

GTGTAGGCGTAGGTTTAGCAGTATCAACTAAATATGCTTTACTTACCCAACCAGATTTACCATTCACAGAGCCATATAAATAAATTGATTTATCAATTTGTTGTTGCTTTGAAGCCTTAAATGTTTGGTTTCCAGAGCCAGACACACTACCAGCAACTTGTTTAGATGTACCCCAAGGTAC  >  USA300TCH1516_ALE/1060145‑1060325

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