Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,787,508 C→T 7.7% G2111D (GGT→GAT)  pvdI ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,787,5080CT 7.7% 149.7 / 2.7 52G2111D (GGT→GAT) pvdInon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (24/24);  new base T (2/2);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

GCCTTCTTCGGTCAGGGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGTTTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGCCCAGATCGCTGCGACCTTGCATTCGGTATCGCTGTACGGCGGCACATAGCGCACCTGAAGCTGGGCCAGATCGGGTTGCGGCAGCG  >  NC_002947/4787360‑4787652
                                                                                                                                                    |                                                                                                                                                
gCCTTCTTCGGTCAGGGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGTTTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGac                                                                                                                                                  >  8:251975/1‑149 (MQ=255)
  cttcttCGGTCAGGGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGTTTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGacca                                                                                                                                                >  5:186198/1‑149 (MQ=255)
    ccggttctgcgggcctttgcccaggcgctggaaggcgcgggccaggcccGCAGAACCGGTGTTTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGaccacc                                                                                                                                              <  1:240671‑M2/100‑1 (MQ=255)
     cttcGGTCAGGGCCTAGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGTTTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGaccac                                                                                                                                               <  5:109378/147‑1 (MQ=255)
        ttctgcgggcctttgcccaggcgctggaaggcgcgggccaggcccGCAGAACCGGTGTTTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGt                                                                                                                                          >  2:48120‑M2/46‑149 (MQ=255)
        cGGTCAGGGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGTTTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGATCACCCAGt                                                                                                                                          >  8:246814/1‑149 (MQ=255)
            gcgggcctttgcccaggcgctggaaggcgcgggccaggcccGCAGAACCGGTGTTTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCa                                                                                                                                       >  6:171580‑M2/42‑148 (MQ=255)
             aGGGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGTTTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGATCACCCAGTTCaaa                                                                                                                                     <  7:246814/149‑1 (MQ=255)
                  CTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCAGCAGAACCGGTGTTTCGAAGATGAGCCGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACGACCCAGTTCAAAGAAGt                                                                                                                                <  7:82897‑M1/149‑1 (MQ=255)
                      cccGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGTTTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTgtc                                                                                                                            <  7:111662/149‑1 (MQ=255)
                             cTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGTTTCGAAGATGAGACCGATCGGCATAACGATGCCCAGTTTGCTGTTGATACGTGAGAGGGCCTGGGCCGCCCCCAGCGAGTGACCCCCCCGTTCAAAGAAGTTGTCGTCCAGc                                                                                                                     >  8:32272/1‑149 (MQ=255)
                             cTTCCAGCGCCTGGGCAAAGGCCCGCAGAACAGGTGTTTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTCGCTGTTGATACGTGAGAGGTCCTGGGCCGCCAGCAGAGAGTGACCACCCAGTTCAAATAAGTTGTCGACCAGc                                                                                                                     >  8:146925/1‑149 (MQ=255)
                                  aGCGCCTGGGCAAAGGCCCGCAGAACCGGTGTTTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGAc                                                                                                                >  4:112159/1‑149 (MQ=255)
                                      ccTGGGCAAAGGCCCGCAGAACCGGTGTTTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGc                                                                                                            <  7:32272/149‑1 (MQ=255)
                                       cTGGGCAAAGGCCCGCAGAACCGGTGTTTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCt                                                                                                           <  5:288367/149‑1 (MQ=255)
                                          ggCAAAGGCCCGCAGAACCGGTGTTTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGa                                                                                                        >  3:74925/1‑149 (MQ=255)
                                           gCAAAGGCCCGCAGAACCGGTGTTTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGAt                                                                                                       <  7:251975/149‑1 (MQ=255)
                                           gCAAAGGCCCGCAGAACCGGTGTTTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGAt                                                                                                       <  7:138350/149‑1 (MQ=255)
                                            cAAAGGCCCGCAGAACCGGTGTTTCGAAGATGAGCCGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGCCCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATg                                                                                                      >  2:54371/1‑149 (MQ=255)
                                            cAAAGGCCCGCAGAACCGGTGTTTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTTTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATg                                                                                                      <  2:114938/149‑1 (MQ=255)
                                                ggcccGCAGAACCGGTGTTTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATgcg                                                                                                    >  2:117953‑M2/6‑147 (MQ=255)
                                                 gcccGCAGAACCGGTGTTTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAg                                                                                                 <  6:19640‑M2/145‑1 (MQ=255)
                                                   ccGCAGAACCGGTGGTTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCGGACCTGCTCGATGCGCAGGa                                                                                               >  2:31822‑M2/3‑149 (MQ=255)
                                                                 tttCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGcc                                                                                       >  3:159783/1‑143 (MQ=255)
                                                                  ttCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGCCCAGATCg                                                                                >  8:204886/1‑149 (MQ=255)
                                                                     gAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGCCCAGATCGCTg                                                                             >  4:115063/1‑149 (MQ=255)
                                                                            agaCGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGATCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGCCCAGATCGCTGCGACCtt                                                                      >  7:53394/1‑149 (MQ=255)
                                                                             gaCGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGCCCAGATCGCTGCGACCTTg                                                                     <  4:159783/149‑1 (MQ=255)
                                                                                   tCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGCCCAGATCGCTGCGACCTTGCATTCg                                                               <  3:100310/149‑1 (MQ=255)
                                                                                      gCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCTTCCAGCCCGACCTGCTCGCTGCGCAGGACATCGGCCCAGATCGCTGCGACCTTGCATTCGGTa                                                            >  8:262462/1‑149 (MQ=255)
                                                                                      gCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGCCCAGATCGCTGCGACCTTGCATTCGGTa                                                            >  2:219064/1‑149 (MQ=255)
                                                                                       cATATCGATGCCCAGTTGTCTGTTGATACGTTAGAGGGCCTGGTCCGCCAGCAGCGAGTGACCAGCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCTGGACATCGTCCCAGATCGCTGCGACCTTGCATTCGGTAt                                                           <  7:146925/149‑1 (MQ=255)
                                                                                       cATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGCCCAGATCGCTGCGACCTTGCATTCGGTAt                                                           <  1:117953/149‑1 (MQ=255)
                                                                                         tatCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGCCCAGATCGCTGCGACCTTGCATTCGGTATCg                                                         <  3:115063/149‑1 (MQ=255)
                                                                                              aTGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGCCCAGATCGCTGCGACCTTGCATTCGGTATCGCTGTa                                                    <  3:112159/149‑1 (MQ=255)
                                                                                               tGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGCCCAGATCGCTGCGACCTTGCATTCGGTATCGCTGTAc                                                   <  5:171580/149‑1 (MQ=255)
                                                                                               tGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGCCCAGATCGCTGCGACCTTGCATTCGGTATCGCTGTAc                                                   <  7:204886/149‑1 (MQ=255)
                                                                                                gCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGCCCAGATCGCTGCGACCTTGCATTCGGTATCGCTGTAcg                                                  >  5:17268/1‑149 (MQ=255)
                                                                                                  ccAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGCCCAGATCGCTGCGACCTTGCATTCGGTATCGCTGTAcggc                                                >  1:65597/1‑149 (MQ=255)
                                                                                                  ccAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGCCCAGATCGCTGCGACCTTGCATTCGGTATCGCTGTAcggc                                                >  7:20348/1‑149 (MQ=255)
                                                                                                   cAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGCCCAGATCGCTGCGACCTTGCATTCGGTATCGCTGTAcggcg                                               <  6:186198/149‑1 (MQ=255)
                                                                                                    aGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGATCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGCCCAGATCGCTGCGACCTTGCATTCGGTATCGCTGTAcggcgg                                              <  8:53394/149‑1 (MQ=255)
                                                                                                    aGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGCCCAGATCGCTGCGACCTTGCATTCGGTATCGCTGTAcggcgg                                              <  2:162050/149‑1 (MQ=255)
                                                                                                         gCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGCCCAGATCGCTGCGACATTGCATTCGGTATCGCTGTACGGCGGCACAt                                         <  1:31822/149‑1 (MQ=255)
                                                                                                            gTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGCCCAGATCACTGCGACCTTGCATTCGGTATCGCTGTACGGCGGCACATAgc                                      >  1:175393/1‑149 (MQ=255)
                                                                                                             ttGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGCCCAGATCGCTGCGACCTTGCATTCGGTATCGCTGTACGGCGGCACATAgc                                      <  8:20348/148‑1 (MQ=255)
                                                                                                                aTACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGCCCAGATCGCTGCGACCTTGCATTCGGTATCGCTGTACGGCGGCACATAGCGCAc                                  >  4:127977/1‑149 (MQ=255)
                                                                                                                     tGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGCCCAGATCGCTGCGACCTTGCATTCGGTATCGCTGTACGGCGGCACATAGCGCACCTGaa                             <  1:48120/149‑1 (MQ=255)
                                                                                                                          gggCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGCCCAGATCGCTGCGACCTTGCATTCGGTATCGCTGTACGGCGGCACATAGCGCACCTGAAGCTgg                        >  6:245083/1‑149 (MQ=255)
                                                                                                                          gggCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGCCCAGATCGCTGCGACCTTGCATTCGGTATCGCTGTACGGCGGCACATAGCGCACCTGAAGCTgg                        <  8:142382/149‑1 (MQ=255)
                                                                                                                                    cgccAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGCCCAGATCGCTGCGACCTTGCATTCGGTATCGCTGTACGGCGGCACATAGCGCACCTGAAGCTGGGCCAGATCgg              <  2:65597/149‑1 (MQ=255)
                                                                                                                                    cgccAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGCCCAGATCGCTGCGACCTTGCATTCGGTATCGCTGTACGGCGGCACATAGCGCACCTGAAGCTGGGCCAGATCgg              >  2:263167/1‑149 (MQ=255)
                                                                                                                                          cagcGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGCCCAGATCGCTGCGACCTTGCATTCGGTATCGCTGTACGGCGGCACATAGCGCACCTGAAGCTGGGCCAGATCGGGTTGCg        >  5:212825/1‑149 (MQ=255)
                                                                                                                                           agcGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGCCCAGATCGCTGCGACCTTGCATTCGGTATCGCTGTACGGCGGCACATAGCGCACCTGAAGCTGGGCCAGATCGGGTTGCgg       <  3:127977/149‑1 (MQ=255)
                                                                                                                                            gcGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGCCCAGATCGCTGCGACCTTGCATTCGGTATCGCTGTACGGCGGCACATAGCGCACCTGAAGCTGGGCCAGATCGGGTTGCgg       >  1:126297/1‑148 (MQ=255)
                                                                                                                                                gTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGCCCAGATCGCTGCGACCTTGCATTCGGTATCGCTGTACGGCGGCACATAGCGCACCTGAAGCTGGGCCAGATCGGGTTGCGGCAgcg  >  6:88360/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GCCTTCTTCGGTCAGGGCCTGGCCCGCGCCTTCCAGCGCCTGGGCAAAGGCCCGCAGAACCGGTGTTTCGAAGATGAGACGGATCGGCATATCGATGCCCAGTTGGCTGTTGATACGTGAGAGGGCCTGGGCCGCCAGCAGCGAGTGACCACCCAGTTCAAAGAAGTTGTCGTCCAGCCCGACCTGCTCGATGCGCAGGACATCGGCCCAGATCGCTGCGACCTTGCATTCGGTATCGCTGTACGGCGGCACATAGCGCACCTGAAGCTGGGCCAGATCGGGTTGCGGCAGCG  >  NC_002947/4787360‑4787652

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