Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,000,765 T→A 6.8% V77V (GTA→GTT PP_2624 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,000,7650TA 6.8% 174.5 / 3.0 60V77V (GTA→GTTPP_2624hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (27/27);  new base A (2/2);  total (30/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.75e-01

CAGGGCGACCAGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACC  >  NC_002947/3000618‑3000912
                                                                                                                                                   |                                                                                                                                                   
cAGGGCGACCAGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGt                                                                                                                                                     <  8:138583/148‑1 (MQ=255)
     cGACCAGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGCACCGTc                                                                                                                                               <  4:132639/149‑1 (MQ=255)
           gTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTc                                                                                                                                          >  2:341184/1‑148 (MQ=255)
              cTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTcc                                                                                                                                         >  2:9156/1‑146 (MQ=255)
              cTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTcc                                                                                                                                         <  1:9156/146‑1 (MQ=255)
              cTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAgg                                                                                                                                      >  2:149281/1‑149 (MQ=255)
                                        aGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCt                                                                                                            <  2:149320/149‑1 (MQ=255)
                                          cTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTcac                                                                                                                                 >  3:295192/1‑126 (MQ=255)
                                          cTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTcac                                                                                                                                 <  4:295192/126‑1 (MQ=255)
                                             gggCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTg                                                                                                        >  5:348652/1‑148 (MQ=255)
                                               gCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAat                                                                                                     <  6:262358/149‑1 (MQ=255)
                                                cACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGTAACGGTACCGTCCAGGc                                                                                                                                     <  8:328818/116‑1 (MQ=255)
                                                cACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGTAACGGTACCGTCCAGGc                                                                                                                                     >  7:328818/1‑116 (MQ=255)
                                                 aCCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGAACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACAGCGTGTTGaacga                                                                                                   >  7:24841/1‑146 (MQ=255)
                                                  ccTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGAACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACAGCGTGTTGaacgac                                                                                                  <  8:24841/149‑5 (MQ=255)
                                                         cgcgATCCGGTCCGGCCTGAATTCGTTGATGGTGCTGAACTGCAGGTGGACGGCCTGGTCTGTATCAGGGGCGGTGCCGAGCCTGGACGGTACCGTCATGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTAt                                                                                           <  1:88591/149‑1 (MQ=255)
                                                          gcgATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGc                                                                                                             >  5:30759/1‑130 (MQ=255)
                                                          gcgATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGc                                                                                                             <  2:276260/130‑1 (MQ=255)
                                                          gcgATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGc                                                                                                             >  1:276260/1‑130 (MQ=255)
                                                          gcgATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGc                                                                                                             <  6:30759/130‑1 (MQ=255)
                                                           cgATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATc                                                                                          >  1:305387/1‑148 (MQ=255)
                                                             aTCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGaa                                                                                                      >  8:97616/1‑134 (MQ=255)
                                                             aTCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGaa                                                                                                      <  7:97616/134‑1 (MQ=255)
                                                             aTCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGAt                                                                                       >  8:85904/1‑149 (MQ=255)
                                                                       gCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGAcc                                                                             >  7:9983/1‑149 (MQ=255)
                                                                          cGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTGATAGATGCGcctgacctga                                                                          >  4:57158/1‑149 (MQ=255)
                                                                          cGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCCTCACCTCCACGCACACATTGGCGTCGTACAACACGCTGTTGAAAGTATGTTTATCGATCCGCATTAcatta                                                                          >  4:61875/1‑145 (MQ=255)
                                                                           gAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACAGCGTGTTGaac                                                                                                     <  7:174585/121‑2 (MQ=255)
                                                                           gAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACAGCGTGTTGaac                                                                                                     >  8:174585/1‑120 (MQ=255)
                                                                           gAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGAc                                                                         <  7:85904/149‑1 (MQ=255)
                                                                               tCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTc                                                                     <  1:149281/149‑1 (MQ=255)
                                                                                 gTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGAc                                                                              >  5:117357/1‑138 (MQ=255)
                                                                                 gTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGAc                                                                              <  6:117357/138‑1 (MQ=255)
                                                                                  ttGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGtt                                                                                                         <  4:276487/110‑1 (MQ=255)
                                                                                  ttGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGtt                                                                                                         >  3:276487/1‑110 (MQ=255)
                                                                                      tGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAgc                                                              <  1:341184/149‑1 (MQ=255)
                                                                                              aaCTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCAc                                                          <  5:257127/145‑1 (MQ=255)
                                                                                              aaCTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCAc                                                          >  6:257127/1‑145 (MQ=255)
                                                                                              aaCTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACgg                                                        <  2:305387/147‑1 (MQ=255)
                                                                                              aaCTGCAGGTGGACGACCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAatgta                                                                                                  <  8:323753/105‑1 (MQ=255)
                                                                                              aaCTGCAGGTGGACGACCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAatgta                                                                                                  >  7:323753/1‑105 (MQ=255)
                                                                                                 tGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGt                                                    >  2:206745/1‑148 (MQ=255)
                                                                                                   cAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGGGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCa                                                  >  6:47147/1‑148 (MQ=255)
                                                                                                    aGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATc                                                >  1:150255/1‑149 (MQ=255)
                                                                                                      gTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGaa                                                                                                      <  1:263952/93‑1 (MQ=255)
                                                                                                      gTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGaa                                                                                                      >  2:263952/1‑93 (MQ=255)
                                                                                                        ggACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCt                                                      <  6:200186/139‑1 (MQ=255)
                                                                                                        ggACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCt                                                      >  5:200186/1‑139 (MQ=255)
                                                                                                                aggtCGGTATGCGGGGCGGTGCCGAGCCTGGACGGTACCTTCTCTTGCAGGCTCACATTGGCGTCGTTCAACTCGCTGTTGTATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTc                                    <  3:61875/146‑1 (MQ=25)
                                                                                                                aGCTCGGTATCCGGGGCGGTGCCGCGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGGTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTc                                    <  3:57158/149‑1 (MQ=255)
                                                                                                                            gtgtcGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGccc                                <  5:59249/137‑1 (MQ=255)
                                                                                                                            ggggCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCccagcaccag                        >  2:307712/1‑149 (MQ=255)
                                                                                                                            ggggCGGTGCCGAGCCAGGACGGTAACGTCACGTCCAGGCTCACATTGGCGTCGTTAAACACGCTGTTGAATGTATGTTTATAGATTCGCATTACCATACCGTCCACAACCGCACGGCTGTCATACAGGCCGAACTCACccaataacg                         >  6:59249/1‑142 (MQ=255)
                                                                                                                               gCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCAt                     <  8:9983/149‑1 (MQ=255)
                                                                                                                                     ccGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTtgtg               >  8:254847/1‑149 (MQ=255)
                                                                                                                                          ccTGGACGTGACGGTACCGTCCAGGCTCGCATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCgg            >  1:318929/1‑147 (MQ=18)
                                                                                                                                           cTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAg         <  2:150255/149‑1 (MQ=255)
                                                                                                                                              gACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGc        <  5:47147/147‑1 (MQ=255)
                                                                                                                                                 ggTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTAt                                                                                           <  4:255471/61‑1 (MQ=255)
                                                                                                                                                 ggTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTAt                                                                                           >  3:255471/1‑61 (MQ=255)
                                                                                                                                                  gTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCAcc  <  7:254847/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CAGGGCGACCAGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACC  >  NC_002947/3000618‑3000912

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: