Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,550,547 A→T 6.6% H316L (CAC→CTC)  groL → chaperonin GroEL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,550,5470AT 6.6% 176.4 / 3.0 61H316L (CAC→CTC) groLchaperonin GroEL
Reads supporting (aligned to +/- strand):  ref base A (31/26);  new base T (2/2);  total (33/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.36e-01

AAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGATCGAAGCACGCGTCAAGCAGATCCGTGCCCAGATCGAAGAAACTTCCTCGGACTACGACCGTGA  >  NC_002947/1550414‑1550688
                                                                                                                                     |                                                                                                                                             
aaGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCAGGTGCTCCAGGGTAGCg                                                                                                                                    <  3:78709/145‑1 (MQ=255)
aaGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCAGGTGCTCCAGGGTAGCg                                                                                                                                    >  4:78709/1‑145 (MQ=255)
   gTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAgcg                                                                                                                             >  1:52025/1‑149 (MQ=255)
   gTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGAAGCACCTGGGTAACGCCAAgcg                                                                                                                             >  3:156148/1‑149 (MQ=255)
       ctgcAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGt                                                                                                                          >  1:185601/1‑148 (MQ=255)
         gcAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCAt                                                                                                                       >  3:250764/1‑149 (MQ=255)
             tCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTg                                                                                                                   >  7:159244/1‑149 (MQ=255)
               aaGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTc                                                                                                                 <  4:250764/149‑1 (MQ=255)
                       gggCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGAAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGaaa                                                                                                         <  4:156148/149‑1 (MQ=255)
                        ggCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGaaaa                                                                                                        <  4:232655/149‑1 (MQ=255)
                         gCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGaaaa                                                                                                        >  7:303943/1‑148 (MQ=255)
                               gCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAAcacca                                                                                                   <  5:84219/147‑1 (MQ=255)
                               gCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAAcacca                                                                                                   >  6:84219/1‑147 (MQ=255)
                                 gATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGaa                                                                                                          <  3:237740/138‑1 (MQ=255)
                                 gATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGaa                                                                                                          >  4:237740/1‑138 (MQ=255)
                                     gccgcAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACcatcat                                                                                            <  2:52025/148‑1 (MQ=255)
                                       cgcAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCAGGTGCTCCAggg                                                                                                                                         >  4:28780/1‑101 (MQ=255)
                                       cgcAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCAGGTGCTCCAggg                                                                                                                                         <  3:28780/101‑1 (MQ=255)
                                               cATGCTGCAGGACATCTCTGTCGTGACTGGCGGCCAGGTCATCTCCGAAGATATCGGCCTGTCCGTTGAAAGCGCTAGCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGAct                                                                                       <  4:187965/143‑2 (MQ=255)
                                               cATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAac                                                                                                       <  2:362865/127‑1 (MQ=255)
                                               cATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAac                                                                                                       >  1:362865/1‑127 (MQ=255)
                                               cATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAAGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACAACACCATCATCGACGGTGCTg                                                                                 >  6:189823/1‑149 (MQ=255)
                                               aaTGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGCTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGAcc                                                                                       >  3:187965/2‑142 (MQ=255)
                                                  cagtcAGGACAGCGATGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAac                                                                                                       >  3:37948/5‑124 (MQ=255)
                                                    tgcAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCt                                                                            <  8:303943/149‑1 (MQ=255)
                                                      cAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTGATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGa                                                                          <  5:234874/149‑1 (MQ=255)
                                                       aGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGa                                                                                                           >  2:206253/1‑115 (MQ=255)
                                                               gCTGTCCTGACTGGCGGCCAAGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACAGCACCATCATCGACGGTGCTGGCGCTGACACCGAGAt                                                                 <  5:189823/149‑1 (MQ=255)
                                                                   tCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGATCGaa                                                             <  8:333774/149‑1 (MQ=255)
                                                                      tGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCt                                                                            >  1:317568/1‑131 (MQ=255)
                                                                      tGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCt                                                                            <  2:317568/131‑1 (MQ=255)
                                                                              ggcCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAgg                                                                                                            >  7:202169/1‑91 (MQ=255)
                                                                              ggcCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAgg                                                                                                            <  8:202169/91‑1 (MQ=255)
                                                                               gcCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGATCGAAGCAcggatctgc                                                 >  1:293303/1‑142 (MQ=255)
                                                                                    gTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGATCGa                                                              >  1:236591/1‑131 (MQ=255)
                                                                                    gTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGATCGa                                                              <  2:236591/131‑1 (MQ=255)
                                                                                       aTCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGATCGAAGCACGCGTCAAGCAGATCCg                                         >  1:254710/1‑149 (MQ=255)
                                                                                          tcCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGATCGAAGCACGCGTCAAGCAGATCCGTGc                                      >  4:295609/1‑149 (MQ=255)
                                                                                             gaagaaATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCgaga                                                                  >  7:270341/1‑118 (MQ=255)
                                                                                             gaagaaATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCgaga                                                                  <  8:270341/118‑1 (MQ=255)
                                                                                             gaagaaATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGATCGAAGCACGCGTc                                                    >  5:199873/1‑132 (MQ=255)
                                                                                             gaagaaATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGATCGAAGCACGCGTc                                                    <  6:199873/132‑1 (MQ=255)
                                                                                             gaagaaATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGATCGAAGCACGCGTCAAGCAGATCCGTGcc                                     >  5:328777/1‑147 (MQ=255)
                                                                                             gaagaaATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGATCGAAGCACGCGTCAAGCAGATCCGTGcc                                     >  1:37540/1‑147 (MQ=255)
                                                                                                     cGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGATCGAAGCACGCGTCAAGCAGATCCGTGCCCAGATCg                              >  4:303877/1‑146 (MQ=255)
                                                                                                      ggCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGATCGAAGCACGCGTCAAGCAGATCCGTGCCCAGATCgaag                           <  6:11157/148‑1 (MQ=255)
                                                                                                      ggCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGATCGAAGCACGCGTCAAGCAGATCCGTGCCCAGATCgaag                           <  6:6006/148‑1 (MQ=255)
                                                                                                              ccTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGATCGAAg                                                            >  5:115344/1‑107 (MQ=255)
                                                                                                              ccTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGATCGAAg                                                            <  6:115344/107‑1 (MQ=255)
                                                                                                              ccTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGATCGAAGCACGCGTCAAGCAGATCCGTGCCCAGATCGAAGAAACTTCCt                  >  2:308067/1‑149 (MQ=255)
                                                                                                               cTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGATCGAAGCACGCGTCAAGCAGATCCGTGCCCAGATCGAAGAAACTTCCTc                 >  6:240346/1‑149 (MQ=255)
                                                                                                                tGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGATCGAAGCACGCGTCAAGCAGATCCGTGCCCAGATCGAAGAAACTTCCTCg                >  5:139078/1‑149 (MQ=255)
                                                                                                                   aaaCCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGATCGAAGCACGCGTCAAGc                                                <  8:162880/114‑1 (MQ=255)
                                                                                                                   aaaCCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGATCGAAGCACGCGTCAAGc                                                >  7:162880/1‑114 (MQ=255)
                                                                                                                        gCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGATCGAAGCACGCGTCAAGCAGATCCGTGCCCAGATCGAAGAAACTTCCTCGGACTACGa        <  6:139078/149‑1 (MQ=255)
                                                                                                                            cccTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGATCGAAGCACGCGTCAAGCAGATCCGTGCCCAGATCGAAGAAACTTCCTCGGACTACGAcc      <  5:240346/147‑1 (MQ=255)
                                                                                                                            cccTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGATCGAAGCACGCGTCAAGCAGATCCGTGCCCAGATCGAAGAAACTTCCTCGGACTACGAcc      <  2:37540/147‑1 (MQ=255)
                                                                                                                             ccTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGATCGAAGCACGCGTCAAGCAGATCCGTGCCCAGATCGAAGAAACTTCCTCGGACTACGACCGt    <  4:117276/148‑1 (MQ=255)
                                                                                                                             ccTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGATCGAAGCACGCGTCAAGCAGATCCGTGCCCAGATCGAAGAAACTTCCTCGGACTACGACCGt    >  3:117276/1‑148 (MQ=255)
                                                                                                                              cTGGAGCACCTGGGAAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGATCGAAGCACGCGTCAAGCAGATCCGTGCCCAGATCGAAGAAACTTCCTCGGACTACGACCGTga  <  2:254710/149‑1 (MQ=255)
                                                                                                                               tGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGATCGAAGCACGCGTCAAGCAGATCCGTGCCCAGATCGAAGAAACTTCCTCGGACTACGACCGt    >  1:345997/1‑146 (MQ=255)
                                                                                                                                     |                                                                                                                                             
AAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGATCGAAGCACGCGTCAAGCAGATCCGTGCCCAGATCGAAGAAACTTCCTCGGACTACGACCGTGA  >  NC_002947/1550414‑1550688

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: