Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 308,924 G→A R137C (CGC→TGC)  ecpC ← putative fimbrial usher protein EcpC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913308,9240GA97.1% 112.6 / NA 34R137C (CGC→TGC) ecpCputative fimbrial usher protein EcpC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (19/14);  minor base T (0/1);  total (19/15)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

TACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAACGGGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCGTTATCTTCACTCTCTTCCAGCTGTATTTTCCGGATGCGTAGCTGTCCATCGTCCAGCCAGATAAAAGCGCTGCCGATTCGCTGATCGTCCTGGCGACCCTGGCTACCGGCGAGATGAATAT  >  NC_000913/308689‑309157
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                         
tACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCCGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCCGCAGCAGTTGGCACAAGCTGAGATcc                                                                                                                                                                                                                              >  1:397434/1‑249 (MQ=255)
          cGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCACAAGCTGAGATCCAGCTGCGCGTTg                                                                                                                                                                                                                  >  2:106697/1‑251 (MQ=255)
            ccgacccccaaatgtgttgcgcgaggtgcagtAACTTTATCAAGGGACAGTTAGCGGGTTGTGTCGCCCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGGAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCACAAGCTGAGATCCAGCTGCGCGTTgtc                                                                                                                                                                                                                <  2:43607/226‑1 (MQ=255)
                 caccacATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCACAAGCTGAGATCCAGCTGCGCGTTgtcagtc                                                                                                                                                                                                            >  1:74548/1‑250 (MQ=255)
                          atgTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGCTAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGATGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCACAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATgg                                                                                                                                                                                                  >  2:198370/1‑251 (MQ=255)
                                                  aTTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCACAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAACGGGGCGt                                                                                                                                                                          >  2:315856/1‑251 (MQ=255)
                                                        cagaagtggggggggttgtttGCCCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACTGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCACAAGCTGAGATCCAGCTGCGCGTTGGCAGTCAGGGGGATGGTAAGGGCCTCATTGAACGGTGCGTTGGCCa                                                                                                                                                                    <  2:392866/232‑1 (MQ=255)
                                                             gATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCACAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAACGGGGCGTTGGCCAGAGcc                                                                                                                                                               >  1:368708/1‑251 (MQ=255)
                                                                                cccGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCACAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAACGGGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGttt                                                                                                                                            <  1:150467/251‑1 (MQ=255)
                                                                                           gCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGTGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCACAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAACGGGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGac                                                                                                                                  >  1:169571/1‑250 (MQ=255)
                                                                                             aaCTGTTTGGTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGGCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTTACCACTAGCTGCAGCAGCAGTTGGCACAAGCTTAGATCCAGCTGCGCGTTGGCAGTCAGGGGGATGGTAAGGGCCTCATTGAACGTGGCGTGGTCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACt                                                                                                                               <  1:395568/251‑1 (MQ=255)
                                                                                             aaCTGGTTGTTAGAGACGCCCAAGTTATAGCTCATAGTACTGCTGCGGGTGTTAACACTGGACTGCCCGAGGTCTTCGCTGCGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCACAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAACGGGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACt                                                                                                                               <  1:213685/251‑1 (MQ=255)
                                                                                               cTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCACAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAACGGGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTgg                                                                                                                             >  1:28865/1‑251 (MQ=255)
                                                                                                 ggTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCTCAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAACGGGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCg                                                                                                                           <  1:44585/251‑1 (MQ=255)
                                                                                                 ggTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCACAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAACGGGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGc                                                                                                                            <  1:139480/250‑1 (MQ=255)
                                                                                                        ataGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCACAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAACGGGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCGTTATCtt                                                                                                                    <  1:66338/251‑1 (MQ=255)
                                                                                                             cGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCACAAGCTGAGGTCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAACGGGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTTACACTGGCGTTATCTTCActc                                                                                                               <  2:383989/251‑1 (MQ=255)
                                                                                                                 cAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCACAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAACGGGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCGTTATCTTCACTCTCtt                                                                                                           >  1:187769/1‑251 (MQ=255)
                                                                                                                    gctataGCTAAGATTACTGCTGAGGGCGTTAACACTGGCATGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCACAAGCTGAGATCCAGCCGCGCGCTGTAAGTCAGGGGGATGGTAAGGGCCTCATTGAACGGGGCGTAGGCCAGAGCCATCAGATGCTGTCGAGTTTGTTCGCTGACACTGGCGTTATCTTCACTCTCTTcc                                                                                                         >  2:276824/3‑250 (MQ=255)
                                                                                                                              aGATTACTGCTGAGGGTGTGAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCACAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAACGGGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCGTTATCTTCACTCTCTTCCAGCTGTAtttt                                                                                              <  1:159667/251‑1 (MQ=255)
                                                                                                                                    ctgctgAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCACAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAACGGGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCGTTATCTTCACTCTCTTCCAGCTGTATTTTCCGGAt                                                                                        >  1:85188/1‑251 (MQ=255)
                                                                                                                                      gctgAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCACAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAACGGGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCGTTATCTTCACTCTCTTCCAGCTGTATTTTCCGGATg                                                                                       <  2:368708/250‑1 (MQ=255)
                                                                                                                                         gAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCACAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAACGGGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCGTTATCTTCACTCTCTTCCAGCTGTATTTTCCGGATGCGt                                                                                    >  1:307039/1‑250 (MQ=255)
                                                                                                                                         gAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCACAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAACGGGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCGTTATCTTCACTCTCTTCCAGCTGTATTTTCCGGATGCGTa                                                                                   <  1:98859/251‑1 (MQ=255)
                                                                                                                                                       tGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCACAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAACGGGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCGTTATCTTCACTCTCTTCCAGCTGTATTTTCCGGATGCGTAGCTGTCCATCGTcc                                                                     >  2:59137/1‑251 (MQ=255)
                                                                                                                                                                 gttGTCTTCGCTACGTGAGCGTACTACGGTTCCCAGCGCTTCGCGCTTGACCGCTAGCTGCAGCAGCAGTTGGCACAAGCGGAGATCCAGCTGCGCGTGGTCAGTCAGGGGGATGGTAGGGGCCTCATTGAACGGGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCGTTATCTTCTCTCTCTTCCAGCTGTATTTTCCGGATGCGTAGCTGTCCATCTTCCAGCCAGATaa                                                           <  2:323582/249‑1 (MQ=255)
                                                                                                                                                                          gCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCACAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAACGGGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCGTTATCTTCACTCTCTTCCAGCTGTATTTTCCGGATGCGTAGCTGTCCATCGTCCAGCCAGATAAAAGCGCTGc                                                  >  2:176429/1‑251 (MQ=255)
                                                                                                                                                                                            ggTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCACAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAACGGGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCGTTATCTTCACTCTCTTCCAGCTGTATTTTCCGGATGCGTAGCTGTCCATCGTCCAGCCAGATAAAAGCGCTGCCGATTCGCTGATCGTCCt                                >  2:95466/1‑251 (MQ=255)
                                                                                                                                                                                                   aGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCACAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAACGGGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCGTTATCTTCACTCTCTTCCAGCTGTATTTTCCGGATGCGTAGCTGTCCATCGTCCAGCCAGATAAAAGCGCTGCCGATTCGCTGATCGTCCTGGCGAcc                         >  1:353201/1‑251 (MQ=255)
                                                                                                                                                                                                            cgcTTGACCACTAGCTGCAGCAGCAGTTGGCACAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAACGGGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCGTTATCTTCACTCTCTTCCAGCTGTATTTTCCGGATGCGTAGCTGTCCATCGTCCAGCCAGATAAAAGCGCTGCCGATTCGCTGATCGTCCTGGCGACCCTGGCTAcc                >  2:397181/1‑251 (MQ=255)
                                                                                                                                                                                                                tGACCACTAGCTGCAGCAGCAGTTGGCACAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAACGGGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCGTTATCTTCACTCTCTTCCAGCTGTATTTTCCGGATGCGTAGCTGTCCATCGTCCAGCCAGATAAAAGCGCTGCCGATTCGCTGATCGTCCTGGCGACCCTGGCTACCGGCg            >  1:73166/1‑251 (MQ=255)
                                                                                                                                                                                                                        aGCTTCAGCAGCAGTTGGCACAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAACGGGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCGGTATCTTCACTCTCTTCCAGCTGTATTTTCCGGATGCGTAGCTGTCCATCGTCCAGCCAGATAAAAGCGCTGCCGATTCGCTGAGCGTCCTGGCGACCCTGGCTACCGGCGATATGaa     <  2:172031/250‑1 (MQ=255)
                                                                                                                                                                                                                            gccgcagcagTTGGCACAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAACGGGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCGTTATCTTCACTCTCTTCCAGCTGTATTTTCCGGATGCGTAGCTGTCCATCGTCCAGCCAGATAAAAGCGCTGCCGATTCGCTGATCGTCCTGGCGACCCTGGCTACCGGCGAGATGAatat  <  2:38161/246‑1 (MQ=255)
                                                                                                                                                                                                                             cagcagcagTTGGCACAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGCAAGGGCCTCATTGAACGGGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCGTTATCTTCACTCTCTTCCAGCTGTATTTTCCGGCTGCGTAGCTGTCCATCGTCCAGCCAGATAAAAGCGCTGCCGATTCGCTGATCGTCCTGGCGACCCTGGCTACCGGCGAGATGAatat  >  1:47843/1‑248 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                         
TACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGACCACTAGCTGCAGCAGCAGTTGGCGCAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAACGGGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCGTTATCTTCACTCTCTTCCAGCTGTATTTTCCGGATGCGTAGCTGTCCATCGTCCAGCCAGATAAAAGCGCTGCCGATTCGCTGATCGTCCTGGCGACCCTGGCTACCGGCGAGATGAATAT  >  NC_000913/308689‑309157

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