Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,696,568 C→T 7.0% intergenic (+20/‑192) pyrC' → / → hslV dihydroorotase‑like protein/ATP‑dependent HslVU protease peptidase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,696,5680CT 7.0% 161.3 / 2.1 57intergenic (+20/‑192)pyrC'/hslVdihydroorotase‑like protein/ATP‑dependent HslVU protease peptidase subunit
Reads supporting (aligned to +/- strand):  ref base C (16/37);  new base T (2/2);  total (18/39)
Fisher's exact test for biased strand distribution p-value = 5.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCAAGCGTCCACCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCCCCCCACCCCCATA  >  NC_002947/5696420‑5696712
                                                                                                                                                    |                                                                                                                                                 
ccGCAAGCGTCCACCGTTGCTGGCGAACAATGGTTCTCCCGTCGCGAGAACAGCCCCTTCATCGGCCACTGCCTGCCGCGCGCGGTGCGTTATACCTAGGTCGATGGCCACATCTGCCACTAGGCCTGAGTACCATGCTTAGCTGGCAc                                                                                                                                                   >  7:259680/1‑149 (MQ=255)
  gCAAGCGTCCACCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGc                                                                                                                                                 >  1:139593/1‑149 (MQ=255)
  gCAAGCGTCCACCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGc                                                                                                                                                 >  7:370073/1‑149 (MQ=255)
      gCGTCCACCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGc                                                                                                                                             <  4:384200/149‑1 (MQ=255)
      gCGTCCACCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGc                                                                                                                                             >  4:414777/1‑149 (MQ=255)
       cGTCCACCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCt                                                                                                                                            <  4:87016/149‑1 (MQ=255)
        gTCCACCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCATGCCCGCTc                                                                                                                                           >  8:352167/1‑149 (MQ=255)
        gTCCACCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTc                                                                                                                                           <  7:115576/149‑1 (MQ=255)
         tCCACCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCGTGCCCGCtaa                                                                                                                                          <  3:161841/149‑3 (MQ=255)
         tCCACCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTcc                                                                                                                                          <  7:196625/149‑1 (MQ=255)
          ccACCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTcc                                                                                                                                          >  5:132344/1‑148 (MQ=255)
               gTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAgg                                                                                                                                    <  6:132344/149‑1 (MQ=255)
                   cTGGCGAACAATGGTTCTCCCGTGGCGACAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATc                                                                                                                                 <  5:387432/148‑1 (MQ=255)
                           cAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAAtt                                                                                                                        >  2:122311/1‑149 (MQ=255)
                           cAATGGTTCTCACGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTCGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCATGCCCGCTCCCACAGGGTTACcacatg                                                                                                                         >  3:109498/1‑146 (MQ=255)
                                 ttCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAgg                                                                                                                  >  5:31719/1‑149 (MQ=255)
                                   ctcCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGc                                                                                                                 <  1:122311/148‑1 (MQ=255)
                                   ctcCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCa                                                                                                                >  7:189581/1‑149 (MQ=255)
                                      ccGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCAACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCACACACGATAACCACAAATACAAGGccagt                                                                                                             >  7:349977/1‑145 (MQ=255)
                                           gCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTaa                                                                                                       >  3:27386/1‑149 (MQ=255)
                                            cGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAAc                                                                                                      >  1:37326/1‑149 (MQ=255)
                                                 aCTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCtgt                                                                                                 <  1:121302/149‑1 (MQ=255)
                                                 aCTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCtgt                                                                                                 <  7:137644/149‑1 (MQ=255)
                                                  cTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCtgtg                                                                                                <  6:31719/149‑1 (MQ=255)
                                                  cTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCtgtg                                                                                                <  3:414777/149‑1 (MQ=255)
                                                   tGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCAt                                                                                                              <  7:88203/134‑1 (MQ=255)
                                                   tGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCAt                                                                                                              >  8:88203/1‑134 (MQ=255)
                                                   tGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAAACCTGTgg                                                                                               >  5:414812/1‑149 (MQ=255)
                                                   tGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGCTGCGTTAAACATTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACAAAAAATCCAAGGCATGTGGTAAACCtgtg                                                                                                >  2:335039/1‑148 (MQ=255)
                                                     cccGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTgg                                                                                               <  4:99131/147‑1 (MQ=255)
                                                     cccGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTgg                                                                                               <  7:211366/147‑1 (MQ=255)
                                                     cccGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTgg                                                                                               >  4:129590/1‑147 (MQ=255)
                                                        gTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCATGCCCGCTCCCACAGGGTTACCACATGCCTTGGAATTGTGGTGATCCTGTGGGAGCg                                                                                          <  7:352167/149‑1 (MQ=21)
                                                         ttCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCgg                                                                                         <  8:347611/149‑1 (MQ=255)
                                                         ttCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCgg                                                                                         <  2:139593/149‑1 (MQ=255)
                                                                   aCTGGCTGCCGGGGGCGTTGCGTTATACCTTGGTCGATGGCCACATCTTCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCcgc                                                                               <  8:349977/149‑1 (MQ=255)
                                                                        ctgccGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATg                                                                          >  3:165737/1‑149 (MQ=255)
                                                                           ccGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGc                                                                         <  2:281558/147‑1 (MQ=255)
                                                                           ccGGGCGCGGGCCGTAATCCCCTGGTCGCTGGCCACATGTGCCACGAGGCCGGAGTCCCATGCTTAGGGGGCTCGCCAGCTCCCACAGGATCCCCACATCTCCAAGGCCGGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCt                                                                       <  6:153511/149‑1 (MQ=255)
                                                                             gggCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTc                                                                     <  1:396861/149‑1 (MQ=255)
                                                                             gggCGCGGTGCGTTATACCTTGGGCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCACCCACAGGATCACCACAATTCCCAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTc                                                                     <  6:54342/149‑1 (MQ=255)
                                                                              ggCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTca                                                                    <  8:189581/149‑1 (MQ=255)
                                                                                     tGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTc                                                             <  4:256070/149‑1 (MQ=255)
                                                                                       cGTTATACCTTGGTCTCTGGCCACATCTGCCTCGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCTCCACAATTCCAACGACTCGGGTAACCCTGGCGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTccc                                                           <  8:238621/149‑1 (MQ=255)
                                                                                        gTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTcccc                                                          <  8:370073/149‑1 (MQ=255)
                                                                                           ttgccTTCGCCGATGGCGACTTCTGCCACGTTGCCGGAGTTCCATCGGTAGCGGGCACGCCCGCTCGCCCATGATCTCCACATTTGCAAGGCCTGTGGTTTCCCTGGGGGTGCGGGCATGCCCGCGAATGCCTTCACATAGTcccccc                                                        <  8:259680/145‑1 (MQ=255)
                                                                                               cTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTCGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGAATGTGGTAACCCTGTGTGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCaa                                                   <  4:360238/149‑1 (MQ=255)
                                                                                               cTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGAt                                                                                                                                  <  3:268697/71‑1 (MQ=255)
                                                                                               cTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGAt                                                                                                                                  >  4:268697/1‑71 (MQ=255)
                                                                                                  ggTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCAAGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAg                                                  <  8:398400/147‑1 (MQ=255)
                                                                                                    tCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGAATGTGGTATCCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGcgcg                                              <  1:335039/149‑1 (MQ=255)
                                                                                                           ccACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTc                                       <  4:27386/149‑1 (MQ=255)
                                                                                                                          gggcTGAGTACCATGCTTTGGGGGTTCGCCAGCTCCTTGTGGATCTCCACATTTCCAACGCATGTTGTTTACCTTTGTGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGa                        <  2:37326/146‑1 (MQ=255)
                                                                                                                                aGTTCCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCt                  <  3:81541/149‑1 (MQ=255)
                                                                                                                                  tACCATGCTTAGCGGGCACGCGCGCTCCCTCAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACTTAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTc                <  1:108233/149‑1 (MQ=255)
                                                                                                                                    ccttGCTTAGCGGGCACGCAAGATCACTCAGGATCACCACAATTCCAAGGATTGTGGTTAACCTTTGGTAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTcc               <  6:414812/145‑1 (MQ=255)
                                                                                                                                             gCGGGCACGCCCGCTCCCACAGGATCCCCACAATTCCAAGGCTGGTGGTAACCCTGGGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCcccccaccccc     <  4:375790/149‑1 (MQ=255)
                                                                                                                                                ggCACGCCCGCTCCCAGAGGATCTCCTCAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCCCCCCACCCCCata  <  3:230634/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CCGCAAGCGTCCACCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCCCCCCACCCCCATA  >  NC_002947/5696420‑5696712

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