Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,340,805 G→A 6.8% A154T (GCC→ACC)  PP_3814 → polyamine ABC transportersubstrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,340,8050GA 6.8% 173.0 / 2.5 59A154T (GCC→ACC) PP_3814polyamine ABC transportersubstrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (28/27);  new base A (2/2);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

GCGCGACAAGATCGACCCACGCTTCGTCTCCGACCATGGCGTCGGCTCATTCTTCTTCTCCTTCGTGCTCGGCTACAACGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGCCGACGGCGTAGCCCCCGACAAGCTGTACCCGCTGGACCTGGACCGCGCCTTCAAGAAGCTCGACACCATCAAGC  >  NC_002947/4340669‑4340946
                                                                                                                                        |                                                                                                                                             
gcgcGACAAGATCGACCCACGCTTCGTCTCCGACCATGGCGTCGGCTCATTCTTCTTCTCCTTCGTGCTCGGCTACAACGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGc                                                                                                                                    >  8:108446/1‑148 (MQ=255)
gcgcGACAAGATCGACCCACGCTACGTCTCCGACCATGGCGTCGGCTCATTCTTCTTCTCCTTCGTGCTCGGCTACAACGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGc                                                                                                                                    >  3:94339/1‑148 (MQ=255)
             gACCCACGCTTCGTCTCCGACCATGGCGTCGGCTCATTCTTCTTCTCCTTCGTGCTCGGCTACAACGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTa                                                                                                                      <  1:76809/149‑1 (MQ=255)
                  aCGCTTCGTCTCCGACCATGGCGTCGGCTCATTCTTCTTCTCCTTCGTGCTCGGCTACAACGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAAt                                                                                                                 >  7:42615/1‑149 (MQ=255)
                            tcCGACCATGGCGTCGGCTCATTCTTCTTCTCCTTCGTGCTCGGCTACAACGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGccc                                                                                                       >  1:235773/1‑149 (MQ=255)
                            tcCGACCATGGCGTCGGCTCATTCTTCTTCTCCTTCGTGATCGGCTACAAAGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTAGACACAAAGGCCTCCGAAGGCAAGACAGCCGGGTACGACTGGCCGAGccc                                                                                                       >  2:239172/1‑149 (MQ=255)
                                 ccATGGCGTCGGCTCATTCTTCTTCTCCTTCGTGCTCGGCTACAACGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCg                                                                                                  <  2:208318/149‑1 (MQ=255)
                                 ccATGGCGTCGGCTCATTCTTCTTCTCCTTCGTGCTCGGCTACAACGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCg                                                                                                  >  7:89125/1‑149 (MQ=255)
                                   aTGGCGTCGGCTCATTCTTCTTCTCCTTCGTGCTCGGCTACAACGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGcc                                                                                                                           <  4:300336/122‑1 (MQ=255)
                                   aTGGCGTCGGCTCATTCTTCTTCTCCTTCGTGCTCGGCTACAACGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGcc                                                                                                                           >  3:300336/1‑122 (MQ=255)
                                      gCGTCGGCTCATTCTTCTTCTCCTTCGTGCTCGGCTACAACGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCt                                                                                              <  6:353845/148‑1 (MQ=255)
                                                  tcttcttctCCTTCGTGCTCGGCTACAACGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGt                                                                                                                       >  4:126506/1‑111 (MQ=255)
                                                  tcttcttctCCTTCGTGCTCGGCTACAACGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGt                                                                                                                       <  3:126506/111‑1 (MQ=255)
                                                   cttcttctCCTTCGTGCTCGGCTACAACGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCctg                                                                                 <  4:94339/148‑1 (MQ=255)
                                                    ttcttctCCTTCGTGCTCGGCTACAACGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCctgct                                                                               >  8:20056/1‑149 (MQ=255)
                                                     tcttctCCTTCGTGCTCGGCTACAACGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCctgctg                                                                              >  6:172720/1‑149 (MQ=255)
                                                     tcttctCCTTCGTACTCGGCTACAACGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCctgctg                                                                              >  1:298105/1‑149 (MQ=255)
                                                         ctcCTTCGTGCTCGGCTACAACGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGc                                                                            <  3:324599/147‑1 (MQ=255)
                                                            cTTCGTGCTCGGCTACAACGAAGGCAAGCTCGGGGCCCACAAGCCGGGGGACTGGCCCGCGCGGTTCGACACCACGGCCTACCCCCGCAAGCGCGCCCGGTACAAATGGCCCCGCCCCGGCGGGCGCGCACTGGCCGAGCGGGccgccg                                                                       >  8:337732/1‑149 (MQ=255)
                                                            cTTCGTGCTCGGCTACAACGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGCCGACg                                                                       >  3:115541/1‑149 (MQ=255)
                                                             ttCGTGCTCGGCTACAACGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGCCGACgg                                                                      >  5:141183/1‑149 (MQ=255)
                                                               cGTGCTCGGCTACAACGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGCCGACgg                                                                      <  1:51031/147‑1 (MQ=255)
                                                                gtacTCGGCTACAACGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGCCGACGGCGt                                                                   <  2:298105/146‑1 (MQ=255)
                                                                 tGCTGGGCTACAACGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGcc                                                                                    <  3:113057/131‑1 (MQ=255)
                                                                 tGCTCGGCTACAACGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACt                                                                                        >  7:200529/1‑127 (MQ=255)
                                                                 tGCTCGGCTACAACGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACt                                                                                        <  8:200529/127‑1 (MQ=255)
                                                                 tGCTCGGCTACAACGAAGGCAAGCTAGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGcc                                                                                    >  4:113057/1‑131 (MQ=255)
                                                                    tCGGCTACAACGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGCCGACGGCGTAGcc                                                               >  2:57824/1‑149 (MQ=255)
                                                                    tCGGCTACAACGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGCCGACGGCGTAGc                                                                >  8:194988/1‑148 (MQ=255)
                                                                      ggCTACAACGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGCCGACGGCGTAGc                                                                <  7:194988/146‑1 (MQ=255)
                                                                          aCAACGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGACCTACCCCGGCAAGCGAGCCCTGTACAAATGTACCAGCCCCGGCCTGCTGGAACTGGCCCTGCTGGAAGAAGGCGTAGCCCCCGa                                                          >  8:255164/1‑148 (MQ=255)
                                                                          aCAACGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGACCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGCCGACGGCGTAGCCCCCGa                                                          >  2:203947/1‑148 (MQ=255)
                                                                            aaCGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGc                                                                            <  8:145397/128‑1 (MQ=255)
                                                                            aaCGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGc                                                                            >  7:145397/1‑128 (MQ=255)
                                                                            aaCGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGCCGACGGCGTAGCCCCCGACaa                                                       >  8:336469/1‑149 (MQ=255)
                                                                                   gCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGCCGACGGCGTAGCCCCCGACAAGCTGTAc                                                <  6:141183/149‑1 (MQ=255)
                                                                                   gCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGCCGACGGCGTAGCCCCCGACAAGCTGTAc                                                >  7:243700/1‑149 (MQ=255)
                                                                                     aaGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGc                                                                                                                                              <  8:2121/53‑1 (MQ=255)
                                                                                     aaGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGc                                                                                                                                              >  7:2121/1‑53 (MQ=255)
                                                                                          cGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTgg                                                                                      <  1:187637/104‑1 (MQ=255)
                                                                                          cGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTgg                                                                                      >  2:187637/1‑104 (MQ=255)
                                                                                           ggggCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGCCGACGGCGTAGCCCCCGACAAGCTGTACCCGctg                                          <  7:20056/147‑1 (MQ=255)
                                                                                           ggggCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGCCGACGGCGTAGCCCCCGACAAGCTGTACCCGctg                                          <  7:336469/147‑1 (MQ=255)
                                                                                            gggCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGCCGACGGCGTAGCCCCCGACAAGCTGTACCCGctggac                                       <  7:108446/149‑1 (MQ=255)
                                                                                              gCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGCCGACGGCGTAGCCCCCGACAAGCTGTACCCGctggacct                                     <  8:89125/149‑1 (MQ=255)
                                                                                               ccAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGCCGACGGCGTAGCCCCCGACAAGCTGTACCCGctggacctg                                    >  3:67931/1‑149 (MQ=255)
                                                                                                  acaaGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGCCGACGGCGTAGCCCCCGACAAGCTGTACCCGctggacctggac                                 <  6:293790/149‑1 (MQ=255)
                                                                                                     aGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGCCGACGGCGTAGCCCCCGACAAGCTGTACCCGCTGGACCTGGACcgc                              <  1:57824/149‑1 (MQ=255)
                                                                                                          gTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGCCGACGGCGTAGCCCCCGACAAGCTGTACCCGCTGGACCTGGACCGCGCCtt                         >  7:109634/1‑149 (MQ=255)
                                                                                                          gTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGCCGACGGCGTAGCCCCCGACAAGCTGTACCCGCTGGACCTGGACCGCGCCtt                         <  6:344419/149‑1 (MQ=255)
                                                                                                             gactggacCGCGCTGTTCGACACCAAGACCTACCCCGGCAAGCTAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGCCGACGGCGTAGCCCCCGACAAGCTGTACCCGCTGGACCTGGACCGCGCCTTCaa                      <  7:255164/149‑1 (MQ=255)
                                                                                                             gactggacCGCGCTGTTCGACACCAAGACCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGCCGACGGCGTAGCCCCCGACAAGCTGTACCCGCTGGACCTGGACCGCGCCTTCaa                      <  1:203947/149‑1 (MQ=255)
                                                                                                                     cgcgCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGCCGACGGCGTAGCCCCCGACAAGCTGTACCCGCTGGACCTGGACCGCGCCTTCAAGAAGc                 <  6:319565/146‑1 (MQ=255)
                                                                                                                     cgcgCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGCCGACGGCGTAGCCCCCGACAAGCTGTACCCGCTGGACCTGGACCGCGCCTTCAAGAAGc                 >  5:319565/1‑146 (MQ=255)
                                                                                                                      gcgcTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGCCGACGGCGTAGCCCCCGACAAGCTGTACCCGCTGGACCTGGACCGCGCCTTCAAGAAGCTCGa             <  5:172720/149‑1 (MQ=255)
                                                                                                                       cgcTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGCCGACGGCGTAGCCCCCGACAAGCTGTACCCGCTGGACCTGGACCGCGCCTTCAAGAAGCTCGa             >  6:204812/1‑148 (MQ=255)
                                                                                                                              cGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGCCGACGGCGTAGCCCCCGACAAGCTGTACCCGCTGGACCTGGACCGCGCCTTCAAGAAGCTCGACACCATc      >  8:258406/1‑148 (MQ=255)
                                                                                                                               gaaacCAACGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCAAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGCCGACGGCGTAGCCCCCGACAAGCTGTACCCGCTGGACCTGGACCGCGCCTTCAAGAAGCTCGACACCATcaa    <  4:115541/146‑1 (MQ=255)
                                                                                                                                 cacCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGCCGACGGCGTAGCCCCCGACAAGCTGTACCCGCTGGACCTGGACCGCGCCTTCAAGAAGCTCGACACCATcaagc  <  2:235773/149‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                             
GCGCGACAAGATCGACCCACGCTTCGTCTCCGACCATGGCGTCGGCTCATTCTTCTTCTCCTTCGTGCTCGGCTACAACGAAGGCAAGCTCGGGGCCAACAAGCCGGTGGACTGGACCGCGCTGTTCGACACCAAGGCCTACCCCGGCAAGCGAGCCCTGTACAAATGGCCCAGCCCCGGCGTGCTGGAACTGGCCCTGCTGGCCGACGGCGTAGCCCCCGACAAGCTGTACCCGCTGGACCTGGACCGCGCCTTCAAGAAGCTCGACACCATCAAGC  >  NC_002947/4340669‑4340946

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