Predicted mutation
evidence seq id position mutation annotation gene description
RA BS168_reference_genome 399,950 A→G A473A (GCA→GCG srfAC → surfactin synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome399,9500AG96.9% 230.6 / ‑1.8 64A473A (GCA→GCGsrfACsurfactin synthetase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (31/31);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-01

CCAGAATCCTGATCAGCCTGTCAGCACGATCAACCTCGTTGACGACAGGGAGAGAGAATTTTTGCTAACCGGCTTAAACCCGCCGGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCAGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCAGACCCTGTCTTATCGCGAATTAGATGAGGAAGCGAACCGCATTGCACGCCGGCTGCAAAAACACGGTGCGGGCAAAGGCTCTGTTGTA  >  BS168_reference_genome/399818‑400082
                                                                                                                                    |                                                                                                                                    
ccAGAATCCTGATCAGCCTGTCAGCACGATCAACCTCGTTGACGACAGGGAGAGAGAATTTTTGCTAACCGGCTTAAACCCGCCGGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAAc                                                                                                                                >  5:122942/1‑139 (MQ=255)
      tCCTGATCAGCCTGTCAGCACGATCAACCTCGTTGACGACAGGGAGAGAGAATTTTTGCTAACCGGCTTAAACCCGCCGGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAAt                                                                                                                          >  3:180095/1‑139 (MQ=255)
      tCCTGATCAGCCTGGCAGCACGATCAACCTCGTTGACGACAGGGAGAGAGAATTTTTGCTAACCGGCTTAAACCCGCCGGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAAt                                                                                                                          <  3:44427/139‑1 (MQ=255)
          gATCAGCCTGTCAGCACGATCAACCTCGTTGACGACAGGGAGAGAGAATTTTTGCTAACCGGCTTAAACCCGCCGGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCg                                                                                                                       >  3:281512/1‑138 (MQ=255)
           aTCAGCCTGTCAGCACGATCAACCTCGTTGACGACAGGGAGAGAGAATTTTTGCTAACCGGCTTAAACCCGCCGGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGa                                                                                                                     >  2:206715/1‑139 (MQ=255)
             cAGCCTGTCAGCACGATCAACCTCGTTGACGACAGGGAGAGAGAATTTTTGCTAACCGGCTTAAACCCGCCGGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCg                                                                                                                       >  7:374332/1‑135 (MQ=255)
              aGCCTGTCAGCACGATCAACCTCGTTGACGACAGGGAGAGAGAATTTTTGCTAACCGGCTTAAACCCGCCGGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGc                                                                                                                  >  6:313620/1‑139 (MQ=255)
                ccTGTCAGCACGATCAACCTCGTTGACGACAGGGAGAGAGAATTTTTGCTAACCGGCTTAAACCCGCCGGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCAc                                                                                                                >  2:369022/1‑139 (MQ=255)
                             tCAACCTCGTTGACGACAGGGAGAGAGAATTTTTGCTAACCGGCTTAAACCCGCCGGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGcgtacgta                                                                                                   >  2:205759/1‑139 (MQ=255)
                               aaCCTCGTTGACGACAGGGAGAGAGAATTTTTGCTAACCGGCTTAAACCCGCCGGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTAtt                                                                                                 <  5:230687/139‑1 (MQ=255)
                               aaCCTCGTTGACGACAGGGAGAGAGAATTTTTGCTAACCGGCTTAAACCCGCCGGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTAt                                                                                                  <  3:369641/138‑1 (MQ=255)
                                aCCTCGTTGACGACAGGGAGAGAGAATTTTTGCTAACCGGCTTAAACCCGCCGGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTa                                                                                                   <  7:187980/136‑1 (MQ=255)
                                     gTTGACGACAGGGAGAGAGAATTTTTGCTAACCGTCTTAAACCCGCCGGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGc                                                                                            <  1:206715/138‑1 (MQ=255)
                                     gTTGACGACAGGGAGAGAGAATTTTTGCTAACCGGCTTAAACCCGCCGGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGc                                                                                            <  4:281512/138‑1 (MQ=255)
                                         acgacAGGGAGAGAGAATTTTTGCTAACCGGCTTAAACCCGCCGGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTGACGTATTACGGcccgca                                                                                       >  2:102471/1‑135 (MQ=255)
                                          cgacAGGGAGAGAGAATTTTTGCTAACCGGCTTAAACCCGCCGGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCAGAcc                                                                                      <  4:347100/139‑1 (MQ=255)
                                                ggAGAGAGAATTTTTGCTAACCGGCTTAAACCCGCCGGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTc                                                                                                <  3:334506/123‑1 (MQ=255)
                                                ggAGAGAGAATTTTTGCTAACCGGCTTAAACCCGCCGGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTc                                                                                                >  4:334506/1‑123 (MQ=255)
                                                     gagaATTTTTGCTAACCGGCTTAAACCCGCCGGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCgg                                                                                                             <  3:352972/105‑1 (MQ=255)
                                                     gagaATTTTTGCTAACCGGCTTAAACCCGCCGGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCgg                                                                                                             >  4:352972/1‑105 (MQ=255)
                                                     gagaATTTTTGCTAACCGGCTTAAACCCGCCGGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCAGACCCTGTCTTATcg                                                                           <  5:3604/139‑1 (MQ=255)
                                                     gagaATTTTTGCTAACCGGCTTAAACCCGCCGGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCAGACCCTGTCTTATcg                                                                           <  4:153970/139‑1 (MQ=255)
                                                        aaTTTTTGCTAACCGGCTTAAACCCGCCGGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCAGTGAACGCCAATCCGGATGCAc                                                                                                                >  2:204786/1‑99 (MQ=255)
                                                        aaTTTTTGCTAACCGGCTTAAACCCGCCGGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCAGTGAACGCCAATCCGGATGCAc                                                                                                                <  1:204786/99‑1 (MQ=255)
                                                          tttttGCTAACCGGCTTAAACCCGCCGGCTTAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGgcg                                                                                                           <  1:125371/102‑1 (MQ=255)
                                                          tttttGCTAACCGGCTTAAACCCGCCGGCTTAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGgcg                                                                                                           >  2:125371/1‑102 (MQ=255)
                                                            tttGCTCACCGGCTTAACCCCGCCGGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGGTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCAGACCCTGTCTTATCGCGAAtt                                                                     <  5:382380/138‑1 (MQ=255)
                                                                  aaCCGGCTTAAACCCGCCGGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGCACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCAGCACCTGTCTTATCGCGCATTAgaggag                                                              >  2:37295/1‑139 (MQ=255)
                                                                           aaaCCCGCCGGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCAGACCCTGTCTTATCGCGAATTAGATGAGGAAGCGaa                                                      <  6:99/138‑1 (MQ=255)
                                                                           aaaCCCGCCGGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCAGACCCTGTCTTATCGCGAATTAGATGAGGAAGCGaa                                                      <  6:313782/138‑1 (MQ=255)
                                                                                  ccgGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCAGAcc                                                                                      <  5:273390/99‑1 (MQ=255)
                                                                                  ccgGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCAGAcc                                                                                      >  6:273390/1‑99 (MQ=255)
                                                                                  ccgGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCAGACCCTGTCTTATCGCGAATTAGATGAGGAAGCGAACCGCATTg                                              >  1:64939/1‑139 (MQ=255)
                                                                                  ccgGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCAGACCCTGTCTTATCGCGAATTAGATGAGGAAGCGAACCGCATTg                                              >  6:74758/1‑139 (MQ=255)
                                                                                     gCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCAGACCCTGTCTTATCGCGAATTAGATGAGGAAGCGAACCGCATTg                                              >  7:172372/1‑136 (MQ=255)
                                                                                               atcaaaCAAAGCCTCTGACTTTTTGTTTCATGTAAGCGGTGAACGACAATCCGGATGCACGGGCGCTTACGTATTCCGGCCAGACCCTGTATTATCGCGAATTAGATGAGGAAGCGAACCGCATTGCACGCCGGCTg                                   <  5:313620/134‑1 (MQ=255)
                                                                                                tGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCAGACGCTGTCTTATGGCTCATTAGATGAGGAAGCGAGCCGCATTGCACGCCGGCTGCa                                 >  6:182599/1‑138 (MQ=255)
                                                                                                  aaaCAAAGCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTACGGCCAGACCCTGTCTTATCGCGAATTAGATGAGGAAGCGAACCGCATTGCACGCCGGCTGc                                  >  7:369758/1‑135 (MQ=255)
                                                                                                         gCCTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCACACCCTGTCATATCGCGAATTAGATGAGGAAGCGAACCGCATTGCCCGCCGGCTGCAAAACAACGGt                       >  2:337329/1‑139 (MQ=255)
                                                                                                          ccTCTGACGTCTTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCAGACCCTGTCTTATCGCGAATTAGATGAGGAAGCGAACCGCATTGCACGCCGGCTGCAAAAACACGGTg                      <  1:102471/139‑1 (MQ=255)
                                                                                                          ccTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCAGCCCCTGTCTTATCGCGAATTTGATGAGGAAGCGAACCGCATTGCACGCCGGCTGCAAAAACCCGGTg                      >  6:349814/1‑139 (MQ=255)
                                                                                                          ccTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCAGACCCTGTCTTATCGCGAATTAGATGAGGAAGCGAACCGCATTGCACGCCGGCTGCAAAAACACGGTg                      <  2:122675/139‑1 (MQ=255)
                                                                                                          ccTCTGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCAGACCCTGTCTTATCGCGAATTAGATGAGGAAGCGAACCGCATTGCACGCCGGCTGCAAAAACACGGTg                      <  5:74758/139‑1 (MQ=255)
                                                                                                             ctGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCAGCCCCTGTCTTATCGCGAATTAGATGAGGAAGCGAACCGCATTGCACGCCGGCTGCAAAAACACGGTGCgg                   >  3:251601/1‑139 (MQ=255)
                                                                                                             ctGACGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCAGACCCTGTCTTATCGCGAATTAGATGAGGAAGCGAACCGCATTGCACGCCGGCTGCAAAAACACGGTGCgg                   >  1:61256/1‑139 (MQ=255)
                                                                                                                 cGTATTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCAGACCCTGTCTTATCGCGAATTAGATGAGGAAGCGAACCGCATTGCACGCCGGCTGCAAAAACACGGTGCggg                  <  8:18490/136‑1 (MQ=255)
                                                                                                                    aTTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCAGACCCTGTCTTATCGCGAATTAGATGAGGAAGCGAACCGCATTGCACGCCGGCTGCAAAAACACGGTGCGGGCAAAgg            >  1:6964/1‑139 (MQ=255)
                                                                                                                    aTTGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCGTTAGTATTACGGCCCGCCCCGGGATTTTAGCGAATTAGATGAGGAATCGAACCGCCTTTCACGCCCGCTGCAAAAACACGGGGCCGGCAAAgg            >  2:192959/1‑139 (MQ=255)
                                                                                                                     ttGGTTCATGGAAGCGGTGAACGCCAATCCGGATGCAGCGGCGCTTACGTATTCCGGCCAGACCCTGTCTTATGGCGAATTAGATGAGGAAGCGAACCGCATTGCACGCCGGCTGCAAAAACACGGTGCGGGCAAAgg            <  3:68555/138‑1 (MQ=255)
                                                                                                                     ttGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCAGACCCTGTCTTATCGCGAATTAGATGAGGAAGCGAACCGCATTGCACGCCGGCTGCAAAAACACGGTGCGGGCAAAGGc           >  1:114840/1‑139 (MQ=255)
                                                                                                                     ttGGTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCAGACCCTGTCTTATCGCGAATTAGATGAGGAAGCGAACCGCATTGCACGCCGGCTGCAAAAACACGGTGCGGGCAAAGGc           >  4:115072/1‑139 (MQ=255)
                                                                                                                       ggTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCGTACGTATTCCGGCCAGACCCTGTCTTAACGCGAATTAGATGAGGAAGCGAACCGCAGTGCACGCCGGCTGCAAAAACACGGTGCGGGCAAAGGctc         <  1:306422/139‑1 (MQ=255)
                                                                                                                        ggtCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCAGACCCTGTCTTATCGCGAATTAGATGAGGAAGCGAACCGCATTGCACGCCGGCTGCAAAAACACGGTGCGGGCAAAgg            <  8:224555/133‑1 (MQ=255)
                                                                                                                        gTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCg                                                                                              <  7:275735/53‑1 (MQ=255)
                                                                                                                        gTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCg                                                                                              >  8:275735/1‑53 (MQ=255)
                                                                                                                        gTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCAGACCCTGTCTTATCGCGAATTAGATGAGGAAGCGAACCGCATTGCACGCCGGCTGCAAAAACACGGTGCGGGCAAAGGctc         <  4:180095/138‑1 (MQ=255)
                                                                                                                        gTTCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCAGACCCTGTCTTATCGCGAATTAGATGAGGAAGCGAACCGCATTGCACGCCGGCTGCAAAAACACGGTGCGGGCAAAGGctc         <  2:61256/138‑1 (MQ=255)
                                                                                                                        gTTCAAGGAAGCGGTGAACGCCAATACGGATGCACCGGCGCTTACGTATTCCGGCCAGACCCTGTCTTATCGCGAATTAGATGAGGAAGCGAACCGCATTGCACGCCGGCTGCAAAAACACGGTGCGGGCAAAgg            <  8:284188/135‑1 (MQ=255)
                                                                                                                          tCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCAGACCCTGTCTTATCGCGAATTAGATGAGGAAGCGAACCGCATTGCACGCCGGCTGCAAAAACACGGTGCGGGCAAAGGCTCtgt      >  1:215438/1‑139 (MQ=255)
                                                                                                                          tCAAGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCGTACGTAATCCGGCCCGCCCCTGGCCTTTAGCGAAATTGAAGTGGAAGCGAACCCCAATGCACGCCGGGGGCAAAAACACCGTGCGGGCAAAGGAGCtgt      >  4:68555/1‑139 (MQ=255)
                                                                                                                             aGGAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTACGGCCAGACCCGGTCTTATAGCGAATTAGATGAGGAAGCGCACCGCCATGCACGCCGGCTGCAAAAACACGGTGCGGGCAAAGGCTCggttg    >  6:304175/1‑138 (MQ=255)
                                                                                                                              ggAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCAGACCCTGTCTTATCGCGAATTAGATGAGGAAGCGAACCGCATTGCACGCCGGCTGCAAAAACACGGTGCGGGCAAAGGCTCtgtt     <  7:253724/136‑1 (MQ=255)
                                                                                                                              ggAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCAGACCCTGTCTTATCGCGAATTAGATGAGGAAGCGAACCGCATTGCACGCCGGCTGCAAAAACACGGTGCGGGCAAAGGCTCTGTTGTa  <  3:191316/139‑1 (MQ=255)
                                                                                                                              ggAAGCGGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCAGACCCTGTCTTATCGCGAATTAGATGAGGAAGCGAACCGCATTGCACGCCGGCTGCAAAAACACGGTGCGGGCAAAGGCTCTGTTGTa  <  2:215438/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                    
CCAGAATCCTGATCAGCCTGTCAGCACGATCAACCTCGTTGACGACAGGGAGAGAGAATTTTTGCTAACCGGCTTAAACCCGCCGGCTCAAGCTCATGAAACAAAGCCTCTGACGTATTGGTTCAAGGAAGCAGTGAACGCCAATCCGGATGCACCGGCGCTTACGTATTCCGGCCAGACCCTGTCTTATCGCGAATTAGATGAGGAAGCGAACCGCATTGCACGCCGGCTGCAAAAACACGGTGCGGGCAAAGGCTCTGTTGTA  >  BS168_reference_genome/399818‑400082

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