BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,704,428	G→C	6.1%	Q81E (CAA→GAA)	crtM ←	Dehydrosqualene synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,704,428	0	G	C	6.1%	752.5 / 28.4	212	Q81E (CAA→GAA)	crtM	Dehydrosqualene synthase
Reads supporting (aligned to +/- strand): ref base G (105/94); new base C (6/7); total (111/101)
Fisher's exact test for biased strand distribution p-value = 7.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 226 total aligned reads displayed.

CAGTATCAATGAGATTATAAAAAGATTGAAAGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTGATTTAAAAATTGAATATCGCCATAAACATCTAT  >  USA300TCH1516_ALE/2704338‑2704518

cAGTATCAATGAGATTATAAAAAGATTGAAAGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTg                                                                                            >  1:3381554/1‑91 (MQ=255)

cAGTATCAATGAGATTATAAAAAGATTGAAAGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTg                                                                                            <  1:580838/91‑1 (MQ=255)

cAGTATCAATGAGATTATAAAAAGATTGAAAGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTg                                                                                            <  2:702005/91‑1 (MQ=255)

cAGTATCAATGAGATTATAAAAAGATTGAAAGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTg                                                                                            >  1:1442730/1‑91 (MQ=255)

cAGTATCAATGAGATTATAAAAAGATTGAAAGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTg                                                                                            <  1:3201228/91‑1 (MQ=255)

cAGTATCAATGAGATTATAAAAAGATTGAAAGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTg                                                                                            <  1:3774215/91‑1 (MQ=255)

  gTATCAATGAGATTATAAAAAGATTGAAAGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGaa                                                                                          <  2:376251/91‑1 (MQ=255)

    aTCAATGAGATTATAAAAAGATTGAAAGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAg                                                                                        <  2:2413479/91‑1 (MQ=255)

    aTCAATGAGATTATAAAAAGATTGAAAGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAg                                                                                        >  2:3706446/1‑91 (MQ=255)

    aTCAATGAGATTATAAAAAGATTGAAAGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAg                                                                                        <  2:352221/91‑1 (MQ=255)

    aTCAATGAGATTATAAAAAGATTGAAAGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAg                                                                                        <  1:2326669/91‑1 (MQ=255)

     tCAATGAGATTATAAAAAGATTGAAAGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGt                                                                                       <  1:394728/91‑1 (MQ=255)

     tCAATGAGATTATAAAAAGATTGAAAGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTCAAAGt                                                                                       <  2:2624314/91‑1 (MQ=255)

         tGAGATTATAAAAAGATTGAAAGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGtgatg                                                                                   >  1:3594174/1‑91 (MQ=255)

         tGAGATTATAAAAAGATTGAAAGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGtgatg                                                                                   >  2:905499/1‑91 (MQ=255)

          gagaTTATAAAAAGATTGAAAGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGtgatga                                                                                  >  1:2169731/1‑91 (MQ=255)

           agaTTATAAAAAGATTGAAAGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTgatgat                                                                                 >  1:2537812/1‑91 (MQ=255)

            gaTTATAAAAAGATTGAAAGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGatgatg                                                                                >  2:3296948/1‑91 (MQ=255)

                  aaaaaGATTGAAAGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCata                                                                          >  1:1015508/1‑91 (MQ=255)

                   aaaaGATTGAAAGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCatat                                                                         <  1:217644/91‑1 (MQ=255)

                    aaaGATTGAAAGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATg                                                                        <  2:2464173/91‑1 (MQ=255)

                     aaGATTGAAATGAGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATgg                                                                       <  1:1824892/91‑1 (MQ=255)

                     aaGATTGAAAGGCGATATTTTTATGTTGTGCGACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATgg                                                                       >  2:3357507/1‑91 (MQ=255)

                     aaGATTGAAAGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATgg                                                                       >  1:555427/1‑91 (MQ=255)

                      aGATTGAAAGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATggg                                                                      <  1:3862679/91‑1 (MQ=255)

                          tGAAAGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTAtt                                                                  >  1:3156253/1‑91 (MQ=255)

                          tGAAAGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTAtt                                                                  >  1:222780/1‑91 (MQ=255)

                           gAAAGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTAttt                                                                 >  1:2281446/1‑91 (MQ=255)

                           gAAAGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTAttt                                                                 >  1:1134810/1‑91 (MQ=255)

                             aaGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTAttttt                                                               >  2:1118937/1‑91 (MQ=255)

                             aaGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTAttttt                                                               >  2:1678754/1‑91 (MQ=255)

                             aaGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTAttttt                                                               <  2:3096129/91‑1 (MQ=255)

                             aaGGCGATAGTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTGGAAAGTGATGATGTTCATATGGGTAttttt                                                               <  2:928641/91‑1 (MQ=255)

                              aGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTc                                                              >  1:3527213/1‑91 (MQ=255)

                              aGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTc                                                              <  1:3469148/91‑1 (MQ=255)

                              aGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTc                                                              >  2:604120/1‑91 (MQ=255)

                              aGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTc                                                              <  1:1666547/91‑1 (MQ=255)

                               ggCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCGTATGGGTATTTTTCa                                                             <  2:2768591/91‑1 (MQ=255)

                                  gATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATa                                                          >  2:2233792/1‑91 (MQ=255)

                                    tatTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGa                                                        <  1:1136973/91‑1 (MQ=255)

                                    tatTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGa                                                        <  1:1139896/91‑1 (MQ=255)

                                    tatTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGa                                                        >  2:1818849/1‑91 (MQ=255)

                                    tatTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGa                                                        <  1:671655/91‑1 (MQ=255)

                                    tatTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGa                                                        >  1:3033729/1‑91 (MQ=255)

                                     atTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGAt                                                       <  2:2692787/91‑1 (MQ=255)

                                       ttttATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTg                                                     >  2:2850253/1‑91 (MQ=255)

                                       ttttATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTg                                                     <  2:1185644/91‑1 (MQ=255)

                                        tttATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGt                                                    >  1:2784722/1‑91 (MQ=255)

                                         ttATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGta                                                   <  1:2705506/91‑1 (MQ=255)

                                          tATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTACTTTTCAATAGATTGtat                                                  >  1:3319867/1‑91 (MQ=255)

                                           aTGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGtata                                                 <  1:990896/91‑1 (MQ=255)

                                             gttgtGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATAtc                                               >  1:1859606/1‑91 (MQ=255)

                                              ttgtGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATAtct                                              >  1:3729723/1‑91 (MQ=255)

                                               tgtgCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATAtctt                                             >  1:1742780/1‑91 (MQ=255)

                                               ggtgCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATAtctt                                             <  2:3849215/90‑1 (MQ=255)

                                                 tgCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATAtcttct                                           >  2:3402039/1‑91 (MQ=255)

                                                 tgCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGAGATGATGTTCATATGGGTATTTTTCAATAGATTGTAAAtcttct                                           >  2:3487492/1‑91 (MQ=255)

                                                 tgCAACATGCTGAAGCGCCATCATGACTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATAtcttct                                           <  1:2094653/91‑1 (MQ=255)

                                                  gCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATcttctt                                          <  2:3065736/91‑1 (MQ=255)

                                                  gCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATcttctt                                          <  2:2156693/91‑1 (MQ=255)

                                                  gCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATcttctt                                          <  2:220346/91‑1 (MQ=255)

                                                  gCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATcttctt                                          <  2:218649/91‑1 (MQ=255)

                                                  acAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATcttctt                                          >  2:1481198/2‑91 (MQ=255)

                                                    aaCATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCtttt                                        <  2:543155/91‑1 (MQ=255)

                                                       aTGCGGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTAtt                                     >  1:3065564/1‑91 (MQ=255)

                                                         gCTGAAGCGCCATCATGATTCTATGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTg                                   >  1:2159323/1‑91 (MQ=255)

                                                         gCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTg                                   >  2:2159179/1‑91 (MQ=255)

                                                         gCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTg                                   >  1:3645990/1‑91 (MQ=255)

                                                         gCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTg                                   >  2:1549310/1‑91 (MQ=255)

                                                          cTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTGa                                  <  1:427248/91‑1 (MQ=255)

                                                          cTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTGa                                  <  1:934067/91‑1 (MQ=255)

                                                           tGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTGAt                                 <  2:3594174/91‑1 (MQ=255)

                                                            gAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTGAtt                                >  2:406818/1‑91 (MQ=255)

                                                            gAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTGAtt                                <  1:566249/91‑1 (MQ=255)

                                                             aaGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTGAttt                               >  2:105892/1‑91 (MQ=255)

                                                             aaGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTGAttt                               >  2:2309842/1‑91 (MQ=255)

                                                               gcgcCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTGATTTaa                             >  2:479719/1‑91 (MQ=255)

                                                                  ccATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAAt                                                           >  1:2349950/1‑58 (MQ=255)

                                                                  ccATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAAt                                                           <  2:2349950/58‑1 (MQ=255)

                                                                  ccATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTGATTTaaaaa                          >  1:2824488/1‑91 (MQ=255)

                                                                  ccATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTGATTTaaaaa                          >  1:899826/1‑91 (MQ=255)

                                                                  ccATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTGATTTaaaaa                          >  2:3085956/1‑91 (MQ=255)

                                                                   catcatGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTGATTTAAAAAt                         <  2:1771401/91‑1 (MQ=255)

                                                                    atcatGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTGATTTAAAAAtt                        >  2:815820/1‑91 (MQ=255)

                                                                    atcatGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTGATTTAAAAAtt                        <  2:152323/91‑1 (MQ=255)

                                                                     tcatGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTGATTTAAAAATTg                       <  1:1495414/91‑1 (MQ=255)

                                                                        tGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTGATTTAAAAATTGAat                    >  1:1471931/1‑91 (MQ=255)

                                                                           ttCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTGATTTAAAAATTGAATATc                 <  1:482178/91‑1 (MQ=255)

                                                                            tCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTTTTTTATTTGATTTAAAAATTGAATATCg                >  2:2017274/1‑91 (MQ=255)

                                                                                cGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTGATTTAAAAATTGAATATCGCCAt            <  2:3688685/91‑1 (MQ=255)

                                                                                cGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTGATTTAAAAATTGAATATCGCCAt            >  2:1294655/1‑91 (MQ=255)

                                                                                cGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTGATTTAAAAATTGAATATCGCCAt            >  1:2392308/1‑91 (MQ=255)

                                                                                  aTCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTGATTTAAAAATTGAATATCGCCATaa          <  1:3253292/91‑1 (MQ=255)

                                                                                  aTCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTGATTTAAAAATTGAATATCGCCATaa          <  1:2559951/91‑1 (MQ=255)

                                                                                   tCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTGATTTAAAAATTGAATATCGCCATaaa         >  2:1944999/1‑91 (MQ=255)

                                                                                    cACTTTCAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTGATTTAAAAATTGAATATCGCCATAAAc        >  1:1249846/1‑91 (MQ=255)

                                                                                        ttGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTGATTTAAAAATTGAATATCGCCATAAACATCt    >  2:3162739/1‑91 (MQ=255)

                                                                                         tGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTGATTTAAAAATTGAATATCGCCATAAACATCta   <  2:1442730/91‑1 (MQ=255)

                                                                                         tGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTGATTTAAAAATTGAATATCGCCATAAACATCta   >  2:2564545/1‑91 (MQ=255)

                                                                                          gAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTGATTTAAAAATTGAATATCGCCATAAACATCtat  >  2:3072762/1‑91 (MQ=255)

CAGTATCAATGAGATTATAAAAAGATTGAAAGGCGATATTTTTATGTTGTGCAACATGCTGAAGCGCCATCATGATTCTACGATCACTTTGAAAGTGATGATGTTCATATGGGTATTTTTCAATAGATTGTATATCTTCTTTTATTTGATTTAAAAATTGAATATCGCCATAAACATCTAT  >  USA300TCH1516_ALE/2704338‑2704518

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