Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,259,589 G→A 12.3% P59P (CCG→CCA)  fimV → motility protein FimV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,259,5890GA12.3% 148.3 / 9.7 57P59P (CCG→CCA) fimVmotility protein FimV
Reads supporting (aligned to +/- strand):  ref base G (23/27);  new base A (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 = 8.38e-01

CATGGCTGCAGCATCCGCGCTGTCGTCTGGCATGGCGAATGCGCTGGGCCTGGGGGAGCTGACCCTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGCAAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGT  >  NC_002947/2259442‑2259735
                                                                                                                                                   |                                                                                                                                                  
cATGGCTGCAGCATCCGCGCTGTCGTCTGGCATGGCGAATGCGCTGGGCCTGGGGGAGCTGACCCTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCgc                                                                                                                                                   >  8:64822/1‑148 (MQ=255)
 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                                                                                                                           <  5:10115/113‑1 (MQ=255)
                                                            gACCCTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCa                                                                                                                           >  6:10115/1‑113 (MQ=255)
                                                             aCCCTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGaa                                                                                        <  6:131224/147‑1 (MQ=255)
                                                             aCCCTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAga                                                                                      >  4:87165/1‑149 (MQ=255)
                                                               ccTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAgacc                                                                                    >  3:207603/1‑149 (MQ=255)
                                                               ccTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAgacc                                                                                    >  6:299199/1‑149 (MQ=255)
                                                               ccTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAgacc                                                                                    >  8:29378/1‑149 (MQ=255)
                                                                 tGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTcc                                                                                                     >  5:244818/1‑130 (MQ=255)
                                                                 tGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTcc                                                                                                     <  6:244818/130‑1 (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                                                                                        >  7:82891/1‑105 (MQ=255)
                                                                                                       gAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGaa                                                                                        <  8:82891/105‑1 (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)
                                                                                                                                                   |                                                                                                                                                  
CATGGCTGCAGCATCCGCGCTGTCGTCTGGCATGGCGAATGCGCTGGGCCTGGGGGAGCTGACCCTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGCAAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGT  >  NC_002947/2259442‑2259735

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