BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,077,629	C→A	6.8%	M171I (ATG→ATT)	yxlF_3 ←	putative ABC transporter ATP‑binding protein YxlF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,077,629	0	C	A	6.8%	360.2 / 14.5	103	M171I (ATG→ATT)	yxlF_3	putative ABC transporter ATP‑binding protein YxlF
Reads supporting (aligned to +/- strand): ref base C (56/40); new base A (5/2); total (61/42)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 107 total aligned reads displayed.

TTGAACAAAATGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGCTCATCCAAGATT  >  USA300TCH1516_ALE/2077542‑2077719

ttGAACAAAATGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCata                                                                                         >  1:1431312/1‑91 (MQ=255)

 tGAACAAAATGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCatat                                                                                        <  1:1621929/91‑1 (MQ=255)

  gAACAAAATGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATc                                                                                       >  1:3513461/1‑91 (MQ=255)

  gAACAAAATGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATc                                                                                       <  1:1567628/91‑1 (MQ=255)

  gAACAAAATGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATGAGAATTCTAATATc                                                                                       <  1:158712/91‑1 (MQ=255)

     cAAAATGTCCGTCTCTTAAAAATACTGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAg                                                                                    <  1:3871935/91‑1 (MQ=255)

      aaaaTGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGt                                                                                   <  1:412038/91‑1 (MQ=255)

      aaaaTGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATGAGAATTCTAATATCAAGt                                                                                   <  1:1683570/91‑1 (MQ=255)

       aaaTGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGtt                                                                                  >  1:509804/1‑91 (MQ=255)

        aaTGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTc                                                                                 <  1:1208479/91‑1 (MQ=255)

        aaTGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTc                                                                                 >  1:2779202/1‑91 (MQ=255)

        aaTGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTc                                                                                 >  1:2937516/1‑91 (MQ=255)

         aTGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGACTTCTCATATCAAGTTCa                                                                                >  1:3947348/1‑91 (MQ=255)

         aTGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCa                                                                                >  1:2187900/1‑91 (MQ=255)

         aTGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCa                                                                                >  1:3947763/1‑91 (MQ=255)

           gTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCAtt                                                                              <  1:452142/91‑1 (MQ=255)

           gTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAACTCAATATCTTCTACCTTATGACTTGATCTTAGAATTCTCATATCAAGCTCAtt                                                                              >  1:3852390/1‑91 (MQ=255)

            tCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCGttt                                                                             >  1:3346022/1‑91 (MQ=255)

             ccGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTa                                                                            >  1:1174055/1‑91 (MQ=255)

              cGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTAc                                                                           >  1:1651237/1‑91 (MQ=255)

                   ctTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAg                                                                      <  1:4105414/91‑1 (MQ=255)

                     tAAAAATACAGCTCTATGACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGAc                                                                    <  1:3551980/91‑1 (MQ=255)

                     tAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGAc                                                                    <  1:3166479/91‑1 (MQ=255)

                     tAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGAc                                                                    <  1:3162336/91‑1 (MQ=255)

                     tAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGAc                                                                    <  1:2839585/91‑1 (MQ=255)

                     tAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGAc                                                                    <  1:3162282/91‑1 (MQ=255)

                     tAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGAc                                                                    <  1:2162501/91‑1 (MQ=255)

                       aaaaTACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACtt                                                                  <  1:2967040/91‑1 (MQ=255)

                       aaaaTACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACtt                                                                  <  1:1357896/91‑1 (MQ=255)

                       aaaaTACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACtt                                                                  <  1:1357897/91‑1 (MQ=255)

                        aaaTACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTa                                                                 >  1:704981/1‑91 (MQ=255)

                          aTACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAAt                                                               <  1:3445281/91‑1 (MQ=255)

                          aTACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAAt                                                               >  1:3701054/1‑91 (MQ=255)

                           tACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAAtt                                                              <  1:480440/91‑1 (MQ=255)

                           tACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAAtt                                                              <  1:2004357/91‑1 (MQ=255)

                           tACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACGTGATATTAGATTTCTCATATCAAGTTCATTTACTAAAGACTTAAtt                                                              <  1:2004358/91‑1 (MQ=255)

                             cAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGt                                                            >  1:3596456/1‑91 (MQ=255)

                             cAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGt                                                            >  1:2280317/1‑91 (MQ=255)

                             cAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGt                                                            <  1:2713467/91‑1 (MQ=255)

                             cAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGt                                                            >  1:146434/1‑91 (MQ=255)

                              aGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTa                                                           >  1:1620891/1‑91 (MQ=255)

                               gCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAg                                                          >  1:842443/1‑91 (MQ=255)

                                      cacaAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGcac                                                   <  1:2877014/91‑1 (MQ=255)

                                       acaAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGcaca                                                  >  1:2776631/1‑91 (MQ=255)

                                       acaAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGcaca                                                  >  1:2776489/1‑91 (MQ=255)

                                       acaAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGcaca                                                  <  1:469589/91‑1 (MQ=255)

                                        caAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACAt                                                 >  1:3979759/1‑91 (MQ=255)

                                         aaaTTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATc                                                >  1:3599326/1‑91 (MQ=255)

                                           attaattCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCag                                              >  1:1492642/1‑90 (MQ=255)

                                              aattCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTg                                           >  1:1360154/1‑91 (MQ=255)

                                                ttCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATCTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAg                                         <  1:3503657/91‑1 (MQ=255)

                                                 tCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCTCATCAATTGAGc                                        <  1:3177590/91‑1 (MQ=255)

                                                 tCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGc                                        <  1:1286192/91‑1 (MQ=255)

                                                 tCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGc                                        >  1:474274/1‑91 (MQ=255)

                                                 tCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGc                                        >  1:505995/1‑91 (MQ=255)

                                                 tCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGc                                        >  1:2134152/1‑91 (MQ=255)

                                                 tCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGc                                        >  1:2397429/1‑91 (MQ=255)

                                                     tatCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATc                                    <  1:726027/91‑1 (MQ=255)

                                                        cttctAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTgg                                 <  1:1115004/91‑1 (MQ=255)

                                                        cttctAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTgg                                 >  1:690011/1‑91 (MQ=255)

                                                         ttctAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGa                                >  1:2225115/1‑91 (MQ=255)

                                                              aCTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCAt                           <  1:1667398/91‑1 (MQ=255)

                                                              aCTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCAt                           >  1:3097064/1‑91 (MQ=255)

                                                              aCTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCAt                           <  1:2240508/91‑1 (MQ=255)

                                                              aCTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCAt                           <  1:1797945/91‑1 (MQ=255)

                                                              aCTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCAt                           >  1:2061582/1‑91 (MQ=255)

                                                              aCTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCAt                           >  1:844009/1‑91 (MQ=255)

                                                                 tATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATAcc                        >  1:422927/1‑91 (MQ=255)

                                                                 tATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATAcc                        >  1:3028901/1‑91 (MQ=255)

                                                                  aTGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCa                       >  1:3480842/1‑91 (MQ=255)

                                                                  aTGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCa                       >  1:3892249/1‑91 (MQ=255)

                                                                      cTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTg                   <  1:4094126/91‑1 (MQ=255)

                                                                      cTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTg                   >  1:2022422/1‑91 (MQ=255)

                                                                       ttGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGt                  >  1:2600228/1‑91 (MQ=255)

                                                                         gATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAg                >  1:2153355/1‑91 (MQ=255)

                                                                         gATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAg                >  1:70053/1‑91 (MQ=255)

                                                                         gATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAg                >  1:4043519/1‑91 (MQ=255)

                                                                          atatTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAgg               >  1:1849738/1‑91 (MQ=255)

                                                                          atatGAGAATTCTAATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAgg               >  1:3825363/1‑91 (MQ=255)

                                                                          atatGAGAATTCTAATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAgg               >  1:3982242/1‑91 (MQ=255)

                                                                           tatTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGc              <  1:3434539/91‑1 (MQ=255)

                                                                           tatTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGc              >  1:220769/1‑91 (MQ=255)

                                                                           tatTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGc              <  1:2286330/91‑1 (MQ=255)

                                                                           tatTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGc              <  1:2286728/91‑1 (MQ=255)

                                                                             ttAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGCTc            >  1:3567113/1‑91 (MQ=255)

                                                                               aGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGCTCAt          >  1:849941/1‑91 (MQ=255)

                                                                               aGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGCTCAt          >  1:1246824/1‑91 (MQ=255)

                                                                               aGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGCTCAt          >  1:4026621/1‑91 (MQ=255)

                                                                               aGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGCTCAt          <  1:146152/91‑1 (MQ=255)

                                                                               aGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGCTCAt          >  1:3747806/1‑91 (MQ=255)

                                                                               aGAATTCTAATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGCTCAt          >  1:1560660/1‑91 (MQ=255)

                                                                               aGAATTCTAATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGCTCAt          >  1:2878024/1‑91 (MQ=255)

                                                                                gAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGCGCCATCTGGATCCATACCATTTGTAGGCTCATc         >  1:3395229/1‑91 (MQ=255)

                                                                                gAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGCTCATc         <  1:3267661/91‑1 (MQ=255)

                                                                                gAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGCAGGCTCATc         <  1:3268774/91‑1 (MQ=255)

                                                                                 aaTTCTAATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGCTCATcc        >  1:3265194/1‑91 (MQ=255)

                                                                                     ctcATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGCTCATCCAAGa    >  1:3612586/1‑91 (MQ=255)

                                                                                     ctcATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGCTCATCCAAGa    >  1:493343/1‑91 (MQ=255)

                                                                                     ctcATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGCTCATCCAAGa    >  1:2950334/1‑91 (MQ=255)

                                                                                       cATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGCTCATCCAAGAtt  <  1:3539483/91‑1 (MQ=255)

TTGAACAAAATGTCCGTCTCTTAAAAATACAGCTCTATCACAAATTAATTCAATATCTTCTAACTTATGACTTGATATTAGAATTCTCATATCAAGTTCATTTACTAAAGACTTAATTGTAGTCAGCACATCAATTGAGCCATCTGGATCCATACCATTTGTAGGCTCATCCAAGATT  >  USA300TCH1516_ALE/2077542‑2077719

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