Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,611,155 T→C 6.8% Y308C (TAC→TGC)  PP_1412 ← 3‑hydroxyisobutyryl‑CoA hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,611,1550TC 6.8% 171.3 / 2.6 59Y308C (TAC→TGC) PP_14123‑hydroxyisobutyryl‑CoA hydrolase
Reads supporting (aligned to +/- strand):  ref base T (31/24);  new base C (2/2);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

CTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGCAACCCTCATGCAGACGGTGGCCGGCATCTGCCAGCAATGGGTCGTCATGCTGCTTGAGCCTTTC  >  NC_002947/1611010‑1611300
                                                                                                                                                 |                                                                                                                                                 
cTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTAtt                                                                                                                                                <  4:113573/149‑1 (MQ=255)
             gcCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAacac                                                                                                                                   >  4:227614/1‑149 (MQ=255)
              ccTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACAcc                                                                                                                                  <  4:110782/149‑1 (MQ=255)
                  gACTACCGCCGTCGGCACTTGCGCCACATCAGGCCATTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTgc                                                                                                                               <  3:9808/148‑1 (MQ=255)
                     tACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGcc                                                                                                                            >  6:37662/1‑148 (MQ=255)
                        cgccgTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTgc                                                                                                                               <  3:128869/142‑1 (MQ=255)
                        cgccgTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTgc                                                                                                                               >  4:128869/1‑142 (MQ=255)
                        cgccgTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAAc                                                                                                                         >  7:182213/1‑148 (MQ=255)
                        cgccgTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAAc                                                                                                                         <  4:301754/148‑1 (MQ=255)
                        cgccgTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAAc                                                                                                                         >  3:301754/1‑148 (MQ=255)
                          ccgTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCCGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGCATTCCATc                                                                                                                                           <  6:165602/128‑1 (MQ=255)
                          ccgTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGCATTCCATc                                                                                                                                           <  6:161426/128‑1 (MQ=255)
                          ccgTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGCATTCCATc                                                                                                                                           >  5:165602/1‑128 (MQ=255)
                          ccgTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGCATTCCATc                                                                                                                                           >  5:161426/1‑128 (MQ=255)
                             tCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAgg                                                                                                                   >  8:263789/1‑149 (MQ=255)
                               ggCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCAGGAACACCTGCGCCAACGACAGGTa                                                                                                                 >  4:177102/1‑149 (MQ=255)
                                       cgcCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGcc                                                                                                          <  7:178923/148‑1 (MQ=255)
                                        gcCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCaa                                                                                                                          <  8:240930/131‑1 (MQ=255)
                                        gcCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCaa                                                                                                                          >  7:240930/1‑131 (MQ=255)
                                             aTCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGa                                                                                                   >  7:207272/1‑149 (MQ=255)
                                             aTCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGa                                                                                                   <  5:54995/149‑1 (MQ=255)
                                                aGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCt                                                                                                >  5:284795/1‑149 (MQ=255)
                                                  gCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGCATCTgg                                                                                              >  6:260499/1‑148 (MQ=255)
                                                  gCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCCGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGc                                                                                              <  3:177102/149‑1 (MQ=255)
                                                      gTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTcc                                                                                           >  8:232582/1‑148 (MQ=255)
                                                       tGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTgcgc                                                                                                                             <  8:200856/113‑1 (MQ=255)
                                                       tGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTgcgc                                                                                                                             >  7:200856/1‑113 (MQ=255)
                                                           aGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCccagacc                                                                                     <  7:232582/149‑1 (MQ=255)
                                                             cgcgGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCccagaccag                                                                                   <  8:182213/149‑1 (MQ=255)
                                                               cgGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGAt                                                                                 <  7:263789/149‑1 (MQ=255)
                                                                     gTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCccagacc                                                                                     >  4:214706/1‑139 (MQ=255)
                                                                     gTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCccagacc                                                                                     <  3:214706/139‑1 (MQ=255)
                                                                     gTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGcc                                                                            <  4:145015/148‑1 (MQ=255)
                                                                     gTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGCATCTGGTCTGGGAGCAGATccgccg                                                                           <  5:260499/149‑1 (MQ=255)
                                                                         tCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAg                                                                       >  6:236158/1‑149 (MQ=255)
                                                                            ttGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCgg                                                                    >  7:75832/1‑149 (MQ=255)
                                                                                 cagcagGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTcc                                                                                           <  6:250502/121‑1 (MQ=255)
                                                                                 cagcagGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTcc                                                                                           >  5:250502/1‑121 (MQ=255)
                                                                                 cagcagGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGCAAc                                                               >  1:241741/1‑149 (MQ=255)
                                                                                   gcagGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGCAAccc                                                             >  7:77144/1‑149 (MQ=255)
                                                                                       gCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGCAACCCTCAt                                                         <  5:236158/149‑1 (MQ=255)
                                                                                       gCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGCAACCCTCAt                                                         <  3:227614/149‑1 (MQ=255)
                                                                                                 aCGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCCGGCAACCCTCATGCAGACCGTg                                               >  1:220353/1‑149 (MQ=255)
                                                                                                 aCGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTAAGGCTCATTCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTAACGCGCACGGCGGATCTGATCCCAGACCAGATGCGCCGTCATCTGGCAACCAAAATGCAGACtgtg                                               >  1:216190/1‑149 (MQ=255)
                                                                                                     cTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGCAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGCAACCCTCATGCAGACGGTGGCCg                                           >  1:248486/1‑149 (MQ=255)
                                                                                                        cACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTGTTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGCaa                                                                >  2:125205/1‑125 (MQ=255)
                                                                                                        cACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTCTTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGCaa                                                                <  1:125205/125‑1 (MQ=255)
                                                                                                           tGAATTCAGGGTGGCGGCAACAGTTCAGGCTAATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGCAACCCTCATGCAGACGGTGGCCGGCATCt                                     >  5:200418/1‑149 (MQ=255)
                                                                                                                     gTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGCAACCCTCATGCAGACGGTGGCCGGCATCTGCCAGCAATg                           <  2:241741/149‑1 (MQ=255)
                                                                                                                          ggcAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGcc                                                                            >  7:274311/1‑95 (MQ=255)
                                                                                                                          ggcAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGcc                                                                            <  8:274311/95‑1 (MQ=255)
                                                                                                                             aaCAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGCAACCCTCAAGCAGACGGTGGCCGGCATCTGCCAGCAATGGGTCGTCa                   >  8:86392/1‑149 (MQ=255)
                                                                                                                               cAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGCAACCCTCATGCAGACGGTGGCCGGCATCTGCCAGCAATGGGTCGTCAtg                 <  7:86392/149‑1 (MQ=255)
                                                                                                                                      ggCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGCAACCCTCATGCAGACGGTGGCCGGCATCTGCCAGCAATGGGTCGTCATGCTGCtt           >  7:278587/1‑148 (MQ=255)
                                                                                                                                         tCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGCAACCCTCATGCAGACGGTGGCCGGCATCTGCCAGCAATGGGTCGTCATGCTGCTTGAGc       >  8:189103/1‑149 (MQ=255)
                                                                                                                                           aTGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGCAACCCTCATGCAGACGGTGGCCGGCATCTGCCAGCAATCGGTCGTCATGCTTCTTGAGCCt     >  6:278621/1‑149 (MQ=255)
                                                                                                                                             gCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGCAACCCTCATGCAGACGGTGGCCGGCATCTGCCAGCAATGGGTCGTCAtgc                >  7:24633/1‑136 (MQ=255)
                                                                                                                                             gCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGCAACCCTCATGCAGACGGTGGCCGGCATCTGCCAGCAATGGGTCGTCAtgc                <  8:24633/136‑1 (MQ=255)
                                                                                                                                                gTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGCAACCCTCATGCAGACGGTGGCCGGCATCTGCCAGCAATGGGTCGTCATGCTGCTTGAGCCTTTc  >  1:122631/1‑147 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGCAACCCTCATGCAGACGGTGGCCGGCATCTGCCAGCAATGGGTCGTCATGCTGCTTGAGCCTTTC  >  NC_002947/1611010‑1611300

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: