BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	665,392	T→G	32.2%	L79* (TTA→TGA)	paiA_1 →	Spermidine/spermine N(1)‑acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	665,392	0	T	G	32.2%	150.9 / 127.2	112	L79* (TTA→TGA)	paiA_1	Spermidine/spermine N(1)‑acetyltransferase
Reads supporting (aligned to +/- strand): ref base T (38/38); new base G (20/16); total (58/54)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01

Only 100 of 116 total aligned reads displayed.

AAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGA  >  USA300TCH1516_ALE/665303‑665482

aaagaaaagaCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATta                                                                                           <  1:3348535/91‑1 (MQ=255)

aaagaaaagaCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATta                                                                                           <  1:445819/91‑1 (MQ=255)

   gaaaagaCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTatat                                                                                        >  1:3824194/1‑91 (MQ=255)

    aaaagaCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATAtt                                                                                       <  1:591796/91‑1 (MQ=255)

     aaagaCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATAttt                                                                                      >  1:2146165/1‑91 (MQ=255)

      aagaCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAGCTTTATTTACGGTGAAGAATTATATTTa                                                                                     <  1:3276369/91‑1 (MQ=255)

      aagaCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTa                                                                                     >  1:1334017/1‑91 (MQ=255)

      aagaCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTa                                                                                     >  1:1333993/1‑91 (MQ=255)

       agaCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTAt                                                                                    <  1:2440277/91‑1 (MQ=255)

        gaCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATc                                                                                   <  1:1141916/91‑1 (MQ=255)

        gaCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATc                                                                                   <  1:2284557/91‑1 (MQ=255)

         aCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCa                                                                                  >  1:1832877/1‑91 (MQ=255)

         aCTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCa                                                                                  <  1:1070807/91‑1 (MQ=255)

          cTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAg                                                                                 >  1:1852655/1‑91 (MQ=255)

           ttCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGc                                                                                >  1:2064369/1‑91 (MQ=255)

           ttCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGc                                                                                <  1:2967815/91‑1 (MQ=255)

            tCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCTTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCt                                                                               <  1:4136656/91‑1 (MQ=255)

            tCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCt                                                                               >  1:2158999/1‑91 (MQ=255)

            tCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCt                                                                               >  1:499896/1‑91 (MQ=255)

            tCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCt                                                                               >  1:3105524/1‑91 (MQ=255)

             cAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTc                                                                              >  1:3589203/1‑91 (MQ=255)

             cAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTc                                                                              >  1:3589530/1‑91 (MQ=255)

             cAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTCTATTTACGGTGAAGTATGATATTTATCAGCTc                                                                              >  1:3556974/1‑91 (MQ=255)

             cAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATAGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTc                                                                              >  1:3590322/1‑91 (MQ=255)

             cAAGAACAATTATTCATAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTGTATTTACGGTGAAGAATGATATTTATCAGCTc                                                                              >  1:3530850/1‑91 (MQ=255)

                    aaTTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGt                                                                       <  1:17810/91‑1 (MQ=255)

                    aaTTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGt                                                                       >  1:3747209/1‑91 (MQ=255)

                    aaTTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGt                                                                       <  1:942832/91‑1 (MQ=255)

                    aaTTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGt                                                                       <  1:2017697/91‑1 (MQ=255)

                     attattCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGtt                                                                      <  1:2411409/91‑1 (MQ=255)

                     attattCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGtt                                                                      >  1:3338302/1‑91 (MQ=255)

                     attattCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGtt                                                                      <  1:1017134/91‑1 (MQ=255)

                     attattCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGtt                                                                      <  1:1017327/91‑1 (MQ=255)

                      ttattCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTa                                                                     <  1:666786/91‑1 (MQ=255)

                      ttattCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTa                                                                     >  1:445702/1‑91 (MQ=255)

                      ttattCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTa                                                                     <  1:2289971/91‑1 (MQ=255)

                        attCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTaaa                                                                   >  1:3099019/1‑91 (MQ=255)

                          tCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAAcc                                                                 <  1:376190/91‑1 (MQ=255)

                           cTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCa                                                                >  1:2611317/1‑91 (MQ=255)

                            ttAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAg                                                               <  1:2708855/91‑1 (MQ=255)

                            ttAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAg                                                               <  1:1009684/91‑1 (MQ=255)

                             tAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACCTTATTTACGGTGAAGCATTATATGTATCAGCTCATTATGTTAAACCAGa                                                              <  1:1940330/91‑1 (MQ=255)

                              aGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGaa                                                             >  1:2225166/1‑91 (MQ=255)

                               gTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAAt                                                            <  1:1117071/91‑1 (MQ=255)

                               gTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAAt                                                            <  1:1117244/91‑1 (MQ=255)

                               gTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAAt                                                            <  1:3364145/91‑1 (MQ=255)

                               gTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAAt                                                            >  1:3980148/1‑91 (MQ=255)

                                tCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATc                                                           <  1:2108060/91‑1 (MQ=255)

                                tCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATc                                                           <  1:2107690/91‑1 (MQ=255)

                                  gTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGc                                                         >  1:510110/1‑91 (MQ=255)

                                    tGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCaa                                                       >  1:1923979/1‑91 (MQ=255)

                                     gaagaaGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAAc                                                      >  1:3599045/1‑91 (MQ=255)

                                     gaagaaGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAAc                                                      >  1:1323521/1‑91 (MQ=255)

                                       agaagGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTACACCAGAATCGCAACAt                                                    >  1:2754268/1‑91 (MQ=255)

                                       agaagGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACAt                                                    >  1:2466448/1‑91 (MQ=255)

                                       agaagGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACAt                                                    <  1:3638736/91‑1 (MQ=255)

                                         aagGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATAc                                                  >  1:7995/1‑91 (MQ=255)

                                           ggTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAg                                                >  1:3630151/1‑91 (MQ=255)

                                            gTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAgg                                               >  1:3225651/1‑91 (MQ=255)

                                            gTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAgg                                               >  1:1164849/1‑91 (MQ=255)

                                             tAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGt                                              <  1:2129038/91‑1 (MQ=255)

                                             tAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGt                                              <  1:2129531/91‑1 (MQ=255)

                                              aaTGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGtt                                             <  1:915264/91‑1 (MQ=255)

                                              aaTGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGtt                                             <  1:1469147/91‑1 (MQ=255)

                                              aaTGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGtt                                             <  1:1469001/91‑1 (MQ=255)

                                               aTGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTa                                            <  1:2957111/91‑1 (MQ=255)

                                               aTGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTa                                            <  1:2957337/91‑1 (MQ=255)

                                                 gACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATg                                          >  1:2286906/1‑91 (MQ=255)

                                                  aCATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATgg                                         <  1:145830/91‑1 (MQ=255)

                                                  aCATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATgg                                         <  1:2624299/91‑1 (MQ=255)

                                                  aCATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATgg                                         <  1:2623455/91‑1 (MQ=255)

                                                  aCATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATgg                                         <  1:3027683/91‑1 (MQ=255)

                                                    aTCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTa                                       >  1:3952411/1‑91 (MQ=255)

                                                     tCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTAc                                      <  1:1413681/91‑1 (MQ=255)

                                                      cGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACa                                     >  1:1543916/1‑91 (MQ=255)

                                                       gTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAg                                    >  1:2067814/1‑91 (MQ=255)

                                                       gTTGGCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAg                                    <  1:2300582/91‑1 (MQ=255)

                                                        ttGGCTTTGCTAACTTTATTTACGGTGAATAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGc                                   <  1:2895931/91‑1 (MQ=255)

                                                        ttGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGc                                   <  1:2466332/91‑1 (MQ=255)

                                                          ggCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCAt                                 >  1:2811497/1‑91 (MQ=255)

                                                          ggCTTTGCTAACTTTATTTACGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCAt                                 <  1:287072/91‑1 (MQ=255)

                                                                      tttatttaCGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAg                     >  1:2396687/1‑91 (MQ=255)

                                                                        gatctaCGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGa                   <  1:184425/87‑1 (MQ=255)

                                                                        gatctaCGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGa                   <  1:182093/87‑1 (MQ=255)

                                                                         atttaCGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGAt                  >  1:3825176/1‑91 (MQ=255)

                                                                         atttaCGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGAt                  >  1:3825493/1‑91 (MQ=255)

                                                                          tttaCGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGAtt                 >  1:1826715/1‑91 (MQ=255)

                                                                            taCGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTAt               >  1:1301155/1‑91 (MQ=255)

                                                                            taCGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTAt               >  1:1301119/1‑91 (MQ=255)

                                                                             aCGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATc              <  1:2026144/91‑1 (MQ=255)

                                                                             aCGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATc              <  1:2749542/91‑1 (MQ=255)

                                                                             aCGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATc              <  1:1944124/91‑1 (MQ=255)

                                                                              cGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCa             <  1:3509822/91‑1 (MQ=255)

                                                                              cGGTGAAGAATGATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCa             >  1:612393/1‑91 (MQ=255)

                                                                               ggTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCAc            >  1:1448509/1‑91 (MQ=255)

                                                                                   aagaaTTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGttt        <  1:387286/91‑1 (MQ=255)

                                                                                   aagaaTTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGttt        >  1:3820604/1‑91 (MQ=255)

                                                                                     gaaTGATATTTATCAGCTCATTATGCTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTg      <  1:3854780/91‑1 (MQ=255)

                                                                                        ttATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAg   <  1:2378017/91‑1 (MQ=255)

                                                                                         tataTTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGa  >  1:2260431/1‑91 (MQ=255)

AAAGAAAAGACTTCAAGAACAATTATTCTTAGTCGTTGAAGAAGGTAATGACATCGTTGGCTTTGCTAACTTTATTTACGGTGAAGAATTATATTTATCAGCTCATTATGTTAAACCAGAATCGCAACATACAGGTTATGGTACAGCATTGTTAAATGAAGGATTATCACGTTTTGAAGA  >  USA300TCH1516_ALE/665303‑665482

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