BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,582,141	Δ1 bp	5.2%	coding (1169/1353 nt)	pnbA →	Para‑nitrobenzyl esterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,582,141	0	T	.	5.2%	708.5 / 26.9	210	coding (1169/1353 nt)	pnbA	Para‑nitrobenzyl esterase
Reads supporting (aligned to +/- strand): ref base T (96/102); new base . (2/9); total (98/112)
Fisher's exact test for biased strand distribution p-value = 6.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

Only 100 of 222 total aligned reads displayed.

TCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACA  >  USA300TCH1516_ALE/2582051‑2582231

 cAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATAtt                                                                                           <  1:2818908/91‑1 (MQ=255)

  aaTATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTc                                                                                          <  2:1819014/91‑1 (MQ=255)

   atatACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCt                                                                                         >  2:1070119/1‑91 (MQ=255)

   atatACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCt                                                                                         >  2:2054212/1‑91 (MQ=255)

   atatACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCt                                                                                         >  1:2778670/1‑91 (MQ=255)

   atatACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCt                                                                                         >  1:1932774/1‑91 (MQ=255)

   atatACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAATGCTTACCatatcta                                                                                         <  1:3367991/91‑4 (MQ=255)

    tataCAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCatatctag                                                                                        <  1:1980938/91‑5 (MQ=255)

     ataCGAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATA‑TCTAGa                                                                                      <  1:3296309/91‑1 (MQ=255)

     ataCAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATA‑TCTAGa                                                                                      <  1:3506810/91‑1 (MQ=255)

      taCAAGATTCAAACGCACAGCTATGGCTTGCTGAGGTTGCATGGCAGGATACTTCTAGCGCACACTATCGCAGTGCTTACCATA‑TCTAGaa                                                                                     <  2:2739183/91‑2 (MQ=255)

      taCAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATA‑TCTAGaa                                                                                     <  1:2555856/91‑2 (MQ=255)

        cAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGata                                                                                    <  1:494057/91‑1 (MQ=255)

        cAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGata                                                                                    >  1:3911398/1‑91 (MQ=255)

        cAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGata                                                                                    <  2:2802753/91‑1 (MQ=255)

         aaGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGatat                                                                                   <  1:3636794/91‑1 (MQ=255)

         aaGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGatat                                                                                   <  1:1189127/91‑1 (MQ=255)

         aaGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGatat                                                                                   <  2:2017265/91‑1 (MQ=255)

         aaGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCAGGATACTTCTAGCGCACACTATCGCAGTGCTTACCATA‑TCTAGaatat                                                                                  <  2:2404136/91‑5 (MQ=255)

          aGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATg                                                                                  >  1:1427882/1‑91 (MQ=255)

           gATTCAAACGCACAACTATGGCTTGTTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATgg                                                                                 >  1:1225812/1‑91 (MQ=255)

            aTTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGt                                                                                <  2:189931/91‑1 (MQ=255)

             ttCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTa                                                                               <  1:6088/91‑1 (MQ=255)

             ttCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTa                                                                               <  1:103442/91‑1 (MQ=255)

              tCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTAt                                                                              <  1:2082224/91‑1 (MQ=255)

              tCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTAt                                                                              <  1:3851937/91‑1 (MQ=255)

              tCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTAACATATTCTAGATATGGTAt                                                                              <  1:413769/91‑1 (MQ=255)

              gcAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTGGCGCACACTATCGCAGTGCTTACCATATTCTAGATGTGGTAt                                                                              <  2:1071522/90‑1 (MQ=255)

               cAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTAtt                                                                             >  1:3061932/1‑91 (MQ=255)

                aaaCGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTAttt                                                                            <  1:1760702/91‑1 (MQ=255)

                 aaCGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTAtttt                                                                           >  1:805375/1‑91 (MQ=255)

                  aCGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTg                                                                          >  1:581881/1‑91 (MQ=255)

                  aCGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTg                                                                          >  1:1169272/1‑91 (MQ=255)

                  aCGCACAACTATGGCTTGCTGAAGTTGCAGGGCATGATAATTCTAGCGAACCCTCTCGCAGTGCTAACCATACTCTAGATATGGTATTTTg                                                                          <  2:30021/91‑1 (MQ=255)

                   cGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTgg                                                                         >  2:3248991/1‑91 (MQ=255)

                      acaACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTc                                                                      <  1:1235948/91‑1 (MQ=255)

                       caACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCg                                                                     >  2:1114067/1‑91 (MQ=255)

                         aCTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGa                                                                   <  1:566523/91‑1 (MQ=255)

                         aCTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGa                                                                   >  2:541718/1‑91 (MQ=255)

                         aCTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTAGTTCGGa                                                                   >  2:3500193/1‑91 (MQ=255)

                          cTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGaa                                                                  <  1:2444095/91‑1 (MQ=255)

                           tATGGCTTGCTGAATTTGCATGGCATGATACTTCTTGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGaaa                                                                 <  1:3068904/91‑1 (MQ=255)

                               gCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTa                                                             <  2:583182/91‑1 (MQ=255)

                                cTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTAc                                                            >  2:163216/1‑91 (MQ=255)

                                 ttGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACa                                                           >  1:3795816/1‑91 (MQ=255)

                                 ttGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTGCGGAAACTTACa                                                           <  1:3034884/91‑1 (MQ=255)

                                   gCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACaaa                                                         >  1:2631768/1‑91 (MQ=255)

                                    cTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAat                                                        <  2:3403428/91‑1 (MQ=255)

                                      gAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAatat                                                      >  1:101309/1‑91 (MQ=255)

                                      gAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAatat                                                      <  2:3925016/91‑1 (MQ=255)

                                         tttGCATGGCATGATTCTTCTAGCGCACACTATCGCTGTGCTTACCATATTCTAGTTATGGTATTTTGGTTCGGAAACTTACAATTATTag                                                   <  2:3412444/91‑1 (MQ=255)

                                         tttGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTag                                                   <  2:3412265/91‑1 (MQ=255)

                                         tttGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTag                                                   >  1:696418/1‑91 (MQ=255)

                                         tttGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTag                                                   <  1:1001227/91‑1 (MQ=255)

                                         tttGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTag                                                   <  2:3723192/91‑1 (MQ=255)

                                          ttGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTagc                                                  >  2:3061795/1‑91 (MQ=255)

                                            gcatggcatgATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTagcag                                                >  1:901953/1‑91 (MQ=255)

                                             catggcatggTACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTagcagc                                               >  1:1500117/1‑91 (MQ=255)

                                                ggcatgATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGcaca                                            <  2:3320836/91‑1 (MQ=255)

                                                ggcatgATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGcaca                                            <  1:3273511/91‑1 (MQ=255)

                                                  catgATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATc                                          >  2:3124145/1‑91 (MQ=255)

                                                  catgATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATc                                          <  1:298369/91‑1 (MQ=255)

                                                   atgATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCa                                         <  2:228605/91‑1 (MQ=255)

                                                    tgATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAATCTTACAAATATTAGCAGCACATCaa                                        >  1:3681677/1‑91 (MQ=255)

                                                    tgATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCAGATCaa                                        >  2:3653096/1‑91 (MQ=255)

                                                     gATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAat                                       <  2:1021439/91‑1 (MQ=255)

                                                         cTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATcc                                   >  1:1892343/1‑91 (MQ=255)

                                                          ttCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCca                                  >  2:1492264/1‑91 (MQ=255)

                                                           tCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCcaa                                 <  1:24009/91‑1 (MQ=255)

                                                            cTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCcaac                                >  1:1123164/1‑91 (MQ=255)

                                                            cTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCcaac                                <  1:939832/91‑1 (MQ=255)

                                                             tAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCcaaca                               <  2:2882404/91‑1 (MQ=255)

                                                              aGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCcaacaa                              >  1:317548/1‑91 (MQ=255)

                                                              aGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCcaacaa                              <  2:1736173/91‑1 (MQ=255)

                                                                cgcgCACTATCGCAGTGCTTACCGTATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAGTATTAGCAGCGCATCAATATCCAACAACt                            <  2:844463/91‑1 (MQ=255)

                                                                 gcACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTg                           >  2:3622426/1‑91 (MQ=255)

                                                                 gcACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTg                           <  2:794864/91‑1 (MQ=255)

                                                                 gcACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTg                           >  1:3528033/1‑91 (MQ=255)

                                                                    cacTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTc                        >  2:1172570/1‑91 (MQ=255)

                                                                      cTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCAc                      <  2:1449923/91‑1 (MQ=255)

                                                                        aTCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCt                    >  2:3357684/1‑91 (MQ=255)

                                                                        aTCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCt                    >  2:2372553/1‑91 (MQ=255)

                                                                        aTCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCt                    <  1:2212304/91‑1 (MQ=255)

                                                                        aTCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCt                    >  2:2840040/1‑91 (MQ=255)

                                                                         tCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCagtga                   <  1:653379/91‑5 (MQ=255)

                                                                         tCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCtt                   <  1:452692/91‑1 (MQ=255)

                                                                          cGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTa                  <  2:2572667/91‑1 (MQ=255)

                                                                          cGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTa                  <  2:2692351/91‑1 (MQ=255)

                                                                          cGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTa                  >  1:3021348/1‑91 (MQ=255)

                                                                          cGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTa                  >  1:974071/1‑91 (MQ=255)

                                                                          cGCAGTGCTGACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTa                  <  1:3777781/91‑1 (MQ=255)

                                                                            cAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTaaa                <  2:3275987/91‑1 (MQ=255)

                                                                              gTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAAtt              >  1:289512/1‑91 (MQ=255)

                                                                              gTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAAtt              >  1:3864121/1‑91 (MQ=255)

                                                                                gCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAAtttt            >  1:2497059/1‑91 (MQ=255)

                                                                                   tACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTaa         >  1:3107430/1‑91 (MQ=255)

                                                                                    aCCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAg        >  1:1730167/1‑91 (MQ=255)

                                                                                      cATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTc      >  1:3894274/1‑91 (MQ=255)

                                                                                        tatTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGa    <  1:2077107/91‑1 (MQ=255)

                                                                                          ttCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACa  >  1:2226559/1‑91 (MQ=255)

TCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACA  >  USA300TCH1516_ALE/2582051‑2582231

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