Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,266,910 G→C 8.6% R360G (CGC→GGC)  quiP ← acyl‑homoserine lactone acylase QuiP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,266,9100GC 8.6% 162.7 / 5.1 56R360G (CGC→GGC) quiPacyl‑homoserine lactone acylase QuiP
Reads supporting (aligned to +/- strand):  ref base G (25/25);  new base C (3/2);  total (28/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 = 4.90e-01

GATCACCTTTGAACTGGGGCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGACCGCGCTGCTGCTCCAGGCCAGCTTGCCATTGTAGCCGGCCAGCA  >  NC_002947/1266772‑1267053
                                                                                                                                          |                                                                                                                                                 
gATCACCTTTGAACTGGGGCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGcc                                                                                                                                         >  8:245655/1‑149 (MQ=255)
gATCACCTTTGAACTGGGGCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGcc                                                                                                                                         >  4:179288/1‑149 (MQ=255)
gATCACCTTTGAACTGGGGCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGcc                                                                                                                                         >  5:344719/1‑149 (MQ=255)
gATCACCTTTGAACTGGGGCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGcc                                                                                                                                         <  2:236281/149‑1 (MQ=255)
      cTTTGAACTGGGGCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGc                                                                                                                                   >  2:187699/1‑149 (MQ=255)
      cTTTGAACTGGGGCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCCGGCCCGCGCCCGCAGc                                                                                                                                   <  4:5347/149‑1 (MQ=255)
             cTGGGGCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTAc                                                                                                                            >  6:216545/1‑149 (MQ=255)
                gggCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTAcc                                                                                                                           >  5:204161/1‑147 (MQ=255)
                gggCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTAcc                                                                                                                           <  6:204161/147‑1 (MQ=255)
                gggCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCAt                                                                                                                         >  2:283638/1‑149 (MQ=255)
                  gCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCg                                                                                                                       >  2:267738/1‑149 (MQ=255)
                  gCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGGGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCg                                                                                                                       <  2:175568/149‑1 (MQ=255)
                  gCAGGTGCCGCGCCAGGCCCAGCCTGGGGTTTTCGGGCGGGCCCAGCAGGGGGCCGGGGCGGGGGTCGTACCTCACTTCGCGCACGGGGCGTTGCCCACGTACAAAGTTGGTTTCGCTGCGGGCGCGGGCGGGCAGCCATTTACCATCg                                                                                                                       <  6:293738/149‑1 (MQ=255)
                   gaGGTTCAGCGCCAGGCACAGGCTTGCGTTTTCTGGCTGGGCCAGCAGGGTGGCGTGGCTGGTGTAGTTCATCACTTAGCGCACGGGGCGTTGGCCACGTACATAGTAGGTTTCGGTTCTGGCGAGGGCCGTCAGCCATTTACGATCgg                                                                                                                      <  6:229804/148‑1 (MQ=255)
                     ggTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGcc                                                                                                                    <  1:283638/149‑1 (MQ=255)
                      gTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGTGTGCCGTGGCGGGAGTAGTACATCACTTCGCGCACGGGGCGTTGGCGACGTGCAAAGGAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATATACCCTGGGc                                                                                                                     >  6:207905/1‑147 (MQ=255)
                           gcgcCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAg                                                                                                             <  1:267738/149‑1 (MQ=255)
                             gcCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCt                                                                                                           <  3:168373/149‑1 (MQ=255)
                                     caggcTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATGTACCATCGGCCAGGTAGCTCAGCTgcc                                                                                                   >  6:303376/1‑147 (MQ=255)
                                     caggcTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCCGGCCCGCGCCCGCa                                                                                                                                     <  5:38004/116‑1 (MQ=255)
                                     caggcTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCCGGCCCGCGCCCGCa                                                                                                                                     >  6:38004/1‑116 (MQ=255)
                                        gcTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCTCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTTACTGc                                                                                                <  4:244608/149‑1 (MQ=255)
                                              gTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTggc                                                                                          >  2:244322/1‑149 (MQ=255)
                                              gTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTggc                                                                                          <  3:250593/149‑1 (MQ=255)
                                               ttttCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTggcg                                                                                         >  5:140113/1‑149 (MQ=255)
                                                          ggCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAgctg                                                                                <  1:244322/147‑1 (MQ=255)
                                                           gCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTc                                                                             >  3:202947/1‑149 (MQ=255)
                                                                    gTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAg                                                                   <  3:196675/149‑1 (MQ=255)
                                                                    gTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAg                                                                   <  2:114031/149‑1 (MQ=255)
                                                                       ccGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTc                                                                <  3:179288/149‑1 (MQ=255)
                                                                       ccGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTc                                                                >  3:340236/1‑149 (MQ=255)
                                                                       ccGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTc                                                                <  4:340236/149‑1 (MQ=255)
                                                                           ggCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGa                                                            <  4:202947/149‑1 (MQ=255)
                                                                            gCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGAt                                                           >  2:273046/1‑149 (MQ=255)
                                                                                 tgtCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCg                                                      <  5:216545/149‑1 (MQ=255)
                                                                                 tgtCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCg                                                      <  6:140113/149‑1 (MQ=255)
                                                                                  gtCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGa                                                            <  2:155180/142‑1 (MQ=255)
                                                                                  gtCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGa                                                            >  1:155180/1‑142 (MQ=255)
                                                                                   tCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGc                                                    >  8:138657/1‑149 (MQ=255)
                                                                                         aTCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGAc                                              <  1:187699/149‑1 (MQ=255)
                                                                                            aCTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGACcgc                                           >  5:197664/1‑149 (MQ=255)
                                                                                                         ggCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGACCGCGCTGCTGCTCCAg                              <  7:245655/149‑1 (MQ=255)
                                                                                                                 ccACGTACAAAGTAGGTTTCGCTGCGTGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTc                                                       <  5:368412/116‑1 (MQ=255)
                                                                                                                 ccACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTc                                                       >  6:368412/1‑116 (MQ=255)
                                                                                                                 ccACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGACCGCGCTGCTGCTCCAGGCCAGCtt                      <  7:138657/149‑1 (MQ=255)
                                                                                                                      tACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTATACTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGACCGCGCTGCTGCTCCAGGCCAGCTTGCCa                  >  1:98363/1‑149 (MQ=255)
                                                                                                                       aCAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGACCGCGCTGCTGCTCCAGGCCAGCTTGCCAtt                >  3:221494/1‑149 (MQ=255)
                                                                                                                          aaGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACACGGTCCTGATTGTCGGCCATGACCGCGCTGCTGCTCCAGGCCAGCTTGCCATTGt              <  1:273046/149‑1 (MQ=255)
                                                                                                                           aGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGACCGCGCTGCTGCTCCAGGCCAGCTTGCCATTGTAg            >  1:217198/1‑149 (MQ=255)
                                                                                                                             tAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGACCGCGCTGCTGCTCCAGGCCAGCTTGCCATTGTAGcc          >  1:106434/1‑149 (MQ=255)
                                                                                                                               ggTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGACCGCGCTGCTGCTCCAGGCCAGCTTGCCATTGTAGCCgg        <  4:221494/149‑1 (MQ=255)
                                                                                                                                 tttCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGACCGCGCTGCTGCTCCAGGCCAGCTTGCCATTGTAGCCGGc       >  8:134645/1‑148 (MQ=255)
                                                                                                                                  ttCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGACCGCGCTGCTGCTCCAGGCCAGCTTGCCATTGTAGCCGGc       >  1:330833/1‑147 (MQ=255)
                                                                                                                                     gCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGACCGCGCTGCTGCTCCAGGCCAGCTTGCCATTGTAGCCGGCCAGca  <  6:197664/149‑1 (MQ=255)
                                                                                                                                     gCTGCCGGCCCGCGCCCGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGACCGCGCTGCTGCTCCAGGCCAGCTTGCCATTGTAGCCGGCCAGca  >  6:13023/1‑149 (MQ=255)
                                                                                                                                      cTGCCGGCCCGCGCCCGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGACCGCGCTGCTGCTCCAGGCCAGCTTGCCATTGTAGCCGGCCAGca  >  3:350549/1‑148 (MQ=255)
                                                                                                                                          |                                                                                                                                                 
GATCACCTTTGAACTGGGGCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGACCGCGCTGCTGCTCCAGGCCAGCTTGCCATTGTAGCCGGCCAGCA  >  NC_002947/1266772‑1267053

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: