Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,968,096 G→T 11.2% E90D (GAG→GAT wecD → TDP‑fucosamine acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,968,0960GT11.2% 103.3 / 2.2 36E90D (GAG→GATwecDTDP‑fucosamine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (18/14);  new base T (2/2);  total (20/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.69e-01

TTCGTGCCCGTCCGCGCCAGTATTGAGCCACTAACCTGGGAAAACGCCTTCTTTGGTGTTAACAGCGCCATCGTGCGCATTACGTCTGAAGCACCGCTCCTGACGCCAGACGCGTTAGCGCCGTGGTCACGGGTGCAGGCCAAAATTGCGGCATCAAATACGGGTGAACTGGATGCCCTGCAACAGCTGGGATTCTCCCTGGTAGAAGGTGAAGTTGATTTGGCGCTACCCGTGAACAATGCCAGTGATAGCGGTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAA  >  CP009273/3967824‑3968364
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                            
ctgggcccgntCCGCGCCAGTAGTGTGCCACTCCCCTCGGATATCGCCTTCTTTGGTGTTAACAGCGCCATCGTGCCGATTACGTCTGAAGCAACGGTCCTGAAGCAGGCCGCGTTAGCGCCGTGGTCACGGGTGCAGGCAAAAATTGCGGCATCAAATACGGGTGAACTGGATGCCCTGCAACAGCTGGGAGTCTCCCTGGTAGAAGGTGGAGTGGATTTGGCGCAACCCGTGATCAATGCTAGGGATAGCGGTGCTGTACTGGCGCAagag                                                                                                                                                                                                                                                                              <  2:44121/263‑1 (MQ=255)
              cgcCAGTATTGAGCCACTAACCTGGGAAAACGCCTTCTTTGGTGTTAACAGCGCCATCGTGCGCATTACGTCTGAAGCACCGCTCCTGACGCCAGACGCGTTAGCGCCGTGGTCACGGGTGCAGGCCAAAATTGCGGCATCAAATACGGGTGAACTGGATGCCCTGCAACAGCTGGGATTCTCCCTGGTAGAAGGAGAAGTTGATTTGGCGCTACCCGTGAACAATGCCAGTGATAGCGGTGCTGTAGTGGCACAAGAGACCGATATTCCCGCAt                                                                                                                                                                                                                                                              >  1:486822/1‑275 (MQ=255)
              cgcCAGTATTGAGCCACTAACCTGGGAAAACGCCTTCTTTGGTGTTAACAGCGCCATCGTGCGCATTACGTCTGAAGCACCGCTCCTGACGCCAGACGCGTTAGCGCCGTGGTCACGGGTGCAGGCCAAAATTGCGGCATCAAATACAGGTGAACTGGNTGCCCTGCAACAGCTGGGATTCTCCCTGGTAGAAGGTGAAGTTGATTTGGCGCTACCCGTGAACCATGCCAGTGCTAGCGGTGCTGTAGTGGCACAAGAGACCGATATTCCCGCa                                                                                                                                                                                                                                                               >  1:444437/1‑274 (MQ=255)
                           ccACTAACCTGGGAAAACGCCTTCTTTGGTGTTAACAGCGCCATCGTGCGCATTACGTCTGAAGCACCGCTCCTGACGCCAGACGCGTTAGCGCCGTGGTCACGGGTGCAGGCCAAAATTGCGGCATCAAATACGGGTGAACTGGATGCCCTGCAACAGCTGGGATTCTCCCTGGTAGAAGGTGAAGTTGATTTGGCGCTACCCGTGAACAATGCCAGTGATAACTGTGCTGTAGGGGCACAAGAGAACAATATTCCAGCATTACGTCAGtttt                                                                                                                                                                                                                                                  >  2:343464/1‑272 (MQ=255)
                                           aCGCCTTCTTTGGTGTTAACAGCGCCATCGTGCGCATTACGTCTGAAGCACCGCTCCTGACGCCAGACGCGTTAGCGCCGTGGTCACGGGTGCAGGCCAAAATTGCGGCATCAAATACGGTTGAACTGGATGCCCTGCAACAGCTGGGATTCTCCCTGGTAGAAGGTGAAGTTGATTTGGCGCTACCCGTGAACAATGCCAGTGATAGCGGTGCTGTAGTGGCACAAGAGACCGATAtt                                                                                                                                                                                                                                                                     >  1:169960/1‑239 (MQ=255)
                                                            aaCAGCGCCATCGTGCGCATTACGTCTGAAGCACCGCTCCTGACGCCAGACGCGTTAGCGCCGTGGTCACGGGTGCAGGCCAAAATTGCGGCATCAAATACGGGTGAACTGGATGCCCTGCAACAGCTGGGATTCTCCCTGGTAGAAGGAGAAGTTGATTTGGCGCTAACCGTGAACACTGCCAGTGATAGCGGTGCTGTAGTGGCACAAGAGACCGATATGCCCGAATTACGTCAGTTCGCCAGGGCAGCATTTGCTCAAGGCCCTTTTctgt                                                                                                                                                                                                                 >  2:434482/1‑271 (MQ=255)
                                                                  gcCAGCGTGCGCATTACGTCTGACGCCCCGCTCCGGACGCCAGACGCGTTAGGGCCGTGGTCACGGGTGCAGGCCAAAATTGCGGCATCAAATACGGGTGAACTGGATGCCCTGCAACAGCTGGGATTCTCCCTGGTAGAAGGTGAAGTTGATTTGGCGCTACCCGTGAACAATGCCAGTGATAGCGGTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGt                                                                                                                                                                                                           <  2:486809/274‑1 (MQ=255)
                                                                                          gCACCGCTCCTGACGCCAGACGCGTTAGCGCCGTGGTCACGGGTGCAGGCCAAAATTGCGGCATCAAATACGGGTGAACTGGATGCCCTGCAACAGCTGGGATTCTCCCTGGTAGAAGGTGTAGTTGCTTTGGCGCTACCCGTGAACAATGCCAGTGATAGCGGTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGACAGCAGTc                                                                                                                                                                                   <  1:434495/274‑1 (MQ=255)
                                                                                             acGCTCCTGACGCCAGACGCGTTAGCGCCGTGGTCACGGGTGCAGGCCAAAATTGCGGCATCAAATACGGGTGAACTGGATGNCCTGCAACNGCTGGGATTCTCCCTGGTAGAAGGTGAAGTTGATTTGGCGCTACCCGTGAACAATGCCAGTGATAGCGGTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTTCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCtt                                                                                                                                                                               <  1:322532/274‑1 (MQ=255)
                                                                                                                               aaCGGGGGCAGGCCAAAAGTGCGGCATCAAATACGGGTTAACTGGATGGCTTGCAACAGGTGGGATTCTCCCTGGTAGAAGGTGAAGTTGAGTTGGCGCTACCCGTGAACAATGCCAGTGATAGCGGTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGcgcg                                                                                                                                               <  2:444424/272‑1 (MQ=255)
                                                                                                                                                                      aaCTGGATGCCCTGCAACAGCTGGGATTCTCCCTGGTAGAAGGTGAAGTTGATTTGGCGCTACCCGTGAACAATGCCAGTGATAGCGGTGCTGTAGTGGCACAAGAGACCGata                                                                                                                                                                                                                                                                       >  1:165477/1‑114 (MQ=255)
                                                                                                                                                                          ggATGCCCTGCAACAGCTGGGATTCTCCCTGGTAGAAGGTGAAGTTGATTTGGCGCTACCCGTGAACAATGCCAGTGATAGCGGTGCTGTAGTGGCACAAGAGACCGATATTcc                                                                                                                                                                                                                                                                   <  2:165471/114‑1 (MQ=255)
                                                                                                                                                                                          cTGGGATTCTCCCTGGTAGAAGGTGAAGTTGATTTGGCGCTACCCGTGAACAATGCCAGTGATAGCGGTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGCGGTATGCGCCTGACGCCAGCAGTCGCTTTTATg                                                                                                                                                                          >  2:6738/1‑187 (MQ=255)
                                                                                                                                                                                           tGGGATTCTCCCTGGTAGAAGGTGAGGTTGATTTGGCGCTACCCGTGAACAATGCCAGTGATAGCGGTGCTGTAGTGGCACAAGAGACCGNTATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCCTCCGGCGATATTCGCGGCTATGt                                                                                  <  1:343476/274‑1 (MQ=255)
                                                                                                                                                                                              gATTCTCCCTGGTAGAAGGTGAAGTTGATTTGGCGCTACCCGTGAACAATGCCAGTGATAGCGGTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGcaca                                                                                                                                                                      <  1:6738/187‑1 (MQ=255)
                                                                                                                                                                                                  ctcCCTGGTAGAAGGTGAAGTTGATTTGGCGCTACCCGTGAACAATGCCAGTGATAGCGGTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCgg                                                                                        >  1:215757/1‑261 (MQ=255)
                                                                                                                                                                                                      cTGGTGGAAGGTGAAGTTGATTTGGCGCTACCCGTGAACAATGCCAGTGATAGCGGTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTAt                                                                                    <  2:215750/261‑1 (MQ=255)
                                                                                                                                                                                                          tAGAAGGTGAAGTTGATTTGGCGCTACCCGTGAACAATGCCAGTGATAGCGGTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTc                                                                  >  1:114219/1‑275 (MQ=255)
                                                                                                                                                                                                                                              aTGCCAGTGATAGCGGTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTcgc                                                                                          >  1:287481/1‑215 (MQ=255)
                                                                                                                                                                                                                                               tGCCAGTGATAGCGGTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGca                                                                                                                                                                        >  2:222599/1‑136 (MQ=255)
                                                                                                                                                                                                                                               tGCCAGTGATAGCGGTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACg                                                                                                                >  2:389463/1‑192 (MQ=255)
                                                                                                                                                                                                                                               tGCCAGTGATAGCGGTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATTCGCGAATTGACcttcttccgttac                             >  2:344789/1‑267 (MQ=255)
                                                                                                                                                                                                                                                 ccAGTGATAGCGGTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCGGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATAATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTATTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGAcgc                           >  1:277259/1‑275 (MQ=255)
                                                                                                                                                                                                                                                  cAGTGATAGCGGTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCt                                                                                      <  2:287473/215‑1 (MQ=255)
                                                                                                                                                                                                                                                   aGTGATAGCGGTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGt                                                                                                                                                                    <  1:222606/136‑1 (MQ=255)
                                                                                                                                                                                                                                                   aGTGATAGCGGTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTgcg                                                                                                            <  1:389476/192‑1 (MQ=255)
                                                                                                                                                                                                                                                         aGCGGTGCTGTAGTGGCACAAGATACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTg                                                                                                                                                   >  2:272787/1‑147 (MQ=255)
                                                                                                                                                                                                                                                          gCGGTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTgc                                                                                                                                                  >  1:200161/1‑147 (MQ=255)
                                                                                                                                                                                                                                                             gTGCTGTAGTGGCACAAGATACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGcgcg                                                                                                                                               <  1:272795/147‑1 (MQ=255)
                                                                                                                                                                                                                                                              tGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCgg                                                                                                                                              <  2:200155/147‑1 (MQ=255)
                                                                                                                                                                                                                                                                cTGTAGTGGCACAATAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCcagca                                                                                                                                                                                      >  1:58136/1‑105 (MQ=255)
                                                                                                                                                                                                                                                                    aGTGGCACAATAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCg                                                                                                                                                                                  <  2:58134/105‑1 (MQ=255)
                                                                                                                                                                                                                                                                      tGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGCCGCGGTGCAGGTgctgatgtaca       >  2:158609/1‑268 (MQ=255)
                                                                                                                                                                                                                                                                        gCACAAGATACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGt                                                                                                                                                    >  1:274615/1‑131 (MQ=255)
                                                                                                                                                                                                                                                                          acaAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCNGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGTACGCGGTGCAGGTACTGAGCTGATGCaaa  <  1:344801/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                           caAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCg                                                                                                                                                                                                                   >  2:84668/1‑65 (MQ=255)
                                                                                                                                                                                                                                                                            aaGATACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTgcgc                                                                                                                                                <  2:274607/131‑1 (MQ=255)
                                                                                                                                                                                                                                                                             agagACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTg                                                                                                                                                                                                                 >  2:463642/1‑65 (MQ=255)
                                                                                                                                                                                                                                                                             agagACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTc                                                                                                                                                                                   >  1:198016/1‑95 (MQ=255)
                                                                                                                                                                                                                                                                               agaCCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTgcg                                                                                                                                                                                                               <  1:84670/65‑1 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                            
TTCGTGCCCGTCCGCGCCAGTATTGAGCCACTAACCTGGGAAAACGCCTTCTTTGGTGTTAACAGCGCCATCGTGCGCATTACGTCTGAAGCACCGCTCCTGACGCCAGACGCGTTAGCGCCGTGGTCACGGGTGCAGGCCAAAATTGCGGCATCAAATACGGGTGAACTGGATGCCCTGCAACAGCTGGGATTCTCCCTGGTAGAAGGTGAAGTTGATTTGGCGCTACCCGTGAACAATGCCAGTGATAGCGGTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAA  >  CP009273/3967824‑3968364

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: