Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,744,476 C→G 21.1% V66L (GTA→CTA)  apt ← Adenine phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,744,4760CG21.1% 242.9 / 47.6 107V66L (GTA→CTA) aptAdenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (36/41);  new base G (14/2);  total (57/50)
Fisher's exact test for biased strand distribution p-value = 4.65e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 108 total aligned reads displayed.

ATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTG  >  USA300TCH1516_ALE/1744386‑1744562
                                                                                          |                                                                                      
atTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCTATCCCCATTGAATAAGCTAc                                                                                        <  1:1170020/91‑1 (MQ=255)
   cTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTATTCTAACAGGTGCAAAGCCAATCCCCATTAAATAAGCTACAgg                                                                                     >  1:4145085/1‑91 (MQ=255)
    tAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGa                                                                                    <  1:3039853/91‑1 (MQ=255)
     aGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCCATCCCCATTGAATAAGCTACAGGAc                                                                                   >  1:3863594/1‑91 (MQ=255)
     aGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTAGAGGAc                                                                                   >  1:1610205/1‑91 (MQ=255)
     aGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTAGAGGAc                                                                                   >  1:1609852/1‑91 (MQ=255)
      ggTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTAGAGGACa                                                                                  >  1:2790098/1‑91 (MQ=255)
      ggTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACa                                                                                  <  1:287890/91‑1 (MQ=255)
      ggTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACa                                                                                  >  1:822519/1‑91 (MQ=255)
        tCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGc                                                                                <  1:1946180/91‑1 (MQ=255)
        tCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGc                                                                                <  1:785427/91‑1 (MQ=255)
        tCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGc                                                                                <  1:785489/91‑1 (MQ=255)
         cATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGcc                                                                               <  1:3370847/91‑1 (MQ=255)
          aTACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTAGAGGACAGCCa                                                                              >  1:917808/1‑91 (MQ=255)
           tACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTAAAGGACAGCCaa                                                                             <  1:84640/91‑1 (MQ=255)
            aCTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAAt                                                                            <  1:3847872/91‑1 (MQ=255)
            aCTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAAt                                                                            >  1:1615540/1‑91 (MQ=255)
            aCTCATAACGAATGACTTCACGAGGCAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAAt                                                                            <  1:3846160/91‑1 (MQ=255)
               cATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGAt                                                                         <  1:812131/91‑1 (MQ=255)
               cATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTAAAGGACAGCCAATGAt                                                                         >  1:1747320/1‑91 (MQ=255)
               cATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTAAAGGACAGCCAATGAt                                                                         >  1:1749193/1‑91 (MQ=255)
                aTAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATa                                                                        >  1:381446/1‑91 (MQ=255)
                aTAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATa                                                                        >  1:3013096/1‑91 (MQ=255)
                 tAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATaa                                                                       <  1:88441/91‑1 (MQ=255)
                  aaCGAATGACTTCACGAGGTAATTTCCCTTCTTGTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATaaa                                                                      <  1:3394995/91‑1 (MQ=255)
                  aaCGAATGACTTCACGAGGTAATTTCCCTTCTTGTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATaaa                                                                      <  1:3393893/91‑1 (MQ=255)
                   aCGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAg                                                                     <  1:3999/91‑1 (MQ=255)
                   aCGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAg                                                                     <  1:3998/91‑1 (MQ=255)
                   aCGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTAAAGGACAGCCAATGATAAAg                                                                     >  1:1456616/1‑91 (MQ=255)
                     gAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGcc                                                                   <  1:1590417/91‑1 (MQ=255)
                     gAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGcc                                                                   <  1:3749126/91‑1 (MQ=255)
                      aaTGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCa                                                                  <  1:2565241/91‑1 (MQ=255)
                       aTGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCAc                                                                 >  1:3901586/1‑91 (MQ=255)
                       aTGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCAc                                                                 >  1:44979/1‑91 (MQ=255)
                       aTGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAACTACAGGACAGCCAATGATAAAGCCAc                                                                 >  1:1452232/1‑91 (MQ=255)
                        tGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCAcg                                                                <  1:3886077/91‑1 (MQ=255)
                        tGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCAcg                                                                <  1:3886113/91‑1 (MQ=255)
                        tGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTAAAGGACAGCCAATGATAAAGCCAcg                                                                <  1:2829982/91‑1 (MQ=255)
                           cTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTAGAGGACAGCCAATGATAAAGCCACgcgc                                                             >  1:2878795/1‑91 (MQ=255)
                           cTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTAGAGGACAGCCAATGATAAAGCCACgcgc                                                             >  1:2879718/1‑91 (MQ=255)
                           cTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACgcgc                                                             >  1:3019323/1‑91 (MQ=255)
                           cTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACgcgc                                                             >  1:3019455/1‑91 (MQ=255)
                           cTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACgcgc                                                             >  1:3019473/1‑91 (MQ=255)
                           cTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACgcgc                                                             <  1:1700212/91‑1 (MQ=255)
                           cTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTAAAGGACAGCCAATGATAAAGCCACgcgc                                                             >  1:3359300/1‑91 (MQ=255)
                            ttCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCt                                                            <  1:429140/91‑1 (MQ=255)
                             tCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCtt                                                           <  1:3956426/91‑1 (MQ=255)
                             tCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCtt                                                           <  1:3242618/91‑1 (MQ=255)
                              cACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTc                                                          >  1:466415/1‑91 (MQ=255)
                              cACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTAAAGGACAGCCAATGATAAAGCCACGCGCTTc                                                          <  1:328147/91‑1 (MQ=255)
                               aCGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTAGAGGACAGCCAATGATAAAGCCACGCGCTTCa                                                         >  1:1957029/1‑91 (MQ=255)
                                 gAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAgg                                                       <  1:2641648/91‑1 (MQ=255)
                                 gAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAgg                                                       <  1:2239231/91‑1 (MQ=255)
                                    gTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTAAAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTcc                                                    <  1:910368/91‑1 (MQ=255)
                                      aaTTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGAGCTTCAGGTCCTa                                                  >  1:4153748/1‑91 (MQ=255)
                                       aTTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTAc                                                 <  1:2932661/91‑1 (MQ=255)
                                          tCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTAGAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAAc                                              >  1:1545738/1‑91 (MQ=255)
                                              ttcttTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGata                                          <  1:613288/91‑1 (MQ=255)
                                              ttcttTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGata                                          <  1:3667892/91‑1 (MQ=255)
                                               tcttTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGatat                                         >  1:1291584/1‑91 (MQ=255)
                                                cttTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATc                                        >  1:1372386/1‑91 (MQ=255)
                                                cttTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTAAAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATc                                        <  1:1142824/91‑1 (MQ=255)
                                                  tttCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCaa                                      <  1:910792/91‑1 (MQ=255)
                                                     cTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAt                                   >  1:928234/1‑91 (MQ=255)
                                                     cTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAt                                   >  1:68120/1‑91 (MQ=255)
                                                      tAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTAGAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtc                                  >  1:4101360/1‑91 (MQ=255)
                                                       aaCAGGTGCAAAGCCAATCCCCATTGAATAAACTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtct                                 >  1:858335/1‑91 (MQ=255)
                                                         cAGGTGCAAAGCCAATCCCCATTGAATAAGCTAAAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATctct                               >  1:1522950/1‑91 (MQ=255)
                                                           ggTGCAAAGCCAATCCCCATTGAATAAGCTAAAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGt                             >  1:1053905/1‑91 (MQ=255)
                                                            gTGCAAAGCCAATCCCCATTGAATAAGCTAGAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTc                            >  1:2810273/1‑91 (MQ=255)
                                                            gTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTc                            >  1:3265309/1‑91 (MQ=255)
                                                               cAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCttt                         <  1:3554073/91‑1 (MQ=255)
                                                               cAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCttt                         <  1:3623058/91‑1 (MQ=255)
                                                                aaaGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCtttt                        >  1:53220/1‑91 (MQ=255)
                                                                aaaGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCtttt                        >  1:53089/1‑91 (MQ=255)
                                                                 aaGCCAATCCCCTTTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTg                       >  1:3179505/1‑91 (MQ=255)
                                                                 aaGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTg                       >  1:3631241/1‑91 (MQ=255)
                                                                  aGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGc                      <  1:2576296/91‑1 (MQ=255)
                                                                   gCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCg                     >  1:3838757/1‑91 (MQ=255)
                                                                    ccAATCCCCATTGAATAAACTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGt                    <  1:33268/91‑1 (MQ=255)
                                                                    ccAATCCCCATTGAATAAACTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGt                    <  1:33506/91‑1 (MQ=255)
                                                                     cAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTa                   >  1:453215/1‑91 (MQ=255)
                                                                      aaTCCCCATTGAAGAAGCTAAAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTAt                  <  1:1184719/91‑1 (MQ=255)
                                                                       aTCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTAtt                 >  1:359633/1‑91 (MQ=255)
                                                                       aTCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTAtt                 >  1:3398625/1‑91 (MQ=255)
                                                                        tCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTc                <  1:1176247/91‑1 (MQ=255)
                                                                        tCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTc                <  1:1459472/91‑1 (MQ=255)
                                                                           ccATTGAATAAGCTAGAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTAc             >  1:39627/1‑91 (MQ=255)
                                                                            cATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACa            <  1:1367596/91‑1 (MQ=255)
                                                                             aTTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACaa           >  1:1834413/1‑91 (MQ=255)
                                                                              ttGAATAAGCTAGAGGACAGCCAATGATAAATCCACGCGCTTCATGTCCTACAGCGAGATCAACATCTCTGTCTTTTGCGTATTCTACAAt          <  1:3307791/91‑1 (MQ=255)
                                                                              ttGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAAt          <  1:1837289/91‑1 (MQ=255)
                                                                                   tAAGCTAAAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTAt     <  1:2851311/91‑1 (MQ=255)
                                                                                    aaacTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATc    <  1:543331/88‑1 (MQ=255)
                                                                                    aaGCTAGAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATc    >  1:524489/1‑91 (MQ=255)
                                                                                    aaGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATc    >  1:849595/1‑91 (MQ=255)
                                                                                    aaGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATc    >  1:1858423/1‑91 (MQ=255)
                                                                                     aacTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCt   >  1:2119836/3‑91 (MQ=255)
                                                                                     aGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCt   >  1:1570043/1‑91 (MQ=255)
                                                                                      gCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTg  >  1:3639027/1‑91 (MQ=255)
                                                                                          |                                                                                      
ATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTG  >  USA300TCH1516_ALE/1744386‑1744562

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