BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,885,369	T→C	5.3%	intergenic (‑151/+175)	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,369	0	T	C	5.3%	441.1 / 10.0	132	intergenic (‑151/+175)	ebh_2/moeZ_2	Extracellular matrix‑binding protein ebh/putative adenylyltransferase/sulfurtransferase MoeZ
Reads supporting (aligned to +/- strand): ref base T (53/72); new base C (5/2); total (58/74)
Fisher's exact test for biased strand distribution p-value = 2.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 147 total aligned reads displayed.

ACGCTCCCTCATATTATATTTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTAACAGAAAGT  >  USA300TCH1516_ALE/1885279‑1885459

 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)

        tCATATTATATTTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGttttatttt                                                                                    <  1:1503864/91‑1 (MQ=255)

        tCATATTATATTTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGttttatttt                                                                                    <  1:2774207/91‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                       aGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGttt                                                                     >  1:2979771/1‑91 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

                                           tatTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCt                                                 >  1:3744945/1‑91 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                             ttttATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGa                               >  1:2107047/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)

                                                                tATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACg                            >  1:3488948/1‑91 (MQ=255)

                                                                tATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACg                            <  1:3857728/91‑1 (MQ=255)

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

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

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

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

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

                                                                     tcatcatTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTa                       <  1:2027754/91‑1 (MQ=255)

                                                                       atcatTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGa                     <  1:2888057/91‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                                                           tgtaataaaataaaACTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGt                 >  1:2580384/5‑91 (MQ=38)

                                                                            tttATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGtt                <  1:1899609/91‑1 (MQ=255)

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

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

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

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

                                                                             ttATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGttt               <  1:3008715/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:1221081/91‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                                    ataaaGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTAACa        >  1:2964273/1‑91 (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)

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

                                                                                          ttttattttaTTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTAACAGAAAGt  >  1:2777042/1‑91 (MQ=255)

                                                                                          ttttattttaTTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTAACAGAAAGt  <  1:254212/91‑1 (MQ=255)

ACGCTCCCTCATATTATATTTATAGCTATATTATTTTGTACAATATTTTAAGTTCACTCAATTTTATCATCATCATTTTATAAAATAAAGTTTTATTTTATTACAATTTAGTTTCCGTTTTGTAATTAAAATCTTAATAAAACATAACAAGAACGCTTTAGATTGTTTTGTAACAGAAAGT  >  USA300TCH1516_ALE/1885279‑1885459

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