Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,122,886 G→T 6.9% Q107H (CAG→CAT)  PP_2740 → MerR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,122,8860GT 6.9% 185.0 / 3.4 58Q107H (CAG→CAT) PP_2740MerR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (25/29);  new base T (2/2);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

ATCTTGCGCTGGACCGCCAACGGCCTGCCCATCAGCAAGGTTGGCGACAGGCTCGCCGGCCAGCCTGAGGCGCCTGTGCAGCAGCCCTCGGCCTTTGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGCTGATGCCGGTTTGGCATGACCTGGCTTCAGGTACTGCGCCAGGCCAGCGTAGCCAATGGCTGTTTCTCGATACGTTCT  >  NC_002947/3122740‑3123031
                                                                                                                                                  |                                                                                                                                                 
aTCTTGCGCTGGACCGCCAACGGCCTGCCCATCAGCAAGGTTGGCGACAGGCTCGCCGGCCAGCCTGAGGCGCCTGTGCAGCAGCCCTCGGCCTTTGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGgc                                                                                                                                                 >  4:135349/1‑149 (MQ=255)
     gcgcTGGACCGCCAACGGCCTGCCCATCAGCAAGGTTGGCGACAGGCTCGCCGGCCAGCCTGAGGCGCCTGTGCAGCAGCCCTCGGCCTTTGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTg                                                                                                                                             <  7:128656/148‑1 (MQ=255)
      cgcTGGACCGCCAACGGCCTGCCCATCAGCAAGGTTGGCGACAGGCTCGCCGGCCAGCCTGAGGCGCCTGTGCAGCAGCCCTCGGCCTTTGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTgg                                                                                                                                            <  3:327367/148‑1 (MQ=255)
       gcTGGACCGCCAACGGCCTGCCCATCAGCAAGGTTGGCGACAGGCTCGCCGGCCAGCCTGAGGCGCCTGTGCAGCAGCCCTCGGCCTTTGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTgg                                                                                                                                            >  6:102906/1‑147 (MQ=255)
           gACCGCCAACGGCCTGCCCATCAGCAAGGTTGGCGACAGGCTCGCCGGCCAGCCTGAGGCGCCTGTGCAGCAGCCCTCGGCCTTTGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCATGCGCTGGAAGCGa                                                                                                                                      >  8:239761/1‑149 (MQ=255)
            aCCGCCAACGGCCTGCCCATCAGCAAGGTTGGCGACAGGCTCGCCGGCCAGCCTGAGGCGCCTGTGCAGCAGCCCTCGGCCTTTGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGAt                                                                                                                                     >  6:82871/1‑149 (MQ=255)
                ccAACGGCCTGCCCATCAGCAAGGTTGGCGACAGGCTCGCCGGCCAGCCTGAGGCGCCTGTGCAGCGGCCCTCGGCCTTTGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCAt                                                                                                                                 >  6:173271/1‑149 (MQ=255)
                ccAACGGCCTGCCCATCAGCAAGGTTGGCGACAGGCTCGCCGGCCAGCCTGAGGCGCCTGTGCAGCAGCCCTCGGCCTTTGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCAt                                                                                                                                 >  1:2407/1‑149 (MQ=255)
                 cAACGGCCTGCCCATCAGCAAGGTTGGCGACAGGCTCGCCGGCCAGCCTGAGGCGCCTGTGCAGCAGCCCTCGGCCTTTGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATg                                                                                                                                <  5:82871/149‑1 (MQ=255)
                 cAACGGCCTGCCCATCAGCAAGGTTGGCGACAGGCTCGCCGGCCAGCCTGAGGCGCCTGTGCAGCAGCCCTCGGCCTTTGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATg                                                                                                                                <  7:79698/149‑1 (MQ=255)
                   aCGGCCTGCCCATCAGCAAGGTTGGCGACAGGCTCGCCGGCCAGCCTGAGGCGCCTGTGCAGCAGCCCTCGGCCTTTGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGa                                                                                                                              <  5:149728/149‑1 (MQ=255)
                        cTGCCCATCAGCAAGGTTGGCGACAGGCTCGCCGGCCAGCCTGAGGCGCCTGTGCAGCAGCCCTCGGCCTTTGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCATGCGCTGGAAGCGATCCATGGACAGCt                                                                                                                         <  7:239761/149‑1 (MQ=255)
                        cTGCCCATCAGCAAGGTTGGCGACAGGCTCGCCGGCCAGCCTGAGGCGCCTGTGCAGCAGCCCTCGGCCTTTGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCt                                                                                                                         >  7:213345/1‑149 (MQ=255)
                           cccATCAGCAAGGTTGGCGACAGGCTCGCCGGCCAGCCTGAGGCGCCTGTGCAGCAGCCCTCGGCCTTTGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGtt                                                                                                                      <  8:18915/149‑1 (MQ=255)
                                 aGCAAGGTTGGCGACAGGCTCGCCGGCCAGCCTGAGGCGCCTGTGCAGCAGCCCTCGGCCTTTGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCacac                                                                                                                 <  8:213345/148‑1 (MQ=255)
                                    aaGGTTGGCGACAGGCTCGCCGGCCAGCCTGAGGCGCCTGTGCAGCAGCCCTCGGCCTTTGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACtttt                                                                                                             <  1:194611/149‑1 (MQ=255)
                                      ggTTGGCGACAGGCTCGCCGGCCAGCCTGAGGCGCCTGTGCAGCAGCCCTCGGCCTTTGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGc                                                                                                           >  4:73387/1‑149 (MQ=255)
                                      ggTTGGCGACAGGCTCGCCGGCCAGCCTGAGGCGCCTGTGCAGCAGCCCTCGGCCTTTGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGc                                                                                                           >  8:57906/1‑149 (MQ=255)
                                           gCGACAGGCTCGCCGGCCAGCCTGAGGCGCCTGTGCAGCAGCCCTCGGCCTTTGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCaaa                                                                                                       <  3:142910/148‑1 (MQ=255)
                                            cGACAGGCTCGCCGGCCAGCCTGAGGCGCCTGTGCAGCGGCCCTCGGCCTTTGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAgg                                                                                                     <  5:173271/149‑1 (MQ=255)
                                                 ggCTCGCCGGCCAGCCTGAGGCGCCTGTGCAGCAGCCCTCGGCCTTTGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACg                                                                                                <  7:57906/149‑1 (MQ=255)
                                                            cAGCCTGAGGCGCCTGTGCAGCAGCCCTCGGCCTTTGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCAATAGGCCACGGTGCTGCgg                                                                                       >  8:234416/1‑147 (MQ=255)
                                                                 tGAGGCGCCTGTGCAGCAGCCCTCGGCCTTTGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGc                                                                                >  8:30614/1‑149 (MQ=255)
                                                                     gcgcCTGTGCAGCAGCCCTCGGCCTTTGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGCTGAt                                                                            >  6:42572/1‑149 (MQ=255)
                                                                       gcCTGTGCAGCAGCCCTCGGCCTTTGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGCTGATGc                                                                          <  2:2407/149‑1 (MQ=255)
                                                                       gcCTGTGCAGCAGCCCTCGGCCTTTGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGCTGATGc                                                                          <  5:102906/149‑1 (MQ=255)
                                                                         cTGTGCAGCAGCCCTCGGCCTTTGACGCCTGGCGCGCCGCCGGTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGCTGATGCCg                                                                        <  8:61942/149‑1 (MQ=255)
                                                                                  agcCCTCGGCCTTTGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGCTGATGCCGGTTTGGCAt                                                               <  3:73387/149‑1 (MQ=255)
                                                                                      cTCGGCCTTTGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGCTGATGCCGGTTTGGCATGAcc                                                           <  6:196338/149‑1 (MQ=255)
                                                                                       tCGGCCTTTGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGCTCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGCTGATGCCGGTTTGGCATGACCt                                                          >  1:245628/1‑149 (MQ=255)
                                                                                       tCGGACTTTGACGACTGGAGCGCCGCCGTTGCTCGCGCTTCCAGGGCATTCGATGCGAAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCAACGGTGCTGCGGGACGTGCTGATGCCTGTTTGGCATGACCt                                                          >  2:90907/1‑149 (MQ=255)
                                                                                        cGGCCTTTGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGCTGATGCCGGTTTGGCATGACCTg                                                         <  7:64980/149‑1 (MQ=255)
                                                                                             tttGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGCTGATGCCGGTTTGGCATGACCTGGCtt                                                     >  1:349294/1‑148 (MQ=255)
                                                                                                 acgacTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGCTGATGCCGGTTTGGCATGACCTGGCtt                                                     <  4:322517/144‑1 (MQ=255)
                                                                                                 acgacTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGCTGATGCCGGTTTGGCATGACCTGGCtt                                                     >  3:322517/1‑144 (MQ=255)
                                                                                                          gcgcCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGCTGATGCCGGTTTGGCATg                                                              <  2:29459/126‑1 (MQ=255)
                                                                                                          gcgcCGCCGTTGCGCGCGCTTCCAGGGCATACGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGCTGATGCCGGTTTGGCATg                                                              >  1:29459/1‑126 (MQ=255)
                                                                                                           cgccgcCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGCTGATGCCGGTTTGGCATGACCTGGCTTCAGGTACTgcgc                                         >  8:124345/1‑146 (MQ=255)
                                                                                                           cgccgcCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGCTGATGCCGGTTTGGCATGACCTGGCTTCAGGTACTgcgc                                         <  7:124345/146‑1 (MQ=255)
                                                                                                             ccgccgTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGATGCTGCGGGACGTGCTGATGCCGGTTTGGCATGACCTGGCTTCAGGTACTGCgccaggc                                    <  7:304372/149‑1 (MQ=255)
                                                                                                                       gcgcgcTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGc                                                                                >  7:95799/1‑95 (MQ=255)
                                                                                                                       gcgcgcTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGc                                                                                <  8:95799/95‑1 (MQ=255)
                                                                                                                        cgcgcTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGCTGATGCCGGTTTGGCATGACCTGGCTTCAGGTACTGCGCCAGGCCAGCGTAGcc                          >  2:361480/1‑148 (MQ=255)
                                                                                                                        atcgcTTCCAGCGCCTGCGCATCGAATGCCCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGCTGATGCCGGTTTGGCATGACCTGGCTTCAGGTACTGCgccaggccag                                 >  6:172830/3‑141 (MQ=37)
                                                                                                                        atcgcTTCCAGCGCCTGCGCATCGAATGCCCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGCTGATGCCGGTTTGGCATGACCTGGCTTCAGGTACTGCgccaggccag                                 <  5:172830/139‑1 (MQ=37)
                                                                                                                            cTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCTCGGTGCTGCGGGACGTGCTGATGCCGGTTTGGCATGACCTGGCTTCAGGTACTGCGCCAGGCCAGCGTAGCCAATgg                     <  1:361480/149‑1 (MQ=255)
                                                                                                                            cTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGCTGATGCCGGTTTGGCATGACCTGGCTTCAGGTACTGCGCCAGGCCAGCGTAGCCAATgg                     <  2:349294/149‑1 (MQ=255)
                                                                                                                             ttCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGCTGATGCCGGTTTGGCATGACCTGGCTTCAGGTACTGCGCCAGGCCAGc                                >  7:290222/1‑137 (MQ=255)
                                                                                                                             ttCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGCTGATGCCGGTTTGGCATGACCTGGCTTCAGGTACTGCGCCAGGCCAGc                                <  8:290222/137‑1 (MQ=255)
                                                                                                                              tCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGCTGATGCCGGTTTGGCATGACCTGGCTTCAGGTACTGCgcca                                       <  4:20577/129‑1 (MQ=255)
                                                                                                                              tCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGCTGATGCCGGTTTGGCATGACCTGGCTTCAGGTACTGCgcca                                       >  3:20577/1‑129 (MQ=255)
                                                                                                                               ccAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGCTGATGCCGGTTTGGCATGACCTGGCTTCAGGTACTGCGCCAGGCCAGCGTAGCCAATGGCTg                  >  3:242204/1‑149 (MQ=255)
                                                                                                                                cAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGCTGATGCCGGTTTGGCATGACCTGGCTTCAGGTACTGCGCCAGGCCAGCGTAGCCAATGGc                    >  7:106696/1‑146 (MQ=255)
                                                                                                                                cAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGCTGATGCCGGTTTGGCATGACCTGGCTTCAGGTACTGCGCCAGGCCAGCGTAGCCAATGGc                    <  8:106696/146‑1 (MQ=255)
                                                                                                                                  gggCATTCGATGCGCAGGCGCTGGAAGCGATCGATGGACAGCTGTTCACACTTTTGCCaaa                                                                                                       <  1:204831/61‑1 (MQ=255)
                                                                                                                                  gggCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCaaa                                                                                                       >  2:204831/1‑61 (MQ=255)
                                                                                                                                       ttCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGCTGATGCCGGTTTGGCATGACCTGGCTTCAGGTACTGCGCCAGGCCAGCGTAGCCAATGGCTGTTTCTCGa          <  4:242204/149‑1 (MQ=255)
                                                                                                                                               gcAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGCTGATGCCGGTTTGGCATGACCTGGCTTCAGGTACTGCGCCAGGCCAGCGTAGCCAATGGCTGTTTCTCGATACGTTCt  >  5:64675/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ATCTTGCGCTGGACCGCCAACGGCCTGCCCATCAGCAAGGTTGGCGACAGGCTCGCCGGCCAGCCTGAGGCGCCTGTGCAGCAGCCCTCGGCCTTTGACGACTGGCGCGCCGCCGTTGCGCGCGCTTCCAGGGCATTCGATGCGCAGGCGCTGGAAGCGATCCATGGACAGCTGTTCACACTTTTGCCAAAGGCCACGGTGCTGCGGGACGTGCTGATGCCGGTTTGGCATGACCTGGCTTCAGGTACTGCGCCAGGCCAGCGTAGCCAATGGCTGTTTCTCGATACGTTCT  >  NC_002947/3122740‑3123031

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: ‑