Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,582,146:1 +A 10.6% coding (1174/1353 nt) pnbA → Para‑nitrobenzyl esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,582,1461.A10.6% 686.2 / 80.7 189coding (1174/1353 nt)pnbAPara‑nitrobenzyl esterase
Reads supporting (aligned to +/- strand):  ref base . (89/80);  new base A (18/2);  total (107/82)
Fisher's exact test for biased strand distribution p-value = 1.39e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

Only 100 of 215 total aligned reads displayed.

ATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACAAATG  >  USA300TCH1516_ALE/2582056‑2582235
                                                                                           |                                                                                         
ataCAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAg                                                                                            <  1:372816/91‑1 (MQ=255)
 taCAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGa                                                                                          >  1:1305982/1‑91 (MQ=255)
 taCAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGa                                                                                          <  2:2310406/91‑1 (MQ=255)
 taCAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATCTAGaa                                                                                         <  2:3107517/91‑2 (MQ=255)
 taCAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATCTAGaa                                                                                         <  2:1644554/91‑2 (MQ=255)
 gaCAAGATTCAAACGCACAACTATGGCTTGCTGAATGTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATCTAGaa                                                                                         <  2:3837186/90‑2 (MQ=255)
  aCAAGATTCAAACGCACAACTATGGCTTACTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGat                                                                                         <  1:742524/91‑1 (MQ=255)
  aCAAGATGCAAACGCACAACCATGGCTTGCTGAAGTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATCTAGaat                                                                                        <  2:1713595/91‑3 (MQ=255)
   cAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGata                                                                                        >  1:1431233/1‑91 (MQ=255)
   cAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGata                                                                                        >  2:2735306/1‑91 (MQ=255)
    aaGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGatat                                                                                       <  1:759800/91‑1 (MQ=255)
    aaGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGatat                                                                                       >  2:1984935/1‑91 (MQ=255)
      gATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATgg                                                                                     >  1:1419767/1‑91 (MQ=255)
      gATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATCTAGAATATgg                                                                                     <  1:861698/91‑1 (MQ=255)
       aTTCAAACGCACAACTATGGCTTGCTGAATTTACATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGt                                                                                    <  2:2393319/91‑1 (MQ=255)
        ttCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTa                                                                                   <  1:711823/91‑1 (MQ=255)
        ttCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTa                                                                                   >  1:574784/1‑91 (MQ=255)
        gtCAAACGCACTACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATCTAGAATATGGTa                                                                                   <  2:927039/90‑1 (MQ=255)
         tCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTAt                                                                                  <  2:324028/91‑1 (MQ=255)
         tCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTAt                                                                                  <  1:2901238/91‑1 (MQ=255)
          cAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTAtt                                                                                 >  1:3537465/1‑91 (MQ=255)
          cAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTAtt                                                                                 <  2:1203652/91‑1 (MQ=255)
           aaaCGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTAttt                                                                                >  1:3856949/1‑91 (MQ=255)
            aaCGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTAtttt                                                                               >  1:2544565/1‑91 (MQ=255)
            aaCGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTAtttt                                                                               >  1:3584477/1‑91 (MQ=255)
            aaCGCACAACTAGGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTAtttt                                                                               >  1:2976289/1‑91 (MQ=255)
             aCGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTg                                                                              >  1:1849654/1‑91 (MQ=255)
               gCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGt                                                                            <  2:518056/91‑1 (MQ=255)
               gCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGt                                                                            <  2:1279749/91‑1 (MQ=255)
               gCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGt                                                                            <  1:1349194/91‑1 (MQ=255)
                cacaACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGtt                                                                           >  2:3067855/1‑91 (MQ=255)
                 acaACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTc                                                                          >  1:2870225/1‑91 (MQ=255)
                  caACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCg                                                                         >  1:255679/1‑91 (MQ=255)
                  caACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCg                                                                         <  1:3500271/91‑1 (MQ=255)
                  caACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCg                                                                         <  1:3347883/91‑1 (MQ=255)
                    aCTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGa                                                                       >  2:1099197/1‑91 (MQ=255)
                    aCTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGa                                                                       >  2:3327527/1‑91 (MQ=255)
                      tATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGaaa                                                                     <  1:1634631/91‑1 (MQ=255)
                        tGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACt                                                                   >  1:24356/1‑91 (MQ=255)
                          gCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTa                                                                 >  1:2034192/1‑91 (MQ=255)
                          gCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTa                                                                 <  1:1207027/91‑1 (MQ=255)
                          gCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTa                                                                 <  2:1738592/91‑1 (MQ=255)
                           cTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTAc                                                                <  2:1451682/91‑1 (MQ=255)
                            ttGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACa                                                               <  1:395470/91‑1 (MQ=255)
                            ttGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACa                                                               <  1:206943/91‑1 (MQ=255)
                              gCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTTCCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACaaa                                                             <  1:3358162/91‑1 (MQ=255)
                              gCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACaaa                                                             <  1:3238585/91‑1 (MQ=255)
                               cTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAat                                                            <  2:3184802/91‑1 (MQ=255)
                                 gAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAatat                                                          <  2:902405/91‑1 (MQ=255)
                                   aTTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAGTATTa                                                        >  2:3588400/1‑91 (MQ=255)
                                   aTTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTa                                                        >  1:1603939/1‑91 (MQ=255)
                                    tttGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTag                                                       >  2:1836225/1‑91 (MQ=255)
                                       gcatggcatgATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTagcag                                                    >  1:1042305/1‑91 (MQ=255)
                                        catggcatgATACTTCTAGCGCACACTATCGCCGTGCTTACCATATTCTGATATGGTATTTTGGTTCGGAAACTTACAAATATTAgcagca                                                  >  2:3804108/1‑91 (MQ=255)
                                        catggcatgATACTTCTAGCGCACACTATCGCAGTGCTTACCATATCCTAGATATGGTATTTTgg                                                                             >  1:3082546/1‑65 (MQ=255)
                                         atggcatgATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAgcagca                                                  >  1:2774816/1‑91 (MQ=255)
                                           ggcatgATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGcaca                                                <  1:2741196/91‑1 (MQ=255)
                                            gcatgATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACAt                                               >  2:585786/1‑91 (MQ=255)
                                              atgATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCa                                             >  1:529513/1‑91 (MQ=255)
                                               tgATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCaa                                            >  2:719206/1‑91 (MQ=255)
                                                gATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAat                                           >  1:68542/1‑91 (MQ=255)
                                                gATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAat                                           <  2:2086868/91‑1 (MQ=255)
                                                 aTACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAata                                          >  2:3760381/1‑91 (MQ=255)
                                                 aTACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAata                                          <  2:1272243/91‑1 (MQ=255)
                                                  tACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAatat                                         >  2:102897/1‑91 (MQ=255)
                                                  tACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAatat                                         >  1:3071996/1‑91 (MQ=255)
                                                  caCTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAatat                                         <  2:1766083/90‑1 (MQ=255)
                                                    cTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATcc                                       <  2:664913/91‑1 (MQ=255)
                                                     ttCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCca                                      >  1:2842473/1‑91 (MQ=255)
                                                     ttCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCca                                      <  1:409961/91‑1 (MQ=255)
                                                        tAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCcaaca                                   <  1:2500812/91‑1 (MQ=255)
                                                         aGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCcaacaa                                  >  2:3574795/1‑91 (MQ=255)
                                                           cgcACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACt                                >  1:1447191/1‑91 (MQ=255)
                                                             cacaCTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGc                              <  1:192218/91‑1 (MQ=255)
                                                              acacTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCt                             <  2:2349497/91‑1 (MQ=255)
                                                              acacTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCt                             >  1:1548252/1‑91 (MQ=255)
                                                              acacTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCt                             >  1:1970079/1‑91 (MQ=255)
                                                               cacTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTc                            >  1:2535281/1‑91 (MQ=255)
                                                                 cTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCAc                          <  1:981577/91‑1 (MQ=255)
                                                                  tATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCAcc                         <  2:2774816/91‑1 (MQ=255)
                                                                   aTCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCt                        <  1:1712705/91‑1 (MQ=255)
                                                                   aTCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCt                        >  2:2876269/1‑91 (MQ=255)
                                                                    tCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGACACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCCCAGCtt                       >  2:875309/1‑91 (MQ=255)
                                                                      gCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTaa                     <  2:2976289/91‑1 (MQ=255)
                                                                        aGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAAt                   <  1:1149196/91‑1 (MQ=255)
                                                                           gCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAAtttt                <  2:2757598/91‑1 (MQ=255)
                                                                           gCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAAtttt                >  1:2363746/1‑91 (MQ=255)
                                                                             atACCATATCTAGAATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTa              >  2:938449/2‑91 (MQ=255)
                                                                              tACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTaa             <  1:3282223/91‑1 (MQ=255)
                                                                              tACCATATCTAGAATATGGTATTTTGGTTCGGAAACTTACAAATATTAGAAGCAGATCAATATCCAACAACTGCTCACCTTAAATTTTTat             >  2:3865110/1‑90 (MQ=255)
                                                                                ccATATCTAGAATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGt           >  1:1909958/1‑91 (MQ=255)
                                                                                  catatCTAGAATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTc          >  1:1724879/5‑91 (MQ=255)
                                                                                  atatTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCg         <  2:2662790/91‑1 (MQ=255)
                                                                                  atatTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCg         >  1:1497628/1‑91 (MQ=255)
                                                                                   atatCTAGAATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCg         >  2:1079294/4‑91 (MQ=255)
                                                                                   atatCTAGAATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCg         >  1:3437389/4‑91 (MQ=255)
                                                                                      tCTAGAATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACa      >  1:2304952/1‑91 (MQ=255)
                                                                                      tCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACaa     >  2:2489948/1‑91 (MQ=255)
                                                                                        tAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACAAAt   <  2:566100/91‑1 (MQ=255)
                                                                                        tAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACAAAt   >  1:2191262/1‑91 (MQ=255)
                                                                                           |                                                                                         
ATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACAAATG  >  USA300TCH1516_ALE/2582056‑2582235

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: