BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,656,547	A→T	5.2%	intergenic (‑42/+36)	mhqA_3 ← / ← mhqR	Putative ring‑cleaving dioxygenase MhqA/HTH‑type transcriptional regulator MhqR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,656,547	0	A	T	5.2%	337.8 / 7.3	96	intergenic (‑42/+36)	mhqA_3/mhqR	Putative ring‑cleaving dioxygenase MhqA/HTH‑type transcriptional regulator MhqR
Reads supporting (aligned to +/- strand): ref base A (49/42); new base T (3/2); total (52/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 101 total aligned reads displayed.

TCATTGCAGTAACATGGTGAATACCTAGTAATTCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGTTTCTTAAACGCTTGTTGTAAGATT  >  USA300TCH1516_ALE/2656457‑2656634

tCATTGCAGTAACATGGTGAATACCTAGTAATTCATGCTTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGaaa                                                                                         >  1:1775058/1‑91 (MQ=255)

tCATTGCAGTAACATGGTGAATACCTAGTAATTCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGaaa                                                                                         >  1:2775483/1‑91 (MQ=255)

tCATTGCAGTAACATGGTGAATACCTAGTAATTCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGaaa                                                                                         >  1:2775567/1‑91 (MQ=255)

tCATTGCAGTAACATGGTGAATACCTAGTAATTCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGaaa                                                                                         >  1:2775169/1‑91 (MQ=255)

tCATTGCAGTAACATGGTGAATACCTAGTAATTCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGaaa                                                                                         >  1:3154490/1‑91 (MQ=255)

   ttGCAGTAACATGGTGAATACCTAGTAATTCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCg                                                                                      <  1:782065/91‑1 (MQ=255)

    tGCAGTAACATGGTGAATACCTAGTAATTCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAGATCgg                                                                                     >  1:1650214/1‑90 (MQ=255)

    tGCAGTAACATGGTGAATACCTAGTAATTCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCga                                                                                     >  1:984783/1‑91 (MQ=255)

    tGCAGTAACATGGTGAATACCTAGTAATTCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCga                                                                                     >  1:225511/1‑91 (MQ=255)

    tGCAGTAACATGGTGAATACCTAGTAATTCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCga                                                                                     <  1:577045/91‑1 (MQ=255)

    tGCAGTAACATGGTGAATACCTAGTAATTCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCga                                                                                     <  1:271103/91‑1 (MQ=255)

      cAGTAACATGGTGAATACCTAGTAATTCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGATTTCgaga                                                                                   <  1:1553185/91‑1 (MQ=255)

         tAACATGGTGAATACCTAGTAATTCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGATTTCGAGATaa                                                                                <  1:1861496/91‑2 (MQ=255)

         tAACATAGTGAATACCTAATAGTTCATGATTATTTGCCATTTTATATACCCTCTTTGT‑CTTTTAGTGTATTTATCTCGAAATCgagaaatt                                                                               >  1:3221074/1‑87 (MQ=255)

           aCATGGTGAATACCTAGTAATTCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACtt                                                                              >  1:3747511/1‑91 (MQ=255)

           aCATGGTGAATACCTAGTAATTCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACtt                                                                              >  1:1977427/1‑91 (MQ=255)

            cATGGTGAATACCTAGTAATTCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACttt                                                                             <  1:4119705/91‑1 (MQ=255)

              tGGTGAATACCTAGTAATTCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTaa                                                                           <  1:1173848/91‑1 (MQ=255)

               ggTGAATACCTAGTAATTCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAAt                                                                          <  1:2215818/91‑1 (MQ=255)

               ggTGAATACCTAGTAATTCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAAt                                                                          <  1:632796/91‑1 (MQ=255)

                gTGAATACCTAGTAATTCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATg                                                                         <  1:3040788/91‑1 (MQ=255)

                gTGAATACCTAGTAATTCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATg                                                                         <  1:3039759/91‑1 (MQ=255)

                 tGAATACCTAGTAATTCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGt                                                                        <  1:3249174/91‑1 (MQ=255)

                  gAATACCTAGTAATTCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTa                                                                       <  1:252394/91‑1 (MQ=255)

                  gAATACCTAGTAATTCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTa                                                                       >  1:1033101/1‑91 (MQ=255)

                   aaTACCTAGTAATTCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTaa                                                                      <  1:3945516/91‑1 (MQ=255)

                   aaTACCTAGTAATTCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTaa                                                                      <  1:3945488/91‑1 (MQ=255)

                    aTACCTAGTAATTCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTaaa                                                                     >  1:2542726/1‑91 (MQ=255)

                     tACCTAGTAATTCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAAt                                                                    >  1:3695504/1‑91 (MQ=255)

                         tAGTAATTCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAAtcttt                                                                >  1:3213753/1‑91 (MQ=255)

                          aGTAATTCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTa                                                               >  1:3960957/1‑91 (MQ=255)

                           gTAATTCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAg                                                              <  1:1514808/91‑1 (MQ=255)

                               ttCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCa                                                          >  1:2652046/1‑91 (MQ=255)

                                 cATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACg                                                        >  1:1473217/1‑91 (MQ=255)

                                 cATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACg                                                        >  1:1473216/1‑91 (MQ=255)

                                 cATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACg                                                        >  1:1310509/1‑91 (MQ=255)

                                 cATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACg                                                        >  1:1270992/1‑91 (MQ=255)

                                   tGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGcc                                                      >  1:3257328/1‑91 (MQ=255)

                                   tGATTATTTATCATTTTATATACCCTCTTTATATTTTCAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGcc                                                      >  1:3498490/1‑91 (MQ=255)

                                    gATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCt                                                     >  1:2017094/1‑91 (MQ=255)

                                     attattTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCtt                                                    <  1:4083330/91‑1 (MQ=255)

                                      ttatttatCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTa                                                   <  1:2524/91‑1 (MQ=255)

                                       tatttatCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTac                                                  >  1:3778425/1‑91 (MQ=255)

                                       tatttatCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTac                                                  >  1:1454812/1‑91 (MQ=255)

                                       tatttatCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTac                                                  <  1:3552949/91‑1 (MQ=255)

                                        atttatCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTaca                                                 >  1:2639257/1‑91 (MQ=255)

                                          ttatCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACt                                               >  1:3218079/1‑91 (MQ=255)

                                             tCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCt                                            >  1:2866679/1‑91 (MQ=255)

                                             tCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCt                                            >  1:1408630/1‑91 (MQ=255)

                                             tCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCt                                            >  1:803849/1‑91 (MQ=255)

                                             tCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCt                                            <  1:1606519/91‑1 (MQ=255)

                                             tCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCt                                            >  1:711650/1‑91 (MQ=255)

                                             tCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCt                                            >  1:1970831/1‑91 (MQ=255)

                                             tCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCt                                            >  1:4056513/1‑91 (MQ=255)

                                             tCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCt                                            >  1:4017040/1‑91 (MQ=255)

                                              cATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTg                                           <  1:1748936/91‑1 (MQ=255)

                                                ttttATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTa                                         >  1:3532757/1‑91 (MQ=255)

                                                ttttATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTa                                         >  1:3532786/1‑91 (MQ=255)

                                                  ttATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGa                                       <  1:2751927/91‑1 (MQ=255)

                                                    atataCCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGAtt                                     >  1:1088656/1‑91 (MQ=255)

                                                      ataCCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATtgt                                   <  1:4003227/91‑1 (MQ=255)

                                                      ataCCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATtgt                                   >  1:2765050/1‑91 (MQ=255)

                                                      ataCCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATtgt                                   <  1:4003395/91‑1 (MQ=255)

                                                        aCCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGc                                 <  1:3515541/91‑1 (MQ=255)

                                                        aCCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGc                                 <  1:3795414/91‑1 (MQ=255)

                                                        aCCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGc                                 <  1:3795377/91‑1 (MQ=255)

                                                        aCCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGc                                 <  1:3756629/91‑1 (MQ=255)

                                                         cccTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCa                                >  1:2287663/1‑91 (MQ=255)

                                                              tttATATTTTTAGTGCAGTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACTCTACTGTAGATTGTGCACTTAg                           <  1:1249921/91‑1 (MQ=255)

                                                               ttATATTTTTAGTGCATTTATCTCGAAATCGTGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGt                          <  1:3697521/91‑1 (MQ=255)

                                                               ttATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGt                          <  1:3697759/91‑1 (MQ=255)

                                                                 atatTTTTAGTGCATTTATCTCGAAATCGAGATACTTTCATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGttt                        >  1:2359500/1‑91 (MQ=255)

                                                                 atatTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGttt                        >  1:2703824/1‑91 (MQ=255)

                                                                 atatTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGttt                        >  1:2360594/1‑91 (MQ=255)

                                                                  tatTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGTTTc                       >  1:2098772/1‑91 (MQ=255)

                                                                   atTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGTTTCt                      <  1:1907487/91‑1 (MQ=255)

                                                                   atTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGTTTCt                      >  1:2785728/1‑91 (MQ=255)

                                                                    tttttAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGTTTCtt                     >  1:789496/1‑91 (MQ=255)

                                                                    tttttAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGTTTCtt                     >  1:889079/1‑91 (MQ=255)

                                                                    tttttAGTGCATATATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGTTTCtt                     >  1:3205594/1‑91 (MQ=255)

                                                                     ttttAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGTTTCTTa                    >  1:797115/1‑91 (MQ=255)

                                                                      tttAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGTTTCTTaa                   <  1:3549168/91‑1 (MQ=255)

                                                                       ttAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGTTTCTTaaa                  >  1:367294/1‑91 (MQ=255)

                                                                          gTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGTTTCTTAAACGc               <  1:2431988/91‑1 (MQ=255)

                                                                           tGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGTTTCTTAAACGCt              <  1:2548979/91‑1 (MQ=255)

                                                                             cATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGTTTCTTAAACGCttg            <  1:492993/91‑1 (MQ=255)

                                                                             aaagtATCTCGATTTCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGTTTCTTAAACGCttg            >  1:3963992/5‑91 (MQ=255)

                                                                              aTTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGTTTCTTAAACGCttgt           <  1:2017662/91‑1 (MQ=255)

                                                                              aTTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGTTTCTTAAACGCttgt           <  1:1702482/91‑1 (MQ=255)

                                                                              aTTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGTTTCTTAAACGCttgt           <  1:1632796/91‑1 (MQ=255)

                                                                              aTTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGTTTCTTAAACGCttgt           <  1:155969/91‑1 (MQ=255)

                                                                               agtATCTCGATTTCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGTTTCTTAAACGCttgtt          >  1:2226721/3‑91 (MQ=255)

                                                                                 tATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGTTTCTTAAACGCTtgttgt        <  1:1931109/91‑1 (MQ=255)

                                                                                 tATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGTTTCTTAAACGCTtgttgt        <  1:3617043/91‑1 (MQ=255)

                                                                                  aTCTCGATTTCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGTTTCTTAAACGCTTGTTGTa       >  1:434849/1‑91 (MQ=255)

                                                                                   tctcGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGTTTCTTAAACGCTTGTTGTaa      >  1:3298631/1‑91 (MQ=255)

                                                                                    ctcGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGTTTCTTAAACGCTTGTTGTAAg     >  1:2464432/1‑91 (MQ=255)

                                                                                    ctcGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGTTTCTTAAACGCTTGTTGTAAg     >  1:2932610/1‑91 (MQ=255)

                                                                                      cGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGTTTCTTAAACGCTTGTTGTAAGAt   >  1:4006984/1‑91 (MQ=255)

                                                                                       gAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGTTTCTTAAACGCTTGTTGTAAGAtt  <  1:440046/91‑1 (MQ=255)

TCATTGCAGTAACATGGTGAATACCTAGTAATTCATGATTATTTATCATTTTATATACCCTCTTTATATTTTTAGTGCATTTATCTCGAAATCGAGATACTTTAATGTAAATTTTTAGTGCACGCCTTACACTTCTGTAGATTGTGCACTTAGTTTCTTAAACGCTTGTTGTAAGATT  >  USA300TCH1516_ALE/2656457‑2656634

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