BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,530,786	A→G	11.1%	intergenic (‑37/+236)	yefM ← / ← bdbD	Antitoxin YefM/Disulfide bond formation protein D

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,530,786	0	A	G	11.1%	308.1 / 24.1	108	intergenic (‑37/+236)	yefM/bdbD	Antitoxin YefM/Disulfide bond formation protein D
Reads supporting (aligned to +/- strand): ref base A (51/45); new base G (4/8); total (55/53)
Fisher's exact test for biased strand distribution p-value = 2.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

Only 100 of 114 total aligned reads displayed.

TGTCATAAGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGATTC  >  USA300TCH1516_ALE/2530696‑2530876

tGTCATAAGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACaa                                                                                            <  1:80874/91‑1 (MQ=255)

tGTCATAAGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACaa                                                                                            <  1:1417711/91‑1 (MQ=255)

     tAAGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATc                                                                                       <  1:1662400/91‑1 (MQ=255)

     tAAGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATc                                                                                       <  1:2945896/91‑1 (MQ=255)

      aaGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCt                                                                                      <  1:3736794/91‑1 (MQ=255)

      aaGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCt                                                                                      >  1:1982401/1‑91 (MQ=255)

      aaGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCt                                                                                      >  1:3152955/1‑91 (MQ=255)

       aGTGCGTTTAAATTCTGTCGAGCGTATGAATAATTTTTATTAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAGatatttg                                                                                     <  1:3856232/91‑4 (MQ=255)

       aGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAGatatttg                                                                                     <  1:3835387/91‑4 (MQ=255)

       aGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAGatatttg                                                                                     <  1:3825360/91‑4 (MQ=255)

       aGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAGatatttg                                                                                     <  1:3838692/91‑4 (MQ=255)

       aGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAGatatttg                                                                                     <  1:3137046/91‑4 (MQ=255)

       aGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTg                                                                                     >  1:1552112/1‑91 (MQ=255)

        gTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAGATATTTGt                                                                                    <  1:1530251/91‑1 (MQ=255)

        gTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAGATATTTGt                                                                                    <  1:1689228/91‑1 (MQ=255)

        gTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAGATATTTGt                                                                                    <  1:4093979/91‑1 (MQ=255)

        gTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGt                                                                                    <  1:786650/91‑1 (MQ=255)

        gTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGt                                                                                    <  1:3665359/91‑1 (MQ=255)

         tGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTa                                                                                   <  1:3404160/91‑1 (MQ=255)

          gCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTAc                                                                                  >  1:1488122/1‑91 (MQ=255)

           ccTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACa                                                                                 >  1:1346496/1‑91 (MQ=255)

           ccTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACa                                                                                 >  1:2090392/1‑91 (MQ=255)

            cTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACaa                                                                                >  1:1804686/1‑91 (MQ=255)

            cTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACaa                                                                                >  1:1883214/1‑91 (MQ=255)

             tttAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAAt                                                                               >  1:293298/1‑91 (MQ=255)

             tttAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAAt                                                                               >  1:165237/1‑91 (MQ=255)

             tttAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAAt                                                                               <  1:1098339/91‑1 (MQ=255)

               tAACTTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATaa                                                                             <  1:3728432/91‑1 (MQ=255)

               tAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATaa                                                                             >  1:3281982/1‑91 (MQ=255)

               tAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATaa                                                                             >  1:3640107/1‑91 (MQ=255)

               tAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATaa                                                                             >  1:3640484/1‑91 (MQ=255)

                     cTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAg                                                                       >  1:2354008/1‑91 (MQ=255)

                     cTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAg                                                                       >  1:246318/1‑91 (MQ=255)

                      tGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGc                                                                      >  1:282723/1‑91 (MQ=255)

                      tGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGc                                                                      >  1:3776534/1‑91 (MQ=255)

                       gTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCg                                                                     >  1:3630208/1‑91 (MQ=255)

                        tCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGa                                                                    <  1:2653188/91‑1 (MQ=255)

                        tCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGa                                                                    <  1:2653172/91‑1 (MQ=255)

                        tCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGa                                                                    <  1:730190/91‑1 (MQ=255)

                           aGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGacac                                                                 <  1:2810450/91‑1 (MQ=255)

                            gCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTCTAACACCCCCCAATTGTACAAATATCTGTACAATAAAACAAGCGAcaca                                                                <  1:1123766/91‑1 (MQ=255)

                             cGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAg                                                               <  1:2446697/91‑1 (MQ=255)

                              gTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGa                                                              <  1:360879/91‑1 (MQ=255)

                                         ttttAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAactact                                                   <  1:938589/91‑1 (MQ=255)

                                          tttAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAACCAAGCGACACAGAAAAAAActacta                                                  >  1:3851404/1‑91 (MQ=255)

                                          tttAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAACCAAGCGACACAGAAAAAAActacta                                                  >  1:3851272/1‑91 (MQ=255)

                                          tttAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAActacta                                                  >  1:2682324/1‑91 (MQ=255)

                                          tttAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAActacta                                                  <  1:3742827/91‑1 (MQ=255)

                                          tttAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAActacta                                                  >  1:3259619/1‑91 (MQ=255)

                                            taataaTCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATACAACAAGCGACACAGAAAAAAACTACTAAt                                                >  1:3330514/1‑91 (MQ=255)

                                                   cATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGAtt                                         >  1:2106037/1‑91 (MQ=255)

                                                   cATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGAtt                                         >  1:1607716/1‑91 (MQ=255)

                                                      atTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCg                                      <  1:516029/91‑1 (MQ=255)

                                                          aTCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACg                                  <  1:236092/91‑1 (MQ=255)

                                                          aTCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACg                                  <  1:1734308/91‑1 (MQ=255)

                                                          aTCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACg                                  >  1:2152655/1‑91 (MQ=255)

                                                          aTCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACg                                  <  1:653371/91‑1 (MQ=255)

                                                          aTCGCCTCCTTTTAACACACCTCAAATGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACg                                  <  1:3871654/91‑1 (MQ=255)

                                                           tCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACgg                                 <  1:1071660/91‑1 (MQ=255)

                                                           tCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACgg                                 >  1:1281438/1‑91 (MQ=255)

                                                             gCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACgggg                               <  1:3716246/91‑1 (MQ=255)

                                                             gCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACgggg                               <  1:2775399/91‑1 (MQ=255)

                                                             gCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACgggg                               <  1:2774697/91‑1 (MQ=255)

                                                             gCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACgggg                               <  1:2435777/91‑1 (MQ=255)

                                                              cctcctTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGga                              <  1:466333/91‑1 (MQ=255)

                                                                tcctTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGgaag                            >  1:1029809/1‑91 (MQ=255)

                                                                tcctTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGgaag                            >  1:1029590/1‑91 (MQ=255)

                                                                  ctTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGgaagaa                          >  1:108251/1‑91 (MQ=255)

                                                                     ttAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCa                       <  1:2531999/91‑1 (MQ=255)

                                                                      tAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCaa                      >  1:1771291/1‑91 (MQ=255)

                                                                       aaCACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAAt                     <  1:693682/91‑1 (MQ=255)

                                                                       aaCACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAAt                     >  1:600213/1‑91 (MQ=255)

                                                                       aaCACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAAt                     <  1:2207132/91‑1 (MQ=255)

                                                                        acacacCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATa                    >  1:2734342/1‑91 (MQ=255)

                                                                         cacacCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATAc                   >  1:4124080/1‑91 (MQ=255)

                                                                           cacCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCt                 >  1:3712670/1‑91 (MQ=255)

                                                                           cacCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCt                 >  1:2634590/1‑91 (MQ=255)

                                                                            acCTCAATTGTACAAATATCTGTGCAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTa                <  1:1722801/91‑1 (MQ=255)

                                                                            acCTCAATTGTACAAATATCTGTACAATAAAACAAGCAACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTa                >  1:3500792/1‑91 (MQ=255)

                                                                                cAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCt            >  1:2628445/1‑91 (MQ=255)

                                                                                cAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCt            <  1:2323569/91‑1 (MQ=255)

                                                                                  aTTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGt          >  1:1006733/1‑91 (MQ=255)

                                                                                   ttGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTa         >  1:2886299/1‑91 (MQ=255)

                                                                                   ttGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTa         <  1:2570479/91‑1 (MQ=255)

                                                                                    tGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTAc        <  1:1758050/91‑1 (MQ=255)

                                                                                    tGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTAc        <  1:4104530/91‑1 (MQ=255)

                                                                                     gTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACt       <  1:3288663/91‑1 (MQ=255)

                                                                                     gTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACt       <  1:1660412/91‑1 (MQ=255)

                                                                                     gTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACt       <  1:3572774/91‑1 (MQ=255)

                                                                                      tACAGATATTTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTg      >  1:3475008/1‑91 (MQ=255)

                                                                                      tACAGATATTTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTg      >  1:1797401/1‑91 (MQ=255)

                                                                                      tACAGATATTTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTg      >  1:1797371/1‑91 (MQ=255)

                                                                                      tACAGATATTTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTg      >  1:3474111/1‑91 (MQ=255)

                                                                                       aCAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGa     >  1:1320461/1‑91 (MQ=255)

                                                                                       aCAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGa     >  1:1010766/1‑91 (MQ=255)

                                                                                        cAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTCTCCTGTACTGAt    <  1:3200122/91‑1 (MQ=255)

                                                                                        cAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGAt    <  1:1982768/91‑1 (MQ=255)

                                                                                          aaTATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGATTc  <  1:1859345/91‑1 (MQ=255)

                                                                                          aaTATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGATTc  <  1:1444017/91‑1 (MQ=255)

                                                                                          aaTATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGATTc  <  1:1382201/91‑1 (MQ=255)

TGTCATAAGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGATTC  >  USA300TCH1516_ALE/2530696‑2530876

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