Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 937,148 T→C 6.4% intergenic (+55/‑76) USA300HOU_RS04665 → / → pepA_1 NADH dehydrogenase‑like protein/Cytosol aminopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE937,1480TC 6.4% 445.6 / 15.9 140intergenic (+55/‑76)USA300HOU_RS04665/pepA_1NADH dehydrogenase‑like protein/Cytosol aminopeptidase
Reads supporting (aligned to +/- strand):  ref base T (60/71);  new base C (2/7);  total (62/78)
Fisher's exact test for biased strand distribution p-value = 2.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01

Only 100 of 144 total aligned reads displayed.

GGTATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTA  >  USA300TCH1516_ALE/937058‑937238
                                                                                          |                                                                                          
ggtATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGt                                                                                            >  1:1307837/1‑91 (MQ=255)
 gtATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTc                                                                                           >  1:518985/1‑91 (MQ=255)
   aTCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCt                                                                                         >  1:56392/1‑91 (MQ=255)
   aTCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCt                                                                                         <  1:3427283/91‑1 (MQ=255)
    tCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCtt                                                                                        <  1:101636/91‑1 (MQ=255)
    tCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCtt                                                                                        <  1:1722888/91‑1 (MQ=255)
    tCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCtt                                                                                        <  1:1240259/91‑1 (MQ=255)
    tCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCtt                                                                                        <  1:37953/91‑1 (MQ=255)
       gTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAg                                                                                     <  1:3260519/91‑1 (MQ=255)
         ttAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAgtg                                                                                   <  1:1611847/91‑1 (MQ=255)
         ttAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAgtg                                                                                   <  1:843424/91‑1 (MQ=255)
          tAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAgtgt                                                                                  >  1:2181116/1‑91 (MQ=255)
              aTTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACt                                                                              <  1:1725624/91‑1 (MQ=255)
               ttCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTg                                                                             <  1:489245/91‑1 (MQ=255)
               ttCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTg                                                                             >  1:182460/1‑91 (MQ=255)
               ttCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACTg                                                                             <  1:4110666/91‑1 (MQ=255)
                tCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGa                                                                            >  1:2817560/1‑91 (MQ=255)
                 cAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGaa                                                                           >  1:1237925/1‑91 (MQ=255)
                 cAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGaa                                                                           >  1:3252586/1‑91 (MQ=255)
                 cAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGaa                                                                           >  1:1933731/1‑91 (MQ=255)
                 cAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATATATAAAGTCCTTCAGTGTCACTGaa                                                                           >  1:3868814/1‑91 (MQ=255)
                  aaaaaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAg                                                                          <  1:1286519/91‑1 (MQ=255)
                   aaaaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAgg                                                                         <  1:2992960/91‑1 (MQ=255)
                    aaaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGACACTGAAgga                                                                        <  1:1850135/91‑2 (MQ=255)
                    aaaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACTGAAgga                                                                        <  1:516045/91‑2 (MQ=255)
                    aaaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACTGAAgga                                                                        <  1:250435/91‑2 (MQ=255)
                    aaaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACTGAAgga                                                                        <  1:4129980/91‑2 (MQ=255)
                     aaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGc                                                                       <  1:3186973/91‑1 (MQ=255)
                     aaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACTGAAggac                                                                       <  1:655312/91‑3 (MQ=255)
                       aGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCtt                                                                     <  1:2406645/91‑1 (MQ=255)
                       aGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCtt                                                                     >  1:1759527/1‑91 (MQ=255)
                        ggTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCttt                                                                    <  1:819151/91‑1 (MQ=255)
                           aaaTTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAtt                                                                 >  1:1289911/1‑91 (MQ=255)
                           aaaTTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAtt                                                                 >  1:4030958/1‑91 (MQ=255)
                           aaaTTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAtt                                                                 >  1:2778795/1‑91 (MQ=255)
                            aaTTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAttg                                                                >  1:3029695/1‑91 (MQ=255)
                            aaTTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACTGAAGGACTTTAtaa                                                                <  1:4013448/91‑3 (MQ=37)
                               tCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAttgttg                                                             <  1:1832324/91‑1 (MQ=255)
                                cTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATtgttgt                                                            >  1:493963/1‑91 (MQ=255)
                                   ggCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACa                                                         >  1:562041/1‑91 (MQ=255)
                                   ggCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACa                                                         >  1:1650970/1‑91 (MQ=255)
                                     cTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAAt                                                       <  1:2821316/91‑1 (MQ=255)
                                     cTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAAt                                                       <  1:2821139/91‑1 (MQ=255)
                                      tttAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAAtt                                                      <  1:3889794/91‑1 (MQ=255)
                                      tttAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAAtt                                                      <  1:1151590/91‑1 (MQ=255)
                                       ttAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAAttt                                                     >  1:3371320/1‑91 (MQ=255)
                                       ttAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAAttt                                                     <  1:682371/91‑1 (MQ=255)
                                       ttAAATAAAGTTTCAGCTAAACACATAAAATAAAATTATTCATTGATAAAGTCGTTCGGTGTCACTGAAGGGCTTTATTGTTGTACAAttt                                                     <  1:3316070/91‑1 (MQ=255)
                                         aaataaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGa                                                   <  1:4024728/91‑1 (MQ=255)
                                          aataaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGaa                                                  <  1:2010426/91‑1 (MQ=255)
                                          aataaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGaa                                                  >  1:1382782/1‑91 (MQ=255)
                                           ataaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAg                                                 <  1:2958821/91‑1 (MQ=255)
                                             aaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTg                                               >  1:589132/1‑91 (MQ=255)
                                               aGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAg                                             >  1:3166221/1‑91 (MQ=255)
                                                 tttCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTa                                           <  1:3642033/91‑1 (MQ=255)
                                                 tttCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTa                                           <  1:1359371/91‑1 (MQ=255)
                                                   tCAGCTAAACACATAAACTTAAATTATTCATTTATAAAGTCCTTCAGTGTCCCTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAAc                                         <  1:849827/91‑1 (MQ=255)
                                                   tCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAAc                                         <  1:3571394/91‑1 (MQ=255)
                                                    cAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACa                                        >  1:487809/1‑91 (MQ=255)
                                                       cTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACatat                                     <  1:3906127/91‑1 (MQ=255)
                                                       cTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACatat                                     <  1:2618526/91‑1 (MQ=255)
                                                       cTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACatat                                     <  1:2618500/91‑1 (MQ=255)
                                                        tAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATc                                    <  1:3168045/91‑1 (MQ=255)
                                                         aaaCACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCa                                   <  1:2865475/91‑1 (MQ=255)
                                                          aaCACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCaa                                  >  1:734800/1‑91 (MQ=255)
                                                           acacATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCaaa                                 >  1:1516722/1‑91 (MQ=255)
                                                           acacATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCaaa                                 <  1:2727726/91‑1 (MQ=255)
                                                               aTAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTa                             >  1:582618/1‑91 (MQ=255)
                                                               aTAAAATTAAATTATTCATTTATAAAGTCCTTCAGGGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTa                             >  1:638407/1‑91 (MQ=255)
                                                                tAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTaa                            >  1:3714407/1‑91 (MQ=255)
                                                                  aaattaaattATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAagg                          <  1:714133/91‑1 (MQ=255)
                                                                    attaaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAaggag                        <  1:3243393/91‑1 (MQ=255)
                                                                    attaaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAaggag                        <  1:684625/91‑1 (MQ=255)
                                                                      taaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGt                      >  1:2406481/1‑91 (MQ=255)
                                                                        aattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTAc                    <  1:120228/91‑1 (MQ=255)
                                                                        aattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTAc                    >  1:3922609/1‑91 (MQ=255)
                                                                        aattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTAc                    <  1:822949/91‑1 (MQ=255)
                                                                          ttattCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACAt                  >  1:2466310/1‑91 (MQ=255)
                                                                          ttattCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACAt                  <  1:717952/91‑1 (MQ=255)
                                                                          ttattCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACAt                  <  1:1396910/91‑1 (MQ=255)
                                                                            attCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTa                <  1:1088817/91‑1 (MQ=255)
                                                                            attCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTa                >  1:3642709/1‑91 (MQ=255)
                                                                            attCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTa                <  1:1113093/91‑1 (MQ=255)
                                                                            attCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTa                <  1:2102628/91‑1 (MQ=255)
                                                                            attCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTa                >  1:3642780/1‑91 (MQ=255)
                                                                                aTTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGaa            <  1:3092632/91‑1 (MQ=255)
                                                                                  ttATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAAtt          <  1:3006547/91‑1 (MQ=255)
                                                                                   tataAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAAttt         <  1:3041562/91‑1 (MQ=255)
                                                                                   tataAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAAttt         <  1:3038624/91‑1 (MQ=255)
                                                                                   tataAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAAttt         <  1:2473091/91‑1 (MQ=255)
                                                                                   tataAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAAttt         >  1:3590268/1‑91 (MQ=255)
                                                                                   tataAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAAttt         >  1:1751882/1‑91 (MQ=255)
                                                                                   aataAAGCCCTTCAGTGACACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTAAGAAttt         >  1:2087014/2‑91 (MQ=25)
                                                                                      aaaGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaa      >  1:3465684/1‑91 (MQ=255)
                                                                                      aaaGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaa      >  1:3019204/1‑91 (MQ=255)
                                                                                      aaaGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaa      >  1:162476/1‑91 (MQ=255)
                                                                                       aaGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaa     <  1:273624/91‑1 (MQ=255)
                                                                                         gTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaatt   >  1:1827575/1‑91 (MQ=255)
                                                                                          tCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaatta  >  1:3129359/1‑91 (MQ=255)
                                                                                          tCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaatta  >  1:3128357/1‑91 (MQ=255)
                                                                                          |                                                                                          
GGTATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTA  >  USA300TCH1516_ALE/937058‑937238

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: