BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,972,372	G→A	6.5%	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	6.5%	709.9 / 24.2	217	A23T (GCT→ACT)	prsA	Foldase protein PrsA
Reads supporting (aligned to +/- strand): ref base G (115/87); new base A (5/9); total (120/97)
Fisher's exact test for biased strand distribution p-value = 1.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

Only 100 of 229 total aligned reads displayed.

TTAACTTACAAAGGAGTTTCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCGGTAAAG  >  USA300TCH1516_ALE/1972282‑1972462

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

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

  aaCTTACAAAGGAGTTTCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCt                                                                                          >  1:79363/1‑91 (MQ=255)

  aaCTTACAAAGGAGTTGCCACTATGAAGATGACAAACAAACTAAGCGTTCCGGTAACAGGTAGAGCGTTATGATTAGGCGCTCGTCgcgca                                                                                          >  2:1060649/1‑90 (MQ=255)

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

    cTTACAAAGGAGTTTCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAg                                                                                        <  1:1691273/91‑1 (MQ=255)

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

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

    cTTACAAAGGAGTTCCAACTATGAAGATGATAAACAAATTAATCCTTCCGGTAACAGCTAGGGCTTTATTATTAGGCGCTTGTGGCACTAg                                                                                        <  1:3084484/91‑1 (MQ=255)

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

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

      tACAAAGGAGTTTCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCt                                                                                          >  2:1257392/1‑87 (MQ=255)

      tACAAAGGAGTTTCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCt                                                                                          <  1:1257392/87‑1 (MQ=255)

        cAAAGGAGTTTCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCACTAGCGcc                                                                                    <  1:2183875/91‑1 (MQ=255)

           aGGAGTTTCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCcaca                                                                                 <  1:3167417/91‑1 (MQ=255)

                tttCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGActc                                                                            <  1:3629345/91‑1 (MQ=255)

                 ttCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGActct                                                                           >  1:2380979/1‑91 (MQ=255)

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

                  tCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTa                                                                          <  2:3725310/91‑1 (MQ=255)

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

                    aaCTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTaaa                                                                        <  2:502146/91‑1 (MQ=255)

                    aaCTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTaaa                                                                        >  1:64807/1‑91 (MQ=255)

                     aCTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAg                                                                       >  2:2531126/1‑91 (MQ=255)

                       gaggAAGATGATAAACTAAGTAATCGTTCCGGGAACAGCTAGTGCTTTATTATTATTCTCTTTTGGCGCTAGTGCCACAGACTCTAAAGaa                                                                     <  1:2524050/88‑1 (MQ=255)

                         tGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGaaaa                                                                   <  2:698007/91‑1 (MQ=255)

                         tGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGCCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGaaaa                                                                   >  2:3503328/1‑91 (MQ=255)

                          gAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAAt                                                                  <  1:1025356/91‑1 (MQ=255)

                          gAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAAt                                                                  >  1:1798958/1‑91 (MQ=255)

                           aaGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATa                                                                 >  2:1281568/1‑91 (MQ=255)

                           aaGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATa                                                                 >  1:266306/1‑91 (MQ=255)

                           aaGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATa                                                                 >  1:2038232/1‑91 (MQ=255)

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

                             gatgatAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGTCTCTAAAGAAAATACa                                                               >  2:2973872/1‑91 (MQ=255)

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

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

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

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

                                 atAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTaa                                                           >  2:1607816/1‑91 (MQ=255)

                                 atAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTaa                                                           <  2:1578789/91‑1 (MQ=255)

                                 atAAACAAATTAAACGTTCCGGTAACAGCTAGTGCTTTATTATGAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTaa                                                           <  1:3874374/91‑1 (MQ=255)

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

                                     aCAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCTCTAGTGCCACAGACTCTAAAGAAAATACATTAATttc                                                       <  1:3259259/91‑1 (MQ=255)

                                      cAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATttct                                                      >  1:3295123/1‑91 (MQ=255)

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

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

                                         aTTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGTTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTtcttct                                                   <  2:3417615/91‑1 (MQ=255)

                                         aTTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTtcttct                                                   >  2:2871358/1‑91 (MQ=255)

                                         aTTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTtcttct                                                   <  2:2598512/91‑1 (MQ=255)

                                          ttAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTa                                                  <  2:60950/91‑1 (MQ=255)

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

                                             aTCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAg                                               <  1:37824/91‑1 (MQ=255)

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

                                              tCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGc                                              >  1:1786550/1‑91 (MQ=255)

                                               cGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCt                                             <  2:1517972/91‑1 (MQ=255)

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

                                                    cGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGgag                                         >  2:2023370/1‑90 (MQ=255)

                                                     ggTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGAc                                       >  1:969037/1‑91 (MQ=255)

                                                     ggTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGAc                                       >  2:683995/1‑91 (MQ=255)

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

                                                       tAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGt                                     <  2:2333885/91‑1 (MQ=255)

                                                       tAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGt                                     <  2:1796986/91‑1 (MQ=255)

                                                       tAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGt                                     <  2:2811438/91‑1 (MQ=255)

                                                        aaCAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTa                                    <  1:3521964/91‑1 (MQ=255)

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

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

                                                           aGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACa                                 <  1:2424860/91‑1 (MQ=255)

                                                           aGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACa                                 <  2:2049435/91‑1 (MQ=255)

                                                             cTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGt                               <  1:2492730/91‑1 (MQ=255)

                                                             cTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGt                               <  1:313631/91‑1 (MQ=255)

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

                                                               aGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTg                             <  2:816641/91‑1 (MQ=255)

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

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

                                                                 tGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAACTACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCa                           <  2:3183851/91‑1 (MQ=255)

                                                                  gCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAg                          <  1:3224200/91‑1 (MQ=255)

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

                                                                   cTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGa                         >  2:2071330/1‑91 (MQ=255)

                                                                   cTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGa                         <  1:292651/91‑1 (MQ=255)

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

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

                                                                            taGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGaa                <  2:1361514/91‑1 (MQ=255)

                                                                             aGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGaaa               >  1:2075888/1‑91 (MQ=255)

                                                                              ggCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGaaaa              <  1:1174829/91‑1 (MQ=255)

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

                                                                               gcgcTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGaaaaa             <  1:3494664/91‑1 (MQ=255)

                                                                                cgcTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGaaaaaa            >  2:2822467/1‑91 (MQ=255)

                                                                                 gcTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGaaaaaaa           <  1:2777615/91‑1 (MQ=255)

                                                                                  cTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAAt          >  2:3065915/1‑91 (MQ=255)

                                                                                   ttGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATc         >  1:3156588/1‑91 (MQ=255)

                                                                                    tgtgGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCg        >  2:1900443/1‑91 (MQ=255)

                                                                                    tgtgGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCg        <  1:1904749/91‑1 (MQ=255)

                                                                                    tgtgGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCg        <  1:800286/91‑1 (MQ=255)

                                                                                    tgtgGCACTAGCGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCg        >  1:2854270/1‑91 (MQ=255)

                                                                                    tgtgGCACTAGCGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCg        >  2:3867390/1‑91 (MQ=255)

                                                                                     gtgGCACTAGCGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCgg       >  2:3046753/1‑91 (MQ=255)

                                                                                     gtgGCACTAGCGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCgg       >  1:115216/1‑91 (MQ=255)

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

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

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

                                                                                          gcTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCGGTAAAg  <  1:214349/91‑1 (MQ=255)

TTAACTTACAAAGGAGTTTCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCGGTAAAG  >  USA300TCH1516_ALE/1972282‑1972462

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