BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	803,539	A→G	20.4%	intergenic (‑239/+31)	USA300HOU_RS03920 ← / ← dtpT	Putative lipid kinase/Di‑/tripeptide transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	803,539	0	A	G	20.4%	315.8 / 77.3	148	intergenic (‑239/+31)	USA300HOU_RS03920/dtpT	Putative lipid kinase/Di‑/tripeptide transporter
Reads supporting (aligned to +/- strand): ref base A (47/70); new base G (16/14); total (64/84)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

Only 100 of 150 total aligned reads displayed.

TCTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACAATTGTTGT  >  USA300TCH1516_ALE/803449‑803629

tCTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGa                                                                                             >  1:69172/1‑91 (MQ=255)

   ataATGAGCATATTAAGTTTTTATTTATTCTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACtt                                                                                          <  1:1935491/91‑1 (MQ=255)

   ataATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGTCATTTAAATGGCtt                                                                                          <  1:1537562/91‑1 (MQ=255)

   ataATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGTCATTTAAATGGCtt                                                                                          <  1:1537454/91‑1 (MQ=255)

   ataATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGTCATTTAAATGGCtt                                                                                          <  1:342509/91‑1 (MQ=255)

    taATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGTCATTTAAATGGCTTa                                                                                         <  1:1507259/91‑1 (MQ=255)

    taATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGTCATTTAAATGGCTTa                                                                                         <  1:3244959/91‑1 (MQ=255)

    taATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTa                                                                                         <  1:1000500/91‑1 (MQ=255)

    taATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTa                                                                                         <  1:1000486/91‑1 (MQ=255)

     aaTGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTAc                                                                                        <  1:3130053/91‑1 (MQ=255)

      aTGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACt                                                                                       >  1:752476/1‑91 (MQ=255)

      aTGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACt                                                                                       <  1:906410/91‑1 (MQ=255)

      aTGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACATATAAAAAAACAGTAAGCCATTTAAATGACTTACt                                                                                       >  1:2930543/1‑91 (MQ=255)

       tGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGTCATTTAAATGGCTTACTg                                                                                      <  1:3080906/91‑1 (MQ=255)

       tGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTg                                                                                      <  1:430098/91‑1 (MQ=255)

       tGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTg                                                                                      <  1:2949563/91‑1 (MQ=255)

        gAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGTCATTTAAATGGCTTACTG‑t                                                                                    <  1:3223031/91‑1 (MQ=255)

        gAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGTCATTTAAATGGCTTACTG‑t                                                                                    <  1:2016292/91‑1 (MQ=255)

        gAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑t                                                                                    <  1:3904771/91‑1 (MQ=255)

        gAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑t                                                                                    <  1:943911/91‑1 (MQ=255)

        gAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑t                                                                                    <  1:1399710/91‑1 (MQ=255)

        gAGCAGATTAAGTTTTTATTTATACAGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGTCATTTAAATGGCTTACTG‑t                                                                                    <  1:3861506/91‑1 (MQ=37)

         aGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑tt                                                                                   <  1:1220105/91‑1 (MQ=255)

         aGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑tt                                                                                   <  1:337237/91‑1 (MQ=255)

            atatTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑ttttt                                                                                <  1:3811944/91‑1 (MQ=255)

            atatTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑ttttt                                                                                <  1:1280891/91‑1 (MQ=255)

            atatTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑ttttt                                                                                >  1:781582/1‑91 (MQ=255)

            atatTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑ttttt                                                                                <  1:3320082/91‑1 (MQ=255)

             tatTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑tttttt                                                                               <  1:222773/91‑1 (MQ=255)

             tatTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑tttttt                                                                               <  1:2473594/91‑1 (MQ=255)

             tatTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑tttttt                                                                               <  1:222311/91‑1 (MQ=255)

              atTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGTCATTTAAATGGCTTACTG‑ttttttt                                                                              <  1:233060/91‑2 (MQ=255)

                 aaGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAAt                                                                           >  1:4146088/1‑91 (MQ=255)

                 aaGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAAt                                                                           >  1:3199616/1‑91 (MQ=255)

                      tttatttaTACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGc                                                                      <  1:2456650/91‑1 (MQ=255)

                      tttatttaTACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGc                                                                      <  1:2456843/91‑1 (MQ=255)

                       ttatttatACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGcc                                                                     <  1:1298571/91‑1 (MQ=255)

                        tatttatACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCca                                                                    <  1:370709/91‑1 (MQ=255)

                          tttatACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCcaac                                                                  <  1:3631953/91‑1 (MQ=255)

                           ttatACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGcctggc                                                                 <  1:2659753/91‑5 (MQ=255)

                               aCTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAatat                                                             >  1:2006652/1‑91 (MQ=255)

                               aCTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAatat                                                             >  1:2006414/1‑91 (MQ=255)

                                cTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATAtt                                                            >  1:1799114/1‑91 (MQ=255)

                                 tGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTa                                                           <  1:1743069/91‑1 (MQ=255)

                                 tGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTa                                                           >  1:943874/1‑91 (MQ=255)

                                    tatCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACg                                                        >  1:1048328/1‑91 (MQ=255)

                                     atCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGt                                                       <  1:4124246/91‑1 (MQ=255)

                                          gAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATAcc                                                  >  1:595555/1‑91 (MQ=255)

                                          gAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATAcc                                                  >  1:1194264/1‑91 (MQ=255)

                                            aTTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCtt                                                >  1:1870769/1‑91 (MQ=255)

                                            aTTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCtt                                                >  1:3037516/1‑91 (MQ=255)

                                            aTTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCtt                                                >  1:1502644/1‑91 (MQ=255)

                                               taattaatAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCa                                             >  1:3950607/1‑91 (MQ=255)

                                               taattaatAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCa                                             <  1:2406444/91‑1 (MQ=255)

                                                aattaatAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCAt                                            <  1:876499/91‑1 (MQ=255)

                                                 attaatAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATc                                           >  1:2091174/1‑91 (MQ=255)

                                                   taatAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATCCCTTTCATCGc                                         >  1:3761683/1‑91 (MQ=255)

                                                   taatAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGc                                         >  1:2877706/1‑91 (MQ=255)

                                                    aatAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCt                                        <  1:2544376/91‑1 (MQ=255)

                                                    aatAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCt                                        <  1:2544475/91‑1 (MQ=255)

                                                    aatAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCt                                        >  1:939640/1‑91 (MQ=255)

                                                     atAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCtt                                       <  1:1035019/91‑1 (MQ=255)

                                                      tAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCttt                                      <  1:526296/91‑1 (MQ=255)

                                                      tAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCttt                                      >  1:936906/1‑91 (MQ=255)

                                                       aGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTtg                                     <  1:911597/91‑1 (MQ=255)

                                                        gAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTtga                                    <  1:1475266/91‑1 (MQ=255)

                                                          aaCCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTtgatg                                  <  1:2298952/91‑1 (MQ=255)

                                                           aCCGATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTtaatga                                 <  1:3268111/91‑1 (MQ=255)

                                                            ccTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTgatgat                                >  1:2018683/1‑91 (MQ=255)

                                                            ccTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTgatgat                                >  1:1778874/1‑91 (MQ=255)

                                                             cTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGAtt                               >  1:1199941/1‑91 (MQ=255)

                                                              tataAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTa                              >  1:2803526/1‑91 (MQ=255)

                                                              tataAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTa                              <  1:2840912/91‑1 (MQ=255)

                                                               ataAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTaa                             <  1:2635050/91‑1 (MQ=255)

                                                                  aaaaaaCAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAgg                          <  1:2491773/91‑1 (MQ=255)

                                                                  aaaaaaCAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAgg                          <  1:2287592/91‑1 (MQ=255)

                                                                   aaaaaCAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGt                         >  1:1174408/1‑91 (MQ=255)

                                                                        cAGTAAGTCATTTAAATGGCTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGAt                    >  1:3607434/1‑91 (MQ=255)

                                                                          gTAAGTCATTTAAATGGCTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTg                  >  1:3488346/1‑91 (MQ=255)

                                                                          gTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTg                  <  1:761257/91‑1 (MQ=255)

                                                                            aagtcATTTAAATGGCTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCt                >  1:789015/5‑91 (MQ=255)

                                                                            aaGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTgag                >  1:261474/1‑89 (MQ=255)

                                                                              gCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTaa              >  1:1999666/1‑91 (MQ=255)

                                                                               tcATTTAAATGGCTTACTGTTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTaa              >  1:2873853/2‑91 (MQ=39)

                                                                               ccATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAAt             <  1:588354/91‑1 (MQ=255)

                                                                                cATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATa            >  1:564216/1‑91 (MQ=255)

                                                                                cATTTAAATAACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATa            <  1:3979440/91‑1 (MQ=255)

                                                                                 aTTTAAATGGCTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATAc           >  1:4024388/1‑91 (MQ=255)

                                                                                 aTTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATAc           <  1:316186/91‑1 (MQ=255)

                                                                                  tttAAATGGCTTACTG‑TTTTTTGAATTATGCCAACAATATTAACGTATACCTCTCATCGCTTTGCTGATTACAGGTGACATTGCTAATACa          >  1:1221316/1‑91 (MQ=255)

                                                                                  tttAAATGGCTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACa          >  1:355324/1‑91 (MQ=255)

                                                                                  tttAAATGGCTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACa          >  1:2784774/1‑91 (MQ=255)

                                                                                  tttAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACa          >  1:2997190/1‑91 (MQ=255)

                                                                                   ttAAATGGCTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACaa         >  1:899645/1‑91 (MQ=255)

                                                                                   ttAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACaa         <  1:2138224/91‑1 (MQ=255)

                                                                                   ttAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACaa         <  1:284989/91‑1 (MQ=255)

                                                                                   ttAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACaa         >  1:2901066/1‑91 (MQ=255)

                                                                                      aaTGGCTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACTAttg      >  1:3460748/1‑91 (MQ=255)

                                                                                      aaTGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACAAttg      >  1:1360263/1‑91 (MQ=255)

                                                                                          aCTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACAATtgttgt  >  1:2693275/1‑91 (MQ=255)

TCTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTG‑TTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACAATTGTTGT  >  USA300TCH1516_ALE/803449‑803629

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