Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,464,198 T→A 9.0% T474S (ACC→TCC) ‡ algX ← alginate biosynthesis protein AlgX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,464,1980TA 9.0% 124.9 / 3.4 52T474S (ACC→TCC) ‡algXalginate biosynthesis protein AlgX
Reads supporting (aligned to +/- strand):  ref base T (30/16);  new base A (2/3);  total (32/20)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.45e-01

TTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCG  >  NC_002947/1464051‑1464346
                                                                                                                                                     |                                                                                                                                                    
ttGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTc                                                                                                                                                     >  3:12131/1‑149 (MQ=255)
     cTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTgcgc                                                                                                                                                <  7:129098/149‑1 (MQ=255)
     cTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTgcg                                                                                                                                                 >  8:129098/1‑148 (MQ=255)
        ttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAg                                                                                                                                              <  6:81444/148‑1 (MQ=255)
         tCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCAGACCGCGCAgg                                                                                                                                             <  1:347325/148‑1 (MQ=255)
          cGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCa                                                                                                                                               <  4:13256/145‑1 (MQ=255)
          cGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCa                                                                                                                                               >  3:13256/1‑145 (MQ=255)
          cGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCa                                                                                                                                           >  7:204552/1‑149 (MQ=255)
           gATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCAGACCGCGCAGGCag                                                                                                                                          <  3:91652/149‑2 (MQ=255)
             tCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCCCCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTCAAAAGGGTCAGGGCTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACgg                                                                                                                                        >  5:277315/1‑149 (MQ=255)
             tCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACgg                                                                                                                                        >  8:190642/1‑149 (MQ=255)
                 cTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGaa                                                                                                                                    <  3:231430/149‑1 (MQ=255)
                    gtagtaACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCCGGCACGGCGAAgtt                                                                                                                                 >  1:348764/1‑149 (MQ=255)
                    gtagtaACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAgtt                                                                                                                                 >  8:198812/1‑149 (MQ=255)
                    gtagtaACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCAGACCGCGCAGGCa                                                                                                                                           <  2:151490/139‑1 (MQ=255)
                    gtagtaACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCAGACCGCGCAGGCa                                                                                                                                           >  1:151490/1‑139 (MQ=255)
                       gtaACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGt                                                                                                                                                      <  4:139173/125‑1 (MQ=255)
                       gtaACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGt                                                                                                                                                      >  3:139173/1‑125 (MQ=255)
                            ccTGCGGCGGTACCCGCGGGGGGTGCGGGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGTAATTTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGaa                                                                                                                                    <  3:238845/138‑1 (MQ=255)
                             cTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACATGATGCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGt                                                                                                                          >  1:123845/1‑149 (MQ=255)
                              ggcggggGCACCTGGGGGGCGGGCGGGCAACCACAAAACACGCCCACCGCGCGCGGGGGGGGGCTGTAATTTCGTTTCAAACGGGTCTGGGGCCGCCTCACGGGCCTGCCTGGGGGGTCGGCGCAGCCACGGCGAAGTTGTTGCGTTTg                                                                                                                       <  8:204552/148‑1 (MQ=255)
                                cgggggAACCTGCGCGGCGTGCGCGGCACCACCAAACAAGGCCACCGCGACCAGGGCGGGGCTGGAAATTCGTTTTAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTg                                                                                                                       <  7:327698/147‑1 (MQ=255)
                                        ccAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTgg                                                                                                             >  3:290879/1‑149 (MQ=255)
                                        ccAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTgg                                                                                                             >  3:285772/1‑149 (MQ=255)
                                            cgcgGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGCGTCGCCTCACTGGCCTGCCCGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTc                                                                                                         <  8:142005/149‑1 (MQ=255)
                                                gcgTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTgc                                                                                                                                                  >  4:114819/1‑104 (MQ=255)
                                                gcgTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTgc                                                                                                                                                  <  3:114819/104‑1 (MQ=255)
                                                        ggCACCACCAACCACGCCCAACGCTCGCTGGGGGGGGCTGGAAATTCGTTTCAAACGGGCCTTGGGGCGCCTCAGGGGCCGGCCGGCGCGGGCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGt                                                                                             <  7:198812/149‑1 (MQ=255)
                                                             caccaAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATGAGTTGGAAAAGGGTCATGGGTCGCCTCAATGGCATGCCTGGGCGGTCGGCGAGGGCACGGCGATGTTGTTGCGTTTGCAGCGCTTGGCTTGGCCGTTCTGGTTGCCGc                                                                                        >  6:206568/1‑149 (MQ=255)
                                                              accaAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCt                                                                                       >  3:53051/1‑149 (MQ=255)
                                                               ccaAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCTTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAgttgtt                                                                                                                              <  2:61272/109‑1 (MQ=255)
                                                                caAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTg                                                                                                                       >  3:208563/1‑115 (MQ=255)
                                                                 aaaCAAGGCCAACGCGACCTGGGCGGTGCTGTAATTTGGTTTGAAAAGGGCCATGGGTGGCCTCACTTGCGTGCCTGCGGGGTTTGCGATGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTc                                                                                    <  6:64519/149‑1 (MQ=255)
                                                                     aaGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGGTTGAAAAGGGTCAAGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGCCCTTCTGGGTGCCGCTTTCCGGc                                                                                >  5:270134/1‑149 (MQ=255)
                                                                               cgAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTgg                                                                                               >  7:94100/1‑124 (MQ=255)
                                                                                            gCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGa                                                             >  8:264863/1‑145 (MQ=255)
                                                                                            gCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGa                                                             <  7:264863/145‑1 (MQ=255)
                                                                                                gAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTgg                                                                                               >  7:243263/1‑107 (MQ=255)
                                                                                                     tcttttGAAACGGGCCTTGGGGCCCCTCAGTGGCCTGCCTGCGGGGTCTGCGCTGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCgctggc                                                 <  4:53051/145‑1 (MQ=255)
                                                                                                             aaaGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGt                                        >  6:340060/1‑149 (MQ=255)
                                                                                                               aGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTcc                                      <  7:190642/149‑1 (MQ=255)
                                                                                                                gggTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAgg                                                        >  7:104298/1‑130 (MQ=255)
                                                                                                                   tCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCa                                          >  5:300532/1‑141 (MQ=255)
                                                                                                                     aTGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCt                                >  7:53289/1‑149 (MQ=255)
                                                                                                                         gTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTcgc                            <  5:340060/149‑1 (MQ=255)
                                                                                                                              cacACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTAGCGCATc                       >  3:41665/3‑149 (MQ=255)
                                                                                                                              cTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTccagatgc                       >  8:181769/1‑142 (MQ=255)
                                                                                                                                 aCTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTg                     >  6:264534/1‑148 (MQ=255)
                                                                                                                                  cTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTgg                    >  7:287710/1‑148 (MQ=255)
                                                                                                                                        cgtgcctgcGCAGACCGAGCACGCAGGCCAAGTTGTTGAGTTTGGAGAGCTTGGATGCGAACTTCTGGGTGCCGCTGTACGGCCCCTGTACCGCGACCGCCCGCCGGCGCTGGCAGGCCCAGTCCTCGTCTTCGCGCATCTGGatgcc               >  8:264074/3‑144 (MQ=255)
                                                                                                                                         ctgcctgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCAGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGAc               >  5:204243/1‑148 (MQ=255)
                                                                                                                                          tgcctgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACg              >  3:107997/1‑148 (MQ=255)
                                                                                                                                                cgcgGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGc       >  4:283549/1‑149 (MQ=255)
                                                                                                                                                  cgGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCg     <  8:287710/149‑1 (MQ=255)
                                                                                                                                                     tCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCg  >  7:216603/1‑149 (MQ=255)
                                                                                                                                                     tCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCg  >  5:29350/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCG  >  NC_002947/1464051‑1464346

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