BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,794,352	A→T	12.4%	Y426N (TAT→AAT)	USA300HOU_RS08960 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,794,352	0	A	T	12.4%	552.1 / 57.9	192	Y426N (TAT→AAT)	USA300HOU_RS08960	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (112/56); new base T (15/9); total (127/65)
Fisher's exact test for biased strand distribution p-value = 8.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 200 total aligned reads displayed.

CAAACGGCAATGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTTATATTTCCG  >  USA300TCH1516_ALE/1794262‑1794442

 aaaCGGCAATGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATaa                                                                                           >  2:2459311/1‑91 (MQ=255)

  aaCGGCAATGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAAt                                                                                          <  2:358689/91‑1 (MQ=255)

  aaCGGCAATGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAAt                                                                                          <  1:1514290/91‑1 (MQ=255)

   aCGGCAATGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATg                                                                                         >  1:2442582/1‑91 (MQ=255)

    cGGCAATGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGa                                                                                        >  1:21705/1‑91 (MQ=255)

     ggCAATGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGaa                                                                                       >  1:1173757/1‑91 (MQ=255)

       cAATGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACa                                                                                     >  1:450672/1‑91 (MQ=255)

       cAATGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACa                                                                                     >  2:2317054/1‑91 (MQ=255)

        aaTGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATTATGAACat                                                                                    <  1:2522481/91‑1 (MQ=255)

        aaTGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACat                                                                                    <  2:1929529/91‑1 (MQ=255)

        aaTGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACat                                                                                    >  1:1449752/1‑91 (MQ=255)

         aTGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACata                                                                                   <  2:2102998/91‑1 (MQ=255)

          tGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACataa                                                                                  >  1:1389200/1‑91 (MQ=255)

           gATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACataaa                                                                                 <  2:2651285/91‑1 (MQ=255)

              tCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACataaataa                                                                              >  1:1548139/1‑91 (MQ=255)

              tCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACataaataa                                                                              <  2:2384371/91‑1 (MQ=255)

                cACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATaaataaat                                                                            >  2:1328994/1‑91 (MQ=255)

                cACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATaaataaat                                                                            <  1:2646882/91‑1 (MQ=255)

                cACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATaaataaat                                                                            >  2:3110963/1‑91 (MQ=255)

                cACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATaaataaat                                                                            >  1:2200589/1‑91 (MQ=255)

                cACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATaaataaat                                                                            <  1:3359455/91‑1 (MQ=255)

                cACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAAGGAACATaaatgaat                                                                            >  1:3499941/1‑91 (MQ=255)

                     tataAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACt                                                                       >  1:990514/1‑91 (MQ=255)

                     tataAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACt                                                                       <  2:376736/91‑1 (MQ=255)

                      ataAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTa                                                                      >  1:1272204/1‑91 (MQ=255)

                       taAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTaa                                                                     <  2:1471532/91‑1 (MQ=255)

                        aaaCCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAg                                                                    >  2:3726482/1‑91 (MQ=255)

                        aaaCCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAg                                                                    <  2:2496427/91‑1 (MQ=255)

                         aaCCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGa                                                                   >  2:1224277/1‑91 (MQ=255)

                         aaCCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGa                                                                   >  2:1922080/1‑91 (MQ=255)

                           ccTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGata                                                                 <  2:1260649/91‑1 (MQ=255)

                           ccTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGata                                                                 >  1:1802408/1‑91 (MQ=255)

                              ctTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGatatat                                                              <  1:2299036/91‑1 (MQ=255)

                              ctTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTCTCATGTTCATAATGAACATAAATAAATGCACTAAGatatat                                                              >  2:3409533/1‑91 (MQ=255)

                               ttGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATc                                                             >  1:536933/1‑91 (MQ=255)

                               ttGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATc                                                             >  1:3332150/1‑91 (MQ=255)

                                tGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCTTTTTTATAATGAACATAAATAAATGCACTAAGATATATca                                                            >  2:2715434/1‑91 (MQ=255)

                                tGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATca                                                            >  1:2988112/1‑91 (MQ=255)

                                tGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATca                                                            <  2:2825266/91‑1 (MQ=255)

                                  tCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATcaca                                                          >  1:2885817/1‑91 (MQ=255)

                                   cGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAg                                                         >  2:3669700/1‑91 (MQ=255)

                                   cGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAg                                                         <  1:2024555/91‑1 (MQ=255)

                                    gTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGc                                                        >  1:3367901/1‑91 (MQ=255)

                                        atgatCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGTACATAAATAAATGCACTAAGATATATCACAGCTaaa                                                    >  2:3507532/1‑91 (MQ=255)

                                        atgatCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTaaa                                                    >  2:3165880/1‑91 (MQ=255)

                                         tgatCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTaaaa                                                   <  2:1762849/91‑1 (MQ=255)

                                           atCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACt                                                 >  2:2748396/1‑91 (MQ=255)

                                            tCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTa                                                >  1:2905442/1‑91 (MQ=255)

                                             cTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAg                                               >  1:1571982/1‑91 (MQ=255)

                                                tAATTTGTTGTTTTAATAACTGCACTCTTTTTTCATGTTCATAATGAACGTAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGa                                            <  2:2946387/91‑1 (MQ=255)

                                                tAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGa                                            >  1:2317133/1‑91 (MQ=255)

                                                  aTTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATg                                          <  1:3554216/91‑1 (MQ=255)

                                                   tttGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGc                                         >  1:1781558/1‑91 (MQ=255)

                                                       ttgttTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCt                                     >  1:738018/1‑91 (MQ=255)

                                                       ttgttTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCt                                     >  1:987264/1‑91 (MQ=255)

                                                           tttAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTtaa                                 >  2:2314282/1‑91 (MQ=255)

                                                            ttAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTtaat                                <  1:1944637/91‑1 (MQ=255)

                                                            ttAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTtaat                                <  1:405065/91‑1 (MQ=255)

                                                            ttAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTtaat                                >  2:2145702/1‑91 (MQ=255)

                                                             taataaCTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTtaata                               >  2:646442/1‑91 (MQ=255)

                                                             taataaCTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTtaata                               >  1:936238/1‑91 (MQ=255)

                                                              aataaCTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTtaataa                              >  1:1007487/1‑91 (MQ=255)

                                                              aataaCTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTtaataa                              >  2:2200138/1‑91 (MQ=255)

                                                                  aCTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAAtt                          <  1:1697322/91‑1 (MQ=255)

                                                                   cTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTa                         >  2:803629/1‑91 (MQ=255)

                                                                   cTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTa                         >  2:1844053/1‑91 (MQ=255)

                                                                     gAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAAc                       >  2:387053/1‑91 (MQ=255)

                                                                       aCTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACat                     >  1:1125242/1‑91 (MQ=255)

                                                                       aCTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACat                     >  1:2584502/1‑91 (MQ=255)

                                                                       aCTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACat                     >  1:3133264/1‑91 (MQ=255)

                                                                         tctTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACattg                   <  2:3416737/91‑1 (MQ=255)

                                                                         tctTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACattg                   >  1:2192496/1‑91 (MQ=255)

                                                                           ttttttCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACattgat                 >  2:16215/1‑91 (MQ=255)

                                                                             ttttCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACattgattg               >  1:1726945/1‑91 (MQ=255)

                                                                             ttttCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACattgattg               >  1:2956096/1‑91 (MQ=255)

                                                                              tttCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACAttgattga              >  1:657715/1‑91 (MQ=255)

                                                                               ttCATTTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGAc             >  2:2692142/1‑91 (MQ=255)

                                                                               ttCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGAc             <  1:2481006/91‑1 (MQ=255)

                                                                                ttaTGTTCATTATGAACATGAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACt            >  1:2061094/3‑91 (MQ=255)

                                                                                tCATTTTCATAATGGACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACt            <  2:506860/91‑1 (MQ=255)

                                                                                  aTGTTCATTATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTta          >  1:1802119/1‑91 (MQ=255)

                                                                                  aTGTTCATTATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTta          >  2:769837/1‑91 (MQ=255)

                                                                                  aTGTTCATTATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTta          >  2:2949364/1‑91 (MQ=255)

                                                                                  aTGTTCATTATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTta          >  1:3084105/1‑91 (MQ=255)

                                                                                  aTGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTta          >  2:3736362/1‑91 (MQ=255)

                                                                                  aTGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTta          <  1:29254/91‑1 (MQ=255)

                                                                                   tGTTCATTATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTtat         >  2:369707/1‑91 (MQ=255)

                                                                                   tGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTtat         >  1:123076/1‑91 (MQ=255)

                                                                                   tGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTtat         >  1:123998/1‑91 (MQ=255)

                                                                                     ttCATTATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTTatat       >  1:3302353/1‑91 (MQ=255)

                                                                                      tCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTTATAtt      >  1:805164/1‑91 (MQ=255)

                                                                                       cATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTTATAttt     <  1:2619514/91‑1 (MQ=255)

                                                                                        aTAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTTATAttt     <  1:1858482/90‑1 (MQ=255)

                                                                                        aTAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTTATATTTc    >  2:2048474/1‑91 (MQ=255)

                                                                                         tAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTTATATTTcc   <  2:3499941/91‑1 (MQ=255)

                                                                                         tAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTTATATTTcc   <  2:326325/91‑1 (MQ=255)

                                                                                         tAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTTATATTTcc   >  1:819119/1‑91 (MQ=255)

                                                                                          aaTGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTTATATTTCCg  <  1:141223/91‑1 (MQ=255)

                                                                                          aaTGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTTATATTTCCg  <  2:3289942/91‑1 (MQ=255)

                                                                                          aaTGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTTATATTTCCg  <  1:592402/91‑1 (MQ=255)

CAAACGGCAATGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTTATATTTCCG  >  USA300TCH1516_ALE/1794262‑1794442

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