BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,449,468	A→T	10.3%	S97T (TCT→ACT)	hutI ←	Imidazolonepropionase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,449,468	0	A	T	10.3%	554.3 / 43.8	186	S97T (TCT→ACT)	hutI	Imidazolonepropionase
Reads supporting (aligned to +/- strand): ref base A (72/93); new base T (3/16); total (76/110)
Fisher's exact test for biased strand distribution p-value = 2.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 188 total aligned reads displayed.

TTTTTAAACAAGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATCTACTAATGCTGGGGACACTAC  >  USA300TCH1516_ALE/2449380‑2449558

tttttAAACAAGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACt                                                                                          >  1:565701/1‑91 (MQ=255)

tttttAAACAAGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACt                                                                                          <  1:1782547/91‑1 (MQ=255)

 ttttAAACAAGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGTCtt                                                                                         <  1:2046844/91‑1 (MQ=255)

 ttttAAACAAGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACtt                                                                                         >  2:877501/1‑91 (MQ=255)

 ttttAAACAAGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACtt                                                                                         <  1:368232/91‑1 (MQ=255)

  tttAAACAAGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGTCTTa                                                                                        <  1:2168702/91‑1 (MQ=255)

  tttAAACAAGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGTCTTa                                                                                        <  1:2841829/91‑1 (MQ=255)

  tttAAACAAGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGTCTTa                                                                                        <  2:1106399/91‑1 (MQ=255)

  tttAAACAAGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGTCTTa                                                                                        <  2:392443/91‑1 (MQ=255)

   ttAAACAAGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGTCTTAc                                                                                       <  2:2417888/91‑1 (MQ=255)

   ttAAACAAGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGTCTTAc                                                                                       <  1:430253/91‑1 (MQ=255)

       aCAAGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGTCTTACCTTg                                                                                   <  1:2493933/91‑1 (MQ=255)

       aCAAGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTg                                                                                   <  2:1355517/91‑1 (MQ=255)

       aCAAGTCAGCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGTCTTACCTTg                                                                                   <  2:3593214/91‑1 (MQ=255)

        cAAGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGa                                                                                  <  2:3026547/91‑1 (MQ=255)

        cAAGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGa                                                                                  <  1:1536085/91‑1 (MQ=255)

        cAAGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGa                                                                                  <  2:3646808/91‑1 (MQ=255)

        cAAGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGa                                                                                  <  2:3880147/91‑1 (MQ=255)

         aaGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGTCTTACCTTGaa                                                                                 <  1:281755/91‑2 (MQ=255)

          aGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACg                                                                                <  1:427825/91‑1 (MQ=255)

          aGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACg                                                                                >  1:644562/1‑91 (MQ=255)

           gTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGt                                                                               <  1:627536/91‑1 (MQ=255)

           gTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGt                                                                               >  1:963579/1‑91 (MQ=255)

            tCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGtt                                                                              <  1:877501/91‑1 (MQ=255)

             cATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGttt                                                                             >  2:2380478/1‑91 (MQ=255)

              aTCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGttgt                                                                            <  2:1129460/91‑3 (MQ=255)

               tCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTa                                                                           <  2:3211897/91‑1 (MQ=255)

               tCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTa                                                                           >  2:3738917/1‑91 (MQ=255)

               tCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATACTTCAAGGTAAGACTTACCTTGACGTTTTa                                                                           >  2:1548706/1‑91 (MQ=255)

                cTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTaa                                                                          >  1:1518899/1‑91 (MQ=255)

                cTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTaa                                                                          <  1:610837/91‑1 (MQ=255)

                 ttCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAAt                                                                         <  1:3057221/91‑1 (MQ=255)

                   cGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGa                                                                       <  1:1736646/91‑1 (MQ=255)

                     gATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACa                                                                     >  1:1500434/1‑91 (MQ=255)

                       tGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACAtc                                                                   >  2:148966/1‑91 (MQ=255)

                        gTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACAtct                                                                  <  2:2659951/91‑1 (MQ=255)

                         tttCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACAtctc                                                                 <  1:1391512/91‑1 (MQ=255)

                          ttCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCa                                                                >  2:1777170/1‑91 (MQ=255)

                           tctcTAGTAGCATTAACGGTAGATAATATGCCAGCACCCATTTCTAATATTTCAAGGTATGACTTACCTTGACGTTTTAATGACATCTCAt                                                               <  1:3443063/91‑1 (MQ=255)

                            ctctAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTGTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATg                                                              <  1:82289/91‑1 (MQ=255)

                            ctctAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATg                                                              <  2:66341/91‑1 (MQ=255)

                             tctAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGCCATCCCATGt                                                             >  1:3642924/1‑91 (MQ=255)

                             tctAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGt                                                             >  1:2529014/1‑91 (MQ=255)

                             tctAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGt                                                             <  1:3742992/91‑1 (MQ=255)

                              ctAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGtt                                                            >  2:2898569/1‑91 (MQ=255)

                               tagtagCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTtc                                                           <  1:3777364/91‑1 (MQ=255)

                               tagtagCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTtc                                                           >  2:3179987/1‑91 (MQ=255)

                                 gtagCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTtctc                                                         <  1:1039427/91‑1 (MQ=255)

                                 gtagCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTtctc                                                         <  1:1039636/91‑1 (MQ=255)

                                       ttAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATc                                                   >  1:1611765/1‑91 (MQ=255)

                                        tAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATcc                                                  <  2:2746025/91‑1 (MQ=255)

                                         aaCGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATcca                                                 <  2:1971671/91‑1 (MQ=255)

                                           cGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATccacc                                               <  1:2631636/91‑1 (MQ=255)

                                           cGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATccacc                                               >  1:2036924/1‑91 (MQ=255)

                                           cGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATccacc                                               >  1:2375660/1‑91 (MQ=255)

                                           cGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATccacc                                               >  1:1802844/1‑91 (MQ=255)

                                            ggTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATcaacca                                              >  1:3229870/1‑91 (MQ=255)

                                             gtacataataTGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGAGCCACCaa                                             >  2:3445191/5‑91 (MQ=255)

                                               aGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAAt                                           <  1:1809503/91‑1 (MQ=255)

                                               aGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAAt                                           <  2:3740941/91‑1 (MQ=255)

                                               aGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAAt                                           >  2:1449315/1‑91 (MQ=255)

                                                gATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATg                                          <  2:3385464/91‑1 (MQ=255)

                                                   aataTGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATg                                                                        >  2:1595292/1‑58 (MQ=255)

                                                     tatGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTaaa                                     <  2:328818/91‑1 (MQ=255)

                                                     tatGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTaaa                                     <  1:2248457/91‑1 (MQ=255)

                                                     tatGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTaaa                                     <  2:1886387/91‑1 (MQ=255)

                                                         ccaccaCCCATTTCTAATATTTCAAGGTAAGACTTATCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATggg                                 <  2:3056133/91‑1 (MQ=255)

                                                         ccaccaCCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATggg                                 <  1:2961759/91‑1 (MQ=255)

                                                         ccaccaCCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATggg                                 >  2:3558353/1‑91 (MQ=255)

                                                         ccaccaCCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATggg                                 <  1:3275055/91‑1 (MQ=255)

                                                          caccacCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGt                                <  2:3243296/91‑1 (MQ=255)

                                                           accaccCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTa                               >  1:880783/1‑91 (MQ=255)

                                                             caccCATTTCTAATATTTCAAGGTAAGACTTTCCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTAtg                             <  2:1327475/91‑1 (MQ=255)

                                                             caccCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTAtg                             >  2:785742/1‑91 (MQ=255)

                                                             caccCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTAtg                             <  1:781846/91‑1 (MQ=255)

                                                                ccATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCCTGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGc                          <  1:3835360/91‑1 (MQ=255)

                                                                ccATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGc                          >  1:3581173/1‑91 (MQ=255)

                                                                  aTTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCAt                        >  2:1114042/1‑91 (MQ=255)

                                                                  aTTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCAt                        >  1:1927254/1‑91 (MQ=255)

                                                                    ttCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTaa                                                                          <  1:389026/39‑1 (MQ=255)

                                                                    ttCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCACATGTTAAATGGGTATGTGCATct                      >  2:2971146/1‑91 (MQ=255)

                                                                       tAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATctact                   <  2:2462264/91‑1 (MQ=255)

                                                                       tAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATctact                   <  1:2314738/91‑1 (MQ=255)

                                                                       tAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATctact                   <  1:1675533/91‑1 (MQ=255)

                                                                       tAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATctact                   <  2:2237782/91‑1 (MQ=255)

                                                                       tAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATctact                   <  1:3158550/91‑1 (MQ=255)

                                                                         atatTTCGAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATCTACTaa                 <  1:1935802/91‑1 (MQ=255)

                                                                               cAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATCTACTAATGCTgg           >  2:1930707/1‑91 (MQ=255)

                                                                                 aGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATCTACTAATGCTgggg         >  2:1216453/1‑91 (MQ=255)

                                                                                  ggTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATCTACTAATGCTGGGGa        >  1:801237/1‑91 (MQ=255)

                                                                                  ggTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCAACTACTAATGCTGGTGa        >  2:518905/1‑91 (MQ=255)

                                                                                   gTAAGTCTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATCTACTAATGCTGGGGac       >  1:1690172/1‑91 (MQ=255)

                                                                                   gTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATCTACTAATGCTGGGGac       <  2:2036924/91‑1 (MQ=255)

                                                                                   gTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATCTACTAATGCTGGGGac       <  2:2168249/91‑1 (MQ=255)

                                                                                   gTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATCTACTAATGCTGGGGac       <  2:1008016/91‑1 (MQ=255)

                                                                                    tAAGTCTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATCTACTAATGCTGGGGaca      >  1:1145868/1‑91 (MQ=255)

                                                                                      aGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATCTACTAATGCTGGGGACACt    <  2:1517413/91‑1 (MQ=255)

                                                                                      aGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATCTACTAATGCTGGGGACACt    <  2:617461/91‑1 (MQ=255)

                                                                                       gACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATCTACTAATGCTGGGGACACTa   <  1:3633621/91‑1 (MQ=255)

                                                                                        aCTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATCTACTAATGCTGGGGACACTAc  <  2:2652519/91‑1 (MQ=255)

TTTTTAAACAAGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATCTACTAATGCTGGGGACACTAC  >  USA300TCH1516_ALE/2449380‑2449558

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