BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	803,528	C→T	16.3%	intergenic (‑228/+42)	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,528	0	C	T	16.3%	359.6 / 54.7	141	intergenic (‑228/+42)	USA300HOU_RS03920/dtpT	Putative lipid kinase/Di‑/tripeptide transporter
Reads supporting (aligned to +/- strand): ref base C (52/66); new base T (6/17); total (58/83)
Fisher's exact test for biased strand distribution p-value = 1.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 147 total aligned reads displayed.

GTTATAGCTGTCTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAAT  >  USA300TCH1516_ALE/803439‑803618

gTTATAGCTGTCTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGcc                                                                                           >  1:1574201/1‑91 (MQ=255)

gTTATAGCTGTCTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGcc                                                                                           >  1:958003/1‑91 (MQ=255)

 ttATAGCTGTCTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCa                                                                                          >  1:1566038/1‑91 (MQ=255)

 ttATAGCTGTCTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCa                                                                                          <  1:2945591/91‑1 (MQ=255)

 ttATAGCTGTCTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCa                                                                                          >  1:1987557/1‑91 (MQ=255)

 ttATAGCTGTCTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCa                                                                                          <  1:759657/91‑1 (MQ=255)

 acataGCTGTCTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCa                                                                                          >  1:911137/3‑91 (MQ=255)

  tataGCTGTCTATAATGAGCATATTCAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCAt                                                                                         <  1:3869535/91‑1 (MQ=255)

  tataGCTGTCTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCAt                                                                                         <  1:3702297/91‑1 (MQ=255)

        tGTCTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAAt                                                                                   >  1:2416395/1‑91 (MQ=255)

         gTCTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATg                                                                                  <  1:1500853/91‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                  gAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGTCATTTAAATGGCTTACTGt                                                                         <  1:1507968/91‑1 (MQ=255)

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

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

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

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

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

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

                   aGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGtt                                                                        <  1:1220121/91‑1 (MQ=255)

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

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

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

                       tatTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGtttttt                                                                    <  1:2473497/91‑1 (MQ=255)

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

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

                       tatTAAGTTTTTATTTATACTGATATCGTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGtttttt                                                                    <  1:222774/91‑1 (MQ=255)

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

                        atTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTg                                                                   <  1:1791628/91‑1 (MQ=255)

                          tAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGaa                                                                 <  1:1282861/91‑1 (MQ=255)

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

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

                            aGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAAtt                                                               >  1:3116255/1‑91 (MQ=255)

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

                                 ttatttatACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGcc                                                          <  1:768314/91‑1 (MQ=255)

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

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

                                  tatttatACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCca                                                         >  1:3036548/1‑91 (MQ=255)

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

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

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

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

                                           tGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTa                                                >  1:943329/1‑91 (MQ=255)

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

                                              tatCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCTATTTAAATGTCTTACTGTTTTTTGAATTAGGCCAACAATATTAACg                                             >  1:1050344/1‑91 (MQ=255)

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

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

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

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

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

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

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

                                                          aattaatAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCAt                                 <  1:1198261/91‑1 (MQ=255)

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

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

                                                            ttaatAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCg                               >  1:2425703/1‑91 (MQ=255)

                                                             taatAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGc                              >  1:2878196/1‑91 (MQ=255)

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

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

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

                                                                 aGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTtg                          <  1:911620/91‑1 (MQ=255)

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

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

                                                                    aaCCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTtgatg                       <  1:2299532/91‑1 (MQ=255)

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

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

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

                                                                       cTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGAtt                    >  1:2149006/1‑91 (MQ=255)

                                                                        tataAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTa                   >  1:2819812/1‑91 (MQ=255)

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

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

                                                                           aaaaaaaCAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAg                <  1:2515993/91‑1 (MQ=255)

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

                                                                               aaaCAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTga            <  1:2658830/91‑1 (MQ=255)

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

                                                                                   aGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGAtt        >  1:847324/1‑91 (MQ=255)

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

                                                                                    gTAAGTCATTTAAATGGCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTg       >  1:3488224/1‑91 (MQ=255)

                                                                                     tAAGTCATTTAAATGGCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGc      >  1:3973502/1‑91 (MQ=255)

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

                                                                                      aaGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCt     >  1:56133/1‑91 (MQ=255)

                                                                                      aaGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCt     >  1:958975/1‑91 (MQ=255)

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

                                                                                         ccATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAAt  <  1:588802/91‑1 (MQ=255)

GTTATAGCTGTCTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAAT  >  USA300TCH1516_ALE/803439‑803618

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