Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,145,549 T→G 5.6% intergenic (+83/‑278) USA300HOU_RS10945 → / → USA300HOU_RS10950 hypothetical protein/2‑oxoglutaramate amidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,145,5490TG 5.6% 636.6 / 21.8 190intergenic (+83/‑278)USA300HOU_RS10945/USA300HOU_RS10950hypothetical protein/2‑oxoglutaramate amidase
Reads supporting (aligned to +/- strand):  ref base T (88/91);  new base G (7/4);  total (95/95)
Fisher's exact test for biased strand distribution p-value = 5.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 211 total aligned reads displayed.

AAAAATAATCAATATATTGTTAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATACTAATTTCTAAAGAAAAAGTATTTCTTTATGTTAG  >  USA300TCH1516_ALE/2145459‑2145639
                                                                                          |                                                                                          
 aaaaTAATCAATATATTGTTAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCgctg                                                                                           >  1:1459520/1‑91 (MQ=255)
   aataatCAATATATTGTTAGGAATGGTACAGAAATGATAAAGGATCATTAGTTATGCATTATGTAGTAGTTCTTACACATTAGCCgctgct                                                                                         <  2:530457/91‑1 (MQ=255)
   aataatCAATATATTGTTAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCgctgct                                                                                         <  1:2111770/91‑1 (MQ=255)
    ataatCAATATATTGTTAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTGCACATTAGCAGCGGCTa                                                                                        <  2:1523034/91‑1 (MQ=37)
    ataatCAATATATTATTAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTa                                                                                        >  2:3294386/1‑91 (MQ=255)
     taatCAATATATTGTTAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTaa                                                                                       >  2:21006/1‑91 (MQ=255)
      aatCAATATATTGTTAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAAt                                                                                      <  1:3837015/91‑1 (MQ=255)
      aatCAATATATTGTTAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCAGCGGCTAAt                                                                                      <  2:2634052/91‑1 (MQ=255)
      aatCAATATATTGTTAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCAGCGGCTAAt                                                                                      <  1:2197738/91‑1 (MQ=255)
       atCAATATATTGTTAGGAATGGGACAGAAATGATAAATAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAAtg                                                                                     >  1:3280380/1‑91 (MQ=255)
       atCAATATATTGTTAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAAtg                                                                                     <  2:1833754/91‑1 (MQ=255)
       atCAATATATTGTTAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAAtg                                                                                     <  2:1931415/91‑1 (MQ=255)
       atCAATATATTGTTAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAAtg                                                                                     >  2:19578/1‑91 (MQ=255)
       atCAATATATTGTTAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAAtg                                                                                     <  2:3543165/91‑1 (MQ=255)
        tCAATATATTGTTAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAAtgt                                                                                    >  1:399099/1‑91 (MQ=255)
         cAATATATTGTTAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAAtgtg                                                                                   >  2:2105093/1‑91 (MQ=255)
         cAATATATTGTTAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAAtgtg                                                                                   >  2:2011164/1‑91 (MQ=255)
         cAATATATTGTTAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAAtgtg                                                                                   >  2:2217555/1‑91 (MQ=255)
           atatatTGTTAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCa                                                                                 >  1:1872839/1‑91 (MQ=255)
           atatatTGTTAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCa                                                                                 >  2:615506/1‑91 (MQ=255)
           atatatTGTTAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCa                                                                                 >  1:1830181/1‑91 (MQ=255)
            tatatTGTTAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCAc                                                                                <  1:2598427/91‑1 (MQ=255)
             atatTGTTAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACt                                                                               <  2:1428914/91‑1 (MQ=255)
             atatTGTTAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACt                                                                               <  1:947560/91‑1 (MQ=255)
             atatTGTTAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACt                                                                               <  1:104405/91‑1 (MQ=255)
               atTGTTAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTa                                                                             <  2:2165996/91‑1 (MQ=255)
                ttGTTAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTaa                                                                            <  2:2034037/91‑1 (MQ=255)
                 tGTTAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAg                                                                           <  2:2661692/91‑1 (MQ=255)
                  gTTAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGa                                                                          <  2:3685337/91‑1 (MQ=255)
                    tAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAAt                                                                        >  2:1056873/1‑91 (MQ=255)
                     aGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATa                                                                       <  2:624766/91‑1 (MQ=255)
                     aGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATa                                                                       <  1:1833961/91‑1 (MQ=255)
                       gAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAgg                                                                     >  1:1680562/1‑91 (MQ=255)
                        aaTGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGa                                                                    >  2:3105031/1‑91 (MQ=255)
                         aTGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACt                                                                               >  1:1269566/1‑79 (MQ=255)
                         aTGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGaa                                                                   <  1:728528/91‑1 (MQ=255)
                            ggACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATAc                                                                >  1:943639/1‑91 (MQ=255)
                            ggACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATAc                                                                >  1:3238498/1‑91 (MQ=255)
                            ggACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATAc                                                                <  1:311914/91‑1 (MQ=255)
                             gACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACa                                                               <  1:1962982/91‑1 (MQ=255)
                              aCAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACAt                                                              <  2:2807120/91‑1 (MQ=255)
                                aGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGa                                                            >  2:2733180/1‑91 (MQ=255)
                                 gAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAg                                                           <  1:1587733/91‑1 (MQ=255)
                                 gAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAg                                                           >  2:346978/1‑91 (MQ=255)
                                  aaaTGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGt                                                          <  1:3857777/91‑1 (MQ=255)
                                  aaaTGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGt                                                          <  1:3072619/91‑1 (MQ=255)
                                   aaTGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTa                                                         <  1:1893201/91‑1 (MQ=255)
                                     tGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGTACTTAAGAATAGGAATACATGAGtttt                                                       <  2:837401/91‑4 (MQ=255)
                                     tGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTaaa                                                       >  2:2695026/1‑91 (MQ=255)
                                      gATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTgaaa                                                      >  1:3129171/1‑91 (MQ=255)
                                      gATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTaaaa                                                      >  2:1880176/1‑91 (MQ=255)
                                       aTAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCgctgct                                                                                         <  1:3285015/55‑1 (MQ=255)
                                       aTAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCgctgct                                                                                         >  2:3284994/1‑55 (MQ=255)
                                           aGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCAt                                                 >  2:743218/1‑91 (MQ=255)
                                           aGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCAt                                                 >  2:1399591/1‑91 (MQ=255)
                                            gAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATg                                                <  2:1563782/91‑1 (MQ=255)
                                              aTCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGta                                              <  2:2344267/91‑1 (MQ=255)
                                               tCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGtat                                             <  2:336011/91‑1 (MQ=255)
                                                 aTTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATaa                                           >  2:17929/1‑91 (MQ=255)
                                                 aTTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATaa                                           >  1:2340502/1‑91 (MQ=255)
                                                  ttAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAg                                          >  1:2792592/1‑91 (MQ=255)
                                                   tAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAGTGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGa                                         <  1:3080282/91‑1 (MQ=255)
                                                   tAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGa                                         <  1:1129988/91‑1 (MQ=255)
                                                   tAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGa                                         >  1:2375701/1‑91 (MQ=255)
                                                     gTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGaaa                                       <  1:3282712/91‑1 (MQ=255)
                                                        tttCATTAGGAGGTAGTTCTTAAACATTAGCCGCGGGTAATGTGCACTTAGAAGTAGGAGTACTGGAGTAAAATTCGTGTATAAGAAATAc                                    <  2:2687225/90‑1 (MQ=9)
                                                        aTTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATAc                                    >  1:1713218/1‑91 (MQ=255)
                                                        aTTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAACGTGCACTTAAGAATAGGAATACATGAGTAAACCTCATGTATAAGAAATAc                                    >  2:3196955/1‑91 (MQ=255)
                                                         ttCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATACt                                   <  2:868042/91‑1 (MQ=255)
                                                              tATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATACTAAttt                              >  1:83519/1‑91 (MQ=255)
                                                              tATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATACTAAttt                              >  1:2086462/1‑91 (MQ=255)
                                                               aTGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGTATAGGAATACATGAGTAAAACTCATGTATAAGAAATACTAATTTc                             >  2:2333595/1‑91 (MQ=255)
                                                                tGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATACTAATTTCt                            <  1:2958132/91‑1 (MQ=255)
                                                                tGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATACTAATTTCt                            >  1:3726014/1‑91 (MQ=255)
                                                                 gtagtaGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATACTAATTTCTa                           >  1:3385733/1‑91 (MQ=255)
                                                                   agtagtTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATACTAATTTCTaaa                         <  1:466403/91‑1 (MQ=255)
                                                                   agtagtTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATAAGTAAAACTCATGTATAAGAAATACTAATTTCTaac                         >  2:3841086/1‑90 (MQ=255)
                                                                    gtagtTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATACTAATTTCTAAAg                        <  1:1782638/91‑1 (MQ=255)
                                                                    gtagtTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATACTAATTTCTAAAg                        <  2:1821608/91‑1 (MQ=255)
                                                                      agtTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATACTAATTTCTAAAGaa                      >  2:1229247/1‑91 (MQ=255)
                                                                        ttCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATACTAATTTCTAAAGaaaa                    <  1:894776/91‑1 (MQ=255)
                                                                          cTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATACTAATTTCTAAAGAAAAAg                  <  1:2570045/91‑1 (MQ=255)
                                                                           ttACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATACTAATTTCTAAAGAAAAAGt                 >  2:2837727/1‑91 (MQ=255)
                                                                            tACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATACTAATTTCTAAAGAAAAAGTa                >  1:1232647/1‑91 (MQ=255)
                                                                             gcacaTTAGCAGCGGCTAATGTGTACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATACTAATTTCTAAAGAAAAAGTAt               >  1:792101/2‑91 (MQ=14)
                                                                               acaTTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATACTAATTTCTAAAGAAAAAGTAttt             <  2:1508535/91‑1 (MQ=255)
                                                                               acaTTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATACTAATTTCTAAAGAAAAAGTAttt             <  1:2858072/91‑1 (MQ=255)
                                                                                caTTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATACTAATTTCTAAAGAAAAAGTATTTc            >  2:721280/1‑91 (MQ=255)
                                                                                caTTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATACTAATTTCTAAAGAAAAAGTATTTc            >  1:979834/1‑91 (MQ=255)
                                                                                 aTTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATACTAATTTCTAAAGAAAAAGTATTTCt           >  2:161477/1‑91 (MQ=255)
                                                                                 aTTAGCAGCGGCTAATGTGTACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATACTAATTTCTAAAGAAAAAGTATTTCt           >  1:3587528/1‑91 (MQ=11)
                                                                                  ttAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATACTAATTTCTAAAGAAAAAGTATTTCtt          >  1:3325883/1‑91 (MQ=255)
                                                                                   tAGCAGCGGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATACTAATTTCTAAAGAAAAAGTATTTCttt         >  2:3188964/1‑91 (MQ=255)
                                                                                      ccGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATACTAATTTCTAAAGAAAAAGTATTTCTTTATg      >  2:1793728/1‑91 (MQ=255)
                                                                                       cGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATACTAATTTCTAAAGAAAAAGTATTTCTTTATGt     >  2:1308800/1‑91 (MQ=255)
                                                                                       cGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATACTAATTTCTAAAGAAAAAGTATTTCTTTATGt     <  1:1297127/91‑1 (MQ=255)
                                                                                        gctgctAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATACTAATTTCTAAAGAAAAAGTATTTCTTTATGtt    >  2:2549577/1‑91 (MQ=255)
                                                                                        gctgctAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATACTAATTTCTAAAGAAAAAGTATTTCTTTATGtt    >  1:321338/1‑91 (MQ=255)
                                                                                        gctgctAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATACTAATTTCTAAAGAAAAAGTATTTCTTTATGtt    >  2:3249714/1‑91 (MQ=255)
                                                                                          tgctAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATACTAATTTCTAAAGAAAAAGTATTTCTTTATGTTag  >  1:430373/1‑91 (MQ=255)
                                                                                          |                                                                                          
AAAAATAATCAATATATTGTTAGGAATGGGACAGAAATGATAAAGAATCATTAGTTATTCATTATGTAGTAGTTCTTACACATTAGCCGCTGCTAATGTGCACTTAAGAATAGGAATACATGAGTAAAACTCATGTATAAGAAATACTAATTTCTAAAGAAAAAGTATTTCTTTATGTTAG  >  USA300TCH1516_ALE/2145459‑2145639

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: