Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 712,946 C→T intergenic (‑144/‑41) ybfF ← / → seqA acyl‑CoA esterase/negative modulator of initiation of replication

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913712,9460CT98.3% 217.5 / ‑4.1 59intergenic (‑144/‑41)ybfF/seqAacyl‑CoA esterase/negative modulator of initiation of replication
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (35/22);  total (35/24)
Fisher's exact test for biased strand distribution p-value = 3.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

CATGGCAAAGTTCTTTTTTTCGCGTTGTCGGGTTAGGGTATTATGTTGACCATTGTGCCACAGGGCTGCAACAAATAAGGTTTATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCCTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTATATTGCCAGCCACACTAAGCATATCGGCGAGAGCGCATCCGACATTTTACGGCGTATG  >  NC_000913/712800‑713082
                                                                                                                                                  |                                                                                                                                        
cATGGCAAAGTTCTTTTTTTCGCGTTGTCGGGTTAGGGTATTATGTTGACCATTGTGCCACAGGGCTGCAACAAATAAGGTTTATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGc                                                                                                                                      <  1:867826/151‑1 (MQ=255)
cATGGCAAAGTTCTTTTTTTCGCGTTGTCGGGTTAGGGTATTATGTTGACCATTGTGCCACAGGGCTGCAACAAATAAGGTTTATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGc                                                                                                                                      >  1:431173/1‑151 (MQ=255)
        aGTTCTTTTTTTCGCGTTGTCGGGTTAGGGTATTATGTTGACCATTGTGCCACAGGGCTGCAACAAATAAGGTTTATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGt                                                                                                                               >  2:26292/1‑150 (MQ=255)
        aGTTCTTTTTTTCGCGTTGTCGGGTTAGGGTATTATGTTGACCATTGTGCCACAGGGCTGCAACAAATAAGGTTTATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGt                                                                                                                               >  1:601164/1‑150 (MQ=255)
                ttttCGCGTTGTCGGGTTAGGGTATTATGTTGACCATTGTGCCACAGGGCTGCAACAAATAAGGTTTATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCt                                                                                                                       >  2:869508/1‑150 (MQ=255)
                 tttCGCGTTGTCGGGTTAGGGTATTATGTTGACCATTGTGCCACAGGGCTGCAACAAATAAGGTTTATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTa                                                                                                                      >  1:57004/1‑150 (MQ=255)
                  ttCGCGTTGTCGGGTTAGGGTATTATGTTGACCATTGTGCCACAGGGCTGCAACAAATAAGGTTTATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTa                                                                                                                      >  1:187029/1‑149 (MQ=255)
                  ttCGCGTTGTCGGGTTAGGGTATTATGTTGACCATTGTGCCACAGGGCTGCAACAAATAAGGTTTATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTa                                                                                                                      >  1:119483/1‑149 (MQ=255)
                      cgTTGTCGGGTTAGGGTATTATGTTGACCATTGTGCCACAGGGCTGCAACAAATAAGGTTTATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGacac                                                                                                                >  2:858596/1‑151 (MQ=255)
                            cGGGTTAGGGTATTATGTTGACCATTGTGCCACAGGGCTGCAACAAATAAGGATTATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGAcactgcact                                                                                                          <  2:208156/151‑1 (MQ=255)
                             gggTTAGGGTATTATGTTGACCATTGTGCCACAGGGCTGCAACAAATAAGGTTTATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGAcactgcactg                                                                                                         >  1:526655/1‑151 (MQ=255)
                                 tAGGGTATTATGTTGACCATTGTGCCACAGGGCTGCAAAAAATAAGGTTTATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGAtt                                                                                                     >  1:416560/1‑151 (MQ=255)
                                  aGGGTATTATGTTGACCATTGTGCCACAGGGCTGCAACAAATAAGGTTTATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGAtt                                                                                                     >  1:425387/1‑150 (MQ=255)
                                  aGGGTATTATGTTGACCATTGTGCCACAGGGCTGCAACAAATAAGGTTTATTCCGAGTTTCTCTGCAAGCCAGGCTTGTCGCTATCCGCTGCCGGGATTTACTCAGATACTCGTGGCGACTTGTATTCAGCTAAGACACTGCACTgggt                                                                                                      <  2:233323/149‑3 (MQ=255)
                                       attatGTTGACCATTGTGCCACAGGGCTGCAACAAATAAGGTTTATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATg                                                                                               >  2:80457/1‑151 (MQ=255)
                                       attatGTTGACCATTGTGCCACAGGGCTGCAACAAATAAGGTTTATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATg                                                                                               >  2:407562/1‑151 (MQ=255)
                                       attatGTTGACCATTGTGCCACAGGGCTGCAACAAATAAGGTTTATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATg                                                                                               >  1:128268/1‑151 (MQ=255)
                                         tatGTTGACCATTGTGCCACAGGGCTGCAACAAATAAGGTTTATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGa                                                                                              >  1:322321/1‑150 (MQ=255)
                                          atGTTGACCATTGTGCCACAGGGCTGCAACAAATAAGGTTTATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGaaa                                                                                            <  1:506660/151‑1 (MQ=255)
                                               gACCATTGTGCCACAGGGCTGCAACAAATAAGGTTTATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGa                                                                                        <  2:652744/150‑1 (MQ=255)
                                                       tgCCACAGGGCTGCAACAAATAAGGTTTATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTtg                                                                               >  2:138655/1‑151 (MQ=255)
                                                         ccACAGGGCTGCAACAAATAAGGTTTATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCCTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTtgat                                                                             <  1:815688/151‑1 (MQ=255)
                                                                   gCAACAAATAAGGTTTATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAActct                                                                   >  1:572517/1‑151 (MQ=255)
                                                                    caacaaATAAGGTTTATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTa                                                                  <  1:13979/151‑1 (MQ=255)
                                                                          aTAAGGTTTATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCta                                                            <  1:191089/151‑1 (MQ=255)
                                                                            aaGGTTTATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCtat                                                           <  2:614877/150‑1 (MQ=255)
                                                                            aaGGTTTATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCtat                                                           >  2:902610/1‑150 (MQ=255)
                                                                             aGGTTTATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTatat                                                         <  2:225722/151‑1 (MQ=255)
                                                                               gcttATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTATATTg                                                       >  2:194092/3‑151 (MQ=255)
                                                                                tttATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTATATTgc                                                      >  2:71164/1‑151 (MQ=255)
                                                                                tttATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTATATTgc                                                      >  2:591909/1‑151 (MQ=255)
                                                                                 ttATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTATATTgcc                                                     >  1:677852/1‑151 (MQ=255)
                                                                                 ttATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTGGATGATGAACTCTACAGCTATATTgc                                                      <  1:754400/150‑1 (MQ=255)
                                                                                     tCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTATATTgccagcc                                                 >  1:464692/1‑151 (MQ=255)
                                                                                      ccGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTATATTgccagcc                                                 >  2:503701/1‑150 (MQ=255)
                                                                                      ccGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTATATTgccagcc                                                 >  2:256250/1‑150 (MQ=255)
                                                                                        gAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTATATTGCCAGCcac                                               >  1:227873/1‑150 (MQ=255)
                                                                                        gAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTATATTGCCAGCcac                                               >  2:857750/1‑150 (MQ=255)
                                                                                          gTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTATATTGCCAGCCACACt                                            <  1:579752/151‑1 (MQ=255)
                                                                                             tttCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTATATTGCCAGCCACACTaa                                          <  2:864808/150‑1 (MQ=255)
                                                                                             tttCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTATATTGCCAGCCACACTaa                                          <  2:772399/150‑1 (MQ=255)
                                                                                               tCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTATATTGCCAGCCACACTAAGc                                        >  1:596684/1‑150 (MQ=255)
                                                                                               tCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTATATTGCCAGCCACACTAAGc                                        >  1:519074/1‑150 (MQ=255)
                                                                                                 tGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTATATTGCCAGCCACACTAAGCat                                      >  1:422165/1‑150 (MQ=255)
                                                                                                  gCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTATATTGCCAGCCACACTAAGCatat                                    <  2:281750/151‑1 (MQ=255)
                                                                                                  gCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTATATTGCCAGCCACACTAAGCat                                      <  2:163254/149‑1 (MQ=255)
                                                                                                      gCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTATATTGCCAGCCACACTAAGCATATCGGc                                <  1:22644/151‑1 (MQ=255)
                                                                                                         aGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTATATTGCCAGCCACACTAAGCATATCGGCg                               >  2:776876/1‑149 (MQ=255)
                                                                                                          ggCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTATATTGCCAGCCACACTAAGCATATCGGCgaga                            <  1:857817/151‑1 (MQ=255)
                                                                                                            cTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTATATTGCCAGCCACACTAAGCATATCGGCGAGAgc                          >  1:309645/1‑151 (MQ=255)
                                                                                                                aCGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTATATTGCCAGCCACACTAAGCATATCGGCGAGAGCGCAt                      >  1:613050/1‑151 (MQ=255)
                                                                                                                 cGCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTATATTGCCAGCCACACTAAGCATATCGGCGAGAGCGCATc                     <  1:455202/151‑1 (MQ=255)
                                                                                                                  gCTATCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTATATTGCCAGCCACACTAAGCATATCGGCGAGAGCGCATcc                    <  1:510381/151‑1 (MQ=255)
                                                                                                                     aTCCGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTATATTGCCAGCCACACTAAGCATATCGGCGAGAGCGCATCCGAc                 <  2:499470/151‑1 (MQ=255)
                                                                                                                       ccGCTGCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTATATTGCCAGCCACACTAAGCATATCGGCGAGAGCGCATCCGACAt               <  2:828493/151‑1 (MQ=255)
                                                                                                                            gCCGGGATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTATATTGCCAGCCACACTAAGCATATCGGCGAGAGCGCATCCGACATTTTAc          >  1:881312/1‑151 (MQ=255)
                                                                                                                                ggATTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTATATTGCCAGCCACACTAAGCATATCGGCGAGAGCGCATCCGACATTTTACGGCg      >  2:908879/1‑151 (MQ=255)
                                                                                                                                  aTTTATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTATATTGCCAGCCACACTAAGCATATCGGCGAGAGCGCATCCGACATTTTACGGCGTa    <  2:309619/151‑1 (MQ=255)
                                                                                                                                     tATTCATATACTCTTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTATATTGCCAGCCACACTAAGCATATCGGCGAGAGCGCATCCGACATTTTACGGCGTATg  <  1:571135/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                        
CATGGCAAAGTTCTTTTTTTCGCGTTGTCGGGTTAGGGTATTATGTTGACCATTGTGCCACAGGGCTGCAACAAATAAGGTTTATTCCGAGTTTTTCTGCAAGCCAGGCTTGACGCTATCCGCTGCCGGGATTTATTCATATACTCCTGGCGACTTGTATTCAGCTAAGACACTGCACTGGATTAAGATGAAAACGATTGAAGTTGATGATGAACTCTACAGCTATATTGCCAGCCACACTAAGCATATCGGCGAGAGCGCATCCGACATTTTACGGCGTATG  >  NC_000913/712800‑713082

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