Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,202,857 G→T 7.6% R229L (CGG→CTG)  xcpY → type II secretion pathway protein XcpY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,202,8570GT 7.6% 152.2 / 4.2 53R229L (CGG→CTG) xcpYtype II secretion pathway protein XcpY
Reads supporting (aligned to +/- strand):  ref base G (22/27);  new base T (2/2);  total (24/29)
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.46e-01

GTGCCGCAGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCGGCACAG  >  NC_002947/1202709‑1202996
                                                                                                                                                    |                                                                                                                                           
gTGCCGCAGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCg                                                                                                                                             >  6:242918/1‑149 (MQ=255)
gTGCCGCAGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCAACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCg                                                                                                                                             >  7:138141/1‑149 (MQ=255)
   ccGCAGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGgct                                                                                                                                          >  6:1901/1‑149 (MQ=255)
    cGCAGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGgctg                                                                                                                                         >  4:243502/1‑149 (MQ=255)
      cAGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGgctgct                                                                                                                                       <  5:298519/149‑1 (MQ=255)
       aGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGctgctg                                                                                                                                      <  6:181023/149‑1 (MQ=255)
       aGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGctgctg                                                                                                                                      <  7:328475/149‑1 (MQ=255)
        gCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTgg                                                                                                                                     <  3:313496/149‑1 (MQ=255)
         ccACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCCGCCGCCGCGGCTGCTgg                                                                                                                                     >  8:173027/1‑148 (MQ=255)
         ccACCGTGGGCAGTGCTGCACAGGACGGCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGGTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTgg                                                                                                                                     <  7:173027/148‑1 (MQ=255)
           aCCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCAAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAg                                                                                                                                  >  2:330481/1‑149 (MQ=255)
              gTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCtt                                                                                                                               >  7:60219/1‑149 (MQ=255)
                gggCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCTGCCAGCAGCCGCTGCTGCTGGcgggaggc                                                                                                                             <  3:287729/149‑8 (MQ=255)
                 ggCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGgctgct                                                                                                                                       >  5:152796/1‑138 (MQ=255)
                 ggCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGgctgct                                                                                                                                       <  6:152796/138‑1 (MQ=255)
                   cAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTg                                                                                                                          <  5:1901/149‑1 (MQ=255)
                         tgcACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACAAACTGCCCGGCCAGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCAGGCAGCTTGCCTgagaacg                                                                                                                    >  5:247905/1‑143 (MQ=255)
                          gcACAGGACGTCGCTGAGCGGCAATTGGGCAACCGTATTCGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCGGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTgg                                                                                                                   >  1:205631/1‑149 (MQ=255)
                                 aCGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGt                                                                                                            >  1:346976/1‑149 (MQ=255)
                                                aaTTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCtgtg                                                                                             <  3:243502/149‑1 (MQ=255)
                                                    ggCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCAGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGtt                                                                                          <  1:260843/148‑1 (MQ=255)
                                                     gCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGa                                                                                        <  2:65938/149‑1 (MQ=255)
                                                       cAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCCGCTGCTGCTGGcggga                                                                                                                                >  3:142078/1‑103 (MQ=255)
                                                       cAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCCGCTGCTGCTGGcggga                                                                                                                                <  4:142078/107‑5 (MQ=255)
                                                        aaCCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGt                                                                                      >  7:246014/1‑148 (MQ=255)
                                                               tttGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATg                                                                              <  2:155270/149‑1 (MQ=255)
                                                                ttGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAAt                                                                               >  1:317173/1‑147 (MQ=255)
                                                                    gccgcTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGt                                                                          >  5:43103/1‑148 (MQ=255)
                                                                    gccgcTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTc                                                                         <  3:54832/149‑1 (MQ=255)
                                                                    gccgcTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGGTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGt                                                                          >  1:231931/1‑148 (MQ=255)
                                                                    gccgcTCGACACACTGCCCGGCCTGTTCAACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTc                                                                         <  8:138141/149‑1 (MQ=255)
                                                                          cGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTgg                                                                                            >  4:281787/1‑124 (MQ=255)
                                                                          cGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTgg                                                                                            <  3:281787/124‑1 (MQ=255)
                                                                           gACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGc                                                                            >  4:68839/1‑139 (MQ=255)
                                                                           gACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGc                                                                            <  3:68839/139‑1 (MQ=255)
                                                                            acacacTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCaa                                                                 >  1:275060/1‑149 (MQ=255)
                                                                                acTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGGTTGGGGTGGCCGGTGGTTGAGTCAGCAATGGCGTCAGGCGCAActct                                                             >  3:303344/1‑149 (MQ=255)
                                                                                 cTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTg                                                            <  2:275060/149‑1 (MQ=255)
                                                                                    cccGGCCTGGTCGACCGCGCCCGCACGGGGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCGACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGgc                                                          <  8:60219/148‑1 (MQ=255)
                                                                                     ccGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGgcg                                                         >  6:378396/1‑148 (MQ=255)
                                                                                     acGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCAGCTGGCAGCTTGCCTGGGGCTGGCTAACGTTTTGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAAATCTGgccc                                                        >  5:96284/2‑147 (MQ=255)
                                                                                                      gcCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGt                                       <  2:346976/149‑1 (MQ=255)
                                                                                                               gTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGa                              <  8:246014/149‑1 (MQ=255)
                                                                                                                  cgcTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAg                             >  7:325747/1‑147 (MQ=255)
                                                                                                                  cgcTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCa                           <  7:160544/149‑1 (MQ=255)
                                                                                                                   gcTCGTTGCCGGCGGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAg                          <  2:205631/149‑1 (MQ=255)
                                                                                                                       gTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCa                      <  5:378396/149‑1 (MQ=255)
                                                                                                                       gTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGGTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCa                      <  2:231931/149‑1 (MQ=255)
                                                                                                                             ggccgCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTc                <  2:317173/149‑1 (MQ=255)
                                                                                                                                  ccTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCAc           <  8:325747/149‑1 (MQ=255)
                                                                                                                                     cccGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGGTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCAcgcg        <  4:303344/149‑1 (MQ=255)
                                                                                                                                      ccGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCgg       >  1:52175/1‑149 (MQ=255)
                                                                                                                                       tgCCAGCAGCCGCTGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCgg       >  4:206684/2‑148 (MQ=255)
                                                                                                                                          cagcagcagCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCGGca     >  7:27459/1‑147 (MQ=255)
                                                                                                                                           agcagcagcGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCGGCACAg  <  8:315431/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                           
GTGCCGCAGCCACCGTGGGCAGTGCTGCACAGGACGTCGCTGAGCGGCAATTGGCCAACCGTATTTGCGCCGCTCGACACACTGCCCGGCCTGTTCGACCGCGCCCGCAAGGTGCGCTCGTTGCCGGCCGCCTCCCGCCAGCAGCAGCGGCTGCTGGCAGCTTGCCTGGTGCTGGCTACCGTTTGGGGTGGCCTGTGGTTGAGTCAGCAATGGCGTCAGGCGCAACTCTGGCGCAGCCAGGTCATTGCCGTGACCGGGGAGCAGGCCAGCCCTCGGCACGCGGCACAG  >  NC_002947/1202709‑1202996

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: