BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	598,495	C→T	11.3%	intergenic (+44/+238)	tuf → / ← yxeP_2	Elongation factor Tu/putative hydrolase YxeP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	598,495	0	C	T	11.3%	706.7 / 55.1	244	intergenic (+44/+238)	tuf/yxeP_2	Elongation factor Tu/putative hydrolase YxeP
Reads supporting (aligned to +/- strand): ref base C (117/97); new base T (9/19); total (128/116)
Fisher's exact test for biased strand distribution p-value = 2.79e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 245 total aligned reads displayed.

GTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTAT  >  USA300TCH1516_ALE/598405‑598584

   gACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTAAATCGAg                                                                                          >  1:1962540/1‑91 (MQ=255)

      gTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAgggg                                                                                       <  1:1118163/91‑1 (MQ=255)

      gTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAgggg                                                                                       <  1:1280133/91‑1 (MQ=255)

      gTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAgggg                                                                                       >  2:2728964/1‑91 (MQ=255)

      gTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGACCCCTCGATTGAggga                                                                                       <  1:2312314/91‑2 (MQ=255)

         cTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑t                                                                                   <  2:2377422/91‑1 (MQ=255)

          tGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑tt                                                                                  >  1:1176416/1‑91 (MQ=255)

          tGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑tt                                                                                  <  2:665289/91‑1 (MQ=255)

          tGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAAAAGACCCCTCGATTGAGGGATC‑tt                                                                                  <  1:2532306/91‑1 (MQ=37)

           gTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑tt‑t                                                                                <  2:1494765/91‑1 (MQ=255)

            tAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑tt‑tt                                                                               >  1:3648582/1‑91 (MQ=255)

               gATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑tt‑ttttt                                                                            >  1:2919314/1‑91 (MQ=255)

               gATCAGGCGTTGTTACTGAAATCATTAAAGAATTTCTAATTTCTTAGATTTTATATAAAAAAAAGACCCCTCGATTGAGGGATC‑tt‑ctttt                                                                            <  1:3173383/91‑1 (MQ=25)

                aTCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAAAAGACCCCTCGATTGAGGGATC‑tt‑cttttt                                                                           <  1:2210710/91‑2 (MQ=25)

                  cAGTCGTTGTTACTGAAATCATTAAAGAATTTCTAATTTCTTAGATTTTATATAAAAAAAAGCCCCCTCTATTGAGGGAGC‑TT‑CTTTTTAt                                                                         <  2:216336/91‑1 (MQ=255)

                  cAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAAt                                                                         >  1:1353313/1‑91 (MQ=255)

                    ggCGTTGTTGCTGAAAGCAGTAAATAATTTCTGATTGCTTAGAGTTTATATAAAAAGAGGATCTCTCAAGCGAGGGGTC‑TT‑TTTTTAAtgt                                                                       <  2:600865/91‑1 (MQ=255)

                    ggCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAAtgt                                                                       <  2:3458275/91‑1 (MQ=255)

                    ggCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAAtgt                                                                       <  1:2264499/91‑1 (MQ=255)

                     gCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAAtgtg                                                                      >  1:2785544/1‑91 (MQ=255)

                     gCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAAAAGACCCCTCGATTGAGGGATC‑TTCTTTTTatata                                                                       <  1:873807/91‑5 (MQ=18)

                        ttgttACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTaa                                                                   >  2:1190560/1‑91 (MQ=255)

                        ttgttACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTaa                                                                   <  2:2685199/91‑1 (MQ=255)

                         tgttACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTaaa                                                                  <  2:267798/91‑1 (MQ=255)

                         tgttACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTaaa                                                                  <  1:511837/91‑1 (MQ=255)

                              tagAAATTATTTAATGATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTg                                                             >  2:2222816/3‑91 (MQ=255)

                              cTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTACATTTTg                                                             >  2:1304006/1‑91 (MQ=255)

                              cTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTg                                                             >  1:1537590/1‑91 (MQ=255)

                               tGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGt                                                            <  1:2978291/91‑1 (MQ=255)

                               tGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGt                                                            <  2:1282293/91‑1 (MQ=255)

                               tGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGt                                                            <  1:1679280/91‑1 (MQ=255)

                                gAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTa                                                           >  2:2473879/1‑91 (MQ=255)

                                   aTCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATg                                                        <  2:293974/91‑1 (MQ=255)

                                    tCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATgg                                                       >  2:1610268/1‑91 (MQ=255)

                                    tCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATgg                                                       >  2:1610347/1‑91 (MQ=255)

                                     cATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGc                                                      <  2:3722693/91‑1 (MQ=255)

                                     cATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGc                                                      <  2:1353357/91‑1 (MQ=255)

                                      aTTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCt                                                     >  2:2518457/1‑91 (MQ=255)

                                       ttAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTCTTT‑TTTTTAATGTGTAAATTTTGTAATGGCt                                                     <  1:1323799/91‑1 (MQ=255)

                                       ttAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTa                                                    >  1:323258/1‑91 (MQ=255)

                                        tAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTAt                                                   <  1:3832777/91‑1 (MQ=255)

                                        tAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATACCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATAGatat                                                   >  2:2619871/1‑91 (MQ=255)

                                          aataatTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTc                                                 >  1:3943022/1‑91 (MQ=255)

                                          aataatTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTc                                                 >  1:3311209/1‑91 (MQ=255)

                                           ataatTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCg                                                >  1:1694693/1‑91 (MQ=255)

                                            taatttctaatttcTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGa                                               >  1:839042/1‑91 (MQ=255)

                                             aatttctaatttctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGAt                                              <  1:766928/91‑1 (MQ=255)

                                               tttctaatttctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGAttt                                            >  1:2126103/1‑91 (MQ=255)

                                                ttctaatttctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTa                                           >  2:3846536/1‑91 (MQ=255)

                                                  ctaatttctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTACCTTTTGTAATCGCTATTAGATTTaga                                         >  2:3427207/1‑91 (MQ=255)

                                                  ctaatttctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTaga                                         >  2:838982/1‑91 (MQ=255)

                                                   taatttctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTagaa                                        >  2:2489885/1‑91 (MQ=255)

                                                   taatttctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTagaa                                        <  1:989452/91‑1 (MQ=255)

                                                   taatttctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTagaa                                        >  2:401654/1‑91 (MQ=255)

                                                   taatttctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTagaa                                        <  2:421387/91‑1 (MQ=255)

                                                    aatttctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGTGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTagaag                                       >  2:1941916/1‑91 (MQ=255)

                                                     atttctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTagaaga                                      >  2:73715/1‑91 (MQ=255)

                                                        tctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACa                                   <  2:3591388/91‑1 (MQ=255)

                                                        tctTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACa                                   >  1:619330/1‑91 (MQ=255)

                                                          ttAGATTTTATATAAAAAGAAGATCCCTCAATTGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACaat                                 >  2:2868330/1‑91 (MQ=255)

                                                          ttAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACaat                                 >  1:1463517/1‑91 (MQ=255)

                                                          ttAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACaat                                 >  2:2342595/1‑91 (MQ=255)

                                                          ttAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACaat                                 >  1:2109839/1‑91 (MQ=255)

                                                            aGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACaataa                               >  1:420385/1‑91 (MQ=255)

                                                            aGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACaataa                               >  2:2636889/1‑91 (MQ=255)

                                                             gcttttATATAAACAGAAGCCCCCTCAACCGAGGGGTC‑CT‑TTTTTAATGTGTAAATTTTGTAAGGGCTATTCGATTTAGACGAACaataat                              >  2:3651771/3‑91 (MQ=255)

                                                              aTTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTGTGTGATGGCTATTCGATTTAGAAGAACAATAAtt                             <  2:1421851/91‑1 (MQ=255)

                                                                  tatatgAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtgatg                         <  2:512565/91‑1 (MQ=255)

                                                                  tatataAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtgatg                         <  1:369813/91‑1 (MQ=255)

                                                                   atataAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtgatga                        >  1:1870711/1‑91 (MQ=255)

                                                                   atataAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtgatga                        >  2:1799851/1‑91 (MQ=255)

                                                                   atataAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtgatga                        >  1:2867861/1‑91 (MQ=255)

                                                                   atataAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtgatga                        >  2:78372/1‑91 (MQ=255)

                                                                   atataAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATtgatga                        >  2:3338266/1‑91 (MQ=255)

                                                                    tataAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGaa                       >  2:2038095/1‑91 (MQ=255)

                                                                     ataAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGaaa                      <  2:105678/91‑1 (MQ=255)

                                                                     ataAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGaaa                      >  1:1018340/1‑91 (MQ=255)

                                                                       aaaaaGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAga                    <  2:1216438/91‑1 (MQ=255)

                                                                       aaaaaGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAga                    <  2:1663539/91‑1 (MQ=255)

                                                                        aaaaGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAgac                   <  1:2589366/91‑1 (MQ=255)

                                                                         aaaGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAgact                  <  1:312429/91‑1 (MQ=255)

                                                                           agaagaTCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAgactga                <  1:241134/91‑1 (MQ=255)

                                                                           agaagaTCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAgactga                <  1:2702459/91‑1 (MQ=255)

                                                                              agaccccTCGATTGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACta             >  1:714759/5‑91 (MQ=255)

                                                                              agaTCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACta             <  2:3233152/91‑1 (MQ=255)

                                                                              agaTCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACta             <  2:3095922/91‑1 (MQ=255)

                                                                              agaTCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACta             >  1:139557/1‑91 (MQ=255)

                                                                               gaTCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaa            >  2:3314474/1‑91 (MQ=255)

                                                                                aTCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaat           >  1:3734842/1‑91 (MQ=255)

                                                                                aTCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaat           >  1:3524402/1‑91 (MQ=255)

                                                                                 tCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaata          >  2:3273171/1‑91 (MQ=255)

                                                                                 tCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaata          >  1:1704519/1‑91 (MQ=255)

                                                                                 tCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATCCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaata          >  2:2895903/1‑91 (MQ=255)

                                                                                  cccTCGATTGAGGGATC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACtaataa         >  1:3439256/1‑91 (MQ=255)

                                                                                   ccTCGATTGAGGGATC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATaaa        >  2:3548375/1‑91 (MQ=255)

                                                                                   ccTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATaaa        >  2:1697200/1‑91 (MQ=255)

                                                                                     tCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAAc      <  2:2834162/91‑1 (MQ=255)

                                                                                      cAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACt     <  1:3612792/91‑1 (MQ=255)

                                                                                      cAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACt     >  1:3346782/1‑91 (MQ=255)

                                                                                        aTCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTta   >  2:3012935/1‑91 (MQ=255)

GTGGACGTACTGTAGGATCAGGCGTTGTTACTGAAATCATTAAATAATTTCTAATTTCTTAGATTTTATATAAAAAGAAGATCCCTCAATCGAGGGGTC‑TT‑TTTTTAATGTGTAAATTTTGTAATGGCTATTCGATTTAGAAGAACAATAATTGATGAAAGACTGACTAATAAAACTTAT  >  USA300TCH1516_ALE/598405‑598584

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