BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	164,590	T→A	13.5%	N29Y (AAT→TAT)	phnD ←	Phosphate‑import protein PhnD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	164,590	0	T	A	13.5%	334.4 / 40.6	111	N29Y (AAT→TAT)	phnD	Phosphate‑import protein PhnD
Reads supporting (aligned to +/- strand): ref base T (48/48); new base A (2/13); total (50/61)
Fisher's exact test for biased strand distribution p-value = 1.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 117 total aligned reads displayed.

CATTTTGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACTTAAAATTTTTCATTTGATA  >  USA300TCH1516_ALE/164500‑164680

cATTTTGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGAtt                                                                                            >  1:2026255/1‑91 (MQ=255)

cATTTTGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGAtt                                                                                            >  1:813364/1‑91 (MQ=255)

cATTTTGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGAtt                                                                                            >  1:1018112/1‑91 (MQ=255)

 aTTTTGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTa                                                                                           <  1:2806634/91‑1 (MQ=255)

 aTTTTGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTa                                                                                           <  1:3013009/91‑1 (MQ=255)

  ttttGCGAAGGTTCAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTAt                                                                                          >  1:3175332/1‑91 (MQ=255)

   tttGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATc                                                                                         >  1:1321779/1‑91 (MQ=255)

   tttGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATc                                                                                         >  1:728370/1‑91 (MQ=255)

   tttGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTAGATTATc                                                                                         <  1:40102/91‑1 (MQ=255)

   tttGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTAGATAATc                                                                                         <  1:3013890/91‑1 (MQ=255)

   tttGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTAGATAATc                                                                                         <  1:3015625/91‑1 (MQ=255)

   tttGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTAGATAATc                                                                                         <  1:2657031/91‑1 (MQ=255)

   tttGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTAGATAATc                                                                                         <  1:2657096/91‑1 (MQ=255)

    ttGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTAGATAATca                                                                                        <  1:3930495/91‑2 (MQ=255)

    ttGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTAGATAATca                                                                                        <  1:1103303/91‑2 (MQ=255)

     tGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTAGATAATcaa                                                                                       <  1:1107046/91‑3 (MQ=255)

      gCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTaa                                                                                      <  1:2806186/91‑1 (MQ=255)

      gCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTaa                                                                                      >  1:2945371/1‑91 (MQ=255)

      gCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTaa                                                                                      <  1:487499/91‑1 (MQ=255)

      gCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTaa                                                                                      >  1:1420922/1‑91 (MQ=255)

       cGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTaaa                                                                                     >  1:739621/1‑91 (MQ=255)

       cGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTaaa                                                                                     >  1:739819/1‑91 (MQ=255)

       cGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTaaa                                                                                     >  1:599396/1‑91 (MQ=255)

       cGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTAGATAATCaaaa                                                                                     <  1:416019/91‑1 (MQ=255)

        gAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAg                                                                                    >  1:3221357/1‑91 (MQ=255)

        gAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAg                                                                                    <  1:2339115/91‑1 (MQ=255)

        gAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAg                                                                                    <  1:2339091/91‑1 (MQ=255)

         aaGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGa                                                                                   >  1:1846697/1‑91 (MQ=255)

          aGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTAGATAATCAAAAGaa                                                                                  <  1:3364201/91‑1 (MQ=255)

          aGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTAGATAATCAAAAGaa                                                                                  <  1:3364400/91‑1 (MQ=255)

           ggTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAAc                                                                                 >  1:1986358/1‑91 (MQ=255)

            gTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACt                                                                                <  1:1456064/91‑1 (MQ=255)

             tACAAATTGGACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACtt                                                                               <  1:2875607/91‑1 (MQ=255)

             tACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACtt                                                                               <  1:3078875/91‑1 (MQ=255)

             tACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACtt                                                                               <  1:1664647/91‑1 (MQ=255)

             tACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACtt                                                                               <  1:3712655/91‑1 (MQ=255)

              aCAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTg                                                                              <  1:1510858/91‑1 (MQ=255)

                 aaTTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGt                                                                           >  1:259500/1‑91 (MQ=255)

                  aTTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGtt                                                                          >  1:3566162/1‑91 (MQ=255)

                  aTTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGtt                                                                          <  1:283273/91‑1 (MQ=255)

                    tGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTc                                                                        <  1:214317/91‑1 (MQ=255)

                    tGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTc                                                                        <  1:3323396/91‑1 (MQ=255)

                      aaCGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCca                                                                      <  1:510807/91‑1 (MQ=255)

                      aaCGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCca                                                                      >  1:1429593/1‑91 (MQ=255)

                      aaCGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCca                                                                      >  1:1429232/1‑91 (MQ=255)

                      aaCGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCca                                                                      <  1:1861410/91‑1 (MQ=255)

                      aaCGGTTAAGTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCca                                                                      <  1:3965829/91‑1 (MQ=255)

                       aCGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCcac                                                                     >  1:1913189/1‑91 (MQ=255)

                       aCGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCcac                                                                     <  1:1246320/91‑1 (MQ=255)

                       aCGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCcac                                                                     <  1:2787220/91‑1 (MQ=255)

                        cGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCcaca                                                                    >  1:1075159/1‑91 (MQ=255)

                        cGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCcaca                                                                    >  1:1075507/1‑91 (MQ=255)

                           ttAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACAtgc                                                                 >  1:3009779/1‑91 (MQ=255)

                             aaTTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACAtgctg                                                               <  1:1606907/91‑1 (MQ=255)

                                tCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACAtgctgctg                                                            <  1:2700502/91‑1 (MQ=255)

                                tCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACAtgctgctg                                                            <  1:2700490/91‑1 (MQ=255)

                                   ttAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCaa                                                         >  1:3000387/1‑91 (MQ=255)

                                   ttAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCaa                                                         >  1:3000367/1‑91 (MQ=255)

                                   ttAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCaa                                                         >  1:4013811/1‑91 (MQ=255)

                                   ttAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCaa                                                         >  1:3979227/1‑91 (MQ=255)

                                    tAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCaaa                                                        >  1:3191511/1‑91 (MQ=255)

                                    tAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCaaa                                                        <  1:3893392/91‑1 (MQ=255)

                                      ggTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCaaaaa                                                      >  1:774747/1‑91 (MQ=255)

                                       gTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAAc                                                     <  1:449233/91‑1 (MQ=255)

                                       gTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAAc                                                     >  1:3580704/1‑91 (MQ=255)

                                           gTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATg                                                 >  1:1428434/1‑91 (MQ=255)

                                              tcctccTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACt                                              >  1:2167445/1‑91 (MQ=255)

                                                 tcctGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCt                                           <  1:3673715/91‑1 (MQ=255)

                                                     gATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACa                                       >  1:544701/1‑91 (MQ=255)

                                                       tttAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACAtt                                     <  1:412072/91‑1 (MQ=255)

                                                       tttAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACAtt                                     <  1:1368626/91‑1 (MQ=255)

                                                       tttAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTCCCACATGCTGCTGCAAAAACAATGACTGCTAACAtt                                     <  1:412426/91‑1 (MQ=255)

                                                        ttAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTa                                    >  1:3041214/1‑91 (MQ=255)

                                                        ttAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTa                                    >  1:3040954/1‑91 (MQ=255)

                                                          aGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAAt                                  >  1:3467227/1‑91 (MQ=255)

                                                          aGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAAt                                  >  1:1903004/1‑91 (MQ=255)

                                                            aaTCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATAc                                <  1:963613/91‑1 (MQ=255)

                                                            aaTCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATAc                                >  1:4743/1‑91 (MQ=255)

                                                            aaTCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATAc                                <  1:963666/91‑1 (MQ=255)

                                                               ccGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACaaa                             >  1:619403/1‑91 (MQ=255)

                                                                 gAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATa                           >  1:2779163/1‑91 (MQ=255)

                                                                      aTTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAacac                      <  1:875653/91‑1 (MQ=255)

                                                                      aTTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAacac                      <  1:647064/91‑1 (MQ=255)

                                                                      aTTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAacac                      <  1:3193575/91‑1 (MQ=255)

                                                                       ttACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACt                     >  1:2135500/1‑91 (MQ=255)

                                                                       ttACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACt                     >  1:578862/1‑91 (MQ=255)

                                                                        tACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACtt                    <  1:1623026/91‑1 (MQ=255)

                                                                        tACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACtt                    <  1:3665124/91‑1 (MQ=255)

                                                                          cTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACTTaa                  >  1:3022351/1‑91 (MQ=255)

                                                                          cTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACTTaa                  >  1:3022330/1‑91 (MQ=255)

                                                                            aGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACTTaaaa                <  1:60709/91‑1 (MQ=255)

                                                                            aGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACTTaaaa                <  1:3098957/91‑1 (MQ=255)

                                                                               gTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACTTAAAAttt             >  1:2635703/1‑91 (MQ=255)

                                                                                ttCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACTTAAAAtttt            >  1:72515/1‑91 (MQ=255)

                                                                                 tCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACTTAAAAttttt           <  1:348425/91‑1 (MQ=255)

                                                                                  cTTTAGATAATCAAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACTTAAAATTTTTc          >  1:3772752/1‑91 (MQ=255)

                                                                                  cTTTAGATAATCAAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACTTAAAATTTTTc          >  1:1674885/1‑91 (MQ=255)

                                                                                   ttttGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACTTAAAATTTTTCa         >  1:2511069/1‑91 (MQ=255)

                                                                                    tttGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACTTAAAATTTTTCAt        >  1:3358869/1‑91 (MQ=255)

                                                                                          tATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACTTAAAATTTTTCATTTGATa  <  1:4117822/91‑1 (MQ=255)

CATTTTGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACTTAAAATTTTTCATTTGATA  >  USA300TCH1516_ALE/164500‑164680

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