BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	993,825	A→G	5.4%	intergenic (+29/‑183)	dppE → / → appA	Dipeptide‑binding protein DppE/Oligopeptide‑binding protein AppA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	993,825	0	A	G	5.4%	777.6 / 21.3	243	intergenic (+29/‑183)	dppE/appA	Dipeptide‑binding protein DppE/Oligopeptide‑binding protein AppA
Reads supporting (aligned to +/- strand): ref base A (136/93); new base G (10/3); total (146/97)
Fisher's exact test for biased strand distribution p-value = 2.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 255 total aligned reads displayed.

CACTTAAACATGTATATATTGATAAATCGATAGATAAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGTTTTTTAAATTACTGAAATTGATTAAATGCATA  >  USA300TCH1516_ALE/993735‑993915

cACTTAAACATGTATATATTGATAAATCGATAGATAAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACata                                                                                            >  2:1070094/1‑91 (MQ=255)

  cTTAAACATGTATATATTGATAAATCGATAGATAAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATAtc                                                                                          >  1:3277654/1‑91 (MQ=255)

  cTTAAACATGTATATATTGATAAATCGATAGATAAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATAtc                                                                                          <  1:3319372/91‑1 (MQ=255)

   ttAAACATGTATATATTGATAAATCGATAGATAAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATAtct                                                                                         <  2:2879284/91‑1 (MQ=255)

     aaaCATGTATATATTGATAAATCGATAGATAAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTcc                                                                                       >  2:3834565/1‑91 (MQ=255)

     aaaCATGTATATATTGATAAATCGATAGATAAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTcc                                                                                       <  1:3911966/91‑1 (MQ=255)

     aaaCATGTATATATTGATAAATCGATAGATAAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTcc                                                                                       <  1:2952064/91‑1 (MQ=255)

        cATGTATATATTGATAAATCGATAGATAAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGACTGGAGACATGTCTCCAgc                                                                                    <  1:2789389/91‑2 (MQ=255)

         aTGTATATATTGATAAATCGATAGATAAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTc                                                                                   <  2:3602407/91‑1 (MQ=255)

          tGTATATATTGATAAATCGATAGATAAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCt                                                                                  >  1:3537238/1‑91 (MQ=255)

          tGTATATATTGATAAATCGATAGATAAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCt                                                                                  >  2:2770835/1‑91 (MQ=255)

           gTATATATTGATAAATCGATAGATAAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCtt                                                                                 <  1:1939916/91‑1 (MQ=255)

           gTATATATTGATAAATCGATAGATAAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCtt                                                                                 >  2:1247756/1‑91 (MQ=255)

              tatatTGATAAATCGATAGATAAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCttttt                                                                              <  2:1390809/91‑1 (MQ=255)

              tatatTGATAAATCGATAGATAAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCttttt                                                                              >  1:2522088/1‑91 (MQ=255)

               atatTGATAAATCGATAGATAAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTtg                                                                             >  2:3091701/1‑91 (MQ=255)

                tatTGATAAATCGATAGATAAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTtgt                                                                            >  2:3193407/1‑91 (MQ=255)

                tatTGATAAATCGATAGATAAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTtgt                                                                            >  1:902142/1‑91 (MQ=255)

                  ttGATAAATCGATAGATAAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTgtgt                                                                          >  2:1662643/1‑91 (MQ=255)

                  ttGATAAATCGATAGATAAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTACAGTCTTTTTgtgt                                                                          >  2:3789613/1‑91 (MQ=255)

                     aTAAATCGATAGATAAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTgg                                                                       <  2:1794491/91‑1 (MQ=255)

                      tAAATCGATAGATAAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGa                                                                      >  2:3154324/1‑91 (MQ=255)

                       aaaTCGATAGATAAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGAt                                                                     >  2:1913323/1‑91 (MQ=255)

                        aaTCGATAGATAAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATa                                                                    <  1:993184/91‑1 (MQ=255)

                        aaTCGATAGATAAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATa                                                                    >  1:2041121/1‑91 (MQ=255)

                        aaTCGATAGATAAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATa                                                                    <  1:1478818/91‑1 (MQ=255)

                        aaTCGATAGATAAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATa                                                                    <  2:644093/91‑1 (MQ=255)

                        aaTCGATAGATAAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTAGAGACATATCTCCAGTCTTTTTGTGTTGGATa                                                                    <  2:3878186/91‑1 (MQ=255)

                         aTCGATAGATAAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATaa                                                                   <  2:3946984/91‑1 (MQ=255)

                          tCGATAGATAAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATaaa                                                                  >  2:1347415/1‑91 (MQ=255)

                            gatagataAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATaaaac                                                                >  2:2387423/1‑90 (MQ=255)

                              tagataAAGAAACTGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAAca                                                              >  2:2551720/1‑90 (MQ=255)

                              tagataAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACt                                                              >  1:2678323/1‑91 (MQ=255)

                               agataAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGCTAAAACCtt                                                             >  1:3552290/1‑91 (MQ=255)

                               agataAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACtt                                                             >  1:2919935/1‑91 (MQ=255)

                               agataAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACtt                                                             <  2:756302/91‑1 (MQ=255)

                                  taAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTgg                                                          >  1:225653/1‑91 (MQ=255)

                                   aaaGAAACAGGTAAGAAGAAGAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTggg                                                         <  1:1268931/91‑1 (MQ=255)

                                   aaaGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTggg                                                         <  1:2464225/91‑1 (MQ=255)

                                     aGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGaa                                                       >  2:3309667/1‑91 (MQ=255)

                                     aGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGaa                                                       >  1:2892747/1‑91 (MQ=255)

                                      gAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAAt                                                      >  1:1183448/1‑91 (MQ=255)

                                       aaaCAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATa                                                     <  2:1834102/91‑1 (MQ=255)

                                         aCAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATaaa                                                   >  2:592695/1‑91 (MQ=255)

                                         aCAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATaaa                                                   >  1:2063308/1‑91 (MQ=255)

                                          cAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATaaaa                                                  >  2:2731054/1‑91 (MQ=255)

                                          cAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATaaaa                                                  >  1:1002279/1‑91 (MQ=255)

                                          cAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATaaaa                                                  >  2:2212727/1‑91 (MQ=255)

                                          cAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATaaaa                                                  >  2:512420/1‑91 (MQ=255)

                                          agggTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATaaaa                                                  <  2:3939703/89‑1 (MQ=255)

                                            ggTAAGAAGAAAAAATAATCTGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAAt                                                >  2:3455365/1‑91 (MQ=255)

                                            ggTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAAt                                                <  2:2833314/91‑1 (MQ=255)

                                            ggTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAAt                                                >  1:47134/1‑91 (MQ=255)

                                            ggTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTGTTGTGTTGGATAAAAACTTTGGGAATAAAAAt                                                >  2:3047856/1‑91 (MQ=255)

                                             gTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAAtt                                               <  1:1877988/91‑1 (MQ=255)

                                                  aagaaAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTaaaa                                          >  1:3178292/1‑91 (MQ=255)

                                                   agaaAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAAt                                         <  2:1288145/91‑1 (MQ=255)

                                                      aaaaaTAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAg                                      >  1:3280542/1‑91 (MQ=255)

                                                      aaaaaTAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAg                                      <  1:3667204/91‑1 (MQ=255)

                                                      aaaaaTAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAg                                      <  2:1134875/91‑1 (MQ=255)

                                                      aaaaaTAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAg                                      <  2:2358002/91‑1 (MQ=255)

                                                       aaaaTAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGt                                     <  2:2948962/91‑1 (MQ=255)

                                                          ataataTGCTTTGTAAATTAGGCTGGCGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGt                                  >  2:2244486/1‑91 (MQ=255)

                                                          ataataTGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGt                                  <  1:2212854/91‑1 (MQ=255)

                                                          ataataTGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGt                                  <  1:3815331/91‑1 (MQ=255)

                                                           taataTGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGtt                                 >  2:3502737/1‑91 (MQ=255)

                                                           taataTGCTTTGTAAATTAGGCGGGAGACATTACTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGtt                                 <  2:911670/91‑1 (MQ=255)

                                                            aataTGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGttt                                >  2:1589557/1‑91 (MQ=255)

                                                             atatGCTTTGTAAATTAGGCTGGAGACGTATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGtttt                               <  2:539395/91‑1 (MQ=255)

                                                             atatGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGtttt                               >  1:1667504/1‑91 (MQ=255)

                                                             atatGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGCCATAAAAATTAAAAATAAGTCGtttt                               >  2:3420593/1‑91 (MQ=255)

                                                               atGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGtttttt                             <  1:829085/91‑1 (MQ=255)

                                                               atGCTTTGTAAAGTAGGCTGGAGACATATCTCCAGTCTTTTTGTGGTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGtttttt                             <  2:1574624/91‑1 (MQ=255)

                                                                tGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGTTTTTTa                            <  1:2648759/91‑1 (MQ=255)

                                                                 gCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTTGTTTTTTaa                           <  1:1982808/91‑1 (MQ=255)

                                                                  cTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGTTTTTTaaa                          <  1:674655/91‑1 (MQ=255)

                                                                       tAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGTTTTTTAAATTACt                     <  1:1584176/91‑1 (MQ=255)

                                                                         aaTTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGTTTTTTAAATTACTGa                   >  2:1331644/1‑91 (MQ=255)

                                                                          aTTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGTTTTTTAAATTACTGaa                  >  1:2283519/1‑91 (MQ=255)

                                                                          aTTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGTTTTTTAAATTACTGaa                  <  1:3386306/91‑1 (MQ=255)

                                                                          aTTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGTTTTTTAAATTACTGaa                  >  2:1049081/1‑91 (MQ=255)

                                                                           ttAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGTTTTTTAAATTACTGaaa                 <  2:2756520/91‑1 (MQ=255)

                                                                           ttAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGTTTTTTAAATTACTGaaa                 >  2:2139695/1‑91 (MQ=255)

                                                                             aGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGTTTTTTAAATTACTGAAAtt               <  1:2707352/91‑1 (MQ=255)

                                                                               acTGGAGATATGTCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGTTTTTTAAATTACTGAAATTGa             >  1:530069/2‑91 (MQ=255)

                                                                                cTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGTTTTTTAAATTACTGAAATTGAt            >  1:1747617/1‑91 (MQ=255)

                                                                                cTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGTTTTTTAAATCACTGAAATTGAt            <  1:3030222/91‑1 (MQ=255)

                                                                                 tGGAGATATGTCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGTTTTTTAAATTACTGAAATTGAtt           >  1:1501401/1‑91 (MQ=255)

                                                                                 tGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGTTTTTTAAATTACTGAAATTGAtt           >  1:707284/1‑91 (MQ=255)

                                                                                 tGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGTTTTTTAAATTACTGAAATTGAtt           >  1:1196131/1‑91 (MQ=255)

                                                                                  ggAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAAGAAAAATTTAAAATAAGTCGTTTTTTAAATTACTGAAATTGATTa          >  2:3545883/1‑91 (MQ=255)

                                                                                   gagaCATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGTTTTTTAAATTACTGAAATTGATTaa         >  1:1021893/1‑91 (MQ=255)

                                                                                   gagaCATATCACCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGTTTTTTAAATTACTGAAATTGATTaa         <  2:3568108/91‑1 (MQ=255)

                                                                                    agaCATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGTTTTTTAAATTACTGAAATTGATTaaa        <  1:253022/91‑1 (MQ=255)

                                                                                    agaCATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGTTTTTTAAATTACTGAAATTGATTaaa        <  1:3597154/91‑1 (MQ=255)

                                                                                    agaCATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGTTTTTTAAATTACTGAAATTGATTaaa        >  2:2497355/1‑91 (MQ=255)

                                                                                        atatCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGTTTTTTAAATTACTGAAATTGATTAAATGCa    >  1:3898389/1‑91 (MQ=255)

                                                                                        atatCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGTTTTTTAAATTACTGAAATTGATTAAATGCa    <  2:2837565/91‑1 (MQ=255)

                                                                                        atatCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGTTTTTTAAATTACTGAAATTGATTAAATGCa    >  2:890331/1‑91 (MQ=255)

                                                                                          atCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGTTTTTTAAATTACTGAAATTGATTAAATGCata  <  1:132056/91‑1 (MQ=255)

CACTTAAACATGTATATATTGATAAATCGATAGATAAAGAAACAGGTAAGAAGAAAAAATAATATGCTTTGTAAATTAGGCTGGAGACATATCTCCAGTCTTTTTGTGTTGGATAAAAACTTTGGGAATAAAAATTTAAAATAAGTCGTTTTTTAAATTACTGAAATTGATTAAATGCATA  >  USA300TCH1516_ALE/993735‑993915

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