BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,805,818	T→A	6.9%	L431F (TTA→TTT)	pyk ←	Pyruvate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,805,818	0	T	A	6.9%	354.1 / 14.0	102	L431F (TTA→TTT)	pyk	Pyruvate kinase
Reads supporting (aligned to +/- strand): ref base T (45/50); new base A (2/5); total (47/55)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

Only 100 of 107 total aligned reads displayed.

AGTTTCACCAGTTGGTACACCAGCAGTAATAATGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACATTGACGTGCAGTTTCTTCACTTGGA  >  USA300TCH1516_ALE/1805728‑1805908

aGTTTCACCAGTTGGTACACCAGCAGTAATAATGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGttt                                                                                            <  1:1279824/91‑1 (MQ=255)

aGTTTCACCAGTTAGTACACCAGCAGTAATAATGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGttt                                                                                            <  1:1279444/91‑1 (MQ=255)

 gTTTCACCAGTTGGTACACCAGCAGTAATAATGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTa                                                                                           <  1:3100085/91‑1 (MQ=255)

 gTTTCACCAGTTGGTACACCAGCAGTAATAATGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTa                                                                                           >  1:3692547/1‑91 (MQ=255)

   ttCACCAGTTGGTACACCAGCAGTAATAATGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTAAAc                                                                                         <  1:2459150/91‑1 (MQ=255)

   ttCACCAGTTGGTACACCAGCAGTAATAATGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTAAAc                                                                                         <  1:2459124/91‑1 (MQ=255)

   ttCACCAGTTGGTACACCAGCAGTAATAATGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTAAAc                                                                                         <  1:2458692/91‑1 (MQ=255)

   ttCACCAGTTGGTACACCAGCAGTAATAATGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTAAAc                                                                                         <  1:2459123/91‑1 (MQ=255)

   ttCACCAGTTGGTACACCAGCAGTAATAATGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTAAAc                                                                                         <  1:2459100/91‑1 (MQ=255)

    tCACCAGTTGGTACACCAGCAGTAAGAATGATTAAATCAGCATTAGATACTCTGCCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACa                                                                                        <  1:3858114/91‑1 (MQ=255)

     cACCAGTTGGTACACCAGCAGTAATAATGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACaa                                                                                       >  1:3479391/1‑91 (MQ=255)

       ccAGTTGGTACACCAGCAGTAATAATGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATg                                                                                     >  1:1076997/1‑91 (MQ=255)

           ttGGTACACCAGCAGTAATAATGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATc                                                                                 <  1:1217850/91‑1 (MQ=255)

            tGGTACACCAGCAGTAATAATGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCt                                                                                >  1:3158553/1‑91 (MQ=255)

            tGGTACACCAGCAGTAATAATGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCt                                                                                <  1:3883989/91‑1 (MQ=255)

             ggTACACCAGCAGTAATAATGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTg                                                                               <  1:2078513/91‑1 (MQ=255)

               tACACCAGCAGTAATAATGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTa                                                                             >  1:1440105/1‑91 (MQ=255)

               tACACCAGCAGTAATAATGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTa                                                                             >  1:1321260/1‑91 (MQ=255)

                acacCAGCAGTAATAATGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTAc                                                                            <  1:156783/91‑1 (MQ=255)

                 cacCAGCAGTAATAATGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTAct                                                                           >  1:3975187/1‑91 (MQ=255)

                    cagcagTAATAATGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCtt                                                                        >  1:751656/1‑91 (MQ=255)

                      gcagTAATAATGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTAc                                                                      <  1:2127209/91‑1 (MQ=255)

                      gcagTAATAATGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTAc                                                                      <  1:2128449/91‑1 (MQ=255)

                       cagTAATAATGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACg                                                                     <  1:3525660/91‑1 (MQ=255)

                       cagTAATAATGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACg                                                                     >  1:900074/1‑91 (MQ=255)

                       cagTAATAATGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACg                                                                     <  1:238465/91‑1 (MQ=255)

                        agTAATAATGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGt                                                                    <  1:2935925/91‑1 (MQ=255)

                         gTAATAATGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTc                                                                   <  1:2637324/91‑1 (MQ=255)

                         gTAATAATGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTc                                                                   <  1:2637253/91‑1 (MQ=255)

                         gTAATAATGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTc                                                                   <  1:2549839/91‑1 (MQ=255)

                          taataaTGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTcc                                                                  <  1:1236911/91‑1 (MQ=255)

                           aataatGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCt                                                                 >  1:1067335/1‑91 (MQ=255)

                             taatGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCttt                                                               >  1:801030/1‑91 (MQ=255)

                              ggtGATTTAATCATTATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCtttt                                                              >  1:552386/3‑91 (MQ=255)

                              aatGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCtttt                                                              >  1:2918583/1‑91 (MQ=255)

                               atGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCttttt                                                             <  1:2234076/91‑1 (MQ=255)

                                 gATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTa                                                           <  1:3887947/91‑1 (MQ=255)

                                   ttAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAAc                                                         <  1:2540246/91‑1 (MQ=255)

                                      aaTCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTAc                                                      >  1:92635/1‑91 (MQ=255)

                                      aaTCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTAc                                                      >  1:3440458/1‑91 (MQ=255)

                                       aTCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACa                                                     <  1:1406075/91‑1 (MQ=255)

                                        tCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAg                                                    >  1:1526894/1‑91 (MQ=255)

                                          aCCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGt                                                  <  1:928348/91‑1 (MQ=255)

                                          aCCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGt                                                  <  1:4020585/91‑1 (MQ=255)

                                           ccATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGtt                                                 >  1:374960/1‑91 (MQ=255)

                                           ccATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGtt                                                 >  1:375021/1‑91 (MQ=255)

                                           ccATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGtt                                                 <  1:702273/91‑1 (MQ=255)

                                           ccATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGtt                                                 >  1:1686126/1‑91 (MQ=255)

                                            cATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTg                                                >  1:1848223/1‑91 (MQ=255)

                                              ttAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGaa                                              <  1:2642429/91‑1 (MQ=255)

                                              ttAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGaa                                              <  1:3964869/91‑1 (MQ=255)

                                               tAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAAc                                             <  1:482597/91‑1 (MQ=255)

                                                aGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACt                                            <  1:4043898/91‑1 (MQ=255)

                                                aGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACt                                            >  1:427273/1‑91 (MQ=255)

                                                aGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACt                                            >  1:3819614/1‑91 (MQ=255)

                                                 gATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTc                                           >  1:951123/1‑91 (MQ=255)

                                                  aTACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTcc                                          <  1:3601145/91‑1 (MQ=255)

                                                    aCTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTcccc                                        >  1:2409844/1‑91 (MQ=255)

                                                     ctctACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCa                                       >  1:2103960/1‑91 (MQ=255)

                                                     ctctACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCa                                       >  1:1338311/1‑91 (MQ=255)

                                                      tctACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCaa                                      <  1:2551316/91‑1 (MQ=255)

                                                      tctACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCaa                                      <  1:1482824/91‑1 (MQ=255)

                                                          ccAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACaa                                  >  1:1979914/1‑91 (MQ=255)

                                                            aGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAAtt                                <  1:574654/91‑1 (MQ=255)

                                                            aGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCC‑TTTTTAACTACAGGTTGAACTCCCCAAACAATTg                               <  1:1732522/91‑1 (MQ=255)

                                                                tCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAAc                            >  1:134221/1‑91 (MQ=255)

                                                                gcAACAGCTGTTGAAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCGAACAATTGAAc                            >  1:3762097/2‑91 (MQ=255)

                                                                gcAACAGCTGTTGAAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAAc                            >  1:1307843/2‑91 (MQ=255)

                                                                gcAACAGCTGTTGAAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAAc                            >  1:1569278/2‑91 (MQ=255)

                                                                 cAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACa                           >  1:1630833/1‑91 (MQ=255)

                                                                 cAACAGCTGTTGAAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACa                           >  1:2399879/1‑91 (MQ=255)

                                                                  aaCAGCTGTTGAAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACAt                          >  1:2852248/1‑91 (MQ=255)

                                                                  aaCAGCTGTTGAAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACAt                          >  1:2032962/1‑91 (MQ=255)

                                                                   aCAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACAtt                         >  1:3431887/1‑91 (MQ=255)

                                                                   aCAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACAtt                         >  1:2442959/1‑91 (MQ=255)

                                                                     aGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACATTGa                       >  1:672090/1‑91 (MQ=255)

                                                                     aGCTGTTGAAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACATTGa                       >  1:394439/1‑91 (MQ=255)

                                                                      gCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACATTGAc                      <  1:3307069/91‑1 (MQ=255)

                                                                       cTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACATTGACg                     <  1:1959156/91‑1 (MQ=255)

                                                                         gTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACATTGACGTg                   <  1:3570559/91‑1 (MQ=255)

                                                                             cAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACATTGACGTGCAGt               >  1:2605864/1‑91 (MQ=255)

                                                                              aaCTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACATTGACGTGCAGtt              <  1:4028337/91‑1 (MQ=255)

                                                                              aaCTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACATTGACGTGCAGtt              <  1:4028519/91‑1 (MQ=255)

                                                                              aaCTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACATTGACGTGCAGtt              <  1:4028091/91‑1 (MQ=255)

                                                                              aaCTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACATTGACGTGCAGtt              <  1:3953308/91‑1 (MQ=255)

                                                                               aCTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACATTGACGTGCAGttt             <  1:687081/91‑1 (MQ=255)

                                                                               aCTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACATTGACGTGCAGttt             <  1:687082/91‑1 (MQ=255)

                                                                                 tGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACATTGACGTGCAGTttct           <  1:3369597/91‑1 (MQ=255)

                                                                                 tGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACATTGACGTGCAGTttct           <  1:218050/91‑1 (MQ=255)

                                                                                 tGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACATTGACGTGCAGTttct           <  1:3963649/91‑1 (MQ=255)

                                                                                 tGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACATTGACGTGCAGTttct           >  1:3843326/1‑91 (MQ=255)

                                                                                     ttGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACATTGACGTGCAGTTTCTTCAc       <  1:2408092/91‑1 (MQ=255)

                                                                                     ttGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACATTGACGTGCAGTTTCTTCAc       <  1:2408247/91‑1 (MQ=255)

                                                                                     ttGTTAAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACATTGACGTGCAGTTTCTTCAc       >  1:2327746/1‑91 (MQ=255)

                                                                                     ttGTTAAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACATTGACGTGCAGTTTCTTCAc       >  1:3613047/1‑91 (MQ=255)

                                                                                      tGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACATTGACGTGCAGTTTCTTCACt      <  1:852071/91‑1 (MQ=255)

                                                                                       gTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACATTGACGTGCAGTTTCTTCACtt     <  1:1025951/91‑1 (MQ=255)

                                                                                        tttAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACATTGACGTGCAGTTTCTTCACTTg    <  1:2880911/91‑1 (MQ=255)

                                                                                         ttAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACATTGACGTGCAGTTTCTTCACTTgg   <  1:415457/91‑1 (MQ=255)

                                                                                          tAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACATTGACGTGCAGTTTCTTCACTTGGa  >  1:2272829/1‑91 (MQ=255)

AGTTTCACCAGTTGGTACACCAGCAGTAATAATGATTAAATCACCATTAGATACTCTACCAGTTTCAACAGCTGTTGCAACTGCATTGTTTAACAATGCATCTGTACTCTTACGTCCTTTTTTAACTACAGGTTGAACTCCCCAAACAATTGAACATTGACGTGCAGTTTCTTCACTTGGA  >  USA300TCH1516_ALE/1805728‑1805908

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