BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	803,539	A→G	5.3%	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	5.3%	1222.4 / 35.4	353	intergenic (‑239/+31)	USA300HOU_RS03920/dtpT	Putative lipid kinase/Di‑/tripeptide transporter
Reads supporting (aligned to +/- strand): ref base A (180/153); new base G (8/11); total (189/164)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

Only 100 of 367 total aligned reads displayed.

TCTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACAATTGTTGT  >  USA300TCH1516_ALE/803449‑803629

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

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

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

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

  tataATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACt                                                                                          >  1:5627074/1‑91 (MQ=255)

   ataATGAGCATATTAAGTTTATATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGTCATTTAAATGGCtt                                                                                         <  2:4115180/91‑1 (MQ=38)

    taATGAGCATATTAAGTTTATATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGTCATTTAAATGGCTTa                                                                                        <  2:4311179/91‑1 (MQ=255)

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

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

      aTGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAAGGACTTACt                                                                                      >  2:4823083/1‑91 (MQ=255)

      aTGAGCATATTAAGTTTTTATTTATACTGATATCTTGAACTTAATTAATCGAAACCTATAAAAAACCAGTACGCCATTTAAATGACTTACt                                                                                      >  2:3171956/1‑91 (MQ=255)

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

       tGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTg                                                                                     >  2:301414/1‑91 (MQ=255)

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

         aGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGtt                                                                                   >  2:173985/1‑91 (MQ=255)

         aGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGtt                                                                                   >  1:3965806/1‑91 (MQ=255)

         aGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGtt                                                                                   >  1:5832322/1‑91 (MQ=255)

           cATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGtttt                                                                                 >  1:5219550/1‑91 (MQ=255)

           cATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGtttt                                                                                 <  1:499857/91‑1 (MQ=255)

           cATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGtttt                                                                                 >  2:2322267/1‑91 (MQ=255)

           cATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATACAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGtttt                                                                                 >  1:4414761/1‑91 (MQ=255)

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

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

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

            atatTAAGTTTTTATTTATACTCATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGttttt                                                                                <  2:5773305/91‑1 (MQ=255)

             tatTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGATACCTATAAAGAAACAGTAAGCCATTTAAATGACTTACTGtttttt                                                                               <  2:5645816/91‑1 (MQ=255)

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

              atTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTg                                                                              >  2:2191576/1‑91 (MQ=255)

                   gTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTa                                                                         >  1:459322/1‑91 (MQ=255)

                     ttttATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAgg                                                                       >  1:3606317/1‑91 (MQ=255)

                       ttatttatACTGATATCTTGTATTTAATTAATAGCAACCTATAAAATAACAGTAAGCCATTTATATAACTTACTGTTTTTTGAATTAGGcc                                                                     >  2:3496109/1‑91 (MQ=255)

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

                         atttatACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCcaa                                                                   >  1:4783344/1‑91 (MQ=255)

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

                          tttatACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCcaac                                                                  >  1:3010655/1‑91 (MQ=255)

                           ttatACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCcaaca                                                                 >  1:4705545/1‑91 (MQ=255)

                             ataCTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAat                                                               >  2:1376196/1‑91 (MQ=255)

                             ataCTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAat                                                               >  1:4215620/1‑91 (MQ=255)

                                cTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATAtt                                                            <  2:2630018/91‑1 (MQ=255)

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

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

                                   atatCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAAc                                                         <  2:91283/91‑1 (MQ=255)

                                      tCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGta                                                      >  1:5071084/1‑91 (MQ=255)

                                      tCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGta                                                      <  2:3752536/91‑1 (MQ=255)

                                       cTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGtat                                                     <  1:957441/91‑1 (MQ=255)

                                           aaTTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCt                                                 >  2:5288230/1‑91 (MQ=255)

                                            aTTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCtt                                                <  1:2322267/91‑1 (MQ=255)

                                            aTTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCtt                                                <  2:5345521/91‑1 (MQ=255)

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

                                             tttAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCttt                                               <  1:4258821/91‑1 (MQ=255)

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

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

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

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

                                                 attaatAGAACCCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATc                                           >  2:2738857/1‑91 (MQ=255)

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

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

                                                   taatAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATATCTTTCATCGc                                         <  1:2639596/91‑1 (MQ=255)

                                                     atAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCtt                                       <  1:3440716/91‑1 (MQ=255)

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

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

                                                      tAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTGTCATCGCttt                                      <  2:4870062/91‑1 (MQ=255)

                                                      aaGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCttt                                      <  2:4214308/90‑1 (MQ=255)

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

                                                       aGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCAAACAATATTAACGTATACCTTTCATTGCTTtg                                     >  1:4281808/1‑91 (MQ=255)

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

                                                         aaaCCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTtgat                                   <  1:1960943/91‑1 (MQ=255)

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

                                                               ataAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTaa                             >  2:55051/1‑91 (MQ=255)

                                                                 aaaaaaaCAGGAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAg                           <  2:4669785/91‑1 (MQ=255)

                                                                  aaaaaaCAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAgg                          >  1:5515748/1‑91 (MQ=255)

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

                                                                      aaCAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTgag                      <  2:584175/91‑1 (MQ=255)

                                                                       aCAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTgaga                     <  1:3823173/91‑1 (MQ=255)

                                                                       aCAGTAAGCCATTTAAATGACTTACAGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTgaga                     >  2:5573687/1‑91 (MQ=255)

                                                                        cAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGAt                    <  1:2039505/91‑1 (MQ=255)

                                                                        cAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGAt                    >  1:805264/1‑91 (MQ=255)

                                                                        cAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGAt                    >  2:2470818/1‑91 (MQ=255)

                                                                        cAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGAt                    >  1:3223426/1‑91 (MQ=255)

                                                                        cAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGAt                    >  1:3534674/1‑91 (MQ=255)

                                                                           tAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGc                 >  1:2946838/1‑91 (MQ=255)

                                                                            aaGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCt                <  1:2821087/91‑1 (MQ=255)

                                                                            aaGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCt                <  2:3659451/91‑1 (MQ=255)

                                                                             aGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTa               <  1:2548951/91‑1 (MQ=255)

                                                                             aGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTa               >  2:4547104/1‑91 (MQ=255)

                                                                               tcATTTAAATGGCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGAGAGATTGCTAAt             >  2:3867833/2‑91 (MQ=255)

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

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

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

                                                                                 aTTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATAc           <  2:1142749/91‑1 (MQ=255)

                                                                                  tttAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACa          >  2:3090721/1‑91 (MQ=255)

                                                                                   ttAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACaa         <  2:3966611/91‑1 (MQ=255)

                                                                                   ttAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACaa         <  1:4367184/91‑1 (MQ=255)

                                                                                   ttAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACaa         <  2:1411233/91‑1 (MQ=255)

                                                                                   ttAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACaa         <  1:692310/91‑1 (MQ=255)

                                                                                   ttAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACaa         <  1:1159844/91‑1 (MQ=255)

                                                                                   ttAAATGACTTACTGTTTCTTGAATTAGGCCAGCAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACaa         <  1:5747441/91‑1 (MQ=255)

                                                                                    tAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACAAt        >  2:1561697/1‑91 (MQ=255)

                                                                                     aaaTGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACAAtt       >  2:5813757/1‑91 (MQ=255)

                                                                                        tGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACAAttgtt    >  1:3278086/1‑91 (MQ=255)

TCTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACAATTGTTGT  >  USA300TCH1516_ALE/803449‑803629

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