Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,414,587 G→C 7.5% intergenic (‑119/‑141) PP_1236 ← / → dapA‑I glycine cleavage system transcriptional repressor/4‑hydroxy‑tetrahydrodipicolinate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,414,5870GC 7.5% 165.0 / 4.4 53intergenic (‑119/‑141)PP_1236/dapA‑Iglycine cleavage system transcriptional repressor/4‑hydroxy‑tetrahydrodipicolinate synthase
Reads supporting (aligned to +/- strand):  ref base G (26/23);  new base C (2/2);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.89e-01

ATTGTTCGCGGACGGTGGGGGTGGACATGCAGCATTCCTAAAGCGATGAGCGGTCGGCACAGGACGGGCAAATGCCTCCTGTCGGCAGGGCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAATAAGCCGCCGAGGGGAATGCTCCACCCGCCTGCTTCGCTTGTGCAAGGTCGATGCCCCCAGTACCATTACCGCTCTCTTTTTCCGGCAGGAGCAGTTACATGATTG  >  NC_002947/1414441‑1414734
                                                                                                                                                  |                                                                                                                                                    
aTTGTTCGCGGACGGTGGGGGTGGACATGCAGCATTCCTAAAGCGATGAGCGGTCGGCACAGGACGGGCAAATGCCTCCTGTCGGCAGGGCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCcggcg                                                                                                                                                    >  5:358117/1‑149 (MQ=255)
   gTTCGCGGACGGTGGGGGTGGACATGCAGCATTCCTAAAGCGATGAGCGGTCGGCACAGGACGGGCAAATGCCTCCTGTCGGCAGGGCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGc                                                                                                                                                 >  2:17604/1‑149 (MQ=255)
   gTTCGCGGACGGTGGGGGTGGACATGCAGCATTCCTAAAGCGATGAGCGGTCGGCACAGGACGGGCAAATGCCTCCTGTCGGCAGGGCCCGGGGCTAGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGc                                                                                                                                                 >  8:272912/1‑149 (MQ=255)
    ttCGCGGACGGTGGGGGTGGACATGCAGCATTCCTAAAGCGATGAGCGGTCGGCACAGGACGGGCAAATGCCTCCTGTCGGCAGGGCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGcc                                                                                                                                                >  1:61598/1‑149 (MQ=255)
                   ggTGGACATGCAGCATTCCTAAAGCGATGAGCGGTCGGCACAGGACGGGCAAATGCCTCCTGTCGGCAGGGCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATc                                                                                                                                  >  2:278282/1‑148 (MQ=255)
                      ggACATGCAGCATTCCTAAAGCGATGAGCGGTCGGCACAGGACGGGCAAATGCCTCCTGTCGGCAGGGCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAg                                                                                                                              <  3:73715/149‑1 (MQ=255)
                         cATGCAGCATTCCTAAAGCGATGAGCGGTCGGCACAGGACGGGCAAATGCCTCCTGTCGGCAGGGCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGc                                                                                                                           >  7:188560/1‑149 (MQ=255)
                          aTGCAGCATTCCTAAAGCGATGAGCGGTCGGCACAGGACGGGCAAATGCCTCCTGTCGGCAGGGCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCaaa                                                                                                                               <  8:83308/144‑1 (MQ=255)
                          aTGCAGCATTCCTAAAGCGATGAGCGGTCGGCACAGGACGGGCAAATGCCTCCTGTCGGCAGGGCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCaaa                                                                                                                               >  7:83308/1‑144 (MQ=255)
                           tGCAGCATTCCTGAAGCGATGAGCGGTCGGCACAGGACGGGCAAATGCCTCCTGTCGGCAGGGCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGcc                                                                                                                          >  3:79142/1‑148 (MQ=255)
                            gcagcaTTCCTAAAGCGATGAGCGGTCGGCACAGGACGGGCAAATGCCTCCTGTCGGCAGGGCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGcc                                                                                                                          >  4:329442/1‑147 (MQ=255)
                            gcagcaTTCCTAAAGCGATGAGCGGTCGGCACAGGACGGGCAAATGCCTCCTGTCGGCAGGGCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGcc                                                                                                                          >  3:165207/1‑147 (MQ=255)
                                     cTAAAGCGATGAGCGGTCGGCACAGGACGGGCAAATGCCTCCTGTCGGCAGGGCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCg                                                                                                                >  6:287656/1‑148 (MQ=255)
                                     cTAAAGAGATGAGCGGTCGGCACAGGACGGACAAATGCCTCCTGTCGGCAGGGCCCGGGGCTCGCAAGCTAGGTCTTGGACGCTGGCGGGCAGCAAGGTCGACGGCCCGGCGGGCCCCGCTGCACCAAACAAACGGCCAAGGGTAGcc                                                                                                                >  8:51708/1‑147 (MQ=255)
                                      tAAAGCGATGAGCGGTCGGCACAGGACGGGCAAATGCCTCCTGTCGGCAGGGCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGa                                                                                                               >  3:85141/1‑148 (MQ=255)
                                      gaaaGCGATGAGCGGTCGGCACAGGACGGGCAAATGCCTCCTGTCGCCAGCGCCCGGGGCTCGCATCCGAGGTCTTGGCCGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGa                                                                                                               >  3:85149/2‑148 (MQ=255)
                                          gCGATGAGCGGTCGGCACAGGACGGGCATATGCGTCCTGTCGGCAGGGCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCTGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAATaa                                                                                                         <  4:85149/149‑1 (MQ=255)
                                          gCGATGAGCGGTCGGCACAGGACGGGCAAATGCCTCCTGTCGGCAGGGCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAATaa                                                                                                         <  4:85141/149‑1 (MQ=255)
                                           ggATGAGCGGTCGGCACAGGACGGGCAAATGCCTCCTGTCGGCAGGGCCCGGGGCTTGCATGCGAGGGCTTGGACGCTGGCGGGGAGCATGGTCGACGGCCCGGCGGGCCGCGCTGGACCGATCAAAGGGCCAAGGGTAGCGAAATAAg                                                                                                        <  7:51708/148‑1 (MQ=255)
                                           cGATGAGCGGTCGGCACAGGACGGGCAAATGCCTCCTGTCGGCAGGGCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAATAAg                                                                                                        >  3:347580/1‑149 (MQ=255)
                                            gATGAGCGGTCGGCACAGGACGGGCAAATGCCTCCTGTCGGCAGGGCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAATAAgc                                                                                                       <  4:347580/149‑1 (MQ=255)
                                             aTGAGCGGTCGGCACAGGACGGGCAAATGCCTCCTGTCGGCAGGGCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAATAAgcc                                                                                                      <  2:61598/149‑1 (MQ=255)
                                              tGAGCGGTCGGCACAGGACGGGCAAATGCCTCCTGTCGGCAGGGCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAATAAgccg                                                                                                     <  8:188560/149‑1 (MQ=255)
                                                aGCGGTCGGCACAGGACGGGCAAATGCCTCCTGTCGGCAGGGCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAATGGCCAAGGGTAGCGAAATAAgccgcc                                                                                                   >  6:204385/1‑149 (MQ=255)
                                                            aGGACGGGCAATTGCCTCCTGTCGGCAGGGCCCGGGGCTCGCATGCGAGGTCTTGGAGGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGaa                                                                                                              <  5:59777/127‑1 (MQ=255)
                                                            aGGACGGGCAAATGCCTCCTGTCGGCAGGGCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAATAAGCCGCCGAGGGGAATGCt                                                                                       >  2:89397/1‑149 (MQ=255)
                                                            aGGACGGGCAAATGCCTCCTGTCGGCAGGGCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAATAAGCCGCCGAGGGGAATGCt                                                                                       <  2:11990/149‑1 (MQ=255)
                                                              gACGGGCAAATGCCTCCTGTCGGCAGGGCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCcgc                                                                                                                                              >  8:165323/1‑93 (MQ=255)
                                                              gACGGGCAAATGCCTCCTGTCGGCAGGGCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAACGTCGACGGCCCGGCGGGCcgc                                                                                                                                              <  7:165323/93‑1 (MQ=255)
                                                                 gggCAAATGCCTCCTGTCGGCAGGGCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAATAAGCCGCCGAGGGGAATGCTCCAcc                                                                                  >  2:100321/1‑149 (MQ=255)
                                                                       aTGCCTCCTGTCGGCAGGGCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAATAAGCCGCCGAGGGGAATGCTCCACCCGCCTg                                                                            >  2:351127/1‑149 (MQ=255)
                                                                          cctcctGTCGGCAGGGCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAATAAGCCGCCGAGGGGAATGCTCCACCCGCCTGCtt                                                                         >  7:100337/1‑149 (MQ=255)
                                                                                gTCGGCAGGGCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGCCAAgg                                                                                                                      <  7:291477/99‑1 (MQ=255)
                                                                                gTCGGCAGGGCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAACGGCCAAgg                                                                                                                      >  8:291477/1‑99 (MQ=255)
                                                                                     cAGGGCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAATAAGCCGCCGAGGGGAATGCTCCACCCGCCTGCTTCGCTTGTGCaa                                                              <  1:17604/149‑1 (MQ=255)
                                                                                        ggCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAATAAGCCGCCGAGGGGAATGCTCCACCCGCCTGCTTCGCTTGTGCAAGGt                                                           <  1:351127/149‑1 (MQ=255)
                                                                                        ggCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAATAAGCCGCCGAGGGGAATGCTCCACCCGCCTGCTTCGCTTGTGCAAGGt                                                           <  1:100321/149‑1 (MQ=255)
                                                                                           ccGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAACTAAGCCGCCGAGGGGAATGCTCCACCCGCCTGCTTCGCTTGTGCAAGGTCg                                                         <  1:89397/149‑1 (MQ=255)
                                                                                              gggCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAATAAGCCGCCGAGGGGAATGCTCCACCCGCCTGCTTCGCTTGTGCAAGGTCGATGc                                                     <  4:79142/149‑1 (MQ=255)
                                                                                                     cATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAATAAGCCGCCGAGGGGAATGCTCCACCCGCCTGCTTCGCTTGTGCAAGGTCGATGCCCCCAGt                                              <  5:287656/149‑1 (MQ=255)
                                                                                                         cGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGTGCCAAGGGTAGCGAAATAAGCCGCCGAGGGGAATGCTCCACCCGCCTGCTTCGCTTGTGCAAGGTCGATGCCCCCAGTAcc                                           <  7:272912/148‑1 (MQ=255)
                                                                                                               cTTGGACGCTGGCGGGGAGGAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAATAAGCCGCCGAGGGGAATGCTCCACCCGCCTGCTTCGCTTGTGCAAGGTCGATGCCCCCAGTACCATTACCg                                    <  6:358117/149‑1 (MQ=255)
                                                                                                                ttGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAATAAGCCGCCGAGGGGAATGCTCCACCCGCCTGCTTCGCTTGTGCAAGGTCGATGCCCCCAGTACCATTACCGc                                   >  1:52310/1‑149 (MQ=255)
                                                                                                                ttGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAATAAGCCGCCGAGGGGAATGCTCCACCCGCCTGCTTCGCTTGTGCAAGGTCGATGCCCCCAGTACCATTACCGc                                   >  3:325228/1‑149 (MQ=255)
                                                                                                                  ggACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAATAAGCCGCCGAGGGGAATGCTCCACCCGCCTGCTTCGCTTGTGCAAGGTCGATGCCCCCAGTACCATTACCg                                    <  4:165207/146‑1 (MQ=255)
                                                                                                                          gCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAATAAGCCGCCGAGGGGAATGCTCCACCCGCCTGCTTCGCTTGTGCAAGGTCGATGCCCCCAGTACCATTACCGCTCTCTTTTTc                         <  3:329442/149‑1 (MQ=255)
                                                                                                                              ggAGCAAGGTCGACGGCCCGGCGGGCCGCGCTTCACCGATCAAAGGGCCAAGGGTAGCGAAATAAGCCGCCGAGGGGAATGCTCCACCCGCCTGCTTCGCTTGTGCAAGGTCGATGCCCCCAGTACCATTACCGCTCTCTTTTTCCGGc                     <  5:204385/149‑1 (MQ=255)
                                                                                                                                     ggTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAATAAGCCGCCGAGGGGAATGCTCCACCCGCCTGCTTCGCTTGTGCAAGGTCGATGCCCCCAGTACCATTACCGCTCTCTTTTTCCGGCAGGAg                >  5:184763/1‑147 (MQ=255)
                                                                                                                                     gacgcgcGGCCCCCCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAATAAGCCGCCGAGGGGAATGCTCCACCCGCCTGCTTCGCTTGTGCAAGGTCGATGCCCCCAGTACCATTACCGCTCTCTTTTTCCGGCAGGAGCa              <  7:318632/143‑1 (MQ=255)
                                                                                                                                      gTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAATAAGCCGCCGAGGGGAATGCTCCACCCGCCTGCTTCGCTTGTGCAAGGTCGATGCCCCCAGTACCATTACCGCTCTCTTTTTCCGGCAGGAGCAg             >  2:166693/1‑149 (MQ=255)
                                                                                                                                       agcgcGGCCCGCCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAATAAGCCGCCGAGGGGAATGCTCCACCCGCCTGCTTCGCTTGTGCAAGGTCGATGCCCCCAGTACCATTACCGCTCTCTTTTTc                         <  4:73098/132‑1 (MQ=255)
                                                                                                                                       agcgcGGCCCGCCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAATAAGCCGCCGAGGGGAATGCTCCACCCGCCTGCTTCGCTTGTGCAAGGTCGATGCCCCCAGTACCATTACCGCTCTCTTTTTc                         >  3:73098/5‑136 (MQ=255)
                                                                                                                                       agcgcGGCCCGCCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAATAAGCCGCCGAGGGGAATGCTCCACCCGCCTGCTTCGCTTGTGCAAGGTCGATGCCCCCAGTACCATTACCGCTCTCTTTTTCCGGCAGGAGCAGt            >  1:329464/5‑149 (MQ=255)
                                                                                                                                            ggCCCGGCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAATAAGCCGCCGAGGGGAATGCTCCACCCGCCTGCTTCGCTTGTGCAAGGTCGATGCCCCCAGTACCATTACCGCTCTCTTTTTCCGGCAGGAGCAGTTACAt       >  2:342068/1‑149 (MQ=255)
                                                                                                                                                cggcggGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAATAAGCCGCCGAGGGGAATGCTCCACCCGCCTGCTTCGCTTGTGCAAGGTCGATGCCCCCAGTACCATTACCGCTCTCTTTTTCCGGCAGGAGCAGTTACATGAtt   <  1:166693/149‑1 (MQ=255)
                                                                                                                                                 ggcggGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAATAAGCCGCCGAGGGGAATGCTCCACCCGCCTGCTTCGCTTGTGCAAGGTCGATGCCCCCAGTACCATTACCGCTCTCTTTTTCCGGCAGGAGCAGTTACATGATTg  >  8:78934/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ATTGTTCGCGGACGGTGGGGGTGGACATGCAGCATTCCTAAAGCGATGAGCGGTCGGCACAGGACGGGCAAATGCCTCCTGTCGGCAGGGCCCGGGGCTCGCATGCGAGGTCTTGGACGCTGGCGGGGAGCAAGGTCGACGGCCCGGCGGGCCGCGCTGCACCGATCAAAGGGCCAAGGGTAGCGAAATAAGCCGCCGAGGGGAATGCTCCACCCGCCTGCTTCGCTTGTGCAAGGTCGATGCCCCCAGTACCATTACCGCTCTCTTTTTCCGGCAGGAGCAGTTACATGATTG  >  NC_002947/1414441‑1414734

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: