Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,223,223 G→T 11.0% intergenic (‑530/+30) PP_2823 ← / ← PP_2824 methyl‑accepting chemotaxis transducer/TetR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,223,2230GT11.0% 152.7 / 8.1 54intergenic (‑530/+30)PP_2823/PP_2824methyl‑accepting chemotaxis transducer/TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (20/28);  new base T (4/2);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.04e-01

GACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGCATCATGAGCCAGCCATCAC  >  NC_002947/3223075‑3223371
                                                                                                                                                    |                                                                                                                                                    
gACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGGTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAg                                                                                                                                                      <  3:64158/149‑1 (MQ=255)
  cGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAggg                                                                                                                                                    <  5:278590/149‑1 (MQ=255)
              agagCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGt                                                                                                                                        <  5:72056/149‑1 (MQ=255)
                  ccTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTAc                                                                                                                                    >  7:108929/1‑149 (MQ=255)
                     gCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCgg                                                                                                                                 >  8:265052/1‑149 (MQ=255)
                     gCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCgg                                                                                                                                 >  3:48782/1‑149 (MQ=255)
                      cAAATAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGTAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCgg                                                                                                                                 <  6:246168/148‑1 (MQ=255)
                      cAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCgg                                                                                                                                 <  5:335851/148‑1 (MQ=255)
                      cAAAAAGCAGGGCAGGCGATATCGGGGCCTGGTGCTTGGGGTGGATCGGCAAGGCGCTGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCgg                                                                                                                                 <  7:300650/148‑1 (MQ=255)
                       aaaaaGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGGTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCgg                                                                                                                                 >  7:301281/1‑147 (MQ=255)
                          aaGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGc                                                                                                                            >  3:162829/1‑149 (MQ=255)
                              agggcaggCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACt                                                                                                                        >  1:35445/1‑149 (MQ=255)
                               gggcaggCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCCTGGAAGgc                                                                                                                                               >  5:315188/1‑125 (MQ=255)
                               gggcaggCGCTATCGGGGCCTGGGGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCCTGGAAGgc                                                                                                                                               <  6:315188/125‑1 (MQ=255)
                                ggcaggCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACt                                                                                                                        >  4:106730/1‑147 (MQ=255)
                                 gcaggCGCTATCGGGGCCGGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAg                                                                                                                     <  5:87543/149‑1 (MQ=255)
                                   aggCGCTATCGGGGCCGGGGGCTTGGGTGTGCGCGGCAAGGCTCGGGGCGGCAGTTGCGCGACGACGGAAGGCGTCGAAGACACGATCGGTTGTGTGGGGGGGCCGCCTTCCAGGGAAGGCGACGGGTCTACCGGGAGGCGACTTAGGa                                                                                                                   <  2:209663/149‑1 (MQ=255)
                                     gcgcTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACt                                                                                                                 >  6:267714/1‑149 (MQ=255)
                                       gcTATCGGGGCGTGGTGCTGGGGGTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACTCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTg                                                                                                               <  4:195290/149‑1 (MQ=255)
                                        cTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGATGCGACTTAGGACTTgg                                                                                                              <  7:265052/149‑1 (MQ=255)
                                              gggCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGttt                                                                                                        >  4:323628/1‑149 (MQ=255)
                                                gCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTg                                                                                                      <  2:314367/149‑1 (MQ=255)
                                                 ccTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTgg                                                                                                     <  8:108929/149‑1 (MQ=255)
                                                      tGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGga                                                                                                <  2:223262/149‑1 (MQ=255)
                                                        cTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGgag                                                                                               >  7:299464/1‑148 (MQ=255)
                                                         ttGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGgaga                                                                                              >  6:213416/1‑148 (MQ=255)
                                                            ggTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATc                                                                                          <  3:323628/149‑1 (MQ=255)
                                                             gTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCa                                                                                         <  4:252774/149‑1 (MQ=255)
                                                              ttGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCAt                                                                                        <  8:75933/149‑1 (MQ=255)
                                                                   tCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTgg                                                                                                     <  2:256723/131‑1 (MQ=255)
                                                                   tCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTgg                                                                                                     >  1:256723/1‑131 (MQ=255)
                                                                      gCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTg                                                                                                               >  2:94432/1‑118 (MQ=255)
                                                                      gCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTg                                                                                                               <  1:94432/118‑1 (MQ=255)
                                                                        aaGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAg                                                                              <  4:48782/149‑1 (MQ=255)
                                                                          ggCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGc                                                                            <  5:267714/149‑1 (MQ=255)
                                                                                    ggcAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAg                                                                  <  2:35445/149‑1 (MQ=255)
                                                                                                     ggAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTgc                                                 >  1:228261/1‑149 (MQ=255)
                                                                                                          gcggcgAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTaa                                            <  8:299464/149‑1 (MQ=255)
                                                                                                             gcgAAGACACGTGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAACGCGCCGGGTCTACCGGGATGCGATTTATGACTTGGCTGTTTTGGGTGGAGAAATCAAGACAATGAACGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACgg                                         <  2:259966/149‑1 (MQ=255)
                                                                                                                  gACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCg                                    <  8:301281/149‑1 (MQ=255)
                                                                                                                    cacGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGa                                  <  2:228261/149‑1 (MQ=255)
                                                                                                                          cGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGa                             >  1:27121/1‑148 (MQ=255)
                                                                                                                            gTTTGGTGGGGGGGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAg                          <  6:132352/149‑1 (MQ=255)
                                                                                                                             tttGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAgg                         >  5:331960/1‑149 (MQ=255)
                                                                                                                               tggtggGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGgc                        <  4:162829/148‑1 (MQ=255)
                                                                                                                                      ggTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGcatcat                <  5:213416/149‑1 (MQ=255)
                                                                                                                                      ggTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGcatcat                >  6:31164/1‑149 (MQ=255)
                                                                                                                                       ccggcGCCTTCCCTGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGcatcat                >  7:269033/5‑148 (MQ=255)
                                                                                                                                       ccggcGCCTTCCCTGGAAGACGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGcatcat                >  2:180478/5‑148 (MQ=255)
                                                                                                                                        tGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGCATCATGa              <  5:252486/149‑1 (MQ=255)
                                                                                                                                            gCCTTCCCTGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGga                                                                                                >  2:261756/1‑63 (MQ=255)
                                                                                                                                            gCCTTCCCTGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGga                                                                                                <  1:261756/63‑1 (MQ=255)
                                                                                                                                               ttCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGCATCATGagccagc        >  5:306690/1‑148 (MQ=255)
                                                                                                                                               ttCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGCATCATGagccagc        >  3:167252/1‑148 (MQ=255)
                                                                                                                                                  cAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCAt                                                                                        <  2:251552/65‑1 (MQ=255)
                                                                                                                                                  cAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCAt                                                                                        >  1:251552/1‑65 (MQ=255)
                                                                                                                                                    gggAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGCATCATGAGCCAGCCATCAc  <  5:31164/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GACGCAATGGCGCCAGAGCCTGCAAAAAGCAGGGCAGGCGCTATCGGGGCCTGGTGCTTGGGTTGGATCGGCAAGGCGCGGGGCGGCAGTTGCGCGACGACGGAAGGCGGCGAAGACACGAGCGGTTTGGTGGGGGTGACGCCTTCCAGGGAAGGCGCCGGGTCTACCGGGAGGCGACTTAGGACTTGGCGGTTTTGGGTGGAGAAATCATGACAATGAAGGCGTCAATCAAGGCTTTCCCGCGCTCTGCGCTAACGGGTTCGGACTCGAAAGGCAGCATCATGAGCCAGCCATCAC  >  NC_002947/3223075‑3223371

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: