BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,516,450	T→A	6.4%	*115Y (TAA→TAT)	gpsB ←	Cell cycle protein GpsB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,516,450	0	T	A	6.4%	384.0 / 13.2	112	*115Y (TAA→TAT)	gpsB	Cell cycle protein GpsB
Reads supporting (aligned to +/- strand): ref base T (70/33); new base A (3/4); total (73/39)
Fisher's exact test for biased strand distribution p-value = 2.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 122 total aligned reads displayed.

TATAATATAGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTTGTTGTTGTATTATTCGA  >  USA300TCH1516_ALE/1516363‑1516539

tataatataGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAattat                                                                                        >  1:1382246/1‑91 (MQ=255)

tataatataGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAattat                                                                                        >  1:2565933/1‑91 (MQ=255)

 ataatataGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAattatt                                                                                       <  1:642629/91‑1 (MQ=255)

  taatataGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTAttt                                                                                      <  1:146813/91‑1 (MQ=255)

    atataGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTAc                                                                                    >  1:76251/1‑91 (MQ=255)

     tataGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTAcc                                                                                   >  1:2471918/1‑91 (MQ=255)

      ataGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTGGTAAATAATTActa                                                                                  <  1:4023304/91‑3 (MQ=255)

      ataGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCa                                                                                  <  1:2391383/91‑1 (MQ=255)

       taGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCaa                                                                                 >  1:1480709/1‑91 (MQ=255)

          cGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATa                                                                              >  1:2503685/1‑91 (MQ=255)

          cGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATa                                                                              >  1:2503760/1‑91 (MQ=255)

           gATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATAc                                                                             <  1:2873226/91‑1 (MQ=255)

             ttACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAg                                                                           >  1:686598/1‑91 (MQ=255)

              tACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAgt                                                                          >  1:3094696/1‑90 (MQ=255)

              tACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAgt                                                                          >  1:3094384/1‑90 (MQ=255)

              tACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGc                                                                          <  1:3698705/91‑1 (MQ=255)

               aCCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTGGTAAATAATTACTAAATACaaaa                                                                         <  1:3926135/91‑4 (MQ=37)

               aCCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTGGTAAATAATTACTAAATACaaaa                                                                         <  1:2849317/91‑4 (MQ=37)

               aCCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCt                                                                         >  1:2680315/1‑91 (MQ=255)

               aCCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCt                                                                         >  1:2627360/1‑91 (MQ=255)

               aCCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCt                                                                         >  1:2627013/1‑91 (MQ=255)

                  cAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCtttt                                                                      <  1:2858557/91‑1 (MQ=255)

                   aaaaTATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCttttt                                                                     <  1:3793298/91‑1 (MQ=255)

                   aaaaTATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCttttt                                                                     >  1:3793414/1‑91 (MQ=255)

                   aaaaTATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCttttt                                                                     >  1:540230/1‑91 (MQ=255)

                   aaaaTATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCttttt                                                                     >  1:3793397/1‑91 (MQ=255)

                         tCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGt                                                               <  1:4066447/91‑1 (MQ=255)

                         tCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGt                                                               >  1:3085176/1‑91 (MQ=255)

                         tCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGt                                                               <  1:2878441/91‑1 (MQ=255)

                          cACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGtt                                                              >  1:1165130/1‑91 (MQ=255)

                          cACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGtt                                                              >  1:2329094/1‑91 (MQ=255)

                          cACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGtt                                                              >  1:1165486/1‑91 (MQ=255)

                          cACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGtt                                                              >  1:3601256/1‑91 (MQ=255)

                          cACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGtt                                                              <  1:3101682/91‑1 (MQ=255)

                          cACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATATACCAAATACAGCTTTTTCTAAGtt                                                              >  1:3853038/1‑91 (MQ=255)

                            cTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTg                                                            >  1:2402491/1‑91 (MQ=255)

                            cTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTg                                                            >  1:1105648/1‑91 (MQ=255)

                            cTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTg                                                            >  1:2089925/1‑91 (MQ=255)

                            cTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTg                                                            >  1:2402229/1‑91 (MQ=255)

                            cTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGATTg                                                            >  1:3341260/1‑91 (MQ=255)

                            cTTTTAAAATTAGAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTg                                                            >  1:3177375/1‑91 (MQ=255)

                             ttttAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGa                                                           >  1:3099133/1‑91 (MQ=255)

                              tttAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAATAATTATTTACCAAATACAGCTTTTTCTAAGTTTGaa                                                          <  1:3489215/91‑1 (MQ=255)

                               ttAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGaaa                                                         <  1:667404/91‑1 (MQ=255)

                               ttAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGaaa                                                         >  1:3156740/1‑91 (MQ=255)

                               ttAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGaaa                                                         <  1:2955089/91‑1 (MQ=255)

                               ttAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGaaa                                                         <  1:2955029/91‑1 (MQ=255)

                                      tataACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGttt                                                  >  1:3255208/1‑91 (MQ=255)

                                        taACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTa                                                >  1:1566510/1‑91 (MQ=255)

                                            aTAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAat                                            >  1:3263589/1‑91 (MQ=255)

                                            aTAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAat                                            >  1:3258037/1‑91 (MQ=255)

                                               gTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATc                                         <  1:1421165/91‑1 (MQ=255)

                                                   ttattaGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACa                                     <  1:4033612/91‑1 (MQ=255)

                                                   ttattaGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACa                                     <  1:4033937/91‑1 (MQ=255)

                                                   ttattaGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACa                                     >  1:1902589/1‑91 (MQ=255)

                                                     attaGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACatt                                   <  1:2103470/91‑1 (MQ=255)

                                                     attaGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACatt                                   <  1:2505824/91‑1 (MQ=255)

                                                      ttaGTACGACAGTTAAACTTTTGTATATAGTCATTATTTACCAAATCCAGATTTTTATAAGTTTGAAAAACGTTATCAAATATTTACctta                                  >  1:3743560/1‑91 (MQ=255)

                                                       taGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACattat                                 >  1:1905928/1‑91 (MQ=255)

                                                       taGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACattat                                 >  1:3919962/1‑91 (MQ=255)

                                                       taGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACattat                                 >  1:348656/1‑91 (MQ=255)

                                                        aGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACattatt                                >  1:4017729/1‑91 (MQ=255)

                                                         gTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTAttt                               >  1:2664597/1‑91 (MQ=255)

                                                         gTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTAttt                               >  1:3285268/1‑91 (MQ=255)

                                                         gTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTAttt                               >  1:2970002/1‑91 (MQ=255)

                                                          tAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTg                              <  1:1747388/91‑1 (MQ=255)

                                                             gACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAg                           >  1:2373826/1‑91 (MQ=255)

                                                               cAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGAt                         >  1:2299078/1‑91 (MQ=255)

                                                                 gTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGt                       >  1:1956514/1‑91 (MQ=255)

                                                                  ttAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTa                      >  1:693210/1‑91 (MQ=255)

                                                                  ttAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTa                      <  1:3367851/91‑1 (MQ=255)

                                                                  ttAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTa                      >  1:2643275/1‑91 (MQ=255)

                                                                     aaCTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAGGATGTAttt                   <  1:2082011/91‑1 (MQ=255)

                                                                     aaCTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTAttt                   <  1:2865823/91‑1 (MQ=255)

                                                                      aCTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATttg                  >  1:1823435/1‑91 (MQ=255)

                                                                      aCTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATttg                  >  1:1938948/1‑91 (MQ=255)

                                                                       cTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATttgt                 <  1:2285515/91‑1 (MQ=255)

                                                                       aagctGTATTTGGTAAATAATTACTAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATttgt                 >  1:1848860/5‑91 (MQ=25)

                                                                          ttGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATttgttgt              >  1:4039368/1‑91 (MQ=255)

                                                                           tGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATttgttgtt             >  1:2156302/1‑91 (MQ=255)

                                                                            gTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATttgttgttg            >  1:1753029/1‑91 (MQ=255)

                                                                            gTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATttgttgttg            <  1:2776788/91‑1 (MQ=255)

                                                                            gTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATttgttgttg            >  1:1753068/1‑91 (MQ=255)

                                                                            gTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATttgttgttg            >  1:1753293/1‑91 (MQ=255)

                                                                             tATTTGGTAAATAATTACTAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTtgttgttgt           >  1:2617811/1‑91 (MQ=38)

                                                                             tATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTtgttgttgt           <  1:2915564/91‑1 (MQ=255)

                                                                             tATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTtgttgttgt           <  1:1642570/91‑1 (MQ=255)

                                                                              aTTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTTGTTGTTGta          >  1:2838633/1‑91 (MQ=255)

                                                                              aTTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTTGTTGTTGta          >  1:1035861/1‑91 (MQ=255)

                                                                              aTTTAGTAAGTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTTGTTGTTGta          <  1:3976705/91‑1 (MQ=255)

                                                                              aTTTAGTAAGTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTTGTTGTTGta          <  1:3976831/91‑1 (MQ=255)

                                                                                 tAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTTGTTGTTGtatta       >  1:3940063/1‑91 (MQ=255)

                                                                                 tAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTTGTTGTTGtatta       >  1:332402/1‑91 (MQ=255)

                                                                                 tAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTTGTTGTTGtatta       >  1:2543468/1‑91 (MQ=255)

                                                                                     aaTTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTTGTTGTTGTATTATTCg   >  1:324874/1‑91 (MQ=255)

                                                                                      attattTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTTGTTGTTGTATTATTCga  >  1:661894/1‑91 (MQ=255)

                                                                                      attattTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTTGTTGTTGTATTATTCga  >  1:661927/1‑91 (MQ=255)

                                                                                      attattTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTTGTTGTTGTATTATTCga  >  1:662044/1‑91 (MQ=255)

                                                                                      attattTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTTGTTGTTGTATTATTCga  >  1:2667571/1‑91 (MQ=255)

                                                                                      attattTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTTGTTGTTGTATTATTCga  >  1:4017369/1‑91 (MQ=255)

TATAATATAGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTTGTTGTTGTATTATTCGA  >  USA300TCH1516_ALE/1516363‑1516539

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