BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	803,539	A→G	21.7%	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	21.7%	565.0 / 167.0	294	intergenic (‑239/+31)	USA300HOU_RS03920/dtpT	Putative lipid kinase/Di‑/tripeptide transporter
Reads supporting (aligned to +/- strand): ref base A (103/126); new base G (25/38); total (128/166)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.30e-01

Only 100 of 295 total aligned reads displayed.

CTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACAATTGTTGT  >  USA300TCH1516_ALE/803450‑803629

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

  ataATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCAATTTAAATGACtt                                                                                         >  2:1161551/1‑91 (MQ=255)

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

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

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

    aaTGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGTCATTTAAATGGCTTAc                                                                                       <  2:189766/91‑1 (MQ=255)

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

    aaTCAGCATGTTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTAc                                                                                       <  1:3554370/91‑1 (MQ=255)

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

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

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

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

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

       gtgCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGTCATTTAAATGGCTTACTGt                                                                                    <  2:1294314/89‑1 (MQ=255)

       gAGCATATTAAGTTTTTATTTATACTGATATGTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGTCATTTAAATGGCTTACTGt                                                                                    <  2:660598/91‑1 (MQ=38)

       gAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAATCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGt                                                                                    >  2:1550186/1‑91 (MQ=255)

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

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

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

       gAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGt                                                                                    >  1:1454280/1‑91 (MQ=255)

       gAGCATATTAAGTTTTTATGTATACTGATATCTTGAATTTAATTAATAGAAACCTATAACAAAACAGTAAGCCATTTAAAGGACTTACTGt                                                                                    >  2:778855/1‑91 (MQ=255)

        ggCATATTAAGTTTTTATTTATACTGAGAGCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGTCATTTAAATGGCTTACTGtt                                                                                   <  2:2232846/90‑1 (MQ=37)

        aGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGTCATTTAAATGGCTTACTGtt                                                                                   <  2:827016/91‑1 (MQ=255)

        aGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGTCATTTAAATGGCTTACTGtt                                                                                   <  1:3813417/91‑1 (MQ=255)

        aGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGTCATTTAAATGGCTTACTGtt                                                                                   <  2:3561213/91‑1 (MQ=255)

         gCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGTCATTTAAATGGCTTACTGttt                                                                                  <  1:1250485/91‑1 (MQ=255)

         gCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGttt                                                                                  <  2:40161/91‑1 (MQ=255)

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

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

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

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

            tatTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAACAACAGTAAGCCATTTAAATGACTTACTGtttttt                                                                               <  2:3081625/91‑1 (MQ=255)

             atTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGttttttc                                                                              >  2:1107534/1‑90 (MQ=255)

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

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

                 aGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAAtt                                                                          <  1:3858438/91‑1 (MQ=255)

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

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

                      ttatttatACTGATATCTTGAATTTAATTAATAGAAACCTAT‑AACAAACAGTAAGCCATTTAAATGACTTACTGCTTTTTGAATTAGGCca                                                                    >  2:3769102/1‑91 (MQ=255)

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

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

                        atttatACTGATATCTTGGATTTAATGAATAGATACTTCTAAAAAAACATTGAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCcaa                                                                   <  2:3907622/91‑1 (MQ=255)

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

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

                         tttatACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGCAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCcaac                                                                  <  1:3421040/91‑1 (MQ=255)

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

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

                                  atatCTTGAATTTAATTAATAGAAACCTATAAGAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAAc                                                         <  2:201279/91‑1 (MQ=255)

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

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

                                     tCTTGAATTTAATTAATAGAAACATATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGta                                                      >  2:3355594/1‑91 (MQ=255)

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

                                               aattaatAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCAt                                            >  2:565969/1‑91 (MQ=255)

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

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

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

                                                    atAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCtt                                       >  1:559475/1‑91 (MQ=255)

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

                                                      aGAAACCTATAAAAAAACAGTAAGCTATTTATATGACTCACTGGATTTTGAATTAGGCAAACAATATTAACGTATACCTTTCATCGCTTtg                                     >  2:3170121/1‑91 (MQ=255)

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

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

                                                         aaCCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTtgatg                                  >  1:2587734/1‑91 (MQ=255)

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

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

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

                                                               taAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTaaa                            >  1:1712414/1‑91 (MQ=255)

                                                                aaaaaaaCAGTAAGTCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAg                           <  1:3902186/91‑1 (MQ=255)

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

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

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

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

                                                                  aaaaaCAGTAAGTCATTTAAATGGCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGt                         >  2:2179306/1‑91 (MQ=255)

                                                                  aaaaaCAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGt                         <  1:3774425/91‑1 (MQ=255)

                                                                  aaaaaCAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGt                         >  2:102754/1‑91 (MQ=255)

                                                                   aaaaCAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTg                        >  2:3185999/1‑91 (MQ=255)

                                                                       cAGTCAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGCTGATTACAGGTGAGAt                    >  2:3170923/1‑91 (MQ=255)

                                                                       cAGTAAGCCATTTAGATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGAt                    <  2:332948/91‑1 (MQ=255)

                                                                        aGTAAGTCATTTAAATGGCTTACTGTTTGTTGACTTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAGGGTGAGAtt                   >  2:2803862/1‑91 (MQ=38)

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

                                                                         gTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTg                  <  1:1152508/91‑1 (MQ=255)

                                                                         gTAACCCATTGAAGTGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTGCATCGCTTTGATGTTTAAAGGTGAGATTg                  <  2:2513803/91‑1 (MQ=255)

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

                                                                           aagtcATTTAAATGGCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCt                >  1:1925422/5‑91 (MQ=255)

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

                                                                           aaGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTAGACCTTTCATCGCTTTGATGATTATAGGTGAGATTGCt                >  2:3229973/1‑91 (MQ=255)

                                                                            agtcATTTAAATGGCTTACTGTATTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTa               >  1:1616679/4‑91 (MQ=255)

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

                                                                             gCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCCCTTTGCTGATTAAACGTGAGATTGCTac              >  2:1420800/1‑90 (MQ=255)

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

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

                                                                               cATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATa            <  2:1809158/91‑1 (MQ=255)

                                                                                aTTTAAATGGCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATAc           >  1:1550178/1‑91 (MQ=255)

                                                                                aTTTAAATGGCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATAc           >  1:728975/1‑91 (MQ=255)

                                                                                 tttAAATGGCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACa          >  1:2356803/1‑91 (MQ=255)

                                                                                 tttAAATGGCTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACa          >  1:313543/1‑91 (MQ=255)

                                                                                 tttAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACa          <  2:523343/91‑1 (MQ=255)

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

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

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

                                                                                         aCTTACTGTTTTTTGTATTATTCCAACAATATTAACCTCTACCTTTCATCACTTTGCTGATTAAAGGTGAGCTTGCTCATACAATtgttgt  >  2:2487224/1‑91 (MQ=255)

CTATAATGAGCATATTAAGTTTTTATTTATACTGATATCTTGAATTTAATTAATAGAAACCTATAAAAAAACAGTAAGCCATTTAAATGACTTACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACAATTGTTGT  >  USA300TCH1516_ALE/803450‑803629

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