BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,904,152	C→A	11.4%	A182S (GCC→TCC)	atl_2 ←	Bifunctional autolysin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,904,152	0	C	A	11.4%	430.3 / 39.8	131	A182S (GCC→TCC)	atl_2	Bifunctional autolysin
Reads supporting (aligned to +/- strand): ref base C (77/39); new base A (7/8); total (84/47)
Fisher's exact test for biased strand distribution p-value = 1.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 135 total aligned reads displayed.

CAGCACCGCCTGAAATAGCTTTTTCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACTTTTAACTGCGCCAGTTTCTAATAATGC  >  USA300TCH1516_ALE/1904062‑1904242

cAGCACCGCCTGAAATAGCTTTTTCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGc                                                                                            <  1:3984495/91‑1 (MQ=255)

cAGCACCGCCTGAAATAGCTTTTTCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGc                                                                                            <  1:1283580/91‑1 (MQ=255)

cAGCACCGCCTGAAATAGCTTTTTCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGc                                                                                            >  1:1331211/1‑91 (MQ=255)

 aGCACCGCCTGAAATAGCTTTTTCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGACt                                                                                           >  1:2007824/1‑91 (MQ=255)

 aGCACCGCCTGAAATAGCTTTTTCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGACt                                                                                           >  1:2008468/1‑91 (MQ=255)

  gCACCGCCTGAAATAGCTTTTTCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTc                                                                                          >  1:2994061/1‑91 (MQ=255)

   cACCGCCTGGAATAGCTTTTTCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTcc                                                                                         <  1:3963391/91‑1 (MQ=255)

    aCCGCCTGAAATAGCTTTTTCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCt                                                                                        >  1:3567212/1‑91 (MQ=255)

    aCCGCCTGAAATAGCTTTTTCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCt                                                                                        >  1:1877526/1‑91 (MQ=255)

     ccGCCTGAAATAGCTTTTTCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTa                                                                                       >  1:2392346/1‑91 (MQ=255)

     ccGCCTGAAATAGCTTTTTCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTa                                                                                       <  1:1593630/91‑1 (MQ=255)

     ccGCCTGAAATAGCTTTTTCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTa                                                                                       >  1:4098026/1‑91 (MQ=255)

     ccGCCTGAAATAGCTTTTTCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTa                                                                                       >  1:1707620/1‑91 (MQ=255)

     ccGCCTGAAATAGCTTTTTCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTa                                                                                       >  1:1707573/1‑91 (MQ=255)

     ccGCCTGAAATAGCTTTTTCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGATCCTa                                                                                       <  1:87769/91‑1 (MQ=255)

     ccGCCTGAAATAGCTTTTTCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGATCCTa                                                                                       <  1:87784/91‑1 (MQ=255)

      cGCCTGAAATAGCTTTTTCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTAc                                                                                      >  1:894799/1‑91 (MQ=255)

       gCCTGAAATAGCTTTTTCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACt                                                                                     <  1:1760888/91‑1 (MQ=255)

          tGAAATAGCTTTTTCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCa                                                                                  <  1:2841372/91‑1 (MQ=255)

          tGAAATAGCTTTTTCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGATCCTACTCCa                                                                                  <  1:2783979/91‑1 (MQ=255)

           gAAATAGCTTTTTCAGGTGTATCCCAACCATGCTTTTTAGCATGTTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCAt                                                                                 <  1:4136431/91‑1 (MQ=255)

           gAAATAGCTTTTTCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCAt                                                                                 >  1:2638526/1‑91 (MQ=255)

           gAAATAGCTTTTTCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCAt                                                                                 >  1:2238706/1‑91 (MQ=255)

           gAAATAGCTTTTTCAGGTGTATCCCAACCATGCTCTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCCCCTACTCCAt                                                                                 >  1:1148694/1‑91 (MQ=255)

            aaaTAGCTTTTTCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATa                                                                                <  1:482453/91‑1 (MQ=255)

              aTAGCTTTTTCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGATCCTACTCCATAGa                                                                              >  1:3572239/1‑91 (MQ=255)

                 gCTTTTTCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAAt                                                                           >  1:484287/1‑91 (MQ=255)

                 gCTTTTTCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAAt                                                                           >  1:1363156/1‑91 (MQ=255)

                    ttttCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTgt                                                                        >  1:835867/1‑91 (MQ=255)

                     tttCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTgta                                                                       <  1:222284/91‑1 (MQ=255)

                     tttCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTgta                                                                       >  1:290598/1‑91 (MQ=255)

                       tCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTgtagt                                                                     <  1:2098985/91‑1 (MQ=255)

                       tCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTgtagt                                                                     >  1:1164906/1‑91 (MQ=255)

                       tCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTgtagt                                                                     >  1:1164935/1‑91 (MQ=255)

                       tCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTgtagt                                                                     <  1:3582774/91‑1 (MQ=255)

                       tCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTgtagt                                                                     >  1:3079609/1‑91 (MQ=255)

                       tCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGATCCTACTCCATAGAAATTgtagt                                                                     >  1:3012023/1‑91 (MQ=255)

                        cAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTgtagta                                                                    <  1:83995/91‑1 (MQ=255)

                         aGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTAc                                                                   >  1:1779040/1‑91 (MQ=255)

                            tgtATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACttt                                                                >  1:1822890/1‑91 (MQ=255)

                            tgtATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACttt                                                                >  1:1822953/1‑91 (MQ=255)

                            tgtATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGATCCTACTCCATAGAAATTGTAGTACttt                                                                >  1:2849248/1‑91 (MQ=255)

                            tgtATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGATCCTACTCCATAGAAATTGTAGTACttt                                                                >  1:3098526/1‑91 (MQ=255)

                            tgtATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGCGTCTTTATCAAGGGCTCCTACTCCATAGACATTGTAGTAAttt                                                                >  1:4141481/1‑91 (MQ=255)

                             gtATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACtttt                                                               >  1:816861/1‑91 (MQ=255)

                                tCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTCGTACTTTTTGc                                                            >  1:3759411/1‑91 (MQ=255)

                                tCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGc                                                            <  1:277396/91‑1 (MQ=255)

                                tCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGc                                                            <  1:595852/91‑1 (MQ=255)

                                   cAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCAt                                                         <  1:2948430/91‑1 (MQ=255)

                                    aaCCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATc                                                        >  1:2903971/1‑91 (MQ=255)

                                    aaCCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATc                                                        >  1:2273384/1‑91 (MQ=255)

                                        aTGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAAtt                                                    >  1:3652120/1‑91 (MQ=255)

                                            tttttAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCgc                                                >  1:685402/1‑90 (MQ=255)

                                            tttttAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGa                                                >  1:1504407/1‑91 (MQ=255)

                                             ttttAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGAc                                               <  1:1988020/91‑1 (MQ=255)

                                                    atatTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTa                                        >  1:2270562/1‑91 (MQ=255)

                                                    atatTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTa                                        >  1:2270521/1‑91 (MQ=255)

                                                     tatTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAg                                       <  1:573637/91‑1 (MQ=255)

                                                     tatTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAg                                       >  1:2047637/1‑91 (MQ=255)

                                                     tatTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAg                                       >  1:2047758/1‑91 (MQ=255)

                                                      atTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGc                                      <  1:3223147/91‑1 (MQ=255)

                                                      atTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGc                                      <  1:3376827/91‑1 (MQ=255)

                                                      atTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGc                                      <  1:3223337/91‑1 (MQ=255)

                                                      atTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGATCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGc                                      >  1:2686270/1‑91 (MQ=255)

                                                           gCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAAtt                                 >  1:2105405/1‑91 (MQ=255)

                                                            cACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTc                                >  1:2226224/1‑91 (MQ=255)

                                                            cACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTc                                >  1:1459690/1‑91 (MQ=255)

                                                              ccTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCAc                              <  1:3947288/91‑1 (MQ=255)

                                                                 gTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACttt                           <  1:987240/91‑1 (MQ=255)

                                                                 gTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACttt                           <  1:767410/91‑1 (MQ=255)

                                                                  ttttAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACtttt                          <  1:3409767/91‑1 (MQ=255)

                                                                     tAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACTTTTAAc                       <  1:1995649/91‑1 (MQ=255)

                                                                      aaTTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACTTTTAACt                      >  1:1686414/1‑91 (MQ=255)

                                                                      aaTTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACTTTTAACt                      <  1:3227630/91‑1 (MQ=255)

                                                                       aTTCGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACTTTTTACTg                     <  1:3268938/91‑1 (MQ=255)

                                                                        ttGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACTTTTAACTgc                    >  1:3108414/1‑91 (MQ=255)

                                                                        ttGGGTCTTTATCAAGGGATCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACTTTTAACTgc                    <  1:449232/91‑1 (MQ=255)

                                                                        ttGGGTCTTTATCAAGGGATCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACTTTTAACTgc                    <  1:448482/91‑1 (MQ=255)

                                                                          gggTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACTTTTAACTgcgc                  <  1:2684858/91‑1 (MQ=255)

                                                                          gggTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACTTTTAACTgcgc                  <  1:2176356/91‑1 (MQ=255)

                                                                          gggTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACTTTTAACTgcgc                  >  1:1466703/1‑91 (MQ=255)

                                                                           ggTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACTTTTAACTGCGcc                 <  1:3521355/91‑1 (MQ=255)

                                                                             tCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACTTTTAACTGCGCCAg               <  1:1624470/91‑1 (MQ=255)

                                                                              cTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACTTTTAACTGCGCCAGt              >  1:1102802/1‑91 (MQ=255)

                                                                              cTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACTTTTAACTGCGCCAGt              >  1:1102766/1‑91 (MQ=255)

                                                                               tttATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACTTTTAACTGCGCCAGtt             <  1:3876282/91‑1 (MQ=255)

                                                                               tttATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACTTTTAACTGCGCCAGtt             >  1:811011/1‑91 (MQ=255)

                                                                                ttATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACTTTTAACTGCGCCAGttt            <  1:1679029/91‑1 (MQ=255)

                                                                                ttATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACTTTTAACTGCGCCAGttt            >  1:610692/1‑91 (MQ=255)

                                                                                 tATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACTTTTAACTGCGCCAGTTTc           >  1:3530485/1‑91 (MQ=255)

                                                                                 tATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACTTTTAACTGCGCCAGTTTc           >  1:755140/1‑91 (MQ=255)

                                                                                 tATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACTTTTAACTGCGCCAGTTTc           >  1:3530527/1‑91 (MQ=255)

                                                                                   tCAAGGGATCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACTTTTAACTGCGCCAGTTTCta         >  1:475942/1‑91 (MQ=255)

                                                                                     aaGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACTTTTAACTGCGCCAGTTTCtaat       <  1:331432/91‑1 (MQ=255)

                                                                                     aaGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACTTTTAACTGCGCCAGTTTCtaat       >  1:3481781/1‑91 (MQ=255)

                                                                                      aGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACTTTTAACTGCGCCAGTTTCtaata      >  1:2314058/1‑91 (MQ=255)

                                                                                      aGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACTTTTAACTGCGCCAGTTTCtaata      >  1:2315423/1‑91 (MQ=255)

                                                                                       gggCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACTTTTAACTGCGCCAGTTTCtaataa     >  1:24181/1‑91 (MQ=255)

                                                                                       gggCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACTTTTAACTGCGCCAGTTTCtaataa     >  1:1215012/1‑91 (MQ=255)

                                                                                         gCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACTTTTAACTGCGCCAGTTTCTAATAATg   <  1:242639/91‑1 (MQ=255)

CAGCACCGCCTGAAATAGCTTTTTCAGGTGTATCCCAACCATGCTTTTTAGCATATTCTGCACCTGTTTTAATTGGGTCTTTATCAAGGGCTCCTACTCCATAGAAATTGTAGTACTTTTTGCCATCAATTTCGACTCCATTAGCTAATTCACTTTTAACTGCGCCAGTTTCTAATAATGC  >  USA300TCH1516_ALE/1904062‑1904242

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