Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,774,855 A→G 17.8% D777D (GAT→GAC)  clfB ← Clumping factor B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,774,8550AG17.8% 107.9 / 2.5 60D777D (GAT→GAC) clfBClumping factor B
Reads supporting (aligned to +/- strand):  ref base A (27/21);  new base G (3/9);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTCACTATCTGAGTCTGAGTCGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTC  >  USA300TCH1516_ALE/2774716‑2774989
                                                                                                                                           |                                                                                                                                      
aGTCACTATCTGAGTCTGAGTCGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAGTCTCTGTCTGAAt                                                                                                                                       >  2:9338/1‑141 (MQ=255)
   cACTATCTGAGTCTGAGTCGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCg                                                                                                                                    <  1:151517/141‑1 (MQ=255)
    aCTATCTGAGTCTGAGTCGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGa                                                                                                                                   <  1:247722/141‑1 (MQ=255)
    aCTATCTGAGTCTGAGTCGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCGTAATCGCTGTATGAATCCGa                                                                                                                                   >  2:336388/1‑141 (MQ=255)
     cTATCTGAGTCTGAGTCGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGaa                                                                                                                                  >  1:42010/1‑141 (MQ=255)
        tctgagtctgagTCGCTATCTGAATCCGAGTCGCTATCCGAGACTGAGTCGGTGTCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCAGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTAAATCCGAATcc                                                                                                                               >  2:134302/1‑140 (MQ=17)
               ttgagtcGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTatctg                                                                                                                        <  2:209396/140‑1 (MQ=255)
                 gagtcGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTGTCTGAGTCCGCATCGCTctctgtg                                                                                                                      >  2:284302/1‑139 (MQ=14)
                 gagtcGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCTGAGTCTGAGGCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGGTatctgaa                                                                                                                      >  1:176561/1‑141 (MQ=255)
                   gtcGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAGTCCGAATCGCTatctgaatc                                                                                                                    >  1:63896/1‑141 (MQ=255)
                    tcGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTatctgaatct                                                                                                                   >  1:30728/1‑141 (MQ=255)
                     cGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTatctgaatctg                                                                                                                  >  1:260357/1‑141 (MQ=255)
                          tCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAATCGCTATCTGAGTCTGAGTCGCTGTCTGAATCTGGATCACTGTCCGAATCATACTCGCTACCTGAATCTGTATCg                                                                                                             >  2:243260/1‑141 (MQ=17)
                             gAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTGTCTGAGTCTGAATCTCCATCTGCACCTGAATCGCtc                                                                                                          >  2:120981/1‑140 (MQ=14)
                                tCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAATCTGAATCGCTATCTGACTCTGAATCGCTATcc                                                                                                       >  2:89426/1‑141 (MQ=255)
                                tCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATcc                                                                                                       <  1:147779/141‑1 (MQ=255)
                                   gAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTCTCCGAg                                                                                                    >  1:293671/1‑141 (MQ=255)
                                   gAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCCGAATCGCTATCCgcg                                                                                                    >  1:25081/1‑139 (MQ=255)
                                         cTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAg                                                                                              <  1:285362/141‑1 (MQ=255)
                                         cTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATACAAATCGCTATCTTAATCTGAATCGCTATCCGAGTCTGAg                                                                                              <  1:372331/141‑1 (MQ=255)
                                            tatgAGCCTGAGTCGCTATTTGAGTCTGATTCGCTGTCTGAGTCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAATCCGAATCACTGTCTGAGTCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCg                                                                                           <  2:57895/138‑1 (MQ=17)
                                                 gTCTGATGCGCCATTTGAGTCTGAGTCGCTGTCGGAGTCTGAGTCCCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAGTCTGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGtc                                                                                      <  2:205287/141‑1 (MQ=14)
                                                      gttCGCTATCGGAGTCTGATTCTCTGTCTGAGTCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCACTATCTGAGTCCGAATCGCTATCTGAATCTGAATCGCTCTCCTAGTCTGAGTCGCTGtctgaat                                                                                 <  2:143048/139‑1 (MQ=14)
                                                       gTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGtctgaatc                                                                                >  1:367991/1‑141 (MQ=255)
                                                        tCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAATCCGAGTCGCTATCTGAATCTGAGTCGCTCGCCGAGTCTGAGTCGCTGtctgaatct                                                                               >  2:368320/1‑141 (MQ=17)
                                                         tggttgcTGTGTCTGGGTCGCTGTCTGAGTCTGAGTGGCTATCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCACTGTCTGAGTCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGtctgaatctg                                                                              <  2:108639/135‑1 (MQ=14)
                                                         cGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGtctgaatctg                                                                              <  1:68942/141‑1 (MQ=255)
                                                                    tCTGAATCGCTGTTTGGGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCt                                                                   <  2:191251/141‑1 (MQ=255)
                                                                        atttgCTGTCTGAGTCTGAGTCGCGATTTGAGGCTGAGTCGCTGTCGGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGt                                                               <  2:391197/137‑1 (MQ=17)
                                                                        attCGTTGTCTGACTCTGGGTCGCTATCCGAGTCTGAGTCGTTGTCTGAAGCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGACTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGt                                                               <  2:291065/139‑1 (MQ=17)
                                                                        aaTCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGt                                                               <  1:42669/141‑1 (MQ=255)
                                                                          tCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTcc                                                             >  1:19292/1‑141 (MQ=255)
                                                                           cGCTGTCTGAGGCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCg                                                            <  1:37194/141‑1 (MQ=255)
                                                                            gCTGTCTGAGTCTGATTCGTTATCTGAGTCTGTGGCGCTTTCTTAATCTGAATCACTGTCTGAATCTGAATCGCTATCTGAATCTGACTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGa                                                           <  2:100859/141‑1 (MQ=17)
                                                                             cTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGaa                                                          >  1:18340/1‑141 (MQ=255)
                                                                                 ctgagtctgagtCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCTGAATCGCTATCTGAATCTGAATCGCTATCTGAGTCTGAGTCGCTATCTGATTCTGTATCGCTATCTGATTCTGAATCGc                                                      >  2:8589/1‑141 (MQ=11)
                                                                                     gtctgagtcGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTatc                                                  <  1:235193/141‑1 (MQ=255)
                                                                                      tctgagtcGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGctctct                                                 >  1:178852/1‑141 (MQ=255)
                                                                                          agtcGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTatctgaat                                             >  1:181724/1‑141 (MQ=255)
                                                                                              gCTATCTGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTatctgaatctga                                         <  2:14729/141‑1 (MQ=255)
                                                                                               ttagtcGAGTCTGAGTTGTTGTCTGAAGCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAg                                        <  2:202133/136‑1 (MQ=255)
                                                                                               cTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCTGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAATCTGACTCGCTatctgagtctgat                                        >  2:185595/1‑140 (MQ=14)
                                                                                                      gatCTGAGTCGTTGTCTGAGTGTGAGTCACTGTCTGAGTGCGATTCTCTAGCTTAATCGGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTgt                                 <  2:293671/139‑1 (MQ=14)
                                                                                                       gTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCt                                                                                                           <  2:247809/66‑1 (MQ=14)
                                                                                                         ttgggTCGTTGTTTGAAGCTGAATCTTTGTCTGAGTCTGAAGCGCTATCTGAATCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTgtctg                              <  2:390542/137‑1 (MQ=11)
                                                                                                           gAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGACTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTgtctgag                            >  1:280364/1‑141 (MQ=255)
                                                                                                            ggTCGCTGTCTGAGTCTGAGTCTCTGTCTGAGTCTGAATCGCTATCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTgtctgagt                           <  2:131756/140‑1 (MQ=11)
                                                                                                            ggTCGCTGTCTGAATCTGAGTCACTTTCTGGATCTGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTgtctgagt                           <  2:30728/140‑1 (MQ=255)
                                                                                                            aGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGCGGTGCTgtctgtgg                           >  1:223423/1‑138 (MQ=255)
                                                                                                             gTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTgtctgagtc                          >  1:361550/1‑141 (MQ=255)
                                                                                                                gctggctgAATCTGAATCACTGGCTGAATCCGAATCTCTATCTGAATCTCCCTTTCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTgtctgagtctga                       <  1:370453/141‑1 (MQ=17)
                                                                                                                 cTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTGTCTGAGTCCGAATCGCTGTCTGTATCTGAGTCGCTGTCTGAGGCTGaa                      >  2:284871/1‑141 (MQ=14)
                                                                                                                    tctgaatctgaaTCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTGTATGTGGAGAATTTGCCGTCTGAGTCTGAATCg                   >  1:71525/1‑141 (MQ=17)
                                                                                                                      tgaatctgaatcACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCGGAGTCTGAATCGCt                 >  1:375975/1‑141 (MQ=255)
                                                                                                                       gaatctgaatcACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCACTGTCTGAGTCCGAGTCGCTTTCGGAGTCTGAGTCGCTTTCTGAGTCTGAGTCGCtg                >  2:98591/1‑140 (MQ=11)
                                                                                                                           ctgaatcACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTatctg            <  1:57583/141‑1 (MQ=255)
                                                                                                                             gaatcACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTatctgaa          >  1:189832/1‑141 (MQ=255)
                                                                                                                              gctcACTGTCTGAGTCTGAGTCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTatctgaat         <  2:221665/139‑1 (MQ=11)
                                                                                                                                   cTGTCTGAATCCGAATCGCTATCGGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAg    <  1:134302/141‑1 (MQ=255)
                                                                                                                                     gTCTGAGTCCGATTCGCTCTCTGAGTCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTc  <  2:367991/141‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                      
AGTCACTATCTGAGTCTGAGTCGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTC  >  USA300TCH1516_ALE/2774716‑2774989

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