BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,530,786	A→G	8.6%	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	8.6%	570.5 / 31.8	194	intergenic (‑37/+236)	yefM/bdbD	Antitoxin YefM/Disulfide bond formation protein D
Reads supporting (aligned to +/- strand): ref base A (77/100); new base G (7/10); total (84/110)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 207 total aligned reads displayed.

TGTCATAAGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGATTC  >  USA300TCH1516_ALE/2530696‑2530876

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

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

 gTCATAAGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACaaa                                                                                           <  1:1029479/91‑1 (MQ=255)

 gTCATAAGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACaaa                                                                                           >  1:3464120/1‑91 (MQ=255)

 gTCATAAGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACaaa                                                                                           <  2:689608/91‑1 (MQ=255)

    aTAAGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAatat                                                                                        <  2:424540/91‑1 (MQ=255)

    aTAAGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAatat                                                                                        >  1:1892697/1‑91 (MQ=255)

     tAAGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATc                                                                                       >  2:3721404/1‑91 (MQ=255)

      aaGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAGatattt                                                                                      <  1:2922540/91‑3 (MQ=255)

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

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

         tGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTa                                                                                   >  1:1217251/1‑91 (MQ=255)

           ccTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTACCACACCTCAATTGTACAACTATCTGTACa                                                                                 >  1:3106891/1‑91 (MQ=255)

           ccTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACa                                                                                 <  2:1474034/91‑1 (MQ=255)

             tttAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTActat                                                                               >  2:3499880/1‑88 (MQ=255)

              ttAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCGTTTAACACACCTCAATTGTACAAATATCTGTACAATa                                                                              <  1:3271889/91‑1 (MQ=255)

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

               tAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATaa                                                                             <  1:2703386/91‑1 (MQ=255)

               tAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATaa                                                                             <  2:2225789/91‑1 (MQ=255)

               tAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACATCAATTGTACAAATATCTGTACAATaa                                                                             <  2:3221258/91‑1 (MQ=255)

                 aaTTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATaaaa                                                                           >  2:281827/1‑91 (MQ=255)

                  aTTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAAc                                                                          <  1:1791287/91‑1 (MQ=255)

                   ttCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACa                                                                         >  2:2590856/1‑91 (MQ=255)

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

                      tGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGc                                                                      <  1:728517/91‑1 (MQ=255)

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

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

                        tCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGa                                                                    >  2:2403636/1‑91 (MQ=255)

                            gCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGAcaca                                                                <  2:1947945/91‑1 (MQ=255)

                            gCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGAcaca                                                                <  2:3879381/91‑1 (MQ=255)

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

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

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

                               tATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGaa                                                             >  2:1693501/1‑91 (MQ=255)

                                 tGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGaaaa                                                           >  2:301582/1‑91 (MQ=255)

                                   aataatTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGaaaaaa                                                         >  1:1672953/1‑91 (MQ=255)

                                      aatTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAact                                                      >  2:1386918/1‑91 (MQ=255)

                                        tttttAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACa                                                                                 >  2:3025735/1‑62 (MQ=255)

                                        tttttAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAactac                                                    <  2:2909639/91‑1 (MQ=255)

                                         ttttAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACTCAGa                                                              <  1:3056996/80‑1 (MQ=255)

                                         ttttAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGa                                                              >  2:3056996/1‑80 (MQ=255)

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

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

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

                                           ttAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTaa                                                 <  1:1519022/91‑1 (MQ=255)

                                           ttAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTaa                                                 <  1:2761508/91‑1 (MQ=255)

                                            taatatTCATATTAATCGCCTCCTTTTAACACACCTCTATTGTACAAATATCTGTACCATAAAACAAGCGACACAGAATAAAACTCCTaaa                                                >  2:573453/1‑90 (MQ=255)

                                              ataatCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCa                                              >  2:724325/1‑91 (MQ=255)

                                              ataatCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCa                                              <  2:1858281/91‑1 (MQ=255)

                                              ataatCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCa                                              >  2:2083854/1‑91 (MQ=255)

                                              ataatCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCa                                              <  1:2796270/91‑1 (MQ=255)

                                                aatCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAg                                            >  1:1122345/1‑91 (MQ=255)

                                                 atCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGa                                           <  2:2201105/91‑1 (MQ=255)

                                                  tCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGAt                                          <  1:714843/91‑1 (MQ=255)

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

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

                                                       tttATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGc                                     >  2:935258/1‑91 (MQ=255)

                                                         tATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCAc                                   >  1:2558207/1‑91 (MQ=255)

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

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

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

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

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

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

                                                             gCCTCCTTTTAACACACCTCACTCGTACAAATATCTGTACAATAAAACACGCGCCACAGCCACAAACTACTAATCAAGATCTCGCACggcg                               >  2:3538707/1‑89 (MQ=255)

                                                             gCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACgggg                               >  1:853389/1‑91 (MQ=255)

                                                             gCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACgggg                               >  2:449831/1‑91 (MQ=255)

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

                                                             gCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACgggg                               >  2:3296639/1‑91 (MQ=255)

                                                              cctcctTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGga                              >  1:3650265/1‑91 (MQ=255)

                                                              cctcctTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGga                              >  2:1391335/1‑91 (MQ=255)

                                                               ctcctTTTAACACACTTCAATTGTAGTAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGgaa                             <  2:3367308/91‑1 (MQ=255)

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

                                                                    tttAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATc                        <  1:2679591/91‑1 (MQ=255)

                                                                    tttAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATc                        >  2:2753056/1‑91 (MQ=255)

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

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

                                                                     ttAACACACCTCAATTGCACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCa                       <  1:2313680/91‑1 (MQ=255)

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

                                                                          acacCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATAcc                  <  1:3739475/91‑1 (MQ=255)

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

                                                                           cacCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCt                 <  1:2862236/91‑1 (MQ=255)

                                                                            acCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTa                >  2:1167949/1‑91 (MQ=255)

                                                                            aacTCAATTGAACAAATATCTGTACAATAAAACAAGCGACACAAAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTa                <  2:1547574/89‑1 (MQ=255)

                                                                              cTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATc              <  1:3238945/91‑1 (MQ=255)

                                                                               tCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGAGTTCGCACGGGGAAGAATCAATACCTATcc             <  2:3385053/91‑1 (MQ=255)

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

                                                                                 aaTTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTg           <  2:2129912/91‑1 (MQ=255)

                                                                                  aTTGTACAGATATTTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGt          >  2:596711/1‑91 (MQ=255)

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

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

                                                                                   ttGTACAGATATTTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTa         >  1:1700512/1‑91 (MQ=255)

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

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

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

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

                                                                                      tACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTg      <  1:2054701/91‑1 (MQ=255)

                                                                                       aCAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGa     <  2:2695198/91‑1 (MQ=255)

                                                                                          aaTATCTGTACAATAAAACAAGCGACACAGAAAAAAATTACTAATCCAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGATTc  >  1:3092560/1‑91 (MQ=255)

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

TGTCATAAGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGATTC  >  USA300TCH1516_ALE/2530696‑2530876

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