BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,847,006	G→C	15.1%	intergenic (‑210/‑170)	salL ← / → yceI	Adenosyl‑chloride synthase/Protein YceI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,847,006	0	G	C	15.1%	281.1 / 42.1	99	intergenic (‑210/‑170)	salL/yceI	Adenosyl‑chloride synthase/Protein YceI
Reads supporting (aligned to +/- strand): ref base G (38/46); new base C (7/8); total (45/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-01

Only 100 of 101 total aligned reads displayed.

GTGACTATAAAAATGATAATAATAGTCAGATTTATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATATTCGAATTCGAAATAAAATGCTTGCATC  >  USA300TCH1516_ALE/2846917‑2847096

gtGACTATAAAAATGATAATAATAGTCAGATTTATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGa                                                                                           <  1:755696/91‑1 (MQ=255)

 tGACTATAAAAATGATAATAATAGTCAGATTTATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGAc                                                                                          >  1:1802038/1‑91 (MQ=255)

  gACTATAGAAATGATAATAATAGTCAGATTTATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTCACt                                                                                         <  1:2077744/91‑1 (MQ=255)

  gACTATAGAAATGATAATAATAGTCAGATTTATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTCACt                                                                                         <  1:2077894/91‑1 (MQ=255)

  gACTATAAAAATGATAATAATAGTCAGATTTATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACt                                                                                         >  1:2740904/1‑91 (MQ=255)

  gACTATAAAAATGATAATAATAGTCAGATTTATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACt                                                                                         >  1:2740915/1‑91 (MQ=255)

  gACTATAAAAATGATAATAATAGTCAGATTTATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACt                                                                                         <  1:2772594/91‑1 (MQ=255)

  gACTATAAAAATGATAATAATAGTCAGATTTATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTCACt                                                                                         <  1:2078262/91‑1 (MQ=255)

  gACTATAAAAATGATAATAATAGTCAGATTTATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTCACt                                                                                         <  1:3996798/91‑1 (MQ=255)

   aCTATAAAAATGATAATAATAGTCAGATTTATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTa                                                                                        <  1:1362061/91‑1 (MQ=255)

     tataAAAATGATAATAATAGTCAGATTTATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAgg                                                                                      <  1:1732372/91‑1 (MQ=255)

       taAAAATGGTAATAATAGTCAGATTTATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGt                                                                                    <  1:1941317/91‑1 (MQ=255)

         aaaaTGATAATAATAGTCAGATTTATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGtat                                                                                  >  1:1082800/1‑91 (MQ=255)

              gATAATAATAGTCAGATTTATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACaa                                                                             >  1:3641317/1‑91 (MQ=255)

               ataataataGTCAGATTTATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAAt                                                                            <  1:1233714/91‑1 (MQ=255)

                 aataataGTCAGATTTATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAg                                                                          >  1:2872357/1‑91 (MQ=255)

                   taataGTCAGATTTATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAgtg                                                                        >  1:69631/1‑91 (MQ=255)

                   taataGTCAGATTTATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTCACTAGGGTATAACAATAgtg                                                                        <  1:2433911/91‑1 (MQ=255)

                   taataGTCAGATTTATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTCACTAGGGTATAACAATAgtg                                                                        >  1:1457953/1‑91 (MQ=255)

                    aataGTCAGATTTATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAgtgt                                                                       >  1:1418271/1‑91 (MQ=255)

                      taGTCAGATTTATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTgtgg                                                                     <  1:2445302/91‑1 (MQ=255)

                      taGTCAGATTTATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTgtgg                                                                     <  1:2444233/91‑1 (MQ=255)

                      taGTCAGATTTATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTgtgg                                                                     <  1:448252/91‑1 (MQ=255)

                      taGTCAGATTTATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCGTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTgtgg                                                                     <  1:3597318/91‑1 (MQ=255)

                        gTCAGATTTATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTgtggtg                                                                   <  1:4016752/91‑1 (MQ=255)

                        gTCAGATTTATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTgtggtg                                                                   >  1:975514/1‑91 (MQ=255)

                        gTCAGATTTATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTCGAATTAATTGACTAGGGTATAACAATAGTgtggtg                                                                   >  1:1505954/1‑91 (MQ=255)

                           aGATTTATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTCACTAGGGTATAACAATAGTGTGGTGttt                                                                >  1:682490/1‑91 (MQ=255)

                             aTTTATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTta                                                              <  1:1058639/91‑1 (MQ=255)

                             aTTTATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTta                                                              >  1:3281726/1‑91 (MQ=255)

                               ttATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAAATAATTGACTAGGGTATAACAATAGTGTGGTGTTTtata                                                            >  1:1690518/1‑91 (MQ=255)

                                 aTCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAgg                                                          <  1:184374/91‑1 (MQ=255)

                                 aTCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAgg                                                          <  1:428693/91‑1 (MQ=255)

                                 aTCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTCACTAGGGTATAACAATAGTGTGGTGTTTTATAgg                                                          <  1:894780/91‑1 (MQ=255)

                                 aTCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTCACTAGGGTATAACAATAGTGTGGTGTTTTATAgg                                                          <  1:894553/91‑1 (MQ=255)

                                 aTCTCATTAACTAAATTTCACGGACTGGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAgg                                                          <  1:3660474/91‑1 (MQ=255)

                                    tcATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCa                                                       >  1:70470/1‑91 (MQ=255)

                                    tcATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCa                                                       <  1:2131574/91‑1 (MQ=255)

                                    tcATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCa                                                       <  1:2131884/91‑1 (MQ=255)

                                    tcATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTCACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCa                                                       >  1:1280126/1‑91 (MQ=255)

                                     cATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCAt                                                      <  1:697418/91‑1 (MQ=255)

                                     atctAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCAt                                                      <  1:273384/88‑1 (MQ=255)

                                      aTTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCAtt                                                     <  1:3914121/91‑1 (MQ=255)

                                      aTTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCAtt                                                     <  1:1277282/91‑1 (MQ=255)

                                        tAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCt                                                   <  1:419068/91‑1 (MQ=255)

                                         aaCTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCtt                                                  <  1:3407461/91‑1 (MQ=255)

                                           cTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTa                                                >  1:1434087/1‑91 (MQ=255)

                                            tAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTaa                                               <  1:82893/91‑1 (MQ=255)

                                             aaaTTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTaaa                                              <  1:2094044/91‑1 (MQ=255)

                                                     cGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATcac                                      <  1:1417778/91‑1 (MQ=255)

                                                     cGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATcac                                      >  1:2661977/1‑91 (MQ=255)

                                                     cGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATcac                                      >  1:1375671/1‑91 (MQ=255)

                                                     cGGACTTGTAATCACTTCTTTATGTTTGAATTAATTCACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATcac                                      >  1:35722/1‑91 (MQ=255)

                                                      ggACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATcaca                                     >  1:2427276/1‑91 (MQ=255)

                                                       gACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAg                                    >  1:1318260/1‑91 (MQ=255)

                                                       gACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAg                                    >  1:1220245/1‑91 (MQ=255)

                                                       gACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAg                                    <  1:3870159/91‑1 (MQ=255)

                                                       gACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAg                                    >  1:1219119/1‑91 (MQ=255)

                                                        aCTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCGCAGa                                   >  1:2828102/1‑91 (MQ=255)

                                                        aCTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGa                                   >  1:1619797/1‑91 (MQ=255)

                                                        aCTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGa                                   >  1:3622692/1‑91 (MQ=255)

                                                        aCTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGa                                   >  1:3622671/1‑91 (MQ=255)

                                                        aCTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGa                                   >  1:1368010/1‑91 (MQ=255)

                                                          ttGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGaaa                                 >  1:2048487/1‑91 (MQ=255)

                                                            gTAATCACTTCTTTATGTTTGAATTAATTCACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAat                               >  1:2819140/1‑91 (MQ=255)

                                                              aaTCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAatat                             <  1:3840664/91‑1 (MQ=255)

                                                              aaTCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAatat                             <  1:2810671/91‑1 (MQ=255)

                                                               aTCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATatt                            <  1:2015612/91‑1 (MQ=255)

                                                               aTCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATatt                            <  1:1290608/91‑1 (MQ=255)

                                                               aTCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATatt                            <  1:1435471/91‑1 (MQ=255)

                                                               aTCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATatt                            >  1:3051512/1‑91 (MQ=255)

                                                               aTCACTTCTTTATGTTTGAATTAATTCACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATatt                            >  1:158501/1‑91 (MQ=255)

                                                                 cACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATattcg                          <  1:3039792/91‑1 (MQ=255)

                                                                  aCTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATattcga                         >  1:708903/1‑91 (MQ=255)

                                                                   cttcttTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTATTATAGGTCATTCTTTAAATACATCACAGAAAATattcgaa                        <  1:3602661/91‑1 (MQ=255)

                                                                    ttcttTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATattcgaat                       >  1:4093519/1‑91 (MQ=255)

                                                                    ttcttTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATattcgaat                       >  1:4090310/1‑91 (MQ=255)

                                                                     tcttCATGTTAGATTTAATTGACTAGGGTATAACAGTAGTGTGGTGTTTTATAGGTCATTCTTTAGATACATCACAGAAAATattcgaatt                      <  1:1730105/91‑1 (MQ=255)

                                                                      cttTATGTTTGACTTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATattcgaattc                     >  1:3179382/1‑91 (MQ=255)

                                                                      cttTATGTTTGAATTAATTCACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATattcgaattc                     <  1:3638858/91‑1 (MQ=255)

                                                                       tttATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATattcgaattcg                    <  1:2787511/91‑1 (MQ=255)

                                                                       tttATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATattcgaattcg                    <  1:2787539/91‑1 (MQ=255)

                                                                       tttATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATattcgaattcg                    <  1:2787389/91‑1 (MQ=255)

                                                                          aTGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATATTCGAATTCGaaa                 >  1:334490/1‑91 (MQ=255)

                                                                           tGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATATTCGAATTCGAAAt                <  1:192845/91‑1 (MQ=255)

                                                                            gTTTGAATTAATTGACTAGGGTATAACAATCGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATATTCGAATTCGAAATa               >  1:674517/1‑91 (MQ=255)

                                                                            gTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATATTCGAATTCGAAATa               >  1:3566278/1‑91 (MQ=255)

                                                                            gTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATATTCGAATTCGAAATa               <  1:1585074/91‑1 (MQ=255)

                                                                            gTTTGAATTAATTCACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATATTCGAATTCGAAATa               >  1:2504554/1‑91 (MQ=255)

                                                                                gAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATATTCGAATTCGAAATAAAAt           <  1:878608/91‑1 (MQ=255)

                                                                                gAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATATTCGAATTCGAAATAAAAt           >  1:394573/1‑91 (MQ=255)

                                                                                   ttaattGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATATTCGAATTCGAAATAAAATGCt        <  1:2488480/91‑1 (MQ=255)

                                                                                   ttaattGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATATTCGAATTCGAAATAAAATGCt        <  1:2732400/91‑1 (MQ=255)

                                                                                     aattGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATATTCGAATTCGAAATAAAATGCTTg      >  1:2365122/1‑91 (MQ=255)

                                                                                      attGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATATTCGAATTCGAAATAAAATGCTTGc     >  1:1215910/1‑91 (MQ=255)

                                                                                       ttGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATATTCGAATTCGAAATAAAATGCTTGCa    <  1:3061104/91‑1 (MQ=255)

                                                                                       ttGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATATTCGAATTCGAAATAAAATGCTTGCa    <  1:1331509/91‑1 (MQ=255)

                                                                                       ttGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATATTCGAATTCGAAATAAAATGCTTGCa    >  1:125591/1‑91 (MQ=255)

                                                                                         gACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATATTCGAATTCGAAATAAAATGCTTGCATc  <  1:2684007/91‑1 (MQ=255)

                                                                                         gACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATATTCGAATTCGAAATAAAATGCTTGCATc  <  1:1304616/91‑1 (MQ=255)

GTGACTATAAAAATGATAATAATAGTCAGATTTATCTCATTAACTAAATTTCACGGACTTGTAATCACTTCTTTATGTTTGAATTAATTGACTAGGGTATAACAATAGTGTGGTGTTTTATAGGTCATTCTTTAAATACATCACAGAAAATATTCGAATTCGAAATAAAATGCTTGCATC  >  USA300TCH1516_ALE/2846917‑2847096

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