Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,949,944 C→G 6.9% D134H (GAT→CAT)  fliM ← flagellar biosynthesis switching/energizing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,949,9440CG 6.9% 151.2 / 3.5 58D134H (GAT→CAT) fliMflagellar biosynthesis switching/energizing protein
Reads supporting (aligned to +/- strand):  ref base C (21/32);  new base G (2/2);  total (24/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.64e-02

CGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATGAACAGCGAAGTGCCGCGCAGCGGCTTGATCTTGACCAGGTTGAGGCTGGTGGGCACGTACAGCGAGTGCACGTATTCGCCGAACTTCATG  >  NC_002947/4949800‑4950086
                                                                                                                                                |                                                                                                                                               
cGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATcgccg                                                                                                                                             >  1:2124/1‑149 (MQ=255)
        aTGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGa                                                                                                                                      >  4:102403/1‑148 (MQ=255)
         tGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTAGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGaa                                                                                                                                     >  2:237265/1‑148 (MQ=255)
           tACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAAtt                                                                                                                                 <  3:140304/149‑1 (MQ=255)
             cTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGt                                                                                                                               <  6:147547/149‑1 (MQ=255)
             cTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGt                                                                                                                               >  3:197581/1‑149 (MQ=255)
              tCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTc                                                                                                                              >  8:84395/1‑149 (MQ=255)
              tCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCATTCGATCTTGGCGTGGCGGCCATCGCCGACAAAGAAATAGTc                                                                                                                              >  2:216337/1‑149 (MQ=255)
                gAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATAGCCGCCAAAGAAATTTTcc                                                                                                                             >  7:129034/1‑148 (MQ=255)
                  aGCTCACCGTAATGATCGGCTGCCTGGCTTCCTTGAGGTCGTCGAAACTCTGGTCCAGCACCATGCGCACCACAGGGATTTCTGTCGGGGTGAACTCACGCCGTTCGATCTTGGCGTGTCGGCCATCGCCGCCAAATAAATTTTccac                                                                                                                           <  8:199771/148‑1 (MQ=255)
                     tGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGc                                                                                                                       >  7:83394/1‑149 (MQ=255)
                              tGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAacac                                                                                                               >  3:206196/1‑148 (MQ=255)
                                 tCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTccac                                                                                                                           <  1:118780/133‑1 (MQ=255)
                                    ccTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTg                                                                                                        >  2:207304/1‑149 (MQ=255)
                                           ggCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGa                                                                                                 <  1:167742/149‑1 (MQ=255)
                                            gCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAg                                                                                                <  2:249010/149‑1 (MQ=255)
                                                  ttGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGc                                                                                                                       <  7:50643/120‑1 (MQ=255)
                                                  ttGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGc                                                                                                                       >  8:50643/1‑120 (MQ=255)
                                                    gAGGTCGACGAACCTCTGGTCCAGGACGATGAGCACCAGACGCATTTCGGTCGTGGTGAACTCAGGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATGAAc                                                                                         <  7:49285/148‑1 (MQ=255)
                                                    gAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATGAAc                                                                                         <  4:246545/148‑1 (MQ=255)
                                                          gacgaAACCCTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGGGCACTCACGCCCTGCGCTCTGGGCGTGGCGGCCCTCGCCCCCAAAGGAATAGTCCCCCACCTGGAACACCAGCGAGGCGTACAGGAGGCCCAGCCAAg                                                                                  >  7:118484/1‑149 (MQ=255)
                                                          gacgaAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATGAACAGCGAAg                                                                                  >  1:69572/1‑149 (MQ=255)
                                                            cgaAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGa                                                                                                 <  6:61859/132‑1 (MQ=255)
                                                            cgaAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGa                                                                                                 >  5:61859/1‑132 (MQ=255)
                                                                acacTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATGAACAGCGAAGTGCcgc                                                                            <  1:237265/149‑1 (MQ=255)
                                                                  acTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCGTGGCGTCGAGGAGGACCAGCGAAGTGCcgcg                                                                           >  2:196564/1‑148 (MQ=255)
                                                                    tGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATGAACAGCGAAGTGCCGCGCAg                                                                        >  7:58748/1‑149 (MQ=255)
                                                                      gTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATg                                                                                            <  5:83816/127‑1 (MQ=255)
                                                                      gTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATg                                                                                            >  6:83816/1‑127 (MQ=255)
                                                                         cAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATGAACAGCGAAGTGCCGCGCAGCGGCt                                                                   <  8:83394/149‑1 (MQ=255)
                                                                             aCCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATGAACAGCGAAGTGCCGCGCAGCGGCTTGAt                                                               <  1:196564/149‑1 (MQ=255)
                                                                              cccttcgcAACACACGCAGTTCTGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATGAACAGCGAAGTGCCGCGCAGCGGCTTGATc                                                              <  6:188133/144‑1 (MQ=255)
                                                                              ccATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATGAACAGCGAAGTGCCGCGCAGCGGCTTGATc                                                              <  7:84395/149‑1 (MQ=255)
                                                                              ccATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATGAACAGCGAAGTGCCGCGCAGCGGCTTGATc                                                              <  8:129034/149‑1 (MQ=255)
                                                                               cATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATGAAc                                                                                         <  4:90581/121‑1 (MQ=255)
                                                                                 tGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATGAACAGc                                                                                      <  6:116046/122‑1 (MQ=255)
                                                                                 tGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATGAACAGc                                                                                      >  5:116046/1‑122 (MQ=255)
                                                                                   cgcACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCTCCAGCTTGAACACCAGCTTGGCGTCGAGGATGAACAGCGAAGTGCCGCGCAGCGGCTTGATCTTGAc                                                         >  7:225633/1‑149 (MQ=255)
                                                                                     caccacACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATGAACAGCGAAGTGCCGCGCAGCGGCTTGATCTTGACCa                                                       <  4:197581/149‑1 (MQ=255)
                                                                                                    ggtcggGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATGAACAGCGAAGTGCCGCGCAGCGGCTTGATCTTGACCAGGTTGAGGCtggtg                                         <  3:102403/148‑1 (MQ=255)
                                                                                                      tcggGGTTAACTCACGCCCTTCGATATTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCCGCTTGAACACCAGCTTGTCGTCGAGGATGAACAGCGAAGTGCCTAGCAGCGGCTTGATCTTGACCAGGTTGAGGCTGGTGGGc                                      <  1:216337/149‑1 (MQ=255)
                                                                                                        ggggTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATGAACAGCGAAGTGCCGCGCAGCGGCTTGATCTTGACCAGGTTGAGGCTGGTGGGCAc                                    <  5:25991/149‑1 (MQ=255)
                                                                                                           gTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCTCCAGCTTGAACACCAGCTTGGCGTCGAGGATGAACAGCGAAGTGCCGCGCAGCGGCTTGATCTTGACCAGGTTGAGGCTGGTGGGCACGTa                                 <  8:225633/149‑1 (MQ=255)
                                                                                                                   aCGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATGAACAGCGAAGTGCCGCGCAGCGGCTTGATCTTGACCAGGTTGAGGCTGGTGGGCACGTACAGCGa                           <  8:58748/147‑1 (MQ=255)
                                                                                                                   aCGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATGAACAGCGAAGTGCCGCGCAGCGGCTTGATCTTGACCAGGTTGAGGCTGGTGGGCACGTACAGCGAg                          <  2:19639/148‑1 (MQ=255)
                                                                                                                       ccTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATGAACAGCGAAGTGCCGCGCAGCGGCTTGATCTTGACCAGGTTGAGGCTGGTGGGCACGTACAGCGAGTGCAc                     <  2:2124/149‑1 (MQ=255)
                                                                                                                         ttCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATGAACAGCGAAGTGCCGCGCAGCGGCTTGATCTTGACCAGGTTGAGGCTGGTGGGCACGTACAGCGAGTGCAc                     >  6:246663/1‑147 (MQ=255)
                                                                                                                         ttCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATGAACAGCGAAGTGCCGCGCAGCGGCTTGATCTTGACCAGGTTGAGGCTGGTGGGCACGTACAGCGAGTGCAc                     <  5:246663/147‑1 (MQ=255)
                                                                                                                         ttCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATGAACAGCGAAGTGCCGCGCAGCGGCTTGATCTTGACCAGGTTGAGGCTGGTGGGCACGTACAGCGAGTGCACg                    >  7:204213/1‑148 (MQ=255)
                                                                                                                         ttCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTAGAGGATGAACAGCGAAGTGCCGCGCAGCGGCTTGATCTTGACCAGGTTGAGGCTGGTGGGCACGTACAGCGAGTGCACg                    >  8:149111/1‑148 (MQ=255)
                                                                                                                                  tctttggcgGCGATGGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATGAACAGCGAAGTGCCGCGCAGCGGCTTGATCTTGACCAGGTTGAGGCTGGTGGGCACGTACAGCGa                           >  5:183083/5‑132 (MQ=255)
                                                                                                                                  tctttggcgGCGATGGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATGAACAGCGAAGTGCCGCGCAGCGGCTTGATCTTGACCAGGTTGAGGCTGGTGGGCACGTACAGCGa                           <  6:183083/128‑1 (MQ=255)
                                                                                                                                  ggcgtggcgGCCATCGCCGCCAAAGAAATGTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATGAACAGCGAAGTGCCGCGCAGCGGCTTGATCTTGACCAGGTTGAGGCtggtgg                                        >  1:20719/1‑120 (MQ=255)
                                                                                                                                  ggcgtggcgGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGGGTCGAGGATGAACAGCGAAGTGCCGCGCAGCGGCTTGATCTTGACCAGGTTGAGGCTGGTGGGCACGTACAGCGAGTGCACGTTTTCGc             <  1:235287/146‑1 (MQ=255)
                                                                                                                                  ggcgtggcgGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATGAACAGCGAAGTGCCGCGCAGCGGCTTGATCTTGACCAGGTTGAGGCtggtgg                                        <  2:20719/119‑1 (MQ=255)
                                                                                                                                  ggcgtggcgGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATGAACAGCGAAGTGCCGCGCAGCGGCTTGATCTTGACCAGGTTGAGGCTGGTGGGCACGTACAGCGAGTGCACGTATTc               <  8:204213/144‑1 (MQ=255)
                                                                                                                                    tttggcgGCGATGGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATGAACAGCGAAGTGCCGCGCAGCGGCTTGATCTTGACCAGGTTGAGGCTGGTGGGCACGTACAGCGAGTGc                       <  2:139282/132‑1 (MQ=255)
                                                                                                                                    tttggcgGCGATGGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATGAACAGCGAAGTGCCGCGCAGCGGCTTGATCTTGACCAGGTTGAGGCTGGTGGGCACGTACAGCGAGTGc                       >  1:139282/3‑134 (MQ=255)
                                                                                                                                          ggcCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATGAACAGCGAAGTGCCGCGCAGCGGCTTGATCTTGACCAGGTTGAGGCTGGTGGGCACGTACAGCGAGTGCACGTATTCGCCGAACTTCATg  <  1:207304/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
CGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATCTTGGCGTGGCGGCCATCGCCGCCAAAGAAATTGTCCACCAGCTTGAACACCAGCTTGGCGTCGAGGATGAACAGCGAAGTGCCGCGCAGCGGCTTGATCTTGACCAGGTTGAGGCTGGTGGGCACGTACAGCGAGTGCACGTATTCGCCGAACTTCATG  >  NC_002947/4949800‑4950086

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: