BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,255,815	A→G	5.5%	intergenic (+31/+213)	rplS → / ← USA300HOU_RS06250	50S ribosomal protein L19/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,255,815	0	A	G	5.5%	354.4 / 8.3	109	intergenic (+31/+213)	rplS/USA300HOU_RS06250	50S ribosomal protein L19/hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (49/53); new base G (4/2); total (54/55)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

Only 100 of 113 total aligned reads displayed.

AAATTATATTACTTACGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGTTTAGTTAATGTAAAGTTGATA  >  USA300TCH1516_ALE/1255725‑1255905

aaaTTATATTACTTACGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCa                                                                                            <  1:3624865/91‑1 (MQ=255)

 aaTTATATTACTTACGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCaa                                                                                           <  1:3617420/91‑1 (MQ=255)

  aTTATATTACTTACGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCaaa                                                                                          >  1:2809926/1‑91 (MQ=255)

    tataTTACTTACGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAAtt                                                                                        <  1:1939670/91‑1 (MQ=255)

    tataTTACTTACGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAAtt                                                                                        >  1:1163482/1‑91 (MQ=255)

    tataTTACTTACGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAAtt                                                                                        <  1:1967907/91‑1 (MQ=255)

     atatTACTTACGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTc                                                                                       <  1:3069078/91‑1 (MQ=255)

     atatTACTTACGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTc                                                                                       <  1:3069037/91‑1 (MQ=255)

     atatTACTTACGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTc                                                                                       >  1:2916188/1‑91 (MQ=255)

       atTACTTACGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTTTAAGTCGAAtttga                                                                                     <  1:3196464/91‑4 (MQ=255)

       atTACTTACGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGa                                                                                     >  1:3554071/1‑91 (MQ=255)

       atTACTTACGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGa                                                                                     >  1:226551/1‑91 (MQ=255)

        ttacttacGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGAc                                                                                    >  1:1857456/1‑91 (MQ=255)

          acttacGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTc                                                                                  >  1:2800318/1‑91 (MQ=255)

          acttacGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTc                                                                                  >  1:710498/1‑91 (MQ=255)

          acttacGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTc                                                                                  <  1:3958149/91‑1 (MQ=255)

           cttacGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCa                                                                                 <  1:809906/91‑1 (MQ=255)

            ttacGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCaa                                                                                <  1:1868080/91‑1 (MQ=255)

             tacGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCaaa                                                                               >  1:1338312/1‑91 (MQ=255)

             tacGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCaaa                                                                               <  1:2187344/91‑1 (MQ=255)

              acGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAAc                                                                              <  1:2667228/91‑1 (MQ=255)

              acGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAAc                                                                              <  1:1760529/91‑1 (MQ=255)

               cGTAGTTTACGTGGTAAAGCTGCTAGAATCCGAGAAATTCGTTAGTCAGCATTTAAATAAAGCTATGTTTGAGTCGAATTTGACTCATACa                                                                             <  1:3309231/91‑1 (MQ=25)

                gTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACaa                                                                            <  1:1683192/91‑1 (MQ=255)

                gTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACaa                                                                            <  1:805117/91‑1 (MQ=255)

                  aGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATa                                                                          <  1:3261912/91‑1 (MQ=255)

                    tttACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATaaa                                                                        <  1:1046051/91‑1 (MQ=255)

                    tttACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATaaa                                                                        >  1:3070992/1‑91 (MQ=255)

                     ttACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAg                                                                       >  1:755196/1‑91 (MQ=255)

                      tACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGc                                                                      <  1:1122602/91‑1 (MQ=255)

                       aCGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGcc                                                                     >  1:213330/1‑91 (MQ=255)

                        cGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCa                                                                    >  1:1791476/1‑91 (MQ=255)

                           ggTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCt                                                                 >  1:893026/1‑91 (MQ=255)

                            gTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTa                                                                >  1:148406/1‑91 (MQ=255)

                              aaaGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAgg                                                              >  1:2492679/1‑91 (MQ=255)

                               aaGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGt                                                             <  1:1814707/91‑1 (MQ=255)

                               aaGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGt                                                             <  1:1814494/91‑1 (MQ=255)

                                aGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTc                                                            <  1:1537981/91‑1 (MQ=255)

                                aGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTc                                                            <  1:2610288/91‑1 (MQ=255)

                                 gctgctAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCa                                                           >  1:706460/1‑91 (MQ=255)

                                  ctgctAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCAc                                                          >  1:3123549/1‑91 (MQ=255)

                                  ctgctAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCAc                                                          >  1:441501/1‑91 (MQ=255)

                                   tgctAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACt                                                         <  1:2822857/91‑1 (MQ=255)

                                    gctAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACtt                                                        <  1:3929862/91‑1 (MQ=255)

                                       aGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTa                                                     <  1:3696788/91‑1 (MQ=255)

                                       aGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTa                                                     <  1:3033312/91‑1 (MQ=255)

                                         aaTCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATg                                                   <  1:2466633/91‑1 (MQ=255)

                                          aTCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATgg                                                  >  1:3431646/1‑91 (MQ=255)

                                           tCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCCAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGa                                                 >  1:4135381/1‑91 (MQ=255)

                                            ccAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAg                                                <  1:2880269/91‑1 (MQ=255)

                                             cAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGt                                               >  1:1242999/1‑91 (MQ=255)

                                              aaGTAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTg                                              <  1:3345580/91‑1 (MQ=255)

                                              aaGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTg                                              <  1:3008383/91‑1 (MQ=255)

                                                gAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGat                                            >  1:1371552/1‑91 (MQ=255)

                                                gAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGat                                            >  1:3604719/1‑91 (MQ=255)

                                                gAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATCAGTCAAATCCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGat                                            >  1:3307295/1‑91 (MQ=255)

                                                 aaaTTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGata                                           <  1:867542/91‑1 (MQ=255)

                                                     tCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGa                                       >  1:2021076/1‑91 (MQ=255)

                                                     tCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGa                                       >  1:2020613/1‑91 (MQ=255)

                                                      cGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGAt                                      >  1:2631372/1‑91 (MQ=255)

                                                       gTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATg                                     >  1:2320800/1‑91 (MQ=255)

                                                       gTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATg                                     >  1:401168/1‑91 (MQ=255)

                                                       gTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATg                                     <  1:1250918/91‑1 (MQ=255)

                                                       gTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATg                                     <  1:2540998/91‑1 (MQ=255)

                                                        ttAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATgg                                    >  1:3477165/1‑91 (MQ=255)

                                                        ttAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATgg                                    <  1:3570111/91‑1 (MQ=255)

                                                        ttAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATgg                                    <  1:302770/91‑1 (MQ=255)

                                                        ttAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATgg                                    >  1:2565992/1‑91 (MQ=255)

                                                        ttAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATgg                                    <  1:2327877/91‑1 (MQ=255)

                                                         tAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGc                                   <  1:2549496/91‑1 (MQ=255)

                                                         tAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGc                                   >  1:1822703/1‑91 (MQ=255)

                                                         tAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGc                                   <  1:108199/91‑1 (MQ=255)

                                                              aGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCttttt                              >  1:4125074/1‑91 (MQ=255)

                                                              aGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCttttt                              <  1:2296503/91‑1 (MQ=255)

                                                              aGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCttttt                              <  1:1395328/91‑1 (MQ=255)

                                                               gCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCtttttt                             >  1:1044196/1‑91 (MQ=255)

                                                                cATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGttttaatt                            >  1:2075882/1‑87 (MQ=255)

                                                                cATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCttttttt                            >  1:3129492/1‑91 (MQ=255)

                                                                  tttAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGc                          <  1:2849458/91‑1 (MQ=255)

                                                                   ttAAGCAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCt                         >  1:3277181/1‑91 (MQ=255)

                                                                     aaaCAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGt                       <  1:3116886/91‑1 (MQ=255)

                                                                        cAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGTTTa                    <  1:578970/91‑1 (MQ=255)

                                                                         aaaGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGTTTAg                   >  1:3080657/1‑91 (MQ=255)

                                                                          aaGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGTTTAGt                  >  1:386429/1‑91 (MQ=255)

                                                                           aGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGTTTAGtt                 >  1:3584632/1‑91 (MQ=255)

                                                                           aGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGTTTAGtt                 <  1:1678460/91‑1 (MQ=255)

                                                                             cTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGTTTAGTTaa               <  1:1197344/91‑1 (MQ=255)

                                                                               atgtatgAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGTTGTTTAGTTAATg             >  1:167155/1‑91 (MQ=255)

                                                                               atgtatgAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGTTTAGTTAATg             >  1:1012887/1‑91 (MQ=255)

                                                                               atgtatgAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGTCTAGTTAATg             >  1:3965943/1‑91 (MQ=255)

                                                                                tgtatgAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGTTTAGTTAATGt            <  1:1099682/91‑1 (MQ=255)

                                                                                tgtatgAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGTTTAGTTAATGt            <  1:1051853/91‑1 (MQ=255)

                                                                                  tttgAGTCGAATTTGACTCATACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGTTTAGTTAATGTaa          >  1:5589/3‑91 (MQ=38)

                                                                                   ttgAGTCGAATTTGACTCATACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGTTTAGTTAATGTaaa         >  1:2056763/2‑91 (MQ=255)

                                                                                    tgAGTCGAATTTGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGTTTAGTTAATGTAAAg        >  1:2239680/1‑91 (MQ=255)

                                                                                     gAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGTTTAGTTAATGTAAAGt       <  1:3396707/91‑1 (MQ=255)

                                                                                      aGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGTTTAGTTAATGTAAAGtt      >  1:521958/1‑91 (MQ=255)

                                                                                       gTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGTTTAGTTAATGTAAAGTTg     <  1:3498152/91‑1 (MQ=255)

                                                                                         cAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGTTTAGTTAATGTAAAGTTGAt   <  1:3059567/91‑1 (MQ=255)

                                                                                          aaaTTCGACTCAAACAATAAAGCCATCTAGGGCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGTTTAGTTAATGTAAAGTTGATa  >  1:638164/1‑91 (MQ=255)

AAATTATATTACTTACGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGTTTAGTTAATGTAAAGTTGATA  >  USA300TCH1516_ALE/1255725‑1255905

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