BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,972,372	G→A	12.8%	A23T (GCT→ACT)	prsA →	Foldase protein PrsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,972,372	0	G	A	12.8%	680.7 / 66.0	246	A23T (GCT→ACT)	prsA	Foldase protein PrsA
Reads supporting (aligned to +/- strand): ref base G (106/107); new base A (2/30); total (108/138)
Fisher's exact test for biased strand distribution p-value = 1.41e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 260 total aligned reads displayed.

TTAACTTACAAAGGAGTTTCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCGGTAAAG  >  USA300TCH1516_ALE/1972282‑1972462

ttAACGTACAAAGGAGTTTCAACTAGGAAGATGATAAACTAATTAATCGTTCCCGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGgcg                                                                                            <  1:3672058/91‑1 (MQ=255)

 tAACTTACAAAGGAGTTTCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGgcgc                                                                                           <  1:416084/91‑1 (MQ=255)

  aaCTTACAAAGGAGTTTCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCt                                                                                          <  2:1671542/91‑1 (MQ=255)

   aCTTACAAAGGAGTTTCAACTATGAAGATGATAGACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCACTa                                                                                         <  1:3240383/91‑1 (MQ=255)

   aCTTACAAAGGAGTTTCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCACTa                                                                                         <  2:2187722/91‑1 (MQ=255)

   aCTTACAAAGGAGTTTCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCACTa                                                                                         <  1:3471005/91‑1 (MQ=255)

    cTTACAAAGGAGTTTCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAg                                                                                        >  1:3611362/1‑91 (MQ=255)

    cTTACAAAGGAGTTTCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCACTAg                                                                                        <  1:57278/91‑1 (MQ=255)

    cTTACAAAGGAGTTTCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCACTAg                                                                                        <  1:1618749/91‑1 (MQ=255)

     ttACAAAGGAGTTTCAAGTATGAAGATGATAATCAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCACTAgc                                                                                       <  1:3381762/91‑2 (MQ=255)

     ttACAAAGGAGTTTCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCACTAgc                                                                                       <  1:1062896/91‑2 (MQ=255)

     ttACAAAGGAGTTTCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCACTAgc                                                                                       <  1:624944/91‑2 (MQ=255)

     ttACAAAGGAGTTTCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCACTAgc                                                                                       <  1:1955843/91‑2 (MQ=255)

     ttACAAAGGAGTTTCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATGAGGCGCTTGTGGCACTAgc                                                                                       <  1:1702123/91‑2 (MQ=255)

       aCAAAGGAGTTTCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCCAGTGc                                                                                     >  1:968945/1‑91 (MQ=255)

        cAAAGGAGTTTCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGcc                                                                                    >  1:2812452/1‑91 (MQ=255)

         aaaGGAGTTTCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCca                                                                                   <  1:1109922/91‑1 (MQ=255)

          aaagAGTTTCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCcac                                                                                  <  1:1023343/88‑1 (MQ=255)

           aGGAGTGTCAACTATGAAGATTGTAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCcaca                                                                                 <  2:2650788/91‑1 (MQ=255)

            ggAGTTTCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAg                                                                                >  1:1918375/1‑91 (MQ=255)

              aGTTTCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGAc                                                                              >  1:295668/1‑91 (MQ=255)

                  tCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTa                                                                          >  2:3735974/1‑91 (MQ=255)

                  tCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTa                                                                          >  2:2983731/1‑91 (MQ=255)

                   cAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTaa                                                                         <  1:3495236/91‑1 (MQ=255)

                      cTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGa                                                                      >  1:3845485/1‑91 (MQ=255)

                      cTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGa                                                                      >  2:1467804/1‑91 (MQ=255)

                         tGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGaaaa                                                                   >  1:2549355/1‑91 (MQ=255)

                           aaGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATa                                                                 <  1:1878948/91‑1 (MQ=255)

                            aGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATAc                                                                <  1:3823964/91‑1 (MQ=255)

                            aGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATAc                                                                >  2:504920/1‑91 (MQ=255)

                            aGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATAc                                                                <  2:525948/91‑1 (MQ=255)

                              atgatAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACAt                                                              >  1:837342/1‑91 (MQ=255)

                               tgatAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACAtt                                                             >  1:1588437/1‑91 (MQ=255)

                               tgatAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACAtt                                                             <  1:131914/91‑1 (MQ=255)

                               tgatAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACAtt                                                             <  1:3458930/91‑1 (MQ=255)

                                gatAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTa                                                            >  1:129654/1‑91 (MQ=255)

                                gatAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTa                                                            <  2:3284684/91‑1 (MQ=255)

                                gatAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTa                                                            >  2:378774/1‑91 (MQ=255)

                                  tAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAAt                                                          <  2:3214631/91‑1 (MQ=255)

                                   aaaCAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAAtt                                                         >  1:2880782/1‑91 (MQ=255)

                                   aaaCAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAAtt                                                         <  1:2417940/91‑1 (MQ=255)

                                    aaCAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAAttt                                                        >  2:1089324/1‑91 (MQ=255)

                                    aaCAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAAttt                                                        <  2:3217460/91‑1 (MQ=255)

                                      cAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATttct                                                      <  1:1200488/91‑1 (MQ=255)

                                      cAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATttct                                                      <  1:716781/91‑1 (MQ=255)

                                       aaaTTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATttctt                                                     <  1:2102071/91‑1 (MQ=255)

                                       aaaTTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATttctt                                                     <  2:809017/91‑1 (MQ=255)

                                         aTTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTgcttct                                                   >  2:3650672/1‑91 (MQ=255)

                                         aTTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGAAAAAAAAAAAAAAACATTAATTtcttct                                                   <  2:1958876/91‑1 (MQ=255)

                                          ttAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTa                                                  >  2:3934037/1‑91 (MQ=255)

                                           tAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTaa                                                 >  2:722563/1‑91 (MQ=255)

                                            aaTCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTaaa                                                <  2:3303875/91‑1 (MQ=255)

                                            aaTCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTaaa                                                >  2:3052680/1‑91 (MQ=255)

                                             aTCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGCCTCTAAAGAAAATACATTAATTTCTTCTCAAg                                               >  1:3284563/1‑91 (MQ=255)

                                             aTCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAg                                               >  1:1448016/1‑91 (MQ=255)

                                                 ttCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTgg                                           <  1:1498130/91‑1 (MQ=255)

                                                  tCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGga                                          <  1:554287/91‑1 (MQ=255)

                                                      gTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACg                                      <  2:2782430/91‑1 (MQ=255)

                                                      gTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACg                                      <  2:3464276/91‑1 (MQ=255)

                                                      gTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACg                                      <  1:2591817/91‑1 (MQ=255)

                                                        aaCAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTa                                    >  1:1176193/1‑91 (MQ=255)

                                                         aCAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTaa                                   >  1:3749459/1‑91 (MQ=255)

                                                         aCAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTaa                                   >  2:688768/1‑91 (MQ=255)

                                                         aCAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCAAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTaa                                   <  2:112665/91‑1 (MQ=255)

                                                          cAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAAc                                  >  2:59105/1‑91 (MQ=255)

                                                          cAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAAc                                  >  1:3184897/1‑91 (MQ=255)

                                                             cTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGt                               >  1:2270232/1‑91 (MQ=255)

                                                             cTAGTGCTTTATTATTAGGCGCTTGGGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGt                               <  1:3690585/91‑1 (MQ=255)

                                                              gaGTGCTTTATTATTAGGATCTTGTGCCGCTAGTGCCGCAGACTCTAAAGAAAATACATTAATTTCATCGAAAGCTGGAGACGTAACAGtt                              <  2:2141109/90‑1 (MQ=255)

                                                               aGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTg                             >  2:3493710/1‑91 (MQ=255)

                                                                gTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGc                            <  1:1818578/91‑1 (MQ=255)

                                                                  gCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAg                          >  2:3847014/1‑91 (MQ=255)

                                                                    tttATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGAt                        <  2:1384553/91‑1 (MQ=255)

                                                                      tattattaGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATAc                      >  2:1622886/1‑91 (MQ=255)

                                                                      tattattaGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATAc                      <  2:1777912/91‑1 (MQ=255)

                                                                       attattaGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACa                     >  1:876017/1‑91 (MQ=255)

                                                                         tattaGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAAt                   <  1:2748557/91‑1 (MQ=255)

                                                                         tattaGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAAt                   <  2:2909143/91‑1 (MQ=255)

                                                                         tattaGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAAt                   <  1:166372/91‑1 (MQ=255)

                                                                         tattaGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAAt                   >  2:1628243/1‑91 (MQ=255)

                                                                         tattaGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATAAAAt                   <  2:2436403/91‑1 (MQ=255)

                                                                          attaGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATg                  >  2:1278890/1‑91 (MQ=255)

                                                                          attaGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATg                  <  2:2944312/91‑1 (MQ=255)

                                                                           ttaGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGa                 <  2:2586062/91‑1 (MQ=255)

                                                                           ttaGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGa                 <  2:2708378/91‑1 (MQ=255)

                                                                            taGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGaa                >  1:3782415/1‑91 (MQ=255)

                                                                            taGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGaa                >  1:2258185/1‑91 (MQ=255)

                                                                            taGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGaa                >  2:2539485/1‑91 (MQ=255)

                                                                             aGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGaaa               <  1:2362258/91‑1 (MQ=255)

                                                                             aGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGaaa               <  1:4669/91‑1 (MQ=255)

                                                                              ggCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGaaaa              >  1:1770349/1‑91 (MQ=255)

                                                                                cgcTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGaaaaaa            <  2:2045826/91‑1 (MQ=255)

                                                                                 gcTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGaaaaaaa           >  1:3707851/1‑91 (MQ=255)

                                                                                   ttGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATc         <  2:1774804/91‑1 (MQ=255)

                                                                                   ttGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAAAAAAAAAAAAAAAAATc         <  2:3364330/91‑1 (MQ=255)

                                                                                      tgGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCGGt      >  1:1469880/1‑91 (MQ=255)

                                                                                       ggCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCGGTa     >  1:1793977/1‑91 (MQ=255)

                                                                                        gcgcTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCGGTaa    >  2:1044947/1‑91 (MQ=255)

                                                                                        gcgcTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCGGTaa    >  1:3029747/1‑91 (MQ=255)

                                                                                          gcTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCGGTAAAg  >  2:2271420/1‑91 (MQ=255)

TTAACTTACAAAGGAGTTTCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCGGTAAAG  >  USA300TCH1516_ALE/1972282‑1972462

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