Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,593,713 T→A 11.1% intergenic (+19/+280) mccA → / ← PP_4068 methylcrotonyl‑CoA carboxylase biotin‑containing subunit alpha/Cro/CI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,593,7130TA11.1% 143.4 / 8.7 54intergenic (+19/+280)mccA/PP_4068methylcrotonyl‑CoA carboxylase biotin‑containing subunit alpha/Cro/CI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (23/25);  new base A (3/3);  total (26/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 = 7.44e-01

TTGGTACTGGAGGCCATGAAGATGGAGCACAGCATTCGTGCGCCGCATGCCGGCACGGTGAAGGCGCTGTTTTGCCTGGAGGGGGATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTA  >  NC_002947/4593566‑4593857
                                                                                                                                                   |                                                                                                                                                  
ttGGTACTGGAGGCCATGAAGATGGAGCACAGCATTCGTGCGCCGCATGCCGGCACGGTGAAGGCGCTGTTTTGCCTGGAGGGGGATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCttt                                                                                                                                                   >  6:266404/1‑149 (MQ=255)
   gTACTGGAGGCCATGAAGATGGAGCACAGCATTCGTGCGCCGCATGCCGGCACGGTGAAGGCGCTGTTTTGCCTGGAGGGGGATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCg                                                                                                                                                <  5:216437/149‑1 (MQ=255)
   gTACTGGAGGCCATGAAGATGGAGCACAGCATTCGTGCGCCGCATGCCGGCACGGTGAAGGCGCTGTTTTGCCTGGAGGGGGATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCg                                                                                                                                                <  3:98665/149‑1 (MQ=255)
   gTACTGGAGGCCATGAAGATGGAGCACAGCATTCGTGCGCCGCATGCCGGCACGGTGAAGGCGCTGTTTTGCCTGGAGGGGGATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCg                                                                                                                                                >  5:69254/1‑149 (MQ=255)
   gTACTGGAGGCCATGAAGATGGAGCACAGCATTCGTGCGCCGCATGCCGGCACGGTGAAGGCGCTGTTTTGCCTGGAGGGGGATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCg                                                                                                                                                <  3:120489/149‑1 (MQ=255)
    tACTGGAGGCCATGAAGATGGAGCACAGCATTCGTGCGCCGCATGCCGGCACGGTGAAGGCGCTGTTTTGCCTGGAGGGGGATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCgg                                                                                                                                               >  5:97361/1‑149 (MQ=255)
       tGGAGGCCATGAAGATGGAGCACAGCATTCGTGCGCCGCATGCCGGCACGGTGAAGGCGCTGTTTTGCCTGGAGGGGGATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCAAAGCGGcc                                                                                                                                             >  4:30847/1‑148 (MQ=255)
         gAGGCCATGAAGATGGAGCACAGCATTCGTGCGCCGCATGCCGGCACGGTGAAGGCGCTGTTTTGCCTGGAGGGGGATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCGGccc                                                                                                                                            <  3:329095/147‑1 (MQ=255)
          aGGCCATGAAGATGGAGCACAGCATTCGTGCGCCGCATGCCGGCACGGTGAAGGCGCTGTTTTGCCTGGAGGGGGATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCAAAGCGGCCCtcc                                                                                                                                         <  8:174848/149‑3 (MQ=255)
               aTGAAGATGGAGCACAGCATTCGTGCGCCGCATGCCGGCACGGTGAAGGCGCTGTTTTGCCTGGAGGGGGATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCaa                                                                                                                                    <  6:340011/149‑1 (MQ=255)
               aTGAAGATGGAGCACAGCATTCGTGCGCCGCATGCCGGCACGGTGAAGGCGCTGTTTTGCCTGGAGGGGGATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCAAAGCGGcc                                                                                                                                             <  6:270547/140‑1 (MQ=255)
               aTGAAGATGGAGCACAGCATTCGTGCGCCGCATGCCGGCACGGTGAAGGCGCTGTTTTGCCTGGAGGGGGATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCAAAGCGGcc                                                                                                                                             >  5:270547/1‑140 (MQ=255)
               aTGAAGATGGAGCACAGCATTCGTGCGCCGCATGCCGGCACGGTGAAGGCGCTGTTTTGCCTGGAGGGGGATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCAAAGCGGCCCtccggcac                                                                                                                                    <  3:30847/149‑8 (MQ=255)
                 gAAGATGGAGCACAGCATTCGTGCGCCGCATGCCGGCACGGTGAAGGCGCTGTTTTGCCTGGAGGGGGATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAaca                                                                                                                                  >  5:376685/1‑149 (MQ=255)
                   aGATGGAGCACAGCATTCGTGCGCCGCATGCCGGCACGGTGAAGGCGCTGTTTTGCCTGGAGGGGGATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAAcaca                                                                                                                                <  5:31295/149‑1 (MQ=255)
                      tGGAGCACAGCATTCGTGCGCCGCATGCCGGCACGGTGAAGGCGCTGTTTTGCCTGGAGGGGGATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAACACAAg                                                                                                                              <  5:63235/148‑1 (MQ=255)
                       ggAGCACAGCATTCGTGCGCCGCATGCCGGCACGGTGAAGGCGCTGTTTTGCCTGGAGGGGGATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAACACAAGgc                                                                                                                            >  7:83469/1‑149 (MQ=255)
                          gCACAGCATTCGTGCGCCGCATGCCGGCACGGTGAAGGCGCTGTTTTGCCTGGAGGGGGATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAACACAAGgccgc                                                                                                                         >  5:373060/1‑149 (MQ=255)
                           cacaGCATTCGTGCGCCGCATGCCGGCACGGTGAAGGCGCTGTTTTGCCTGGAGGGGGATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTTCGGCCCTAACGCAACACAAGgccgcc                                                                                                                        >  7:147201/1‑149 (MQ=255)
                                cATTCGTGCGCCGCATGCCGGCACGGTGAAGGCGCTGTTTTGCCTGGAGGGGGATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCAGAGGGCCGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTg                                                                                                                    >  5:47169/1‑148 (MQ=255)
                                          cctcATGACTGCTCGGTGAAGGCGCTGTTTTGCCTGGAGGGGGTTATGGTCAGCGAAGGGAGTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGt                                                                                                        <  2:236720/146‑1 (MQ=255)
                                          ccgcATGCCGGCACGGTGAAGGCGCTGTTTTGCCTGGAGGGGGATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGt                                                                                                        <  8:83469/149‑1 (MQ=255)
                                                    gctcGGTGAAGGCGCTGTTTTGCCTGGAGGGGGATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGt                                                                                              <  5:266404/146‑1 (MQ=255)
                                                    gCACGGTGAAGGCGCTGTTTTGCCTGGAGGGGGATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGt                                                                                              >  6:42529/1‑149 (MQ=255)
                                                     cACGGTGAAGGCGCTGTTTTGCCTGGAGGGGGATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGtt                                                                                             >  4:19183/1‑149 (MQ=255)
                                                      aCGGTGAAGGCGCTGTTTTGCCTGGAGGGGGATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGttt                                                                                            <  3:52058/149‑1 (MQ=255)
                                                      aCGGTGAAGGCGCTGTTTTGCCTGGAGGGGGAAATGGACCGCGAAGGGAATTTGCTGGGGGAAGTGGTGTATTAACGTTCCCGGAGGCCCGATTTCCGCCCAAAACCAAACACAAGGCCGCCCCGCAACACGCCCTAAGGACCGGGtt                                                                                             >  4:312549/1‑148 (MQ=255)
                                                                 gcTGTTTTGCCTGGAGGGGGATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAg                                                                                 >  4:339301/1‑149 (MQ=255)
                                                                  cTGTTTTGCCTGGAGGGGGATATCGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAgg                                                                                <  1:338195/149‑1 (MQ=255)
                                                                     ttttGCCTGGAGGGGGATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAgggc                                                                              <  3:339301/148‑1 (MQ=255)
                                                                       ttGCCTGGAGGGGGATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAgggccgg                                                                           >  7:183640/1‑149 (MQ=255)
                                                                           cTGGAGGGGGATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCct                                                                       >  2:58097/1‑149 (MQ=255)
                                                                             gcaGGGGGATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCAGGAGGGCCGATTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTAACCGCGAAAGGGCCGGGCctgc                                                                     >  7:171504/3‑149 (MQ=255)
                                                                                 ggggATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGc                                                                 <  6:69254/149‑1 (MQ=255)
                                                                                     atatGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAAt                                                             <  5:42529/149‑1 (MQ=255)
                                                                                     atatGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAAt                                                             <  3:19183/149‑1 (MQ=255)
                                                                                           tCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAAtctc                                                          <  4:111457/146‑1 (MQ=255)
                                                                                           tCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAAtctc                                                          >  3:111457/1‑146 (MQ=255)
                                                                                             aGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTAccc                                                    >  7:154591/1‑149 (MQ=255)
                                                                                               cGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACACGCGAAAGGGCCGGGCCTGCTGGCCAAACTCTACCCGt                                                  >  3:287867/1‑149 (MQ=255)
                                                                                                     gACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCa                                            >  3:25388/1‑149 (MQ=255)
                                                                                                     gACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCa                                            >  7:202810/1‑149 (MQ=255)
                                                                                                          tgCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGt                                        <  4:256761/148‑1 (MQ=255)
                                                                                                          tgCTGGTCGAACTGGTGGAATAACGGTGCCAGAGGGCCGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGtt                                       <  6:47169/149‑1 (MQ=255)
                                                                                                           gCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTg                                      <  8:147201/149‑1 (MQ=255)
                                                                                                                   aaCTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCCAAAGGGCCGGGCCTGCTGGCCAAGCGCTACCCTGCGCTCACCGTTgccgcgcc                               >  6:122525/1‑145 (MQ=37)
                                                                                                                              tAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTcgc                   <  6:97361/149‑1 (MQ=255)
                                                                                                                              tAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTcgc                   <  8:183640/149‑1 (MQ=255)
                                                                                                                                aCGGTGCCGGAGGGCCGCTTTGCGGCCCTAAGGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCgg                 <  8:171504/149‑1 (MQ=255)
                                                                                                                                 cGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCgg                 <  1:58097/149‑1 (MQ=255)
                                                                                                                                       gttaGGGCCGCAAAGCGGCCCTCACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAAcc          >  2:120450/4‑149 (MQ=21)
                                                                                                                                        ggAGGGCCGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAAccc         >  2:78158/1‑149 (MQ=255)
                                                                                                                                            ggCCGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCt      <  8:202810/148‑1 (MQ=255)
                                                                                                                                              ccGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTcc    >  2:41686/1‑148 (MQ=255)
                                                                                                                                               cGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTa  >  8:34930/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTGGTACTGGAGGCCATGAAGATGGAGCACAGCATTCGTGCGCCGCATGCCGGCACGGTGAAGGCGCTGTTTTGCCTGGAGGGGGATATGGTCAGCGAAGGGACTGTGCTGGTCGAACTGGTGGAATAACGGTGCCGGAGGGCCGCTTTGCGGCCCTAACGCAACACAAGGCCGCCCCTGCAACAGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCTACCCGTCGCTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTA  >  NC_002947/4593566‑4593857

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: