Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 428,600 T→C 8.0% T231A (ACT→GCT)  ahpF ← Alkyl hydroperoxide reductase subunit F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE428,6000TC 8.0% 546.0 / 24.2 183T231A (ACT→GCT) ahpFAlkyl hydroperoxide reductase subunit F
Reads supporting (aligned to +/- strand):  ref base T (89/79);  new base C (2/13);  total (91/92)
Fisher's exact test for biased strand distribution p-value = 5.30e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 187 total aligned reads displayed.

CAGAACCAGTTGTTTCTTTAACAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCTTTATTTTC  >  USA300TCH1516_ALE/428510‑428690
                                                                                          |                                                                                          
cAGAACCAGTTGTTTCTTTAACAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGt                                                                                            >  1:2635762/1‑91 (MQ=255)
cAGAACCAGTTGTTTCTTTAACAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGt                                                                                            >  2:1253382/1‑91 (MQ=255)
  gAACCAGTTGTTTCTTTAACAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTAc                                                                                          >  1:1810270/1‑91 (MQ=255)
   aaCCAGTTGTTTCTTTAACAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACg                                                                                         <  1:3026745/91‑1 (MQ=255)
   aaCCAGTTGTTTCTTTAACAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACg                                                                                         <  2:2135886/91‑1 (MQ=255)
   aaCCAGTTGTTTCTTTAACAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACg                                                                                         >  2:1617200/1‑91 (MQ=255)
     ccAGTTGTTTCTTTAACAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTa                                                                                       <  2:2096054/91‑1 (MQ=255)
     ccAGTTGTTTCTTTAACAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTa                                                                                       <  2:1474918/91‑1 (MQ=255)
          tgttTCTTTAACAGTAATGAAGTTCTCAATACCAGCAGTATCATTAGCTTGGCCACCGATACGAGCAGCAACTATACCAGCACGTaaaggt                                                                                  <  1:1484424/91‑4 (MQ=255)
          tgttTCTTTAACAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCt                                                                                  >  2:1247742/1‑91 (MQ=255)
           gttTCTTTAACAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCtt                                                                                 >  1:2454632/1‑91 (MQ=255)
           gttTCTTTAACAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCtt                                                                                 >  2:2923442/1‑91 (MQ=255)
           gttTCTTTAACAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGCACGTaaaggtt                                                                                 <  2:652114/91‑5 (MQ=255)
              tCTTTAACAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGAGACGATCAGCAACTATACCAGCACGTAAAGGTTTAc                                                                              <  2:2810170/91‑1 (MQ=255)
                 ttAACAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTg                                                                           >  1:860521/1‑91 (MQ=255)
                 ttAACAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTg                                                                           >  1:2452689/1‑91 (MQ=255)
                 ttAACAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTg                                                                           >  1:2452541/1‑91 (MQ=255)
                 ttAACAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGCACGTAAAGGTTTACGta                                                                           <  2:1887953/91‑2 (MQ=255)
                  tAACAGTAGTGGAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGCACGTAAAGGTTTACGtac                                                                          <  2:2752381/91‑3 (MQ=255)
                  tAACAGTAATGAAGTTCTCAATACCAGCGGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGCACGTAAAGGTTTACGtac                                                                          <  2:2807069/91‑3 (MQ=255)
                   aaCAGTAATGAAGTTCTCAATGCCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCt                                                                         <  2:1341974/91‑1 (MQ=255)
                    aCAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCtg                                                                        >  1:1048641/1‑91 (MQ=255)
                    aCAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGCACGTAAAGGTTTACGta                                                                           >  2:3279780/1‑87 (MQ=255)
                     cAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCtgt                                                                       >  1:409768/1‑91 (MQ=255)
                      aGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCtgtg                                                                      <  1:1358139/91‑1 (MQ=255)
                      aGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCtgtg                                                                      <  1:1407645/91‑1 (MQ=255)
                       gTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTgtgt                                                                     <  1:2337856/91‑1 (MQ=255)
                       gTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTgtgt                                                                     <  1:3528770/91‑1 (MQ=255)
                        tAATGAAGTTCTCAATACCATCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGCACGTAAAGGTTTACGTACtggtat                                                                    <  1:1766991/91‑5 (MQ=37)
                        tAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTa                                                                    >  2:2507221/1‑91 (MQ=255)
                           tGAAGTTCTCAATACCAGCAGTATCATTAACCTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAAt                                                                 <  1:2048445/91‑1 (MQ=255)
                             aaGTTCTCAATACCAGCAGTATCATTAACTTGGTCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCg                                                               <  1:3172212/91‑1 (MQ=255)
                             aaGTTCTCAATACCAGCAGTATCATAAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAACCCTTTACGTGCTGTCTCAATCg                                                               >  2:35380/1‑91 (MQ=255)
                              aGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGc                                                              >  1:1508524/1‑91 (MQ=255)
                              aGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGc                                                              <  2:3390862/91‑1 (MQ=255)
                               gTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCt                                                             >  1:1619160/1‑91 (MQ=255)
                               gTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCt                                                             >  1:1622959/1‑91 (MQ=255)
                                 tctcAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGc                                                           >  1:386918/1‑91 (MQ=255)
                                   tcAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCAc                                                         >  1:2903174/1‑91 (MQ=255)
                                   tcAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCAc                                                         <  1:592227/91‑1 (MQ=255)
                                   tcAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAAACGCTGCAc                                                         >  1:1552102/1‑91 (MQ=255)
                                   gcAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCAc                                                         <  2:348913/90‑1 (MQ=255)
                                    cAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACt                                                        <  1:3354063/91‑1 (MQ=255)
                                    cAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACt                                                        >  1:413731/1‑91 (MQ=255)
                                    cAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACt                                                        <  2:3297149/91‑1 (MQ=255)
                                        aCCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCa                                                    <  2:3154836/91‑1 (MQ=255)
                                        aCCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCa                                                    <  2:3669726/91‑1 (MQ=255)
                                          cagcagTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACt                                                  <  2:3529835/91‑1 (MQ=255)
                                          cagcagTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACt                                                  <  1:1981815/91‑1 (MQ=255)
                                           agcagTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTa                                                 <  1:767156/91‑1 (MQ=255)
                                             cagTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTAGGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGc                                               <  2:2989376/91‑1 (MQ=255)
                                             cagTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTCACGCGCTGTGTAAATCGCTGCACTACCACTAGc                                               >  2:3532917/1‑91 (MQ=255)
                                              agTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTgtgt                                                                     >  2:3790146/1‑68 (MQ=255)
                                              agTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCa                                              >  1:3805279/1‑91 (MQ=255)
                                              agTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCa                                              >  1:3689324/1‑91 (MQ=255)
                                              agTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTAACGTGCTGTGTAAATCGATGCACTACCCCTAGcc                                              >  2:3712823/1‑90 (MQ=255)
                                               gTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAg                                             <  1:1848393/91‑1 (MQ=255)
                                                  tCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGac                                          <  1:1491468/91‑1 (MQ=255)
                                                   cATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGacc                                         <  2:692057/91‑1 (MQ=255)
                                                   cATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGacc                                         >  1:2282409/1‑91 (MQ=255)
                                                     ttAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGaccac                                       >  2:831018/1‑91 (MQ=255)
                                                      tAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGaccacc                                      >  2:921248/1‑91 (MQ=255)
                                                      tAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGaccacc                                      <  2:434370/91‑1 (MQ=255)
                                                      tAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGaccacc                                      >  1:554569/1‑91 (MQ=255)
                                                        aCTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGaccaccac                                    >  2:1428154/1‑91 (MQ=255)
                                                         cTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGaccaccacc                                   <  1:2846295/91‑1 (MQ=255)
                                                           tGGCCACCGATACGATCAGCAACTATATCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCaa                                 <  2:1131898/91‑1 (MQ=255)
                                                           tGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCaa                                 >  2:1557625/1‑91 (MQ=255)
                                                            ggCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAAc                                <  2:2399645/91‑1 (MQ=255)
                                                             gCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACg                               >  1:2090663/1‑91 (MQ=255)
                                                             gCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACg                               >  2:2615863/1‑91 (MQ=255)
                                                                aCCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGCGGAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGAtt                            >  2:3661236/1‑91 (MQ=255)
                                                                   gATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAg                         >  1:2790031/1‑91 (MQ=255)
                                                                   gATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAg                         >  2:3052617/1‑91 (MQ=255)
                                                                   gATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAg                         >  1:485519/1‑91 (MQ=255)
                                                                     tACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCCCTACCACTAGCAGGACCACCACCAACGATTAAGAc                       >  2:5988/1‑91 (MQ=255)
                                                                     tACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGAc                       <  1:660230/91‑1 (MQ=255)
                                                                       cGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGCCCACCACCAACGATTAAGAcat                     >  2:789386/1‑91 (MQ=255)
                                                                       cGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGAcat                     >  2:3207584/1‑91 (MQ=255)
                                                                        gATCAGCAACTATACCATCACGTAAACCTTTACGGGCTGTGTAAATCGCTGCACTACCACGAGCAGAACCACCACCAACGATTAAGAcatc                    <  2:3858081/91‑1 (MQ=255)
                                                                        gATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGAcatc                    >  1:1043974/1‑91 (MQ=255)
                                                                         aTCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGAcatca                   >  1:1734174/1‑91 (MQ=255)
                                                                         aTCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGAcatca                   >  1:3576259/1‑91 (MQ=255)
                                                                          tCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCTCTACCACTAGCAGGACCACCACCAACGATTAAGAcatcat                  <  2:3182944/91‑1 (MQ=255)
                                                                          tCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGAcatcat                  >  1:2489913/1‑91 (MQ=255)
                                                                           cAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATa                 <  2:84736/91‑1 (MQ=255)
                                                                           cAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATa                 <  1:3402924/91‑1 (MQ=255)
                                                                             gCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAg               <  2:3752021/91‑1 (MQ=255)
                                                                              cAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAgg              <  1:1776270/91‑1 (MQ=255)
                                                                                 cTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTc           <  1:1253382/91‑1 (MQ=255)
                                                                                    taCCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCttt        <  1:1386811/91‑1 (MQ=255)
                                                                                     aCCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCTTTa       >  1:3231677/1‑91 (MQ=255)
                                                                                     aCCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCTTTa       >  1:3638793/1‑91 (MQ=255)
                                                                                      ccAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCTTTAt      >  2:2189736/1‑91 (MQ=255)
                                                                                      ccAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGAACACCACCAACGATTAAGCCATCATAAGGTACTTTAt      >  2:3037495/1‑91 (MQ=255)
                                                                                       cAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCTTTAtt     >  2:1737266/1‑91 (MQ=255)
                                                                                       cAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCTTTAtt     <  1:2390183/91‑1 (MQ=255)
                                                                                        aGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCTTTAttt    <  1:3755958/91‑1 (MQ=255)
                                                                                         gTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCTTTAtttt   >  2:769919/1‑91 (MQ=255)
                                                                                          tACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCTTTATTTTc  <  2:1028402/91‑1 (MQ=255)
                                                                                          |                                                                                          
CAGAACCAGTTGTTTCTTTAACAGTAATGAAGTTCTCAATACCAGCAGTATCATTAACTTGGCCACCGATACGATCAGCAACTATACCAGTACGTAAACCTTTACGTGCTGTGTAAATCGCTGCACTACCACTAGCAGGACCACCACCAACGATTAAGACATCATAAGGTTCTTTATTTTC  >  USA300TCH1516_ALE/428510‑428690

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