BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,833,387	G→C	9.9%	intergenic (‑723/+367)	lipA_2 ← / ← hisI	Lipase 1/Phosphoribosyl‑AMP cyclohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,833,387	0	G	C	9.9%	439.7 / 34.9	131	intergenic (‑723/+367)	lipA_2/hisI	Lipase 1/Phosphoribosyl‑AMP cyclohydrolase
Reads supporting (aligned to +/- strand): ref base G (59/59); new base C (3/10); total (62/69)
Fisher's exact test for biased strand distribution p-value = 8.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

Only 100 of 133 total aligned reads displayed.

CGGGATGCATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGCCATTTTTTCCCTT  >  USA300TCH1516_ALE/2833300‑2833477

cGGGATGCATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGtat                                                                                         >  1:3078703/1‑91 (MQ=255)

cGGGATGCATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGtat                                                                                         >  1:3947497/1‑91 (MQ=255)

cGGGATGCATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATACtat                                                                                         <  1:694153/91‑1 (MQ=255)

cGGGATGCATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATACtat                                                                                         <  1:1193652/91‑1 (MQ=255)

  ggatgCATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATaa                                                                                       <  1:4120878‑M2/86‑1 (MQ=255)

  ggatgCATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATaa                                                                                       <  1:4121938‑M2/86‑1 (MQ=255)

  ggatgCATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATaa                                                                                       >  1:413113‑M2/6‑91 (MQ=255)

  ggatgCATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATaa                                                                                       >  1:3680529‑M2/6‑91 (MQ=255)

   gatgCATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATACTATaaa                                                                                      <  1:1601141‑M2/87‑1 (MQ=255)

   gatgCATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATACTATaaa                                                                                      <  1:1601057‑M2/87‑1 (MQ=255)

    atgCATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAAt                                                                                     <  1:1079872‑M2/88‑1 (MQ=255)

    atgCATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATACTATAAAt                                                                                     <  1:685773‑M2/88‑1 (MQ=255)

      gCATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATAc                                                                                   >  1:1297589‑M2/2‑91 (MQ=255)

       cATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACa                                                                                  <  1:3653189/91‑1 (MQ=255)

       cATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACa                                                                                  <  1:556130/91‑1 (MQ=255)

        atatCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACaa                                                                                 >  1:2072423/1‑91 (MQ=255)

        atatCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATACTATAAATACaa                                                                                 <  1:584322/91‑1 (MQ=255)

        atatCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATACTATAAATACaa                                                                                 <  1:1182721/91‑1 (MQ=255)

         tatCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATACTATAAATACaaa                                                                                <  1:1350754/91‑1 (MQ=255)

           tCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGa                                                                              <  1:1381072/91‑1 (MQ=255)

            cTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGaa                                                                             <  1:1989801/91‑1 (MQ=255)

              tAGTTCTAATGGGGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATg                                                                           <  1:1565803/91‑1 (MQ=255)

              tAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATg                                                                           <  1:3821342/91‑1 (MQ=255)

              tAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATg                                                                           >  1:1609008/1‑91 (MQ=255)

              tAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATg                                                                           >  1:1608954/1‑91 (MQ=255)

              tAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATg                                                                           <  1:91828/91‑1 (MQ=255)

               aGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGt                                                                          >  1:4019000/1‑91 (MQ=255)

               aGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGt                                                                          >  1:3057582/1‑91 (MQ=255)

                 ttCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTaa                                                                        >  1:3776937/1‑91 (MQ=255)

                 ttCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTaa                                                                        >  1:1396753/1‑91 (MQ=255)

                 ttCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTaa                                                                        >  1:1396335/1‑91 (MQ=255)

                   cTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAg                                                                      >  1:3246270/1‑91 (MQ=255)

                     aaTGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGAc                                                                    >  1:3958929/1‑91 (MQ=255)

                     aaTGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGAc                                                                    >  1:3149512/1‑91 (MQ=255)

                     aaTGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGAc                                                                    <  1:282423/91‑1 (MQ=255)

                     aaTGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGAc                                                                    >  1:1170078/1‑91 (MQ=255)

                     aaTGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGAc                                                                    >  1:2317098/1‑91 (MQ=255)

                      aTGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACa                                                                   >  1:1115207/1‑91 (MQ=255)

                      aTGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACa                                                                   >  1:2088806/1‑91 (MQ=255)

                       tGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACaa                                                                  >  1:3469898/1‑91 (MQ=255)

                       tGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACaa                                                                  <  1:348330/91‑1 (MQ=255)

                        ggAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAAc                                                                 <  1:166390/91‑1 (MQ=255)

                          agaTTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTa                                                               >  1:3482346/1‑91 (MQ=255)

                          agaTTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTa                                                               <  1:1345642/91‑1 (MQ=255)

                          agaTTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTa                                                               >  1:2176234/1‑91 (MQ=255)

                          agaTTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTa                                                               >  1:2874625/1‑91 (MQ=255)

                          agaTTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTa                                                               >  1:2174895/1‑91 (MQ=255)

                          agaTTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTa                                                               <  1:2359235/91‑1 (MQ=255)

                           gaTTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTAt                                                              <  1:2504403/91‑1 (MQ=255)

                           gaTTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTAt                                                              <  1:2504332/91‑1 (MQ=255)

                            aTTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATg                                                             <  1:4089425/91‑1 (MQ=255)

                            aTTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATg                                                             <  1:1088949/91‑1 (MQ=255)

                             ttGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGa                                                            <  1:916432/91‑1 (MQ=255)

                              tGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGat                                                           <  1:2104148/91‑1 (MQ=255)

                              tGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGat                                                           <  1:45436/91‑1 (MQ=255)

                               gTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGata                                                          <  1:4063437/91‑1 (MQ=255)

                               gTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGata                                                          <  1:4063764/91‑1 (MQ=255)

                                tCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGatat                                                         <  1:2510864/91‑1 (MQ=255)

                                 cTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATAtt                                                        >  1:3241811/1‑91 (MQ=255)

                                 cTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATAtt                                                        >  1:2593837/1‑91 (MQ=255)

                                  ttttAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTa                                                       <  1:2031579/91‑1 (MQ=255)

                                   tttAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTaa                                                      >  1:3933498/1‑91 (MQ=255)

                                     tAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATa                                                    <  1:228225/91‑1 (MQ=255)

                                       aaTAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAgg                                                  >  1:590257/1‑91 (MQ=255)

                                       aaTAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAgg                                                  >  1:3437951/1‑91 (MQ=255)

                                       aaTAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAgg                                                  >  1:3437880/1‑91 (MQ=255)

                                          aGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGgtag                                               <  1:2334004/91‑1 (MQ=255)

                                              tATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAAc                                           <  1:33284/91‑1 (MQ=255)

                                                   cATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACaa                                      >  1:749169/1‑91 (MQ=255)

                                                    aTTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACaaa                                     >  1:3029549/1‑91 (MQ=255)

                                                      ttCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACa                                   <  1:2600073/91‑1 (MQ=255)

                                                       tCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACAt                                  <  1:3713603/91‑1 (MQ=255)

                                                       tCTATTCCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACAt                                  >  1:3225164/1‑91 (MQ=255)

                                                        cTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATa                                 <  1:1813418/91‑1 (MQ=255)

                                                          aTTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAAc                               <  1:832142/91‑1 (MQ=255)

                                                           tttCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAAct                              >  1:719450/1‑91 (MQ=255)

                                                               aGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCttt                          >  1:80401/1‑91 (MQ=255)

                                                               aGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCttt                          >  1:3823921/1‑91 (MQ=255)

                                                                 gCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAg                        >  1:501992/1‑91 (MQ=255)

                                                                   tAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCt                      <  1:1970602/91‑1 (MQ=255)

                                                                   tAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCt                      <  1:262271/91‑1 (MQ=255)

                                                                    aGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCtt                     >  1:2556789/1‑91 (MQ=255)

                                                                     gATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCttt                    <  1:843784/91‑1 (MQ=255)

                                                                        atatTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTc                 >  1:170961/1‑91 (MQ=255)

                                                                        atatTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTc                 >  1:3425159/1‑91 (MQ=255)

                                                                           tttGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGcc              <  1:4111014/91‑1 (MQ=255)

                                                                           tttGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGcc              >  1:4110566/1‑91 (MQ=255)

                                                                           tttGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGcc              <  1:311842/91‑1 (MQ=255)

                                                                            ttGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGCCa             >  1:4004899/1‑91 (MQ=255)

                                                                            ttGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGCCa             <  1:3396540/91‑1 (MQ=255)

                                                                               tATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGCCAttt          <  1:2704510/91‑1 (MQ=255)

                                                                                 tttATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGCCAttttt        <  1:1827833/91‑1 (MQ=255)

                                                                                 tttATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGCCAttttt        <  1:771552/91‑1 (MQ=255)

                                                                                 tttATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGCCAttttt        <  1:3277176/91‑1 (MQ=255)

                                                                                 tttATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGCCAttttt        <  1:932698/91‑1 (MQ=255)

                                                                                  ttATACTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGCCAtttttt       >  1:1068563/1‑91 (MQ=255)

                                                                                    ataGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGCCATTTTTTcc     >  1:2766957/1‑91 (MQ=255)

                                                                                      aGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGCCATTTTTTCCCt   <  1:3436687/91‑1 (MQ=255)

                                                                                       gTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGCCATTTTTTCCCtt  <  1:2857810/91‑1 (MQ=255)

                                                                                       gTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGCCATTTTTTCCCtt  >  1:751107/1‑91 (MQ=255)

CGGGATGCATATCTTAGTTCTAATGGAGATTGTCTTTTAAATAGGCTATGACATTTCTATTTCAGGCTAGATATATTTGTATTTATAGTATAAATACAAAGAATGTAAAGACAACTATGATATTAATAGGTAGTAACGACAAACATAACTCTTTAGCTTTTTCGCCATTTTTTCCCTT  >  USA300TCH1516_ALE/2833300‑2833477

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