Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,885,730 C→A 10.9% R50L (CGC→CTC)  PP_2541 ← transcriptional factor‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,885,7300CA10.9% 161.6 / 9.7 55R50L (CGC→CTC) PP_2541transcriptional factor‑like protein
Reads supporting (aligned to +/- strand):  ref base C (24/25);  new base A (3/3);  total (27/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 = 5.36e-01

TACGACGTGTGTCGGTGGGGTAAGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGGTGTGGCGGGTCGGGTGGCTGGGC  >  NC_002947/2885585‑2885876
                                                                                                                                                 |                                                                                                                                                  
tACGACGTGTGTCGGTGGGGTAAGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCt                                                                                                                                                 <  1:160487/149‑1 (MQ=255)
tACGACGTGTGTCGGTGGGGTAAGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCt                                                                                                                                                 >  2:144274/1‑149 (MQ=255)
    acgTGTGTCGGTGGGGTAAGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGAc                                                                                                                                             <  6:112824/149‑1 (MQ=255)
    acgTGTGTCGGTGGGGTAAGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCGCGACCTGGTCGAGCTCgcc                                                                                                                                             <  6:221159/149‑3 (MQ=255)
           tCGGTGGGGTAAGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGt                                                                                                                                      <  7:175982/149‑1 (MQ=255)
              gTGGGGTAAGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTg                                                                                                                                   <  2:89093/149‑1 (MQ=255)
                ggggTAAGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCt                                                                                                                                 <  1:191648/149‑1 (MQ=255)
                   gTAAGTGCGTCGGGTTGCAGCTGCTTCGCTTCGACCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCtggt                                                                                                                              >  6:348332/1‑149 (MQ=255)
                   gTAAGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCtggt                                                                                                                              >  5:33408/1‑149 (MQ=255)
                    tAAGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCtggtg                                                                                                                             <  7:341075/149‑1 (MQ=255)
                      aGTGCGTCTGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTggtggt                                                                                                                           <  6:320623/149‑1 (MQ=255)
                                 ttGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACgg                                                                                                                >  4:384905/1‑149 (MQ=255)
                                           tgcGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAggg                                                                                                       >  3:358400/1‑148 (MQ=255)
                                                    cccAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATTCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCaa                                                                                             <  1:144274/149‑1 (MQ=255)
                                                    cccAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCaa                                                                                             >  5:234332/1‑149 (MQ=255)
                                                     ccAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAAc                                                                                            >  5:373031/1‑149 (MQ=255)
                                                     ccAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAAc                                                                                            <  1:396846/149‑1 (MQ=255)
                                                         gCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTggg                                                                                        <  6:33408/149‑1 (MQ=255)
                                                          cTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGt                                                                                       >  1:403139/1‑149 (MQ=255)
                                                           tCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTg                                                                                      >  2:344909/1‑149 (MQ=255)
                                                           tCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTg                                                                                      >  2:420288/1‑149 (MQ=255)
                                                           tCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTg                                                                                      >  3:107971/1‑149 (MQ=255)
                                                              aGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAg                                                                                   <  2:99561/149‑1 (MQ=255)
                                                                              ccAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTgg                                                                    <  2:403139/148‑1 (MQ=255)
                                                                                 gCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTggcggc                                                                >  3:380194/1‑149 (MQ=255)
                                                                                 gCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGcc                                                               <  5:348332/149‑1 (MQ=255)
                                                                                     tAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCggg                                                            <  3:384905/149‑1 (MQ=255)
                                                                                      aaGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCaa                                                                                             >  1:250679/1‑115 (MQ=255)
                                                                                      aaGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCaa                                                                                             <  2:250679/115‑1 (MQ=255)
                                                                                       aGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGgca                                                          <  1:420288/149‑1 (MQ=255)
                                                                                          cATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGgcagc                                                        >  8:224094/1‑148 (MQ=255)
                                                                                           aTGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCaa                                                      >  5:101951/1‑149 (MQ=255)
                                                                                               cTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACtt                                                   >  4:335561/1‑148 (MQ=255)
                                                                                               cTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTg                                                  <  2:105196/149‑1 (MQ=255)
                                                                                                tttACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACtt                                                   >  3:128782/1‑147 (MQ=255)
                                                                                                tttACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACtt                                                   >  6:254287/1‑147 (MQ=255)
                                                                                                    cTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAg                                             <  1:362933/149‑1 (MQ=255)
                                                                                                       ttGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGtt                                           <  3:335561/148‑1 (MQ=255)
                                                                                                       ttGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTg                                          >  2:16085/1‑149 (MQ=255)
                                                                                                           aCTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTg                                                  <  4:337098/137‑1 (MQ=255)
                                                                                                           aCTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTg                                                  >  3:337098/1‑137 (MQ=255)
                                                                                                                 gCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGc                                               <  3:342255/134‑1 (MQ=255)
                                                                                                                 gCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGc                                               >  4:342255/1‑134 (MQ=255)
                                                                                                                       aCAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCGGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATg                          >  3:70576/1‑149 (MQ=255)
                                                                                                                       aCAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATg                          <  4:358400/149‑1 (MQ=255)
                                                                                                                         aGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCACCTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCCCTGGCGGGCGATGgt                        >  7:330938/1‑149 (MQ=255)
                                                                                                                         aGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGgt                        <  5:254287/149‑1 (MQ=255)
                                                                                                                             cgtcgAGCGCGACCTGGTCGAGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGGTGTgg                    >  4:405935/4‑149 (MQ=255)
                                                                                                                               tcgAGCGCGACCTGGTCGAGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGGTGTGGcg                  >  4:122474/2‑149 (MQ=255)
                                                                                                                                  aGCGCGACCTGGTCGAGCTCGCCGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCggg                                                            >  6:285721/1‑104 (MQ=255)
                                                                                                                                  aGCGCGACCTGGTCGAGCTCGCCGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCggg                                                            <  5:285721/104‑1 (MQ=255)
                                                                                                                                       gACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGGTGTGGcgggtcgggt          <  4:128782/149‑1 (MQ=255)
                                                                                                                                       gACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGGTGTGGcgg                 >  6:139863/1‑142 (MQ=255)
                                                                                                                                       gACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGGTGTGGcgg                 <  5:139863/142‑1 (MQ=255)
                                                                                                                                             gTCGAGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGGTGTGGCGGGTCGGGTGGCt      <  3:122474/147‑1 (MQ=255)
                                                                                                                                                gcgcTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGGTGTGGCGGGTCGGGTGGCTGGGc  <  7:224094/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TACGACGTGTGTCGGTGGGGTAAGTGCGTCGGGTTGCAGCTGCTGCGCTTCGCCCAGGCTCAAGCCCGGTTGCCAGTACCAGCCATAAGGCATGCCTTTACTGTTGCACTGAAGCGGGAACAGTGCGGCGAGCTCGACCAGGTCGCGCTCGACGGTGCGTTTGCTGGTGGTATGGCCCACGGTTGTCAGGGCAGCCTGCAACTGGGTGGAGCTCATTCCCGGGTGGCGGCCGGGCAGCAACTTGAGCAGTTGCCACTGGCGGGCGATGGTGTGGCGGGTCGGGTGGCTGGGC  >  NC_002947/2885585‑2885876

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: