Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 390,006 A→C 19.2% T755P (ACA→CCA)  srfAB → surfactin synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome390,0060AC19.2% 89.5 / 7.0 60T755P (ACA→CCA) srfABsurfactin synthetase
Reads supporting (aligned to +/- strand):  ref base A (17/31);  new base C (10/2);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 3.76e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTGCTGCAAAGGATGCGGGTATGTTCGGCACACTGCGCCATTTAATCATCGGCGGAGATGCCCTTGTCCCGCATATTGTCAGCAAAGTAAAACAGGCATCGCCGTCGCTTTCATTGTGGAACGGATACGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAGCAGCAATGCCTGCAGCCAATCGGCGCGCCTG  >  BS168_reference_genome/389871‑390144
                                                                                                                                        |                                                                                                                                           
cTTGCTGCAAAGGATGCGGGTATGTTCGGCACACTGCGCCATTTAATCATCGGCGGAGATGCCCTTGTCCCGCATATTGTCAGCAAAGTAAAACAGGCATCGCCGTCGCTTTCATTGTGGAACGGATACGGCCCGa                                                                                                                                             >  2:220213/1‑136 (MQ=255)
               gCGGGTATGTTCGGCACACTGCGCCATTTAATCATCGGCGGAGATGCCCTTGTCCCGCATATTGTCAGCAAAGTAAAACAGGCATCGCCGTCGCTTTCATTGTGGCACGGCTACGGGCCGCCAGCAAACACGACGtttt                                                                                                                          >  3:312254/1‑139 (MQ=255)
                   gTATGTTGGTCACACTGCGCCATTTAATCATCGGCGGAGATGCCCTTGGCCCGATTATTGTCAGAAAAGTAAAACAGGCATCGCCGTCGCTTTCTTTGGGGACCGTATACGGCCCGACAGAAAACACGACGTTTTCAAc                                                                                                                      <  5:285588/139‑1 (MQ=255)
                      tGTTCGGCACACTGCGCCATTTAATCATCGGCGGAGATGCCCTTGTCCCGCATATTGTCAGCAAAGTAAAACAGGCATCGCCGTCGCGTTAATTGTGGAACGGATAAGGCCCCCCCGAAAACACGCCGGTTTGAACCAg                                                                                                                   >  4:316105/1‑139 (MQ=255)
                             cacaCTGCGCCATTTAATCATCGGCGGAGATGCCCTTGTCCCGCATATTGTCAGCAAAGTAAAACAGGCATCGCCGTCGCTTTGATTGTGGTACGGTTAAGGCCCCACAGAAAACACCCCCTTTTAAACCAGTTTTATg                                                                                                            >  3:13404/1‑139 (MQ=255)
                             cacaCTGCGCCATTTAATCATCGGCGGAGATGCCCTTGTCCCGCATATTGTCAGCAAAGTAAAACAGGCATCGCCGTCGCTTTCATTGTGGAACGGATACGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCt                                                                                                             >  7:268491/1‑138 (MQ=255)
                             cacaCTGCGCCATTTAATCATCGGCGGAGATGCCCTTGTCCCGCATATTGTCAGCAAAGTAAAACAGGCATCGCCGTCGCGTTCATTGTGGTACGGGTACGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTg                                                                                                            >  3:67270/1‑139 (MQ=255)
                             cacaCTGCGCCATTTAATCATCGGCGGAGATGCCCTTGTCACGCATATTGTCAGCAAAGTAAAACAGGCATCGCCGTCGCTTTAATTGTGGAACGGATACGGCCCGACAGGAAAAACGAAGTGTTGCACCAGttttat                                                                                                             >  4:348307/1‑136 (MQ=255)
                                  tGCGCCATTTAATCATCGGCGGAGATGCCCTTGTCCCGCATATTGTCAGCAAAGTAAAACAGGCATCGCCGTCGCTTTCATTGTGGAACGGATACGGCCCGACAGAAAACACCACCGTTTTCACCAGTTTTGTGATCGa                                                                                                       >  2:219423/1‑139 (MQ=255)
                                  tGCGCCATTTAATCATCGGCGGAGATGCCCTTGTCCCGCATATTGTCAGCAAAGTAAAACAGGCATCGCCGTCGCGTTGATTTTGGGACCGGTACGGCCCGCCCGCAAACACGCCGGTTTGACCCAGGTTTGTGGTCGa                                                                                                       >  3:269367/1‑139 (MQ=255)
                                  tGCGCCATTTAATCATCGGCGGAGATGCCCTTGTCCCGCATATTGTCAGCAAAGTAAAACAGGCATCGCCGTCGCGTTGATTGTGGAACGGATACGGCCCGCCACAAAACCCGACGTTTTCAACCAGTTTTCTGATCGa                                                                                                       >  4:150932/1‑139 (MQ=255)
                                  tGCGCCATTTAATCATCGGCGGAGATGCCCTTGTCCCGCATATTGTCAGCAAAGTAAAACAGGCATCGCCGTCGCGTTCATTGTGGGAAGGATACGGCCCGCCCGCAAACACGACCTTTTGAACCCGTTTTCTGATCgg                                                                                                       >  7:235014/1‑138 (MQ=255)
                                     gcaaTTTAACCTGCGGCGGAGATGCCCTTGTCCCGAATATTGTCAGCAAAGAAAACCCGGCATCGCCGTCGCTTTCATTGTGTACCGGATACGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACcg                                                                                                    <  4:63836/136‑1 (MQ=255)
                                         tttAATCATCGGCGGATATGCCCTTGCCCCGAATTTTGTCACCAAAGAAAAACCGGCATCGCCGTCGCTTTCATTGTGGAACGGAAACGGCCCGACACAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAg                                                                                                <  4:67270/139‑1 (MQ=255)
                                              tcatcGGCGGAGATGCCCTTGTCCCGCATATTGTCAGCAAAGTAAAACCGCCATCGCCGGCGCTTTCATTGTGGAACGGATACGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATgg                                                                                           <  7:96256/139‑1 (MQ=255)
                                              tcatcGGCGGAGATGCCCTTGTCCCGCATATTGTCAGCAAAGTAAAACAGGCATCGCCGTCGCTTTCATTGTGGAACGGATCCGGCCCGCCAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATgg                                                                                           <  3:132334/139‑1 (MQ=255)
                                                 tcGGCGGAGATGCCCTTGTCCCGCATATTGTCAGCAAAGTAAAACAGGCATCGCCGTCGCTTTCATTGTGGAACGGATACGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATgg                                                                                           <  7:227035/136‑1 (MQ=255)
                                                 tcGGCGGAGATGCCCTTGTCCCGCATATTGTCAGCAAAGTAAAACAGGCATCGCCGTCGCTTTCATTGTGGAACGGATACGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATgg                                                                                           >  8:227035/1‑136 (MQ=255)
                                                   ggcggAGATGCCCTTGTCCCGCATATTGTCAGCAAAGTAAAACAGGCATCGCCGTCGCTTTCATTGTTGAACGGATACGGCCCGACAGAAAACACGACGTTTTGAACCAGTTTTCTGATCGACCGGGCGTATggcggcg                                                                                      >  8:213979/1‑138 (MQ=255)
                                                   ggcggAGATGCCCTTGTCCCGCATATTGTCAGCAAAGTAAAACAGGCATCGCCGTCGCGTTGATTGTGGGACGGATACGGCCCGCCAGAAAACACCACGGTTTAAACCAGTTTTGTGATCGACCGCGAGTATggcggcg                                                                                      >  4:408824/1‑138 (MQ=255)
                                                   ggcggAGATGCCCTTGTCCCGAATATTGTCAGCAAAGTAAAACCGGCAGCGCCGTCGCTTTCATTGGGGACCGGATACGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGct                                                                                      <  6:241313/139‑1 (MQ=255)
                                                   cagacagagGCCCTTGTCCCGCTTATTGGCAGAAAATTAAAACAGGCAGCGCGGTCGCTTTCTTTGTGCACCGGATACGGCCCGACATAAACCACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGct                                                                                      <  3:135887/134‑1 (MQ=255)
                                                       acaccgCCCTGGCCCCGCATATTGTCAGCAAAGTAAAACAGGCATCGCCGTCGCTTTAATTGGGGACCGTATACGGCCCGACATAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTAt                                                                                  <  4:368247/134‑1 (MQ=255)
                                                          aTGCCCTTGTCCCGCATATTGTCAGCAAAGTAAAACAGGCATCGCCGTCGCTTTGATTGTGGAACGGATACGGCCCGACAGAAAACACGACCTTTTCAACCAGTTTTATGATCGACCCCCAGTATGGCGGCTCTATacc                                                                               >  2:100384/1‑139 (MQ=255)
                                                          aTGCCCTTGTCCCGCATATTGTCAGCAAAGTAAAACAGGCATCGCCGTCGCTTTAATTTTGGGAACGTTTCGGCCCCCCAGCAAAAAAGACCTTTTAAACCATGTTTGAGTTAGCCCCCGAGTATTGCCGGTCTaaacc                                                                               >  6:28295/1‑135 (MQ=255)
                                                               cTTGTCCCGCATATTGTCAGCAAAGTAAAACAGGCATCGCCGTCGCGTTAAATGTGGGACGGATACGGCCCCCCACAAAAACCGACGGTTTCAACCACTTTTCTGATAGACCGCGCCGATTGCCGCGCGTTCCCAATCg                                                                          >  1:306557/1‑139 (MQ=255)
                                                                        cATATTGTCAGCAAAGTAAAACAGGCATCGCCGTCGCTTTCATTTTGGAACGGATACGGCCCGACAGAAAACACGCCCTTTTGAACCAGTTTTCTGATCGACCGCCAGTATGGCGGCTCTATACCCATCCGGactgtgg                                                                 >  3:368247/1‑133 (MQ=255)
                                                                        cATATTGTCAGCAAAGTAAAACAGGCATCGCCGTCGCGGTCATTGTGGGTCGGCTACGGCCCGCCACAGAACACGCCGGTTTCAACCAGTTTTATGATAGCCCGCGAGTATGGCGGCGCTATCCCAAACGGGGCGCCgc                                                                 >  7:168170/1‑138 (MQ=255)
                                                                          aatTGTCAGCAACTTAAACCAGCCATCGCCGCCGCTTTCATTGTGGAACGGATACGGCCCCACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATc                                                               <  1:219423/138‑1 (MQ=255)
                                                                             tGTCAGCAAAGTAAAACAGGCATCGCCGTCGCTTTAATTGTGGGAAGGATACGGCCCGACAGAAAACACGCCGTTTTCAACCACTTTTATGATCGACCGCGAGTTTGGCGGCGCTATCCCAAACGGGAAGCCGATCgg                                                             >  4:185928/1‑138 (MQ=255)
                                                                                aatcAAAGTAAAACAGGATTGGCCGTCGCTTGATTTGTGGAACGTATACGGCCCGACATAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGaaa                                                          <  3:348307/135‑1 (MQ=255)
                                                                                   aaacgTAAAACCGCCAGCCCCGGCGCTTTCATTGGGCACCGGATCCGCCCCCACACAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTcc                                                      <  4:269367/135‑1 (MQ=255)
                                                                                    gaaGTAAAACAGGCATCGCCGTCGCTTTCATTGTGGAACGGATACGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCa                                                     <  3:209236/138‑1 (MQ=255)
                                                                                    aacgTAAACCAGGCATCGCCGTCGCTTTCATTGGGGAACGGATACGGCCCGACACAAACCACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCa                                                     <  4:312254/136‑1 (MQ=255)
                                                                                    aaaTTAAAACCGGCAGCGCGGGCGTTTTCATTGTGCAACGTAACGGCCCCGACACAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCa                                                     <  4:13404/139‑1 (MQ=255)
                                                                                    aaaGTAAAACAGGCTTCGCCGACGATTTCTTTTTGGACGGTAGCCGCCCCCACACAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCa                                                     <  5:28295/139‑1 (MQ=255)
                                                                                    aaaGTAAAACAGGCATCGCCGTCGCTTTCATTGTGGAACGGATACGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCa                                                     <  7:213979/139‑1 (MQ=255)
                                                                                    aaaGAAAAACAGGCCTCGCCGTCGCTTTCATTGTGGACCGTTTACGGCCCCCCAGAAAACACGACGTTTTCAACTAGTTTTCTGATCGACTGCGAGTATGGCGTCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCa                                                     <  3:295049/139‑1 (MQ=255)
                                                                                     aaGTAAAACAGGCATCGCCGTCGCTTTCATTGTGGAACGGATACGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCAc                                                    <  1:220213/139‑1 (MQ=255)
                                                                                      aGTAAAACAGGCATCGCCGTCGCGTTGATTTTGGGACGGATACGGCCCGCCCGAAAACACCCCGGTTTCAACCAGTTTTGTGGTCGACCGCGCGGATGGCGGGGCTATCCCAATCGGGGAGCCGCTCGGGAAttgcccg                                                   >  4:189050/1‑134 (MQ=255)
                                                                                        tAAACCCGGCATCGCCGTCGCTTTCATTGGGCAACGTACAGGGCCCGACAGAAACCACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGc                                                 <  2:228613/139‑1 (MQ=255)
                                                                                         aaacccggCATCGCCGTCGCTTTCATTGTGGACCGGATACGGCCCGACATAAACCACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGc                                                 <  3:316105/132‑1 (MQ=255)
                                                                                         aaaaCAGGCATCGCCGGCGCTTTCATTGTGGAACGGATACGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGc                                                 <  3:150932/138‑1 (MQ=255)
                                                                                                gCATCGCCGTCGCTTTGATTTTGGAACGGATACGGCCCGACAGAAAACACGACCTTTTAAACCAGTTTTCTGATAGACCGCGAGTATTGCGGCTCTAACCCAAACGGGAAGCCGATCGGAAATTACACTGCCAAtatta                                         >  4:209236/1‑139 (MQ=255)
                                                                                                         gcGCTTTCATTGTGGAACGGATACGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGagc                                <  8:260966/138‑1 (MQ=255)
                                                                                                              agagttGTGTAATGTATACGGCCGGAAACAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCTAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAgcagca                            <  4:321217/134‑1 (MQ=255)
                                                                                                               tCATTGTGCAACGGATACGGCCCCACAGAAAACACGACGTTTTAAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAgcagca                            <  7:132693/137‑1 (MQ=255)
                                                                                                                 aTTGTGGAACGGATACGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGCTCGGAAATTACACTGCCTATATTATGGATGAGCAGCAATGc                        >  7:29529/1‑139 (MQ=255)
                                                                                                                          cGGATACGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAGCAGCAATGCCTGCAGcc                >  3:82652/1‑138 (MQ=255)
                                                                                                                             aTACGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGAATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAATTTCAACTGCCTATATAATGGATGAGCAGCAATGCCTGCAGCCAATc            >  3:35655/1‑139 (MQ=255)
                                                                                                                             aTACGGCCCGACAGAAAACACGACGTTTTCAACCAGTATTCTGATCGACAGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAGCAGCAATGCCTGCAGcctgtc            >  3:321217/1‑135 (MQ=255)
                                                                                                                               aCGGCCCGACAGAAACCACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAGCAGCAATGCCTGCAGCCAATCgg          <  1:100384/139‑1 (MQ=255)
                                                                                                                               aCGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAGCAGCaa                           >  4:156392/1‑122 (MQ=255)
                                                                                                                               aCGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAGCAGCaa                           <  3:156392/122‑1 (MQ=255)
                                                                                                                               aCGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAGCAGCAATGCCTGCAGCCAATCgg          <  7:312760/139‑1 (MQ=255)
                                                                                                                                  gCCCTACATAAAACAGGACGTTTTCAACGAGTTTTGTGATCGACCGCGAGTATGGCGGCTCTATCCCACTCGGGAAGCCGATCGGAAAATCCACTGCCTTTATTATGGATGAGCAGCAATGGCTGCAGCCAATCGgcgc       <  3:364678/139‑1 (MQ=255)
                                                                                                                                  gCCCGACACAAACCACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAGCAGCAATGCCTGCAGCCAATCGgcgc       <  3:185928/139‑1 (MQ=255)
                                                                                                                                  gCCCGACACAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAGCAGCAATGGCTGCAGCCAATCGgcgc       <  3:370877/139‑1 (MQ=255)
                                                                                                                                   cccccACAGAAACCACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAGCAGCAATGCCTGCAGCCAATCGgcgc       <  5:296940/138‑1 (MQ=255)
                                                                                                                                       gACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTcc                                                      <  2:274588/88‑1 (MQ=255)
                                                                                                                                       gACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTcc                                                      >  1:274588/1‑88 (MQ=255)
                                                                                                                                        aCAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAGCAGCAATGCCTGCCGCCAATCGGCGCGCCTg  >  8:187424/1‑139 (MQ=255)
                                                                                                                                        aCAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAGCAGCAATGCCTGCAGCCAATCGGCGCGCCTg  >  6:20439/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                           
CTTGCTGCAAAGGATGCGGGTATGTTCGGCACACTGCGCCATTTAATCATCGGCGGAGATGCCCTTGTCCCGCATATTGTCAGCAAAGTAAAACAGGCATCGCCGTCGCTTTCATTGTGGAACGGATACGGCCCGACAGAAAACACGACGTTTTCAACCAGTTTTCTGATCGACCGCGAGTATGGCGGCTCTATCCCAATCGGGAAGCCGATCGGAAATTCCACTGCCTATATTATGGATGAGCAGCAATGCCTGCAGCCAATCGGCGCGCCTG  >  BS168_reference_genome/389871‑390144

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