BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,530,791	C→T	8.5%	intergenic (‑42/+231)	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,791	0	C	T	8.5%	321.0 / 16.0	106	intergenic (‑42/+231)	yefM/bdbD	Antitoxin YefM/Disulfide bond formation protein D
Reads supporting (aligned to +/- strand): ref base C (50/47); new base T (6/3); total (56/50)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-01

Only 100 of 110 total aligned reads displayed.

TAAGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGATTCCCTAC  >  USA300TCH1516_ALE/2530701‑2530881

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

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

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

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

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

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

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

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

   gTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAGATATTTGt                                                                                         <  1:1530251/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:1883214/1‑91 (MQ=255)

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

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

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

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

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

          tAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATaa                                                                                  >  1:3281982/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)

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

                   tCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGa                                                                         <  1:2653188/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:3259619/1‑91 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                           tcctTGTAACACACCTCAATTGTACAAATATCTGTACAAGAAAACAAGCGACAGAGAAAAAAACTACTAATCAAGATTTCGCACGGGgaag                                 <  1:3854805/91‑1 (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:1512874/1‑91 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                                gTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACt            >  1:1716839/1‑91 (MQ=255)

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

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

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

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

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

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

                                                                                  acagatatTTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGa          >  1:122705/5‑91 (MQ=255)

                                                                                  acagatatTTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGa          >  1:21514/5‑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:1382201/91‑1 (MQ=255)

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

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

                                                                                       tatCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGATTccc     <  1:3483997/91‑1 (MQ=255)

                                                                                          cTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGATTCCctac  >  1:1390140/1‑91 (MQ=255)

                                                                                          cTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGATTCCctac  >  1:854110/1‑91 (MQ=255)

TAAGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGATTCCCTAC  >  USA300TCH1516_ALE/2530701‑2530881

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