Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 953,557 G→T 7.4% E122D (GAG→GAT cyoC → cytochrome bo terminal oxidase subunit III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947953,5570GT 7.4% 166.1 / 3.3 54E122D (GAG→GATcyoCcytochrome bo terminal oxidase subunit III
Reads supporting (aligned to +/- strand):  ref base G (22/28);  new base T (2/2);  total (24/30)
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.87e-01

TGTTGAAACGCTGTTCCTGCTGCTGTCCTCGATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGATCTGGATGGCAATCATGATGTACCAGATCAACAAGCACG  >  NC_002947/953410‑953702
                                                                                                                                                    |                                                                                                                                                 
tGTTGAAACGCTGTTCCTGCTGCTGTCCTCGATCACCTTCGGCTTCGCCATGTGTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAg                                                                                                                                                   >  1:101885/1‑149 (MQ=255)
 gTTGAAACGCTGTTCCTGCTGCTGTCCTCGATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGtt                                                                                                                                                 >  6:219438/1‑149 (MQ=255)
      aaCGCTTTTCCTGCTGCTGTCCTCGATCACCTTCGGCTTCGCCATGCTGAAGTTGTTCGATGGCAAGAAAGCAGGTGTACTTGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAACTCATAGATTTCCATc                                                                                                                                            <  7:357089/149‑1 (MQ=255)
      aaCGCTGTTCCTGCTGCTGTCCTCGATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCAt                                                                                                                                             <  3:233926/148‑1 (MQ=255)
        cGCTGTTCCTGCTGCTGTCCTCGATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCAc                                                                                                                                          <  6:182996/149‑1 (MQ=255)
        cGCTGTTCCTGCTGCTGTCCTCGATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCAc                                                                                                                                          <  7:107527/149‑1 (MQ=255)
         gCTGTTCCTGCTGCTGTCCTCGATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAACTCATAGATTTCCATcgcg                                                                                                                                         <  2:42224/149‑4 (MQ=255)
                 tgctgctgTCCTCGATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGc                                                                                                                                 >  4:220246/1‑149 (MQ=255)
                   ctgctgTCCTCGATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTACAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGa                                                                                                                               <  3:377168/149‑1 (MQ=255)
                      ctgTCCTCGATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAg                                                                                                                              <  6:162453/147‑1 (MQ=255)
                         tCCTCGATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCtt                                                                                                                         >  3:110969/1‑149 (MQ=255)
                              gATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAg                                                                                                                              <  6:296404/139‑1 (MQ=255)
                              gATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAg                                                                                                                              >  5:296404/1‑139 (MQ=255)
                              gATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGcc                                                                                                                    >  8:344205/1‑149 (MQ=255)
                               aTCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGccc                                                                                                                   <  4:110969/149‑1 (MQ=255)
                               aTCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGccc                                                                                                                   >  1:71146/1‑149 (MQ=255)
                                  aCCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCAGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCcgca                                                                                                                <  5:230714/149‑1 (MQ=255)
                                         tcctCGCCATGCTGTAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGTAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGt                                                                                                         <  8:262396/146‑1 (MQ=255)
                                         gCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGt                                                                                                         >  4:59935/1‑149 (MQ=255)
                                               ccATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTc                                                                                                   <  3:105045/149‑1 (MQ=255)
                                                        aaGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCgg                                                                                            <  3:207163/148‑1 (MQ=255)
                                                                      cAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTgg                                                                             >  7:160237/1‑149 (MQ=255)
                                                                       aaGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTggt                                                                            <  7:344205/149‑1 (MQ=255)
                                                                             gCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTAc                                                                      <  2:71146/149‑1 (MQ=255)
                                                                                 gtgtACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCAc                                                                  >  5:72651/1‑149 (MQ=255)
                                                                                 gtgtACTGGGCTGGTTGGCTGTGACCTTCCTGTACGGTGCACGCTTCATCGCTATGGAAATCGAGGAGTTCCATCACCTGATCAGCTAGGGCTTAGGACCGCAGCGCAGTGGCTTCATGTCGGGCTTCTTCGCCCTGGTAGTTACCCAc                                                                  >  5:399675/1‑149 (MQ=255)
                                                                                    tACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACgg                                                                <  3:220246/148‑1 (MQ=255)
                                                                                    tACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGt                                                               >  6:70150/1‑149 (MQ=255)
                                                                                      cTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTc                                                                                                   <  1:343098/111‑1 (MQ=255)
                                                                                      cTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTc                                                                                                   >  2:343098/1‑111 (MQ=255)
                                                                                         ggCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTg                                                            <  8:272807/147‑1 (MQ=255)
                                                                                         ggCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCa                                                          >  4:76159/1‑149 (MQ=255)
                                                                                         ggCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCa                                                          <  5:219438/149‑1 (MQ=255)
                                                                                            tGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGt                                                       <  3:59935/149‑1 (MQ=255)
                                                                                            gggTTGGCGGTGACCTTCCGGTTCGGTGCAGGCTTCATCGCTATGGAAATCTCTGAGTGCCATAACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTTCACGt                                                       <  4:385172/148‑1 (MQ=255)
                                                                                                   cTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCTCTGGTAGGTACCCACGGTCTGCACGTGActgct                                                >  7:29033/1‑149 (MQ=255)
                                                                                                     gtgACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTgg                                              >  2:291014/1‑149 (MQ=255)
                                                                                                          cTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGa                                         >  7:61410/1‑149 (MQ=255)
                                                                                                          cTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGa                                         >  1:138575/1‑149 (MQ=255)
                                                                                                             ccTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGAt                                        <  3:76159/147‑1 (MQ=255)
                                                                                                                 ttCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGATCTGGa                                   <  2:111420/148‑1 (MQ=255)
                                                                                                                 ttCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGATCTGGAt                                  <  8:160237/149‑1 (MQ=255)
                                                                                                                   cGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGATCTGGATg                                 <  6:98600/148‑1 (MQ=255)
                                                                                                                    ggTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGATCTGGATgg                                <  8:110153/148‑1 (MQ=255)
                                                                                                                        cAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGATCTGGATGGCAAt                            >  4:365559/1‑148 (MQ=255)
                                                                                                                        cAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGATCTGGATGGCAATc                           <  2:259089/149‑1 (MQ=255)
                                                                                                                               cATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCTCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGATCTGGATGGCAATCatgat                      <  8:29033/147‑1 (MQ=255)
                                                                                                                                aTCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGATCTGGATGGCAATCATGATGTa                   >  5:357649/1‑149 (MQ=255)
                                                                                                                                aTCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACCGCTGGTCTGATCTGGATGGCAATCATGATGTa                   >  2:166422/1‑149 (MQ=255)
                                                                                                                                  ggtgATGGAACTCATAGATTTACATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGATCTGGATGGCAATCATGATGTAc                  >  2:305044/4‑148 (MQ=255)
                                                                                                                                   gtgATGGAACTCATAGATTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGTTCTGCACGTGACTGCTGGTCTGATCTGGATGGCAATCATGATGTAc                  >  4:213885/3‑147 (MQ=255)
                                                                                                                                        ggAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGATCTGGATGGCAATCATGATGTACCAGATc            >  7:327976/1‑148 (MQ=255)
                                                                                                                                                 tGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGATCTGGATGGCAATCATGATGTACCAGATCAACAAGCAc   <  2:101885/148‑1 (MQ=255)
                                                                                                                                                 tGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGATCTGGATGGCAATCATGATGTACCAGATCAACAAGCACg  <  1:291014/149‑1 (MQ=255)
                                                                                                                                                 tGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGATCTGGATGGCAATCATGATGTACCAGATCAACAAGCACg  >  1:271714/1‑149 (MQ=255)
                                                                                                                                                    gTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGt                                                       >  2:287641/1‑93 (MQ=255)
                                                                                                                                                    gTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGt                                                       <  1:287641/93‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TGTTGAAACGCTGTTCCTGCTGCTGTCCTCGATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGATCTGGATGGCAATCATGATGTACCAGATCAACAAGCACG  >  NC_002947/953410‑953702

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: