BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,516,453	T→A	6.3%	K114N (AAA→AAT)	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,453	0	T	A	6.3%	387.0 / 13.6	111	K114N (AAA→AAT)	gpsB	Cell cycle protein GpsB
Reads supporting (aligned to +/- strand): ref base T (70/34); new base A (3/4); total (73/38)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

Only 100 of 127 total aligned reads displayed.

TATAATATAGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTTGTTGTTGTATTATTCGAAGAA  >  USA300TCH1516_ALE/1516363‑1516543

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

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

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

     tataGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTAcc                                                                                       >  1:2471919/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: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:2627013/1‑91 (MQ=255)

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

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

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

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

                   aaaaTATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCctttt                                                                         >  1:1399833/1‑91 (MQ=255)

                      atatCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTa                                                                      <  1:2650456/91‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                       ataACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGtttt                                                     >  1:3134177/1‑91 (MQ=255)

                                       ataACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGtttt                                                     <  1:172590/91‑1 (MQ=255)

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

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

                                                 cATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTCATTATTTACCAAATACAGCTGTTTCTAAGTTTGAAATACGTTTTAAAATATCTa                                           >  1:3634251/1‑91 (MQ=255)

                                                  attattAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTAc                                          <  1:3802999/91‑1 (MQ=255)

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

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

                                                   ttattaGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACa                                         <  1:4033612/91‑1 (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:348656/1‑91 (MQ=255)

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

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

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

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

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

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

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

                                                           aaGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGa                                 >  1:3255635/1‑91 (MQ=255)

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

                                                              aCAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGa                              <  1:2606198/91‑1 (MQ=255)

                                                               cAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGAt                             >  1:2299078/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:4091816/91‑1 (MQ=255)

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

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

                                                                        agctGTATTTGGTAAATAATTACTAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATttgtt                    >  1:3954567/4‑91 (MQ=37)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                                      attattTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTTGTTGTTGTATTATTCga      >  1:1805336/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:662103/1‑91 (MQ=255)

                                                                                        tattTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTTGTTGTTGTATTATTCgaag    <  1:2672373/91‑1 (MQ=255)

                                                                                         attTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAATATATCTACATTATTTGAAGATGTATTTGTTGTTGTATTATTCgaaga   <  1:4136958/91‑1 (MQ=255)

                                                                                         attTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTTGTTGTTGTATTATTCgaaga   <  1:378319/91‑1 (MQ=255)

                                                                                         attTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTTGTTGTTGTATTATTCgaaga   <  1:3572753/91‑1 (MQ=255)

                                                                                          tttACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTTGTTGTTGTATTATTCgaagaa  >  1:529116/1‑91 (MQ=255)

TATAATATAGCGATTACCCAAAATATCACTTTTAAAATTATAACATAGTCATTATTAGTAAGACAGTTAAACTTTTGTATTTAGTAATTATTTACCAAATACAGCTTTTTCTAAGTTTGAAATACGTTTTAAAATATCTACATTATTTGAAGATGTATTTGTTGTTGTATTATTCGAAGAA  >  USA300TCH1516_ALE/1516363‑1516543

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