Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,259,590 A→G 12.3% S60G (AGC→GGC)  fimV → motility protein FimV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,259,5900AG12.3% 139.9 / 10.6 57S60G (AGC→GGC) fimVmotility protein FimV
Reads supporting (aligned to +/- strand):  ref base A (23/27);  new base G (3/4);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

ATGGCTGCAGCATCCGCGCTGTCGTCTGGCATGGCGAATGCGCTGGGCCTGGGGGAGCTGACCCTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGCAAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAA  >  NC_002947/2259443‑2259738
                                                                                                                                                   |                                                                                                                                                    
aTGGCTGCAGCATCCGCGCTGTCGTCTGGCATGGCGAATGCGCTGGGCCTGGGGGAGCTGACCCTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAg                                                                                                                                                     >  3:175682/1‑149 (MQ=255)
    cTGCAGCATCCGCGCTGTCGTCTGGCATGGCGAATGCGCTGGGCCTGGGGGAGCTGACCCTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTg                                                                                                                                                 <  2:35207/149‑1 (MQ=255)
            tCCGCGCTGTCGTCTGGCATGGCGAATGCGCTGGGCCTGGGGGAGCTGACCCTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCAcc                                                                                                                                         <  7:127090/149‑1 (MQ=255)
             ccGCGCTGTCGTCTGGCATGGCGAATGCGCTGGGCCTGGGGGAGCTGACCCTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCg                                                                                                                                        >  2:101062/1‑149 (MQ=255)
              cgcgCTGTCGTCTGGCATGGCGAATGCGCTGGGCCTGGGGGAGCTGACCCTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCAGGCTCGGCGCCActtc                                                                                                                                       <  3:9634/149‑4 (MQ=255)
              cgcgCTGTCGTCTGGCATGGCGAATGCGCTGGGCCTGGGGGAGCTGACCCTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCAGGCTCGGCGCCActtc                                                                                                                                       <  3:13889/149‑4 (MQ=255)
                cgcTGTCGTCTGGCATGGCGAATGCGCTGGGCCTGGGGGAGCTGACCCTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGaa                                                                                                                                     <  2:118492/149‑1 (MQ=255)
                  cTGTCGTCTGGCATGGCGAATGCGCTGGGCCTGGGGGAGCTGACCCTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAaga                                                                                                                                   <  1:101062/149‑1 (MQ=255)
                     tcgtcTGGCATGGCGAATGCGCTGGGCCTGGGGGAGCTGACCCTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGtt                                                                                                                                <  3:186399/149‑1 (MQ=255)
                       gtcTGGCATGGCGAATGCGCTGGGCCTGGGGGAGCTGACCCTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTc                                                                                                                               <  1:208983/148‑1 (MQ=255)
                        tcTGGCATGGCGAATGCGCTGGGCCTGGGGGAGCTGACCCTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAg                                                                                                                             >  2:180038/1‑149 (MQ=255)
                            gCATGGCGAATGCGCTGGGCCTGGGGGAGCTGACCCTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGc                                                                                                                            <  1:180038/146‑1 (MQ=255)
                                  cGAATGCGCTGGGCCTGGGGGAGCTGACCCTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCggg                                                                                                                   >  2:244513/1‑149 (MQ=255)
                                               ccTGGGGGAGCTGACCCTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGa                                                                                                      >  2:71521/1‑149 (MQ=255)
                                                    gggAGCTGACCCTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGacctac                                                                                                 <  8:180344/149‑1 (MQ=255)
                                                        gCTGACCCTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATGGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTc                                                                                              >  1:216444/1‑148 (MQ=255)
                                                           gACCCTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCa                                                                                                                              >  6:10115/1‑113 (MQ=255)
                                                           gACCCTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCa                                                                                                                              <  5:10115/113‑1 (MQ=255)
                                                            aCCCTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGaa                                                                                           <  6:131224/147‑1 (MQ=255)
                                                            aCCCTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAga                                                                                         >  4:87165/1‑149 (MQ=255)
                                                              ccTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAgacc                                                                                       >  8:29378/1‑149 (MQ=255)
                                                              ccTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAgacc                                                                                       >  6:299199/1‑149 (MQ=255)
                                                              ccTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAgacc                                                                                       >  3:207603/1‑149 (MQ=255)
                                                                tGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTcc                                                                                                        <  6:244818/130‑1 (MQ=255)
                                                                tGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTcc                                                                                                        >  5:244818/1‑130 (MQ=255)
                                                                tGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGa                                                                                                      >  4:169699/1‑132 (MQ=255)
                                                                tGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGa                                                                                                      <  3:169699/132‑1 (MQ=255)
                                                                 gAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGaa                                                                                                                                     <  8:130860/100‑1 (MQ=255)
                                                                  aaGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAgacctgac                                                                                   <  7:29378/149‑1 (MQ=255)
                                                                  aaGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACATCGAAgacctgac                                                                                   <  1:244513/149‑1 (MQ=255)
                                                                       ggaacaGAAGCAGCCGCTGGACGCCTAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTca                                                                              <  4:218770/146‑1 (MQ=255)
                                                                        gCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTca                                                                              >  4:75131/1‑148 (MQ=255)
                                                                                    ccGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCaa                                                                 <  7:64822/149‑1 (MQ=255)
                                                                                     cGCTGGACGCCGAGATCTAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAAc                                                                <  4:22410/149‑1 (MQ=255)
                                                                                     cGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAAc                                                                <  2:216444/149‑1 (MQ=255)
                                                                                       cTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTcc                                                                                                        <  5:141818/107‑1 (MQ=255)
                                                                                       cTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTcc                                                                                                        >  6:141818/1‑107 (MQ=255)
                                                                                       cTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAAccc                                                              <  4:129710/149‑1 (MQ=255)
                                                                                          gACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCaa                                                           >  8:10216/1‑149 (MQ=255)
                                                                                                    tCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGCAAAAgt                                                 <  7:10216/149‑1 (MQ=255)
                                                                                                      gAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGaa                                                                                           <  8:82891/105‑1 (MQ=255)
                                                                                                      gAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGaa                                                                                           >  7:82891/1‑105 (MQ=255)
                                                                                                      gAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGCCCTGACCTTCACACCGGTGCTCACCCCCAACGGCAAAAgtgt                                               >  4:8240/1‑149 (MQ=255)
                                                                                                       aaTTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGCAAAAGtgtg                                              <  4:175682/149‑1 (MQ=255)
                                                                                                            cTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACgg                                                        >  7:190478/1‑134 (MQ=255)
                                                                                                            cTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACgg                                                        <  8:190478/134‑1 (MQ=255)
                                                                                                                         aCCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGCAAAAGTGTGCTGCGCGTGACCTCCAGt                            <  3:75131/149‑1 (MQ=255)
                                                                                                                            tCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGCAAAAGTGTGCTGCGCGTGACCTCCAGTCAg                         >  4:283818/1‑149 (MQ=255)
                                                                                                                                     ttccggTGGCGCCAGGCTCGGCGCCACTGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGCAAAAGTGTGCTGCGCGTGACCTCCa                              >  2:251069/6‑135 (MQ=255)
                                                                                                                                     ttccggTGGCGCCAGGCTCGGCGCCACTGGAAGAGTTCACCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGCAAAAGTGTGCTGCGCGTGACCTCCa                              <  1:251069/130‑1 (MQ=255)
                                                                                                                                     ccgAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGCAAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGc                  >  6:80070/1‑147 (MQ=255)
                                                                                                                                       gAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGCAAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAgg              >  7:237972/1‑149 (MQ=255)
                                                                                                                                         ggTGGCGCCAGGCTCGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGCAAAAGTGTGCTGCGCGTGa                                    <  6:91624/124‑1 (MQ=255)
                                                                                                                                         ggTGGCGCCAGGCTCGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGCAAAAGTGTGCTGCGCGTGa                                    >  5:91624/2‑125 (MQ=255)
                                                                                                                                          gTGGCGCCAGGCTCGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGCAAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTcc           >  1:305065/1‑149 (MQ=255)
                                                                                                                                           tGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGCAAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCg          <  8:237972/149‑1 (MQ=255)
                                                                                                                                                 cGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCATCCCCAACGGCAACAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCggtggt     >  7:65045/1‑148 (MQ=255)
                                                                                                                                                   aGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGCAAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCaa  <  5:80070/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ATGGCTGCAGCATCCGCGCTGTCGTCTGGCATGGCGAATGCGCTGGGCCTGGGGGAGCTGACCCTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGCAAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAA  >  NC_002947/2259443‑2259738

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