Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,567,469 A→T 8.8% intergenic (‑329/+35) dosC ← / ← yddW diguanylate cyclase, cold‑ and stationary phase‑induced oxygen‑dependent biofilm regulator/liprotein, glycosyl hydrolase homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,567,4690AT 8.8% 188.6 / 7.5 57intergenic (‑329/+35)dosC/yddWdiguanylate cyclase, cold‑ and stationary phase‑induced oxygen‑dependent biofilm regulator/liprotein, glycosyl hydrolase homolog
Reads supporting (aligned to +/- strand):  ref base A (27/25);  new base T (3/2);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

GTGACAGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTAA  >  NC_000913/1567328‑1567616
                                                                                                                                             |                                                                                                                                                   
gTGACAGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCAcc                                                                                                                                              >  1:725463/1‑149 (MQ=255)
   aCAGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGc                                                                                                                                           <  3:464952/149‑1 (MQ=255)
   aCAGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGc                                                                                                                                           >  1:444440/1‑149 (MQ=255)
     aGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCtt                                                                                                                                         <  1:443529/149‑1 (MQ=255)
     aGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGCTTTACCCACCGGCtt                                                                                                                                         <  4:476748/149‑1 (MQ=255)
       cGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGt                                                                                                                                       <  2:157140/149‑1 (MQ=255)
             tCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACaaa                                                                                                                                 <  1:274332/149‑1 (MQ=255)
               gTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGt                                                                                                                               >  3:527341/1‑149 (MQ=255)
                tATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTaaa                                                                                                                                                   <  1:325708/128‑1 (MQ=255)
                tATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTaaa                                                                                                                                                   >  2:325708/1‑128 (MQ=255)
                tATGGTTTTTCATATTCATACATTTTTATGAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGACGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTa                                                                                                                              >  2:598016/1‑149 (MQ=255)
                  tGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAg                                                                                                                            >  2:481326/1‑149 (MQ=255)
                   ggTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGa                                                                                                                           <  1:752848/149‑1 (MQ=255)
                     tttttCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAAt                                                                                                                         >  4:433558/1‑149 (MQ=255)
                                 ataCATTTGTATTAGGGGTTTATGGCGGATTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACt                                                                                                             >  4:114719/1‑149 (MQ=255)
                                  taCATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTg                                                                                                            >  3:631305/1‑149 (MQ=255)
                                    cATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGt                                                                                                                                       >  2:387783/1‑120 (MQ=255)
                                    cATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGt                                                                                                                                       <  1:387783/120‑1 (MQ=255)
                                     aTTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGccc                                                                                                         <  1:654277/149‑1 (MQ=255)
                                     aTTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGccc                                                                                                         >  3:695086/1‑149 (MQ=255)
                                           attaGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGa                                                                                                   <  1:481326/149‑1 (MQ=255)
                                               gggATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCAGCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACttt                                                                                               <  4:723103/149‑1 (MQ=255)
                                                  aTTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCa                                                                                            <  3:499445/149‑1 (MQ=255)
                                                           tGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGa                                                                                   >  2:142816/1‑149 (MQ=255)
                                                           tGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGCTTTACCCACCGGCTTGTTAc                                                                                                                                    >  4:498070/1‑100 (MQ=37)
                                                           tGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGCTTTACCCACCGGCTTGTTAc                                                                                                                                    <  3:498070/100‑1 (MQ=37)
                                                             tttAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGAcc                                                                                 >  4:358499/1‑149 (MQ=255)
                                                                      gtgtGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGa                                                                        >  1:258557/1‑149 (MQ=255)
                                                                      gtgtGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGa                                                                        <  3:114719/149‑1 (MQ=255)
                                                                      gtgtGGTTAATTTGACCTCAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGa                                                                        <  1:598016/149‑1 (MQ=255)
                                                                      atctgGTTAATTTGGATTTAGTAAACATGATTATTAGGTGGATAGTTTAAGGGGGTTACAAACCGGTTGGGAAAGCCCCGGCTTGTTAAAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGa                                                                        <  1:28574/146‑1 (MQ=255)
                                                                           gTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCtt                                                                                                                                         <  4:463172/79‑1 (MQ=255)
                                                                           gTTAATTTGACTTAAGTAAGCATGATTAGTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCtt                                                                                                                                         >  3:463172/1‑79 (MQ=255)
                                                                            ttAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGc                                                                  >  3:690046/1‑149 (MQ=255)
                                                                             tAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCg                                                                 <  1:534810/149‑1 (MQ=255)
                                                                             tAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCg                                                                 <  3:120732/149‑1 (MQ=255)
                                                                                tttGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGtt                                                              >  3:460216/1‑149 (MQ=255)
                                                                                         agtaagCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGAt                                                     <  4:105194/149‑1 (MQ=255)
                                                                                              gCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTc                                                         >  3:438983/1‑140 (MQ=255)
                                                                                              gCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTcc                                                >  3:14001/1‑149 (MQ=255)
                                                                                              gCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACGTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTc                                                         <  4:438983/140‑1 (MQ=255)
                                                                                                aTGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTc                                              >  3:211179/1‑149 (MQ=255)
                                                                                                      attaGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCTAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGaa                                        >  3:716656/1‑149 (MQ=255)
                                                                                                           tGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGa                                   >  3:499143/1‑149 (MQ=255)
                                                                                                               aTAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGt                               >  4:701659/1‑149 (MQ=255)
                                                                                                                   tttAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCg                           >  1:387535/1‑149 (MQ=255)
                                                                                                                       agagGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTaa                       >  3:520978/1‑149 (MQ=255)
                                                                                                                        gagGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAAt                      <  3:165606/149‑1 (MQ=255)
                                                                                                                             tAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCgg                 <  1:142816/149‑1 (MQ=255)
                                                                                                                             tAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCgg                 <  2:258557/149‑1 (MQ=255)
                                                                                                                             tAACAAGCCGGTGCTTTACCCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGAATTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCgg                 >  4:410556/1‑149 (MQ=255)
                                                                                                                              aaCAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCggg                <  4:520978/149‑1 (MQ=255)
                                                                                                                              aaCAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTAACGGAACAAGATGGTGCCGCTAATTTCggg                >  1:426706/1‑149 (MQ=255)
                                                                                                                                 aaGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGcac             <  4:460216/149‑1 (MQ=255)
                                                                                                                                 aaGCCGGTGCTTTACCCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGcac             >  3:418472/1‑149 (MQ=255)
                                                                                                                                  aGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGcaca            <  4:14001/149‑1 (MQ=255)
                                                                                                                                         gggTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAg                                                                       >  3:557678/1‑83 (MQ=255)
                                                                                                                                         gggTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAg                                                                       <  4:557678/83‑1 (MQ=255)
                                                                                                                                            tAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTaa  >  3:144332/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
GTGACAGCGAATATCGTATGGTTTTTCATATTCATACATTTTTATTAGGGATTTATGGCTGTTTAACTAAGTGTGGTTAATTTGACTTAAGTAAGCATGATTATTAGTGGGATAGTTTAAGAGGGTAACAAGCCGGTGGGTAAAGCACCGGCTTGTTACAAAGTAAGAATGGGAGTTTAACTGCCCCAGCGACTTTGCAGATAGCTGACCGCTTGTTGAGTCTGCGGTTTATTCAGATAGTCCTCACGGAACAAGATGGTGCCGCTAATTTCGGGCACAGCATCGTTAA  >  NC_000913/1567328‑1567616

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