BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,255,820	C→T	5.4%	intergenic (+36/+208)	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,820	0	C	T	5.4%	739.2 / 18.7	222	intergenic (+36/+208)	rplS/USA300HOU_RS06250	50S ribosomal protein L19/hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (112/97); new base T (10/2); total (122/100)
Fisher's exact test for biased strand distribution p-value = 7.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

Only 100 of 228 total aligned reads displayed.

ATATTACTTACGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGTTTAGTTAATGTAAAGTTGATAATGTA  >  USA300TCH1516_ALE/1255730‑1255910

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

 tatTACTTACGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCg                                                                                           >  2:2732464/1‑91 (MQ=255)

 tatTACTTACGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCg                                                                                           <  1:2152859/91‑1 (MQ=255)

 tatTACTTACGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCg                                                                                           <  2:505937/91‑1 (MQ=255)

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

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

   ttacttacGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGAc                                                                                         <  2:2712531/91‑1 (MQ=255)

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

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

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

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

      cttacGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCa                                                                                      >  2:3059541/1‑91 (MQ=255)

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

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

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

      ctaacGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCa                                                                                      <  2:1327365/88‑1 (MQ=255)

      cgtacGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCa                                                                                      <  2:3346883/89‑1 (MQ=255)

       ttacGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCaa                                                                                     >  2:2329928/1‑91 (MQ=255)

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

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

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

             aGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATa                                                                               >  2:1814726/1‑91 (MQ=255)

             aGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCATTAATCAGCATTTGAACAAAGCTATGTATGAGTCAAATTCGACTCAATCAATa                                                                               >  1:3679609/1‑91 (MQ=255)

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

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

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

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

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

                   cGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCa                                                                         <  1:1047988/91‑1 (MQ=255)

                    gTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCAt                                                                        <  2:188379/91‑1 (MQ=255)

                     tGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATc                                                                       >  1:2943128/1‑91 (MQ=255)

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

                      ggTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCt                                                                      <  1:396642/91‑1 (MQ=255)

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

                       gTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTa                                                                     <  2:1951112/91‑1 (MQ=255)

                       gTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTa                                                                     <  2:1477248/91‑1 (MQ=255)

                       gTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTa                                                                     <  2:32214/91‑1 (MQ=255)

                       gTAAAGCTGCTAGAATCAAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTa                                                                     >  2:1822910/1‑91 (MQ=255)

                        tAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAg                                                                    <  2:2462182/91‑1 (MQ=255)

                        tAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAg                                                                    >  1:3797639/1‑91 (MQ=255)

                         aaaGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAgg                                                                   <  1:199930/91‑1 (MQ=255)

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

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

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

                              tgctAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACt                                                              >  1:1186578/1‑91 (MQ=255)

                               gctAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACtt                                                             >  2:2173884/1‑91 (MQ=255)

                                ctAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACttt                                                            >  2:2749608/1‑91 (MQ=255)

                                      tCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGa                                                      >  2:2604331/1‑91 (MQ=255)

                                      tCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGa                                                      <  2:2542341/91‑1 (MQ=255)

                                      tCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGa                                                      <  2:3899103/91‑1 (MQ=255)

                                       ccAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAg                                                     >  1:2850437/1‑91 (MQ=255)

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

                                          aGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGa                                                  <  1:219846/91‑1 (MQ=255)

                                          aGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGa                                                  <  1:918885/91‑1 (MQ=255)

                                          aGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGa                                                  >  1:2637526/1‑91 (MQ=255)

                                             aaTTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGatat                                               >  2:2912621/1‑91 (MQ=255)

                                              aTTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGAtata                                              <  2:1865114/91‑1 (MQ=255)

                                                 cGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGATTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGAt                                           <  2:3544729/91‑1 (MQ=255)

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

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

                                                  gTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATACATg                                          >  1:3188741/1‑91 (MQ=255)

                                                  gTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTCATATAGATg                                          >  2:682995/1‑91 (MQ=255)

                                                  gTTAATCAGCATATAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATg                                          >  1:3043798/1‑91 (MQ=255)

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

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

                                                     aaTCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCt                                       <  2:422863/91‑1 (MQ=255)

                                                     aaTCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCt                                       >  2:904241/1‑91 (MQ=255)

                                                        cAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCtttt                                    >  1:3092803/1‑91 (MQ=255)

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

                                                          gCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCtttttt                                  <  2:3913004/91‑1 (MQ=255)

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

                                                           cATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCttttttt                                 <  2:3301123/91‑1 (MQ=255)

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

                                                            aactAAACAAAGATATGTATGAGTCTAATTGAAAACAAAAAAAAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTg                                <  2:613582/88‑1 (MQ=255)

                                                            aTTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTg                                >  1:2558109/1‑91 (MQ=255)

                                                              ttAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCt                              <  1:2605519/91‑1 (MQ=255)

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

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

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

                                                                  aCAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGttt                          >  1:591624/1‑91 (MQ=255)

                                                                  aCAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGttt                          >  2:1131800/1‑91 (MQ=255)

                                                                   cAAAGCTATCTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTTCTGTTTa                         >  1:1967923/1‑91 (MQ=255)

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

                                                                    aaaGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGATTTTTTTGCTGTTTAg                        >  2:3773681/1‑91 (MQ=255)

                                                                       gCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGTTTAGTTa                     >  2:575706/1‑91 (MQ=255)

                                                                        cTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGTTTAGTTaa                    >  1:3610389/1‑91 (MQ=255)

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

                                                                             tatgAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGTTTAGTTAATGTaa               >  2:824491/1‑91 (MQ=255)

                                                                               tgAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGTTTAGTTAATGTAAAg             >  1:1314744/1‑91 (MQ=255)

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

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

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

                                                                                     gaaTTTGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGTGTGATATAGATGGCTTTTTTTGCTGTTTAGCTAATGTACAGTTGATa       >  2:2227984/2‑91 (MQ=255)

                                                                                      aaTTTGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGTTTAGTTAATGTAAAGTTGATaa      >  1:584295/1‑91 (MQ=255)

                                                                                      aaTTTGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGTTTAGTTAATGTAAAGTTGATaa      >  2:2621213/1‑91 (MQ=255)

                                                                                      aaTTTGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGTTTAGCTAATGTAAAGCTGATaa      >  2:660711/1‑91 (MQ=255)

                                                                                      aaTTTGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTATTTTTGCTGTTTAGTTAATGTAAAGTTGATaa      >  1:3939355/1‑91 (MQ=255)

                                                                                      aaTTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGTTTAGTTAATGTAAAGTTGATaa      <  2:1015493/91‑1 (MQ=255)

                                                                                        ttCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTTCTGTTTTTTTAATGTAAAGTTGATAATg    >  2:1056690/1‑91 (MQ=255)

                                                                                          cGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGTTTAGTTAATGTAAAGTTGATAATGTa  <  2:2273878/91‑1 (MQ=255)

ATATTACTTACGTAGTTTACGTGGTAAAGCTGCTAGAATCCAAGAAATTCGTTAATCAGCATTTAAACAAAGCTATGTATGAGTCAAATTCGACTCAAACAATAAAGCCATCTAGGTCACTTTTATGGAGTGATATAGATGGCTTTTTTTGCTGTTTAGTTAATGTAAAGTTGATAATGTA  >  USA300TCH1516_ALE/1255730‑1255910

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