Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,489,475 C→T 8.9% intergenic (+32/‑44) PP_4824 → / → PP_4825 two‑component system sensor histidine kinase/response regulator/MarC family membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,489,4750CT 8.9% 152.4 / 3.2 56intergenic (+32/‑44)PP_4824/PP_4825two‑component system sensor histidine kinase/response regulator/MarC family membrane protein
Reads supporting (aligned to +/- strand):  ref base C (33/17);  new base T (2/3);  total (35/21)
Fisher's exact test for biased strand distribution p-value = 3.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GCAGGCATGGACGGACACATGGCCAAACCGGTGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGCCACTCCAGCAAG  >  NC_002947/5489330‑5489623
                                                                                                                                                 |                                                                                                                                                     
gcaggcaTGGACGGACACATGGCCAAACCGGTGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTcgc                                                                                                                                                     >  2:360380/1‑148 (MQ=255)
    gcaTGGACGGACACATGGCCAAACCCGAGGAGTAGTCGCAGTTGCGTGACCTGATTCAATACTGGCCCAATAACCGGGAAGCCAAGCTGGACCCGGCACAAACGTCCTGACCGATATCCAATCCTCTTGAAGGAGCGGCCACGAGacac                                                                                                                                                >  6:53056/1‑149 (MQ=255)
      aTGGACGGACACATGGCCAAACCGGTGGAGTTGTCGCAGTTGTGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACaa                                                                                                                                              <  6:363064/149‑1 (MQ=255)
         gacggacACATGGCCAAACCGGTGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAgg                                                                                                                                            <  5:240871/148‑1 (MQ=255)
            ggacACATGGCCAAACCGGTGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCcgc                                                                                                                                        >  2:309571/1‑149 (MQ=255)
              acacATGGCCAAACCGGTGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGCGCGGCCTTGTGTCGCGAGGCcgctc                                                                                                                                      <  8:54064/149‑3 (MQ=255)
                  aTGGCCAAACCGGTGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGc                                                                                                                                  <  1:360380/149‑1 (MQ=255)
                   tGGCCAAACCGGTGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCa                                                                                                                                 >  2:117711/1‑149 (MQ=255)
                      ccAAACCGGTGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATa                                                                                                                               >  3:319716/1‑148 (MQ=255)
                       cAAACCGGTGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATaaa                                                                                                                             >  3:300529/1‑149 (MQ=255)
                            cGGTGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAAccc                                                                                                                          >  1:378996/1‑147 (MQ=255)
                               tGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTTACCGATATCCAATGCTATTGCAGGAGCGGCCTCGCGACACAAGTACGCGCCGGCATAAAACAAAgt                                                                                                                      >  5:111986/1‑147 (MQ=255)
                                     tGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTg                                                                                                               >  6:200467/1‑149 (MQ=255)
                                       tCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGc                                                                                                              <  5:304406/148‑1 (MQ=255)
                                               gCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATCCTc                                                                                                      >  4:255541/1‑148 (MQ=255)
                                                 gTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCAt                                                                                                   >  2:382582/1‑149 (MQ=255)
                                                       tGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAgttg                                                                                             >  5:381256/1‑149 (MQ=255)
                                                          ttCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCCAGGAACACTGCTCATGCTCCATGAgttgtt                                                                                           >  2:297408/1‑148 (MQ=255)
                                                          ttCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAgttgtt                                                                                           <  6:196958/148‑1 (MQ=255)
                                                           tCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTTGTGTCGCGAGGCcgc                                                                                                                                        >  8:174256/1‑102 (MQ=255)
                                                           tCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTTGTGTCGCGAGGCcgc                                                                                                                                        <  7:174256/102‑1 (MQ=255)
                                                                  tGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTa                                                                                  >  6:387710/1‑149 (MQ=255)
                                                                   gggCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTAc                                                                                 >  4:217387/1‑149 (MQ=255)
                                                                      ccAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCt                                                                               >  8:205570/1‑148 (MQ=255)
                                                                        aaTCACCGGGAAGCCAAGGTGGACCCGGCACTAACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACCAGGCCGCGCCGGCATAAACCCAAGGAACACAGCTGATGCTCCATGAGTAGTTCAGCGTCTACCTGa                                                                             >  7:18221/1‑148 (MQ=255)
                                                                              cGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCt                                                                      <  3:217387/149‑1 (MQ=255)
                                                                              cGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCt                                                                      >  1:60805/1‑149 (MQ=255)
                                                                               gggAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTc                                                                     <  1:382582/149‑1 (MQ=255)
                                                                                    gCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGc                                                                       >  3:301819/1‑142 (MQ=255)
                                                                                    gCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGc                                                                       <  4:301819/142‑1 (MQ=255)
                                                                                     ccAAGGTGGACCCGGCACATAAGTCCTGACCGATATACAATGCTCTTGCAGGAGCGGCCTATCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGctc                                                              <  5:53056/149‑1 (MQ=255)
                                                                                       aaGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTa                                                            >  8:41846/1‑149 (MQ=255)
                                                                                             gACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGccc                                                      >  8:87496/1‑149 (MQ=255)
                                                                                              aCCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCg                                                     <  5:200467/149‑1 (MQ=255)
                                                                                                  ggCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTg                                                                 >  3:212176/1‑133 (MQ=255)
                                                                                                  ggCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTg                                                                 <  4:212176/133‑1 (MQ=255)
                                                                                                  ggCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGttct                                                 >  1:263517/1‑149 (MQ=255)
                                                                                                    cacaTACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGttcttc                                               <  5:387710/149‑1 (MQ=255)
                                                                                                        tACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGc                                           >  6:220679/1‑149 (MQ=255)
                                                                                                         aCGTCCTGACCGATATCCAATGATCTTGCAGGAGCGGACTCGCGACACAAGGCCGAGCCGGGATAACCCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCATCGTCTACCTGTAAATGCTCGTGCTCTACAGAACGTTCTTCGTGCt                                          >  6:259521/1‑149 (MQ=255)
                                                                                                            tCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTc                                       >  3:20875/1‑149 (MQ=255)
                                                                                                              cTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCgt                                     >  7:157742/1‑149 (MQ=255)
                                                                                                                gACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCgtgt                                   >  8:34072/1‑149 (MQ=255)
                                                                                                                gACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCgt                                     <  5:308096/147‑1 (MQ=255)
                                                                                                                gACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCgt                                     >  6:308096/1‑147 (MQ=255)
                                                                                                                         ccAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGc                          >  3:333578/1‑149 (MQ=255)
                                                                                                                          cAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGcc                         >  6:66635/1‑149 (MQ=255)
                                                                                                                                  ttGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGt                 >  8:332071/1‑149 (MQ=255)
                                                                                                                                    gCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTgg               >  5:19658/1‑149 (MQ=255)
                                                                                                                                     cAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGc              <  7:34072/149‑1 (MQ=255)
                                                                                                                                     cAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGc              <  2:60805/149‑1 (MQ=255)
                                                                                                                                      aGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGcc             >  6:314019/1‑149 (MQ=255)
                                                                                                                                         cgCGGCCTTGTGTCGCGAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCg                                              >  4:324719/2‑113 (MQ=255)
                                                                                                                                         cgCGGCCTTGTGTCGCGAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCg                                              <  3:324719/112‑1 (MQ=255)
                                                                                                                                             gCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGCCACTCCAg      <  1:138334/149‑1 (MQ=255)
                                                                                                                                              ccTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTCGCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGCCACTCCAg      <  1:309571/149‑1 (MQ=255)
                                                                                                                                              ccTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGCCACTCCAGc     <  5:220679/149‑1 (MQ=255)
                                                                                                                                                 cgcgACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGCCACTCCAGCAAg  <  4:333578/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GCAGGCATGGACGGACACATGGCCAAACCGGTGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGCCACTCCAGCAAG  >  NC_002947/5489330‑5489623

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: