BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,449,468	A→T	9.5%	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	9.5%	535.3 / 30.6	184	S97T (TCT→ACT)	hutI	Imidazolonepropionase
Reads supporting (aligned to +/- strand): ref base A (75/90); new base T (2/16); total (78/106)
Fisher's exact test for biased strand distribution p-value = 4.97e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 186 total aligned reads displayed.

TTTTTAAACAAGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATCTACTAATGCTGGGGACACTAC  >  USA300TCH1516_ALE/2449380‑2449558

tttttAAACAAGTCATCTTCGGATGTTTCTCTGGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACt                                                                                          <  1:3379697/91‑1 (MQ=255)

tttttAAACAAGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGTCt                                                                                          <  1:766271/91‑1 (MQ=255)

tttttAAACAAGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCAAATGTTTCAAGGTAAGTCt                                                                                          <  1:3086469/91‑1 (MQ=255)

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

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

   ttAAACAAGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAAGATGCCACCACCCATGTCTAATATTTCAAGGTAAGACTTAc                                                                                       <  2:1830483/91‑1 (MQ=255)

    tAAACAAGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTAcc                                                                                      <  2:506612/91‑1 (MQ=255)

     aaaCGAGTCATCGTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCt                                                                                     <  1:3291801/91‑1 (MQ=255)

      aaCAAGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGTCTTACCtt                                                                                    <  1:1258615/91‑1 (MQ=255)

      aaCAAGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGTCTTACCtt                                                                                    <  2:3582792/91‑1 (MQ=255)

      aaCAAGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCtt                                                                                    >  1:3704268/1‑91 (MQ=255)

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

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

        cAAGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGa                                                                                  >  1:366028/1‑91 (MQ=255)

         aaGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGAc                                                                                 >  1:1509542/1‑91 (MQ=255)

          attCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACg                                                                                >  2:2152203/3‑91 (MQ=255)

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

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

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

           gTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTCATATTTCAAGGTAAGACTTACCTTGACGt                                                                               >  2:1621740/1‑91 (MQ=255)

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

            tCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGtt                                                                              >  2:877058/1‑91 (MQ=255)

            tCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGtt                                                                              >  2:3861660/1‑91 (MQ=255)

             cATCTTCGGATGTTGCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGttt                                                                             <  2:1279452/91‑1 (MQ=255)

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

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

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

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

                 ttCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAAAATTTCAAGGTAAGACTTACCTTGACGTTTTAAt                                                                         <  2:3693870/91‑1 (MQ=255)

                  tCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATg                                                                        <  2:1270756/91‑1 (MQ=255)

                  tCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATg                                                                        <  1:833348/91‑1 (MQ=255)

                   cGGATGTTTCTCTAGTTGCATTAGCGGTAGATAATATGCCAGCACCCCTTTCTAATATTTCAAGGTAGGACTTACCTTGACGTTTTAATGa                                                                       <  2:2521115/91‑1 (MQ=255)

                   cGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGa                                                                       >  2:1281056/1‑91 (MQ=255)

                    ggATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGAc                                                                      <  2:652918/91‑1 (MQ=255)

                     gATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACa                                                                     <  2:3461078/91‑1 (MQ=255)

                     gATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACa                                                                     <  2:1217737/91‑1 (MQ=255)

                      aTGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACAt                                                                    <  1:3260436/91‑1 (MQ=255)

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

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

                       tGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACAtc                                                                   <  1:249520/91‑1 (MQ=255)

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

                        gTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACAtct                                                                  >  2:1839309/1‑91 (MQ=255)

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

                          ttCTCTAGTAGCATTAACGGTAGATAGTATGCCACCACCCATTTCTGAGATTTCGAGGTAAGACTTACCTTGACGTTTTAATGCCGtctcc                                                                >  2:3388855/1‑90 (MQ=255)

                          ttCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCa                                                                <  1:593795/91‑1 (MQ=255)

                          ttCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCa                                                                <  1:2553824/91‑1 (MQ=255)

                          ttCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCa                                                                <  2:2493358/91‑1 (MQ=255)

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

                              ctAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGtt                                                            <  1:3433074/91‑1 (MQ=255)

                                agtagCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTtct                                                          <  2:2319598/91‑1 (MQ=255)

                                  tagCATTAACGGTGGATAATAGGCCACCACCCATTACTGATAGTACAAGGGAAGACTTACCTTGCCGTTTTAATGACATCTCAGGTTCTCg                                                        <  2:2631993/91‑1 (MQ=255)

                                    gCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCgag                                                      >  1:3861600/1‑91 (MQ=255)

                                       ttAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATc                                                   <  1:3749444/91‑1 (MQ=255)

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

                                          aCGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACCTACCTTGACGTTTTAATGACATCTCCTGTTCTCGAGCTcccc                                                >  2:727314/1‑89 (MQ=255)

                                            ggtaggtaATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTGGGCGTTTTAATGACATCGCATGTTCTCGAGATccacca                                              <  2:1862102/91‑1 (MQ=255)

                                             gTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCaa                                             <  1:2117390/91‑1 (MQ=255)

                                               ggATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAAt                                           <  2:1398064/90‑1 (MQ=255)

                                               aGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCaaaa                                           >  1:1468170/1‑90 (MQ=255)

                                                  taataTGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGtt                                        >  2:947120/1‑91 (MQ=255)

                                                  taataTGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGtt                                        <  2:2564893/91‑1 (MQ=255)

                                                   aataTGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTa                                       >  1:3767090/1‑91 (MQ=255)

                                                   aataTGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTa                                       <  1:2448208/91‑1 (MQ=255)

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

                                                      atGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAAt                                    >  1:869301/1‑91 (MQ=255)

                                                      atGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAAt                                    >  1:3355391/1‑91 (MQ=255)

                                                       tGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATg                                   >  1:286572/1‑91 (MQ=255)

                                                        gCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATgg                                  >  1:3311136/1‑91 (MQ=255)

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

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

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

                                                            ccaccCGTTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTAt                              <  1:262519/91‑1 (MQ=255)

                                                            ccaccCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTAt                              <  1:1163201/91‑1 (MQ=255)

                                                            ccaccCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTAt                              >  2:2826286/1‑91 (MQ=255)

                                                            ccaccCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTAt                              <  2:126100/91‑1 (MQ=255)

                                                            ccaccCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTAt                              >  1:1336886/1‑91 (MQ=255)

                                                              accCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTACTGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTAtgt                            <  1:3137753/91‑1 (MQ=255)

                                                              accCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTAtgt                            >  1:3094223/1‑91 (MQ=255)

                                                                 cATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCa                         <  2:3355019/91‑1 (MQ=255)

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

                                                                  aTTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCAt                        <  2:1093014/91‑1 (MQ=255)

                                                                  aTTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCAt                        <  1:2752256/91‑1 (MQ=255)

                                                                   tttCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATc                       >  2:2173479/1‑91 (MQ=255)

                                                                   tttCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATc                       >  2:378234/1‑91 (MQ=255)

                                                                      cTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATctac                    >  1:1288880/1‑91 (MQ=255)

                                                                          tatTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATCTACTAAt                <  1:2692207/91‑1 (MQ=255)

                                                                             ttCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATCTACTAATGCt             <  1:3822518/91‑1 (MQ=255)

                                                                              tCAAGGTACGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATCTCCTAATGCTg            >  2:39374/1‑91 (MQ=255)

                                                                                aaGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATCTACTAATGCTggg          <  1:2197461/91‑1 (MQ=255)

                                                                                aaGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATCTACTAATGCTggg          >  2:2503304/1‑91 (MQ=255)

                                                                                 aGGTAAGTCTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGCAAGTGCATCTACTAATGCTgggg         >  1:3593127/1‑91 (MQ=255)

                                                                                  ggTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATCTACTAATGCTGGGGa        <  2:1202908/91‑1 (MQ=255)

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

                                                                                    tAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGCTCCACCAAAAGTTAAATGGGTATGTGCATCTACTAATGCTGGGGaca      <  2:1036686/91‑1 (MQ=255)

                                                                                    tAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATCTACTAATGCTGGGGaca      <  1:2149756/91‑1 (MQ=255)

                                                                                     aaGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATCTACTAATGCTGGGGacac     >  2:3630689/1‑91 (MQ=255)

                                                                                      aGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAGATGTTAAATGGGTATGTGCATCTACTAATGCTGGGGACACt    >  2:3469016/1‑91 (MQ=255)

                                                                                      aGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGTGTATGTGCATCTACTAATGCTGGGGACACt    >  2:2888028/1‑91 (MQ=255)

                                                                                      aGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATCTACTAATGCTGGGGACACt    >  1:1930233/1‑91 (MQ=255)

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

TTTTTAAACAAGTCATCTTCGGATGTTTCTCTAGTAGCATTAACGGTAGATAATATGCCACCACCCATTTCTAATATTTCAAGGTAAGACTTACCTTGACGTTTTAATGACATCTCATGTTCTCGAGATCCACCAAATGTTAAATGGGTATGTGCATCTACTAATGCTGGGGACACTAC  >  USA300TCH1516_ALE/2449380‑2449558

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