Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,259,594 T→C 12.3% L61P (CTG→CCG) ‡ fimV → motility protein FimV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,259,5940TC12.3% 141.9 / 10.2 57L61P (CTG→CCG) ‡fimVmotility protein FimV
Reads supporting (aligned to +/- strand):  ref base T (22/28);  new base C (3/4);  total (25/32)
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.84e-01

CTGCAGCATCCGCGCTGTCGTCTGGCATGGCGAATGCGCTGGGCCTGGGGGAGCTGACCCTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGCAAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTC  >  NC_002947/2259447‑2259742
                                                                                                                                                   |                                                                                                                                                    
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                                                                                           >  6:299199/1‑149 (MQ=255)
                                                          ccTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAgacc                                                                                           >  8:29378/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                                                                                               >  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                                                            <  8:190478/134‑1 (MQ=255)
                                                                                                        cTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACgg                                                            >  7:190478/1‑134 (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)
                                                                                                                                                   tGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGCAAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTc  >  3:291617/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTGCAGCATCCGCGCTGTCGTCTGGCATGGCGAATGCGCTGGGCCTGGGGGAGCTGACCCTGAAGTCGGCACAGAATCAGCCGCTGGACGCCGAGATCGAATTGCTCGACGTGCGCGACCTCAAAGCCGCCGAAGTGGCGCCGAGCCTGGCGCCACCGGAAGAGTTCAGCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGCAAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTC  >  NC_002947/2259447‑2259742

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: ‑