Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,823,386 C→G 9.1% G2617R (GGC→CGC) ‡ pvdL ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,823,3860CG 9.1% 172.1 / 8.4 55G2617R (GGC→CGC) ‡pvdLnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (25/25);  new base G (2/3);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

AGCAGCACCAGATGCGCAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTCCACCAGGCCATCGGTAC  >  NC_002947/4823244‑4823532
                                                                                                                                              |                                                                                                                                                  
agcagcACCAGATGCGCAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCACTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGGGCGGCAAg                                                                                                                                             <  2:100988/150‑1 (MQ=255)
   agcaCCAGATGCGCAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTcc                                                                                                                                          >  2:867336/1‑150 (MQ=255)
     caCCAGATGCGCAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAg                                                                                                                                        >  1:481268/1‑150 (MQ=255)
     caCCAGATGCACAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAg                                                                                                                                        <  2:292903/150‑1 (MQ=255)
               gcAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCacg                                                                                                                              >  2:702725/1‑150 (MQ=255)
                 aGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCacgac                                                                                                                            >  1:767221/1‑150 (MQ=255)
                  ggCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCacgacg                                                                                                                           <  1:158475/150‑1 (MQ=255)
                    cACCATGTAGTCGGGCAGGTGCTTGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGc                                                                                                                         <  2:327811/150‑1 (MQ=255)
                    cACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACGGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGc                                                                                                                         <  1:559862/150‑1 (MQ=255)
                     aCCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGcc                                                                                                                        >  2:1042408/1‑150 (MQ=255)
                           tAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCtcgtcg                                                                                                                  >  1:1066474/1‑150 (MQ=255)
                              tCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGGGCGGCAAgca                                                                                                                                           <  2:715170/122‑3 (MQ=255)
                              tCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGGGCGGCAAgca                                                                                                                                           >  1:715170/1‑120 (MQ=255)
                                gggCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATc                                                                                                             >  1:7637/1‑150 (MQ=255)
                                     ggTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGt                                                                                                        <  2:879485/150‑1 (MQ=255)
                                      gTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTg                                                                                                       <  1:581714/150‑1 (MQ=255)
                                       tGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTgc                                                                                                      <  2:1066474/150‑1 (MQ=255)
                                            gCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAg                                                                                                 >  1:5551/1‑150 (MQ=255)
                                            gCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAg                                                                                                 >  1:987746/1‑150 (MQ=255)
                                                  aTCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGGGCGGCaatcag                                                                                                                                          >  1:1197717/1‑99 (MQ=255)
                                                  aTCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGGGCGGCAAgcag                                                                                                                                          <  2:1197717/103‑4 (MQ=255)
                                                        tttCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCtctc                                                                                     <  2:1169026/150‑1 (MQ=255)
                                                          tCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGa                                                                                   <  2:481268/150‑1 (MQ=255)
                                                                cTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGc                                                                             >  1:179383/1‑150 (MQ=255)
                                                                cTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGc                                                                             >  1:480605/1‑150 (MQ=255)
                                                                       cagcagTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCCGTTc                                                                      >  2:134488/1‑150 (MQ=255)
                                                                          cagTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTc                                                                               <  1:541286/138‑1 (MQ=255)
                                                                          cagTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTc                                                                               >  2:541286/1‑138 (MQ=255)
                                                                                 tgCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAg                                                            <  1:702725/150‑1 (MQ=255)
                                                                                  gCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGc                                                           >  1:1201281/1‑150 (MQ=255)
                                                                                    aCGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCcg                                                         >  2:227801/1‑150 (MQ=255)
                                                                                       tcgtTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCgg                                                      <  2:1150522/150‑1 (MQ=255)
                                                                                           tACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCt                                                  <  1:689860/150‑1 (MQ=255)
                                                                                                   gTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGTGAGCTTGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACctgc                                          <  2:767221/150‑1 (MQ=255)
                                                                                                     ggCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACctgctg                                        <  1:562003/150‑1 (MQ=255)
                                                                                                       cctgcctgCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGt                                      >  1:618465/1‑150 (MQ=255)
                                                                                                        ctgcctgcACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGAt                                                                   >  1:754247/1‑120 (MQ=255)
                                                                                                        ctgcctgcACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGAt                                                                   <  2:754247/120‑1 (MQ=255)
                                                                                                        ctgcctgcACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTc                                     <  2:5551/150‑1 (MQ=255)
                                                                                                        ctgcctgcACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTc                                     >  1:352200/1‑150 (MQ=255)
                                                                                                         tgcctgcACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCg                                    >  1:500796/1‑150 (MQ=255)
                                                                                                             tgcACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCa                                               >  1:98998/1‑135 (MQ=255)
                                                                                                             tgcACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCa                                               <  2:98998/135‑1 (MQ=255)
                                                                                                             tgcACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATGTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCAc                                >  1:378896/1‑150 (MQ=255)
                                                                                                               cacaCCGGAAAGCCTACCAGCTGCGTGCCGCCCGGCAAGGCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCACGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACgg                              <  2:273665/150‑1 (MQ=255)
                                                                                                                    cAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGAc                         <  2:7637/150‑1 (MQ=255)
                                                                                                                       aTAGCCCACCATCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATa                      <  2:513840/150‑1 (MQ=255)
                                                                                                                           cccACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAACACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTcc                  >  1:675122/1‑150 (MQ=255)
                                                                                                                             caccaGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTc                   >  2:254871/1‑147 (MQ=255)
                                                                                                                             caccaGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTc                   <  1:254871/147‑1 (MQ=255)
                                                                                                                               ccaGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTccacca              <  1:1042408/150‑1 (MQ=255)
                                                                                                                                caGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTCCACCAg             <  2:500796/150‑1 (MQ=255)
                                                                                                                                   ctgctTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTCCACCAGGcc          >  1:1030212/1‑150 (MQ=255)
                                                                                                                                    tgctTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTCCACCAGGCCa         <  1:227801/150‑1 (MQ=255)
                                                                                                                                      ctTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACCCCTTAGTCGGCAAACAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTCCACCAGGCCATc       >  2:296470/1‑150 (MQ=38)
                                                                                                                                           cgccCGGCAAGTCCAGCGCCAGCACGATGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTCCACCAGGCCATCGGTAc  >  2:518250/1‑150 (MQ=255)
                                                                                                                                           cgccCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTCCACCAGGCCATCGGTAc  >  2:208727/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
AGCAGCACCAGATGCGCAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTCCACCAGGCCATCGGTAC  >  NC_002947/4823244‑4823532

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