BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	722,992	A→G	17.0%	I44V (ATA→GTA)	feuB →	Iron‑uptake system permease protein FeuB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	722,992	0	A	G	17.0%	315.9 / 54.2	129	I44V (ATA→GTA)	feuB	Iron‑uptake system permease protein FeuB
Reads supporting (aligned to +/- strand): ref base A (48/59); new base G (17/5); total (65/64)
Fisher's exact test for biased strand distribution p-value = 8.93e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 131 total aligned reads displayed.

TCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATAGCAGCAGTAA  >  USA300TCH1516_ALE/722902‑723079

tCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATta                                                                                         >  1:2187967/1‑91 (MQ=255)

   tataTTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATag                                                                                      <  1:1068627/91‑1 (MQ=255)

   tataTTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATag                                                                                      <  1:3723812/91‑1 (MQ=255)

   tataTTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATag                                                                                      <  1:4050987/91‑1 (MQ=255)

    atatTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAATTGTaga                                                                                     <  1:1524201/91‑1 (MQ=255)

    atatTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAATTGTaga                                                                                     <  1:2214291/91‑1 (MQ=255)

    atatTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAATTGTaga                                                                                     <  1:2212094/91‑1 (MQ=255)

     tattatAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATagag                                                                                    >  1:2997564/1‑91 (MQ=255)

      attatAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAATTGTAGAgg                                                                                   <  1:4095201/91‑1 (MQ=255)

        tataGGTTTTAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCt                                                                                 <  1:3722443/91‑1 (MQ=255)

        tataGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAATTGTAGAGGCt                                                                                 <  1:3056313/91‑1 (MQ=255)

        tataGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCt                                                                                 <  1:2311754/91‑1 (MQ=255)

         ataGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTa                                                                                <  1:3654602/91‑1 (MQ=255)

            ggTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTAttt                                                                             >  1:1206431/1‑91 (MQ=255)

            ggTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTAttt                                                                             >  1:3950486/1‑91 (MQ=255)

            ggTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTAttt                                                                             >  1:1977750/1‑91 (MQ=255)

             gTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTAtttt                                                                            >  1:2471747/1‑91 (MQ=255)

             gTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTAttgt                                                                            >  1:2471770/1‑89 (MQ=255)

             gTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAACAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTAtttt                                                                            <  1:3904831/91‑1 (MQ=255)

                tAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTaa                                                                         >  1:3669548/1‑91 (MQ=255)

                tAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTaa                                                                         >  1:3669530/1‑91 (MQ=255)

                tAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTaa                                                                         >  1:3670837/1‑91 (MQ=255)

                tAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTaa                                                                         <  1:544637/91‑1 (MQ=255)

                tAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTaa                                                                         <  1:544903/91‑1 (MQ=255)

                  aGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAAtt                                                                       <  1:3474233/91‑1 (MQ=255)

                  aGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAAtt                                                                       <  1:1961246/91‑1 (MQ=255)

                     aTGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATtata                                                                    >  1:94209/1‑91 (MQ=255)

                      tGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATaa                                                                   >  1:2605223/1‑91 (MQ=255)

                       gATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAAt                                                                  >  1:780885/1‑91 (MQ=255)

                       gATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAAt                                                                  >  1:605594/1‑91 (MQ=255)

                           ctactaATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTa                                                              >  1:3108259/1‑91 (MQ=255)

                             actaATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGc                                                            >  1:1223124/1‑91 (MQ=255)

                              ctaATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCa                                                           >  1:4004359/1‑91 (MQ=255)

                                 aTCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATc                                                        >  1:1573286/1‑91 (MQ=255)

                                  tCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCa                                                       <  1:2118742/91‑1 (MQ=255)

                                  tCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCa                                                       >  1:514394/1‑91 (MQ=255)

                                  tCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCa                                                       >  1:415980/1‑91 (MQ=255)

                                   cacaCTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCaa                                                      <  1:2722073/91‑1 (MQ=255)

                                     cacTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACa                                                    <  1:2142210/91‑1 (MQ=255)

                                      acTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAg                                                   <  1:2739560/91‑1 (MQ=255)

                                       cTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGc                                                  <  1:3894535/91‑1 (MQ=255)

                                        tatttattTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCa                                                 <  1:3641660/91‑1 (MQ=255)

                                         atttatttCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCaa                                                <  1:3614078/91‑1 (MQ=255)

                                         atttatttCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCaa                                                <  1:1022090/91‑1 (MQ=255)

                                          tttatttCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCaaa                                               >  1:3011349/1‑91 (MQ=255)

                                             atttCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAAca                                            >  1:2674953/1‑91 (MQ=255)

                                                tCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAAcatca                                         <  1:854797/91‑1 (MQ=255)

                                                   aCATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACAtcatca                                      >  1:3257152/1‑91 (MQ=255)

                                                     aTTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAAt                                    <  1:1589197/91‑1 (MQ=255)

                                                     aTTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAAt                                    <  1:1588852/91‑1 (MQ=255)

                                                     aTTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAAt                                    <  1:1588837/91‑1 (MQ=255)

                                                       tAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATga                                  <  1:3813793/91‑1 (MQ=255)

                                                           aGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGAtt                              <  1:2971244/91‑1 (MQ=255)

                                                           aGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGAtt                              <  1:1356985/91‑1 (MQ=255)

                                                           aGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGAtt                              >  1:3430064/1‑91 (MQ=255)

                                                            ggTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTa                             >  1:3816363/1‑91 (MQ=255)

                                                              tGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAgg                           >  1:3223616/1‑91 (MQ=255)

                                                               gACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGa                          <  1:2134699/91‑1 (MQ=255)

                                                                aCGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGAt                         >  1:3591744/1‑91 (MQ=255)

                                                                aCGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGAt                         <  1:478139/91‑1 (MQ=255)

                                                                 cGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGAtt                        <  1:3153898/91‑1 (MQ=255)

                                                                     aaaaTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCa                    <  1:760785/91‑1 (MQ=255)

                                                                      aaaTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAg                   >  1:2485333/1‑91 (MQ=255)

                                                                      aaaTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAg                   >  1:2485695/1‑91 (MQ=255)

                                                                      aaaTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAg                   <  1:443170/91‑1 (MQ=255)

                                                                      aaaTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAg                   <  1:103323/91‑1 (MQ=255)

                                                                       aaTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGa                  >  1:656995/1‑91 (MQ=255)

                                                                       aaTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGa                  <  1:1451935/91‑1 (MQ=255)

                                                                       aaTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGa                  >  1:1609163/1‑91 (MQ=255)

                                                                       aaTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGa                  >  1:1608987/1‑91 (MQ=255)

                                                                        aTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGaa                 <  1:1534489/91‑1 (MQ=255)

                                                                           cAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAata              >  1:2281139/1‑91 (MQ=255)

                                                                            aaGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAatat             <  1:4064041/91‑1 (MQ=255)

                                                                            aaGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAatat             <  1:4063541/91‑1 (MQ=255)

                                                                             aGCCTCTATAATTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAAtata            >  1:2943787/1‑91 (MQ=255)

                                                                              gCCTCTATAATTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATag           >  1:1121480/1‑91 (MQ=255)

                                                                              gCCTCTATAATTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATag           >  1:1830515/1‑91 (MQ=255)

                                                                              gCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATag           >  1:398674/1‑91 (MQ=255)

                                                                              gCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATag           >  1:398782/1‑91 (MQ=255)

                                                                               ccTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagc          >  1:940797/1‑91 (MQ=255)

                                                                               ccTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagc          <  1:3357659/91‑1 (MQ=255)

                                                                               ccTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagc          <  1:1189519/91‑1 (MQ=255)

                                                                                ctctATAATTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagca         >  1:59010/1‑91 (MQ=255)

                                                                                ctctATAATTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagca         >  1:600054/1‑91 (MQ=255)

                                                                                 tctATAATTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagcag        >  1:2380271/1‑91 (MQ=255)

                                                                                 tctATAATTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagcag        >  1:2638698/1‑91 (MQ=255)

                                                                                 tctATAATTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagcag        >  1:1343654/1‑91 (MQ=255)

                                                                                 tctATAATTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagcag        >  1:2638603/1‑91 (MQ=255)

                                                                                 tctACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagcag        <  1:2396670/91‑1 (MQ=255)

                                                                                  ctataaTTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagcagc       >  1:3235340/5‑91 (MQ=255)

                                                                                  ctataaTTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagcagc       >  1:3364088/5‑91 (MQ=255)

                                                                                  ctACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagcagc       <  1:1635603/91‑1 (MQ=255)

                                                                                  ctACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagcagc       <  1:159744/91‑1 (MQ=255)

                                                                                    ataaTTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATAGcagcag     >  1:2949414/3‑91 (MQ=255)

                                                                                    aCAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATAGcagcag     <  1:2682296/91‑1 (MQ=255)

                                                                                    aCAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATAGcagcag     <  1:2682855/91‑1 (MQ=255)

                                                                                     cAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATAGCAGCAGt    >  1:3474544/1‑91 (MQ=255)

                                                                                      aaTTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATAGCAGCAGTa   >  1:3244530/1‑91 (MQ=255)

                                                                                      aaTTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATAGCAGCAGTa   <  1:2201236/91‑1 (MQ=255)

                                                                                       aTTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATAGCAGCAGTaa  >  1:2651436/1‑91 (MQ=255)

TCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATAGCAGCAGTAA  >  USA300TCH1516_ALE/722902‑723079

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