Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,582,141 Δ1 bp 10.3% coding (1169/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,1410T.10.3% 649.2 / 63.1 203coding (1169/1353 nt)pnbAPara‑nitrobenzyl esterase
Reads supporting (aligned to +/- strand):  ref base T (93/88);  new base . (9/12);  total (103/100)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 216 total aligned reads displayed.

TCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACA  >  USA300TCH1516_ALE/2582051‑2582231
                                                                                          |                                                                                           
tCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCatat                                                                                             <  1:1675878/91‑1 (MQ=255)
 cAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATAtt                                                                                            <  1:2565683/91‑1 (MQ=255)
 cAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATAtt                                                                                            >  1:3381736/1‑91 (MQ=255)
  aaTATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTc                                                                                           >  1:3008331/1‑91 (MQ=255)
  aaTATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTc                                                                                           >  1:3624028/1‑91 (MQ=255)
   atatACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCt                                                                                          >  2:1423989/1‑91 (MQ=255)
    tataCAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCatatctag                                                                                         <  1:2344521/91‑5 (MQ=255)
      taCAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGa                                                                                      <  2:2310406/91‑1 (MQ=255)
      taCAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGa                                                                                      >  1:2063110/1‑91 (MQ=255)
      taCAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGa                                                                                      >  1:1305981/1‑91 (MQ=255)
      taCAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATCTAGaa                                                                                     <  1:3879057/91‑2 (MQ=255)
      taCAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATCTAGaa                                                                                     <  1:1923327/91‑2 (MQ=255)
       aCAAGATTCAAACGCACAACTATGGCTTACTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGat                                                                                     <  1:742315/91‑1 (MQ=255)
       aCAAGATGCAAACGCACAACCATGGCTTGCTGAAGTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATCTAGaat                                                                                    <  2:1713595/91‑3 (MQ=255)
        cAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGata                                                                                    >  1:1431233/1‑91 (MQ=255)
        cAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGata                                                                                    <  2:3737924/91‑1 (MQ=255)
        cAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATCTAGaata                                                                                   <  2:3871647/91‑4 (MQ=255)
         aaGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATCTAGaatat                                                                                  <  2:1787821/91‑5 (MQ=255)
           gATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATCTAGAATATgg                                                                                 <  1:861698/91‑1 (MQ=255)
            atgcAAACGCACAACTATGGCTTGCTGAATTGGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGt                                                                                <  1:3548038/88‑1 (MQ=255)
            aTTCAAACGCACAACTATGGCTTGCTGAATTTACATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGt                                                                                <  2:2393319/91‑1 (MQ=255)
             ttCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTa                                                                               <  1:2914681/91‑1 (MQ=255)
             ttCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTa                                                                               >  1:574784/1‑91 (MQ=255)
              tCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTAt                                                                              <  2:324028/91‑1 (MQ=255)
               cAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTAtt                                                                             >  1:134545/1‑91 (MQ=255)
                aaaCGTACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTAttt                                                                            <  2:3581325/91‑1 (MQ=255)
                    gCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGt                                                                        <  2:1279749/91‑1 (MQ=255)
                      acaACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTc                                                                      <  1:1972996/91‑1 (MQ=255)
                      acaACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTc                                                                      >  1:2870225/1‑91 (MQ=255)
                       caACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCg                                                                     >  1:255504/1‑91 (MQ=255)
                       caACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCg                                                                     >  1:255679/1‑91 (MQ=255)
                       caACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCg                                                                     <  1:3347883/91‑1 (MQ=255)
                         aCTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGa                                                                   >  2:3327527/1‑91 (MQ=255)
                         aCTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGa                                                                   >  2:1099197/1‑91 (MQ=255)
                         aCTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGa                                                                   >  1:1886344/1‑91 (MQ=255)
                             tGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACt                                                               >  1:24356/1‑91 (MQ=255)
                               gCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTa                                                             >  1:3032772/1‑91 (MQ=255)
                               gCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTa                                                             <  2:1738592/91‑1 (MQ=255)
                               gCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTa                                                             >  1:1561326/1‑91 (MQ=255)
                                cTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTAc                                                            <  1:3248888/91‑1 (MQ=255)
                                cTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTAc                                                            <  2:1451682/91‑1 (MQ=255)
                                   gCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACaaa                                                         >  1:3738991/1‑91 (MQ=255)
                                   gCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACaaa                                                         <  2:888867/91‑1 (MQ=255)
                                    cTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAat                                                        <  1:3515621/91‑1 (MQ=255)
                                    cTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAat                                                        <  2:3184802/91‑1 (MQ=255)
                                     tGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAata                                                       <  2:255732/91‑1 (MQ=255)
                                     tGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAata                                                       <  2:3798362/91‑1 (MQ=255)
                                         tttGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTag                                                   >  1:1738615/1‑91 (MQ=255)
                                            gcatggcatgATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTagcag                                                >  1:1042305/1‑91 (MQ=255)
                                             catggcatgATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTgg                                                                         <  2:3082546/65‑1 (MQ=255)
                                              atggcatgATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAgcagca                                              <  2:1870467/91‑1 (MQ=255)
                                              atggcatgATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAgcagca                                              <  2:3830417/91‑1 (MQ=255)
                                              atggcatgATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAgcagca                                              >  2:932851/1‑91 (MQ=255)
                                              atggcatgATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAgcagca                                              >  1:2774816/1‑91 (MQ=255)
                                                ggcatgATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGcaca                                            <  1:2741196/91‑1 (MQ=255)
                                                ggcatgATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGcaca                                            >  1:1608629/1‑91 (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:3187753/91‑1 (MQ=255)
                                                     gATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAat                                       <  2:2086868/91‑1 (MQ=255)
                                                      aTACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAata                                      <  2:1272243/91‑1 (MQ=255)
                                                       tACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAatat                                     >  1:3072388/1‑91 (MQ=255)
                                                       tACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAatat                                     >  2:102897/1‑91 (MQ=255)
                                                       caCTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAatat                                     <  2:1766083/90‑1 (MQ=255)
                                                        aCTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATc                                    >  2:338043/1‑91 (MQ=255)
                                                         cTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATcc                                   <  2:664913/91‑1 (MQ=255)
                                                         cTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATcc                                   <  1:754461/91‑1 (MQ=255)
                                                          ttCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCca                                  >  2:3851754/1‑91 (MQ=255)
                                                          ttCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCca                                  <  1:409961/91‑1 (MQ=255)
                                                          ttCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCca                                  >  1:2842473/1‑91 (MQ=255)
                                                            cTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCcaac                                <  1:151953/91‑1 (MQ=255)
                                                             tAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCcaaca                               <  1:2500812/91‑1 (MQ=255)
                                                                cgcACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACt                            >  1:1447191/1‑91 (MQ=255)
                                                                   acacTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCt                         <  1:78117/91‑1 (MQ=255)
                                                                   acacTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCt                         <  2:2349497/91‑1 (MQ=255)
                                                                    cacTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTc                        >  1:2535281/1‑91 (MQ=255)
                                                                      cTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCAc                      <  1:981577/91‑1 (MQ=255)
                                                                      cTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCAc                      <  2:3117251/91‑1 (MQ=255)
                                                                       tATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTTGCAGCACATCAATATCCAACAACTGCTCAcc                     <  1:3450464/91‑1 (MQ=255)
                                                                       tATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCAcc                     <  2:2774816/91‑1 (MQ=255)
                                                                        aTCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCt                    <  1:1712705/91‑1 (MQ=255)
                                                                        aTCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCt                    <  1:408164/91‑1 (MQ=255)
                                                                        aTCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCt                    >  1:2073082/1‑91 (MQ=255)
                                                                        aTCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCt                    <  2:1209048/91‑1 (MQ=255)
                                                                        aTCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCt                    >  2:2876269/1‑91 (MQ=255)
                                                                         tCGCAGTGCTTACCATACCCTAGTGATGTTATTTTGGCTCGGGAACTTGTAGATATTAGCATCAAATCAATATCCTACAACTGCTCAcata                   >  2:875291/1‑88 (MQ=255)
                                                                          cGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTa                  <  2:909563/91‑1 (MQ=255)
                                                                           gCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTaa                 <  2:2976289/91‑1 (MQ=255)
                                                                           gCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTaa                 <  2:1660682/91‑1 (MQ=255)
                                                                             aGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAAt               <  1:1149196/91‑1 (MQ=255)
                                                                              gTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAAtt              >  1:693060/1‑91 (MQ=255)
                                                                                gCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAAtttt            >  2:547795/1‑91 (MQ=255)
                                                                                gCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAAtttt            <  2:2757472/91‑1 (MQ=255)
                                                                                gCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAAtttt            >  1:2363746/1‑91 (MQ=255)
                                                                                   tACCATATCTAGAATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTaa         >  1:1356986/1‑91 (MQ=255)
                                                                                     ccATATCTAGAATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGt       >  1:2305277/1‑91 (MQ=255)
                                                                                     ccATATCTAGAATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGt       >  1:1910212/1‑91 (MQ=255)
                                                                                       atatTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCg     <  2:2662790/91‑1 (MQ=255)
                                                                                          ttCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATCTCCAACAACTGCTCACCTTAAATTTTTAAGTCGACa  >  2:3059216/1‑91 (MQ=255)
                                                                                          |                                                                                           
TCAATATACAAGATTCAAACGCACAACTATGGCTTGCTGAATTTGCATGGCATGATACTTCTAGCGCACACTATCGCAGTGCTTACCATATTCTAGATATGGTATTTTGGTTCGGAAACTTACAAATATTAGCAGCACATCAATATCCAACAACTGCTCACCTTAAATTTTTAAGTCGACA  >  USA300TCH1516_ALE/2582051‑2582231

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: