BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,885,368	G→A	7.0%	intergenic (‑150/+176)	ebh_2 ← / ← moeZ_2	Extracellular matrix‑binding protein ebh/putative adenylyltransferase/sulfurtransferase MoeZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,885,368	0	G	A	7.0%	464.5 / 17.0	144	intergenic (‑150/+176)	ebh_2/moeZ_2	Extracellular matrix‑binding protein ebh/putative adenylyltransferase/sulfurtransferase MoeZ
Reads supporting (aligned to +/- strand): ref base G (59/74); new base A (5/5); total (64/80)
Fisher's exact test for biased strand distribution p-value = 7.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 146 total aligned reads displayed.

TACGCTCCCTCATATTATATTTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTAACAGAAAG  >  USA300TCH1516_ALE/1885278‑1885458

 aCGCTCCCTCATATTATATTTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGt                                                                                           >  1:2164493/1‑91 (MQ=255)

  cGCTCCCTCATATTATATTTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGtt                                                                                          >  1:2662100/1‑91 (MQ=255)

  cGCTCCCTCATATTATATTTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGtt                                                                                          <  1:618571/91‑1 (MQ=255)

   gCTCCCTCATATTATATTTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGttt                                                                                         >  1:2282461/1‑91 (MQ=255)

   gCTCCCTCATATTATATTTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGttt                                                                                         <  1:596849/91‑1 (MQ=255)

   gCTCCCTCATATTATATTTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAAGAAAGttt                                                                                         >  1:2083048/1‑91 (MQ=255)

        cTCATATTATATTTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGttttattt                                                                                    >  1:1441696/1‑91 (MQ=255)

         tCATATTATATTTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAActttatttt                                                                                   <  1:2743273/91‑1 (MQ=255)

         tCATATTACATTTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAActttatttt                                                                                   <  1:2743414/91‑1 (MQ=255)

            tattatattTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTttattttatt                                                                                >  1:1708383/1‑91 (MQ=255)

            tattatattTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTttattttatt                                                                                <  1:1088852/91‑1 (MQ=255)

            tattatattTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAACtttattttata                                                                                <  1:674779/91‑2 (MQ=255)

                 tattTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAAt                                                                           <  1:1475249/91‑1 (MQ=255)

                  attTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAAtt                                                                          <  1:269959/91‑1 (MQ=255)

                   tttATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAAttt                                                                         <  1:2623587/91‑1 (MQ=255)

                   tttATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAAttt                                                                         >  1:3160799/1‑91 (MQ=255)

                    ttATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTa                                                                        >  1:87449/1‑91 (MQ=255)

                     tataGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAg                                                                       <  1:2375988/91‑1 (MQ=255)

                     tataGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAg                                                                       >  1:3081780/1‑91 (MQ=255)

                       taGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGtt                                                                     <  1:1879503/91‑1 (MQ=255)

                         gCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTc                                                                   >  1:4086707/1‑91 (MQ=255)

                         gCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTc                                                                   <  1:703156/91‑1 (MQ=255)

                             tattatTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGtt                                                               <  1:3578489/91‑1 (MQ=255)

                             tattatTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGtt                                                               >  1:285745/1‑91 (MQ=255)

                              attattTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGttt                                                              <  1:2488201/91‑1 (MQ=255)

                               ttattTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGtttt                                                             >  1:397510/1‑91 (MQ=255)

                               ttattTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGtttt                                                             <  1:2912517/91‑1 (MQ=255)

                                 attTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGt                                                           <  1:2976807/91‑1 (MQ=255)

                                     tGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAACTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAAtt                                                       <  1:3490688/91‑1 (MQ=255)

                                           atatTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATc                                                 <  1:2895876/91‑1 (MQ=255)

                                           atatTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATc                                                 <  1:950582/91‑1 (MQ=255)

                                           atatTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATc                                                 <  1:3604941/91‑1 (MQ=255)

                                           atatTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATc                                                 <  1:3792776/91‑1 (MQ=255)

                                           atatTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATc                                                 <  1:982856/91‑1 (MQ=255)

                                           atatTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATc                                                 <  1:1950084/91‑1 (MQ=255)

                                            tatTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCt                                                <  1:805967/91‑1 (MQ=255)

                                              ttttAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTta                                              <  1:1474729/91‑1 (MQ=255)

                                                ttAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTCATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTtaat                                            >  1:3307589/1‑91 (MQ=255)

                                                  aaGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTtaataa                                          <  1:4131143/91‑1 (MQ=255)

                                                    gTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATaaaa                                        >  1:1588896/1‑91 (MQ=255)

                                                    gTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATaaaa                                        <  1:1801958/91‑1 (MQ=255)

                                                     ttCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAAc                                       >  1:2216456/1‑91 (MQ=255)

                                                      tCACGCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACa                                      <  1:3459777/91‑1 (MQ=255)

                                                       cACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACAt                                     >  1:2230479/1‑91 (MQ=255)

                                                         cTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATaa                                   >  1:3782483/1‑91 (MQ=255)

                                                         cTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATaa                                   >  1:4149510/1‑91 (MQ=255)

                                                          tCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAAc                                  >  1:2374632/1‑91 (MQ=255)

                                                           cAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACa                                 <  1:973990/91‑1 (MQ=255)

                                                            aaTTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACaa                                <  1:1594386/91‑1 (MQ=255)

                                                              ttttATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGa                              >  1:2107046/1‑91 (MQ=255)

                                                              ttttATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGa                              >  1:1777595/1‑91 (MQ=255)

                                                               tttATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGaa                             <  1:1759177/91‑1 (MQ=255)

                                                                ttATCATCATCATTTTATAAAATAAAGTTTTATTTTACTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAAc                            >  1:3507239/1‑91 (MQ=255)

                                                                   tcatcatcaTTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCt                         >  1:3494264/1‑91 (MQ=255)

                                                                   tcatcatcaTTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCt                         <  1:212522/91‑1 (MQ=255)

                                                                   tcatcatcaTTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCt                         >  1:257389/1‑91 (MQ=255)

                                                                    catcatcatTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCtt                        <  1:3958938/91‑1 (MQ=255)

                                                                    catcatcatTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCtt                        <  1:1648887/91‑1 (MQ=255)

                                                                      tcatcgtTTTATAAAATAAAGTTTTGTTTTATTACTATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTa                      <  1:3863414/91‑1 (MQ=255)

                                                                      tcatcatTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTa                      >  1:1955852/1‑91 (MQ=255)

                                                                          catTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGAtt                  >  1:4012483/1‑91 (MQ=255)

                                                                          catTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGAtt                  <  1:1236898/91‑1 (MQ=255)

                                                                          catTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGAtt                  <  1:1236552/91‑1 (MQ=255)

                                                                          catTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGAtt                  >  1:4013101/1‑91 (MQ=255)

                                                                            ttttATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGt                <  1:1835054/91‑1 (MQ=255)

                                                                            ttttATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGt                <  1:1382173/91‑1 (MQ=255)

                                                                            ttttATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGt                <  1:1834981/91‑1 (MQ=255)

                                                                             tttATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGtt               >  1:2214725/1‑91 (MQ=255)

                                                                              ttATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGttt              <  1:3008715/91‑1 (MQ=255)

                                                                              ttATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGttt              >  1:1051760/1‑91 (MQ=255)

                                                                              ttATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGttt              <  1:954238/91‑1 (MQ=255)

                                                                              taataaaataaaACTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGttt              >  1:4052347/3‑91 (MQ=39)

                                                                              taataaaataaaACTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGttt              >  1:398739/3‑91 (MQ=39)

                                                                               tataAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGtttt             <  1:1417440/91‑1 (MQ=255)

                                                                               aataaaataaaACTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTTGATTGtttt             >  1:3979512/2‑91 (MQ=255)

                                                                               aataaaataaaACTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGtttt             >  1:2561823/2‑91 (MQ=255)

                                                                                ataaaataaaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTg            <  1:3643362/91‑1 (MQ=255)

                                                                                ataaaataaaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTg            <  1:1221081/91‑1 (MQ=255)

                                                                                ataaaataaaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTCTTg            >  1:3740071/1‑91 (MQ=255)

                                                                                 taaaataaaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGt           <  1:2638636/91‑1 (MQ=255)

                                                                                 taaaataaaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGt           <  1:3800446/91‑1 (MQ=255)

                                                                                 taaaataaaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGt           <  1:3800036/91‑1 (MQ=255)

                                                                                  aaaataaaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTa          >  1:146082/1‑91 (MQ=255)

                                                                                  aaaataaaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTa          >  1:128148/1‑91 (MQ=255)

                                                                                  aaaataaaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTa          >  1:3001449/1‑91 (MQ=255)

                                                                                   aaataaaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTaa         >  1:4128162/1‑91 (MQ=255)

                                                                                   aaataaaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTaa         >  1:4128132/1‑91 (MQ=255)

                                                                                   aaataaaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTaa         >  1:1158641/1‑91 (MQ=255)

                                                                                    aataaaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTAAc        <  1:181003/91‑1 (MQ=255)

                                                                                    aataaaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTAAc        >  1:106364/1‑91 (MQ=255)

                                                                                     ataaaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTAACa       >  1:2964372/1‑91 (MQ=255)

                                                                                     ataaaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTAACa       >  1:838773/1‑91 (MQ=255)

                                                                                     ataaaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTAACa       >  1:839058/1‑91 (MQ=255)

                                                                                     ataaaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTAACa       <  1:1889167/91‑1 (MQ=255)

                                                                                     ataaaGTTTGATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAGCGCTTTAGATTGTTTTGTAACa       <  1:3977590/91‑1 (MQ=255)

                                                                                       aaaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTAACAGa     <  1:2880695/91‑1 (MQ=255)

                                                                                       aaaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTAACAGa     <  1:1887589/91‑1 (MQ=255)

                                                                                        aaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTAACAGaa    <  1:3884018/91‑1 (MQ=255)

                                                                                         aGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTAACAGaaa   <  1:492147/91‑1 (MQ=255)

                                                                                          gTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTAACAGAAAg  <  1:2925389/91‑1 (MQ=255)

TACGCTCCCTCATATTATATTTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTAACAGAAAG  >  USA300TCH1516_ALE/1885278‑1885458

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