BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,794,352	A→T	8.6%	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	8.6%	706.6 / 25.9	226	Y426N (TAT→AAT)	USA300HOU_RS08960	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (116/87); new base T (16/6); total (133/93)
Fisher's exact test for biased strand distribution p-value = 1.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 234 total aligned reads displayed.

CAAACGGCAATGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTTATATTTCCG  >  USA300TCH1516_ALE/1794262‑1794442

  aaCGGCAATGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAAt                                                                                          >  1:2218110/1‑91 (MQ=255)

  aaCGGCAATGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAAt                                                                                          >  1:2862652/1‑91 (MQ=255)

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

      gCAATGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAAc                                                                                      >  2:3350480/1‑91 (MQ=255)

      gCAATGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATATTCATAATGAAc                                                                                      >  2:3282679/1‑91 (MQ=255)

      gCAATGATTCCACTGTATAAACCCCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAAc                                                                                      >  1:106050/1‑91 (MQ=255)

       cAATGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACa                                                                                     <  2:606254/91‑1 (MQ=255)

        aaTGATTCCACTGTATAATCCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACat                                                                                    <  2:2499018/91‑1 (MQ=255)

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

        aaTGATTCCACTCTATAATCCTCTTGGCGTGGATGATCTTTAATTTGTTGTTTTAGTAACTGAACTCTTTTTTTATGTTCATTATGAACat                                                                                    <  2:3446653/91‑1 (MQ=255)

          tGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACataa                                                                                  <  1:2767492/91‑1 (MQ=255)

            aTTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACataaat                                                                                >  2:1380898/1‑91 (MQ=255)

            aTTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACataaat                                                                                >  2:2067204/1‑91 (MQ=255)

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

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

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

               ccACTGTATAAGCCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACAtaaataaa                                                                             <  2:3572395/91‑1 (MQ=255)

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

                 aCTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAAGTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATg                                                                           <  2:3256889/91‑1 (MQ=255)

                 aCTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATg                                                                           >  2:3151818/1‑91 (MQ=255)

                    gTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCAc                                                                        >  1:1834447/1‑91 (MQ=255)

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

                       taAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTTTTTTTTATAATGAACATAAATAAATGCACTaa                                                                     >  2:1751010/1‑91 (MQ=255)

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

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

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

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

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

                          aCCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGat                                                                  >  2:213866/1‑91 (MQ=255)

                          aCCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGat                                                                  <  2:1039414/91‑1 (MQ=255)

                          aCCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGat                                                                  >  1:1084319/1‑91 (MQ=255)

                           ccTCTTGTCGTGGATGATCTTTAATATGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGata                                                                 <  2:2849523/91‑1 (MQ=255)

                              ctTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGatatat                                                              >  2:1439839/1‑91 (MQ=255)

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

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

                                 gTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTTATAATGAACATAAATAAATGCACTAAGATATATcac                                                           >  2:3380956/1‑91 (MQ=255)

                                 gTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATcac                                                           >  2:3380698/1‑91 (MQ=255)

                                 gTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATcac                                                           <  1:800912/91‑1 (MQ=255)

                                  tCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATcaca                                                          <  2:399389/91‑1 (MQ=255)

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

                                    gTGGATGACCTTTAATTTGTTGTCTTAATAACTGAACTCGTCTTTCAAGTTCATGATGAACAAAAATAAAAGCACAAAGACATTACACAGc                                                        >  2:3791165/1‑91 (MQ=255)

                                     tGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCt                                                       <  1:539394/91‑1 (MQ=255)

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

                                          gatCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAAc                                                  <  2:3030873/91‑1 (MQ=255)

                                          gatCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAAc                                                  <  1:1615989/91‑1 (MQ=255)

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

                                              tttAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATTATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTCGt                                              >  1:3278840/1‑91 (MQ=255)

                                              tttAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACCAAGATATATCACAGCTAAAACTAGt                                              >  1:3669622/1‑91 (MQ=255)

                                               ttAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAgcg                                             >  2:2481849/1‑89 (MQ=255)

                                               ttAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTg                                             <  2:1128120/91‑1 (MQ=255)

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

                                                    ttgttgTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCa                                        <  2:589881/91‑1 (MQ=255)

                                                    ttgttgTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCa                                        >  1:508016/1‑91 (MQ=255)

                                                     tgttgtTTTAATAACTGAACTCTTTTTTCATTTTTATATTGAACATAAATAAATGCACTAAGATATATCACAGCTAAAATTAGTGATGCAc                                       >  2:1288801/1‑91 (MQ=255)

                                                     tgttgtTTTAATAACTGAACTCTTTTTTCATATTCATAATGAACATAAATAAATGCAGTAAGATATATCACAGCTAAAACTAGTGATGCAc                                       >  1:3512013/1‑91 (MQ=255)

                                                        tgttTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCtt                                    <  1:3476081/91‑1 (MQ=255)

                                                         gttTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCttt                                   >  1:1381833/1‑91 (MQ=255)

                                                           tttAATAACTGAACTCTTTTTTCTTGTTCATAATGAACATAAAGAGATGCACTATGATATATCACAGCTAAAACTAGAGATGCACCTTtaa                                 >  2:195875/1‑91 (MQ=255)

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

                                                               ataaCTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATaaa                             >  2:3316170/1‑91 (MQ=255)

                                                               ataaCTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATaaa                             >  1:2675473/1‑91 (MQ=255)

                                                                taaCTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATaaaa                            <  1:2259722/91‑1 (MQ=255)

                                                                  aCTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAAtt                          >  2:3831200/1‑91 (MQ=255)

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

                                                                   cTGAACTCTTTTTTCATTTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAACTTa                         >  2:1114042/1‑91 (MQ=255)

                                                                   cTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTa                         <  2:2482738/91‑1 (MQ=255)

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

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

                                                                   cTGAACTCTTTCTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTa                         >  2:1814651/1‑91 (MQ=255)

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

                                                                      aaCTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACa                      <  1:240755/91‑1 (MQ=255)

                                                                      aaCTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACa                      <  2:3801543/91‑1 (MQ=255)

                                                                      aaCTCTTTTTTCATGTTCATAATGAACATAAACAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACa                      <  1:3378710/91‑1 (MQ=255)

                                                                        ctctTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACatt                    >  1:1547939/1‑91 (MQ=255)

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

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

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

                                                                          ctTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACattga                  <  1:39493/91‑1 (MQ=255)

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

                                                                               tttaTGTTCATTATGAACATGAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGAc             >  2:1193020/4‑91 (MQ=255)

                                                                               ttCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGCTTGAc             >  2:3758621/1‑91 (MQ=255)

                                                                               ttCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGAc             >  2:1619704/1‑91 (MQ=255)

                                                                               ttCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGAc             >  1:2055410/1‑91 (MQ=255)

                                                                                 taTGTTCGTTATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACtt           >  2:282864/2‑91 (MQ=255)

                                                                                 cATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACtt           >  1:1662601/1‑91 (MQ=255)

                                                                                    gTTCATTATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTtata        >  1:186241/1‑91 (MQ=255)

                                                                                    gTTCATTATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTtata        >  2:3729130/1‑91 (MQ=255)

                                                                                     ttCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTTatat       <  1:1911181/91‑1 (MQ=255)

                                                                                      tCATTATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTTATAtt      >  2:3699581/1‑91 (MQ=255)

                                                                                       tataataAACATAAATAAATGCACTAAGATATATCACAGCTAAGACTAGTTATTCAACTTTAATAAAATTAACATTGATTGACTTATAttt     >  2:816805/2‑91 (MQ=255)

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

                                                                                        aTAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTTATATTTc    <  1:1487774/91‑1 (MQ=255)

                                                                                        aTAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTTATATTTc    <  2:2741868/91‑1 (MQ=255)

                                                                                        aTAATGAACATAAATAAATGCACCAAGACATAGCACAGCTAAAACTAGTGATGCACCCTTAATAAAAGTAACATTGATTGACTTATATTTc    <  2:354121/91‑1 (MQ=255)

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

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

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

                                                                                          aaTGAACATAAATAAATGCACTAAGATATCTCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTTATATTTCCg  <  1:3169927/91‑1 (MQ=255)

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

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

CAAACGGCAATGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTTATATTTCCG  >  USA300TCH1516_ALE/1794262‑1794442

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