BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	631,535	G→A	12.2%	intergenic (+161/‑341)	gph_1 → / → proP	Phosphoglycolate phosphatase/Proline/betaine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	631,535	0	G	A	12.2%	377.1 / 34.9	131	intergenic (+161/‑341)	gph_1/proP	Phosphoglycolate phosphatase/Proline/betaine transporter
Reads supporting (aligned to +/- strand): ref base G (56/59); new base A (8/8); total (64/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 131 total aligned reads displayed.

TCGCTTTAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCA  >  USA300TCH1516_ALE/631446‑631624

tCGCTTTAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGt                                                                                          >  1:4013779/1‑91 (MQ=255)

tCGCTTTAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGt                                                                                          >  1:586754/1‑91 (MQ=255)

 cGCTTTAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTa                                                                                         >  1:1905659/1‑91 (MQ=255)

  gCTTTAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTaa                                                                                        >  1:2263400/1‑91 (MQ=255)

      tAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAaaat                                                                                    >  1:906583/1‑91 (MQ=255)

      tAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTACAAaaat                                                                                    <  1:1737361/91‑1 (MQ=255)

       aGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAaaata                                                                                   >  1:3552725/1‑91 (MQ=255)

       aGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAaaata                                                                                   <  1:3415072/91‑1 (MQ=255)

       aGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTACAAaaata                                                                                   <  1:92409/91‑1 (MQ=255)

       aGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTACAAaaata                                                                                   <  1:2209248/91‑1 (MQ=255)

       aGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTACAAaaata                                                                                   <  1:2209085/91‑1 (MQ=255)

        gATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAaaataa                                                                                  >  1:2504153/1‑91 (MQ=255)

          ttAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAaaataaat                                                                                <  1:3738551/91‑1 (MQ=255)

              gAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTaaa                                                                            <  1:2916317/91‑1 (MQ=255)

              gAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTaaa                                                                            >  1:66689/1‑91 (MQ=255)

               aaTTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAg                                                                           <  1:159222/91‑1 (MQ=255)

               aaTTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAg                                                                           <  1:157842/91‑1 (MQ=255)

                aTTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGt                                                                          >  1:202507/1‑91 (MQ=255)

                aTTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGt                                                                          <  1:2213726/91‑1 (MQ=255)

                 tttAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTa                                                                         <  1:1608232/91‑1 (MQ=255)

                   tAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAAt                                                                       >  1:1730088/1‑91 (MQ=255)

                    aaaaTAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATg                                                                      >  1:2330357/1‑91 (MQ=255)

                       ataaataTTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACa                                                                   <  1:415622/91‑1 (MQ=255)

                        taaataTTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACaa                                                                  <  1:3424864/91‑1 (MQ=255)

                         aaataTTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACaaa                                                                 >  1:3889229/1‑91 (MQ=255)

                         aaataTTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACaaa                                                                 >  1:3889147/1‑91 (MQ=255)

                         aaataTTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACaaa                                                                 <  1:801952/91‑1 (MQ=255)

                           atatTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGt                                                               <  1:1981789/91‑1 (MQ=255)

                             atttattAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTa                                                             >  1:3402030/1‑91 (MQ=255)

                              tttattAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTAt                                                            >  1:2331215/1‑91 (MQ=255)

                              tttattAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTAt                                                            >  1:2331787/1‑91 (MQ=255)

                               ttattaAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTAtt                                                           >  1:3152959/1‑91 (MQ=255)

                                  ttaAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGaa                                                        <  1:2472297/91‑1 (MQ=255)

                                  ttaAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGaa                                                        <  1:2786331/91‑1 (MQ=255)

                                  ttaAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGaa                                                        <  1:2472405/91‑1 (MQ=255)

                                  ttaAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGaa                                                        >  1:261941/1‑91 (MQ=255)

                                  ttaAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGaa                                                        <  1:2784944/91‑1 (MQ=255)

                                   taAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAg                                                       <  1:39423/91‑1 (MQ=255)

                                   taAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAg                                                       <  1:660157/91‑1 (MQ=255)

                                     aaCATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTa                                                     >  1:746549/1‑91 (MQ=255)

                                      aCATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTaa                                                    >  1:1060081/1‑91 (MQ=255)

                                         ttATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAAtt                                                 >  1:2496218/1‑91 (MQ=255)

                                          tATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTg                                                >  1:741000/1‑91 (MQ=255)

                                          tATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTg                                                >  1:741162/1‑91 (MQ=255)

                                          tATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTg                                                >  1:2317200/1‑91 (MQ=255)

                                          tATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTg                                                <  1:777472/91‑1 (MQ=255)

                                          tATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTg                                                >  1:1259092/1‑91 (MQ=255)

                                             gAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGt                                             <  1:1186433/91‑1 (MQ=255)

                                              aaTTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGta                                            <  1:1793941/91‑1 (MQ=255)

                                                 ttttAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATaa                                         <  1:540708/91‑1 (MQ=255)

                                                 ttttAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATaa                                         <  1:597333/91‑1 (MQ=255)

                                                 ttttAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATaa                                         <  1:2684820/91‑1 (MQ=255)

                                                   ttAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAAc                                       >  1:1480765/1‑91 (MQ=255)

                                                   ttAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAAc                                       <  1:3072763/91‑1 (MQ=255)

                                                   ttAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAAc                                       >  1:1480258/1‑91 (MQ=255)

                                                     aaaGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACAt                                     >  1:3406382/1‑91 (MQ=255)

                                                     aaaGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACAt                                     >  1:108904/1‑91 (MQ=255)

                                                     aaaGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACAt                                     >  1:109046/1‑91 (MQ=255)

                                                      aaGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACAtt                                    >  1:767878/1‑91 (MQ=255)

                                                      aaGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACAtt                                    >  1:4125530/1‑91 (MQ=255)

                                                      aaGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACAtt                                    >  1:2102003/1‑91 (MQ=255)

                                                      aaGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACAtt                                    >  1:3559666/1‑91 (MQ=255)

                                                       agagTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACAttt                                   >  1:1722043/1‑91 (MQ=255)

                                                       agagTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACAttt                                   >  1:2543151/1‑91 (MQ=255)

                                                       agagTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACAttt                                   <  1:2384454/91‑1 (MQ=255)

                                                        gagTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTa                                  <  1:1817571/91‑1 (MQ=255)

                                                        gagTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTa                                  >  1:483633/1‑91 (MQ=255)

                                                        gagTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTa                                  <  1:2089829/91‑1 (MQ=255)

                                                          gTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTaaa                                <  1:2770002/91‑1 (MQ=255)

                                                          gTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTaaa                                <  1:2067123/91‑1 (MQ=255)

                                                          gTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTaaa                                <  1:2769930/91‑1 (MQ=255)

                                                           tAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAAt                               <  1:1745765/91‑1 (MQ=255)

                                                            aaTGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATa                              <  1:1146015/91‑1 (MQ=255)

                                                             aTGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATAc                             <  1:3623591/91‑1 (MQ=255)

                                                               gTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGa                           <  1:555056/91‑1 (MQ=255)

                                                               gTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGa                           <  1:2477394/91‑1 (MQ=255)

                                                                 cTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATg                         <  1:2484010/91‑1 (MQ=255)

                                                                    aCTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCg                      <  1:2645583/91‑1 (MQ=255)

                                                                        gTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGaaaa                  >  1:2422992/1‑91 (MQ=255)

                                                                        gTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGaaaa                  <  1:3079059/91‑1 (MQ=255)

                                                                         ttGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAAt                 <  1:3685038/91‑1 (MQ=255)

                                                                         ttGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAAt                 <  1:213783/91‑1 (MQ=255)

                                                                           tttAATTTATTTTTACAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATgg               >  1:2655737/3‑91 (MQ=38)

                                                                           gATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATgg               >  1:3728475/1‑91 (MQ=255)

                                                                           gATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATgg               >  1:3727758/1‑91 (MQ=255)

                                                                              aaTTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGa            <  1:1269521/91‑1 (MQ=255)

                                                                               atttatttTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGAt           <  1:1893041/91‑1 (MQ=255)

                                                                               atttatttTTACAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGAt           >  1:3412804/1‑91 (MQ=255)

                                                                               atttatttTTACAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGAt           >  1:1917581/1‑91 (MQ=255)

                                                                                tttatttTTACAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATa          >  1:3331128/1‑91 (MQ=255)

                                                                                tttatttTTACAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATa          >  1:1945249/1‑91 (MQ=255)

                                                                                tttatttTTACAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATa          >  1:4134213/1‑91 (MQ=255)

                                                                                 ttatttTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAg         >  1:1845958/1‑91 (MQ=255)

                                                                                  tatttTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGc        <  1:4000302/91‑1 (MQ=255)

                                                                                  tatttTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGc        >  1:3873696/1‑91 (MQ=255)

                                                                                    tttttGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCat      <  1:2791435/91‑1 (MQ=255)

                                                                                     ttttGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCata     <  1:2449993/91‑1 (MQ=255)

                                                                                      tttGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCatat    <  1:1439805/91‑1 (MQ=255)

                                                                                       ttGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATc   >  1:883247/1‑91 (MQ=255)

                                                                                        tGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCa  >  1:2244604/1‑91 (MQ=255)

TCGCTTTAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCA  >  USA300TCH1516_ALE/631446‑631624

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