Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,454,776 G→A 6.9% intergenic (‑174/+31) nicS ← / ← PP_3950 HTH‑type transcriptional repressor NicS/methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,454,7760GA 6.9% 167.9 / 2.8 58intergenic (‑174/+31)nicS/PP_3950HTH‑type transcriptional repressor NicS/methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base G (27/27);  new base A (2/2);  total (29/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 = 6.92e-01

CCAACACTGCCACAGACATTGAGCCATTGCGCGGCTTTTTCGGCTGGGCAAACCCGCTCTCCGTCCGCTACCGGCCTGGAAGTCAGCAGTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGACTTCGAACTGCCCCACCAGCTTCTGCAACTTGCCAGCATCCACCTTCACCTTCTGCGCACTGCTTTGCAGCACCAGGGCGGTC  >  NC_002947/4454631‑4454924
                                                                                                                                                 |                                                                                                                                                    
ccAACACTGCCACAGACATTGAGCCATTGCGCGGCTTTTTCGGCTGGGCAAACCCGCTCTCCGTCCGCTACCGGCCTGGAAGTCAGCAGTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGc                                                                                                                                                   >  8:101456/1‑149 (MQ=255)
  aaCACTGCCACAGACATTGAGCCATTGCGCGGCTTTTTCGGCTGGGCAAACCCGCTCTCCGTCCGCTACCGGCCTGGAAGTCAGCAGTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTc                                                                                                                                                 >  7:121284/1‑149 (MQ=255)
         ccACAGACATTGAGCCATTGCGCGGCTTTTTCGGCTGGGCAAACCCGCTCTCCGTCCGCTACCGGCCTGGAAGTCAGCAGTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGt                                                                                                                                           <  5:289299/148‑1 (MQ=255)
             aGACATTGAGCCATTGCGCGGCTTTTTCGGCTGGGCAAACCCGCTCTCCGTCCGCTACCGGCCTGGAAGTCAGCAGTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCg                                                                                                                                      <  5:8317/149‑1 (MQ=255)
               aCATTGAGCCATTGCGCGGCTTTTTCGGCTGGGCAAACCCGCTCTCCGTCCGCTACCGGCCTGGAAGTCAGCAGTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGa                                                                                                                                    >  7:85179/1‑149 (MQ=255)
                cATTGAGCCATTGCGCGGCTTTTTCGGCTGGGCAAACCCGCTCTCCGTCCGCTACCGGCCTGGAAGTCAGCAGTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGAc                                                                                                                                   <  8:91254/149‑1 (MQ=255)
                 aTTGAGCCATTGCGCGGCTTTTTCGGCTGGGCAAACCCGCTCTCCGTCCGCTACCGGCCTGGAAGTCAGCAGTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGa                                                                                                                                    >  4:262280/1‑147 (MQ=255)
                    gAGCCATTGCGCGGCTTTTTCGGCTGGGCAAACCCGCTCTCCGTCCGCTACCGGCCTGGAAGTCAGCAGTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACAGGt                                                                                                                               >  7:257167/1‑149 (MQ=255)
                         aTTGCGCGGCTTTTTCGGCTGGGCAAACCCGCTCTCCGTCCGCTACCGGCCTGGAAGTCAGCAGTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACAGGCTCCCGGGTTACCGGACAGGTGCGt                                                                                                                           >  4:301387/1‑148 (MQ=255)
                          ttGCGCGGCTTTTTCGGCTGGGCAAACCCGCTCTCCGTCCGCTACCGGCCTGGAAGTCAGCAGTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACAGGTGCGTTa                                                                                                                         <  8:216648/149‑1 (MQ=255)
                                 gCTTTTTCGGCTGGGCAAACCCGCTCTCCGTCCGCTACCGGCCTGGAAGTCAGCAGTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACa                                                                                                                                  <  2:185795/133‑1 (MQ=255)
                                 gCTTTTTCGGCTGGGCAAACCCGCTCTCCGTCCGCTACCGGCCTGGAAGTCAGCAGTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACa                                                                                                                                  >  1:185795/1‑133 (MQ=255)
                                        cGGCTGGGCAAACCCGCTCTCCGTCCGCTACCGGCCTGGAAGTCAGCAGTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGAc                                                                                                            <  8:257167/148‑1 (MQ=255)
                                                cAATCCCGCTCTCCGTCCGCTACCTGCCTTGATGTGAGGAGTCATCGGGGCTGGCTTTGAACTGATCATTTGAAGGTACAGCCTGCGGGATGATGACGGGCTCCGGGGTTTCCGGACAGTTGCGTTACTCAAGGAGCGACGCCTCGTTc                                                                                                   <  4:224939/149‑1 (MQ=255)
                                                cAAACCCGCTCTCCGTCCGCTACCGGCCTGGAAGTCAGCAGTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCg                                                                                                                                              <  6:36277/106‑1 (MQ=255)
                                                cAAACCCGCTCTCCGTCCGCTACCGGCCTGGAAGTCAGCAGTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCg                                                                                                                                              >  5:36277/1‑106 (MQ=255)
                                                cAAACCCGCTCTCCGTCCGCTACCGGCCTGGAAGTCAGCAGTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTc                                                                                                   <  7:101456/149‑1 (MQ=255)
                                                 aaaCCCGCTCTCCGTCCGCTACCGGCCTGGAAGTCAGCAGTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTATAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACAGGt                                                                                                                               <  1:220641/120‑1 (MQ=255)
                                                 aaaCCCGCTCTCCGTCCGCTACCGGCCTGGAAGTCAGCAGTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGt                                                                                                                                           >  2:91303/1‑108 (MQ=255)
                                                 aaaCCCGCTCTCCGTCCGCTACCGGCCTGGAAGTCAGCAGTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGt                                                                                                                                           <  1:91303/108‑1 (MQ=255)
                                                 aaaCCCGCTCTCCGTCCGCTACCGGCCTGGAAGTCAGCAGTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACAGGt                                                                                                                               >  2:220641/1‑120 (MQ=255)
                                                 aaaCCCGCTCTCCGTCCGCTACCGGCCTGGAAGTCAGCAGTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCg                                                                                                  >  4:139299/1‑149 (MQ=255)
                                                     ccGCTCTCCGTCCGCTACCGGCCTGGAAGTCAGCAGTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCca                                                                                               >  8:171584/1‑148 (MQ=255)
                                                           tccgtccgCTACCGGCCTGGAAGTCAGCAGTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCAc                                                                                        >  7:152972/1‑149 (MQ=255)
                                                             cgtccgCTACCGGCCTGGAAGTCAGCAGTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACt                                                                                       <  8:121284/148‑1 (MQ=255)
                                                              gtccgCTACCGGCCTGGAAGTCAGCAGTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGa                                                                                     >  7:129462/1‑149 (MQ=255)
                                                                 cgCTACCGGCCTGGAAGTCAGCAGTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGACtt                                                                                  >  4:27729/1‑149 (MQ=255)
                                                                 cgCTACCGGCCTGGAAGTCAGCAGTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACAGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGACtt                                                                                  >  2:246282/1‑149 (MQ=255)
                                                                      ccGGCCTGGAAGTCAGCAGTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGACTTCGAAc                                                                             >  6:291511/1‑149 (MQ=255)
                                                                                  tCAGCAGTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGACTTCGAACTGCCCCACCAGc                                                                 >  4:44488/1‑149 (MQ=255)
                                                                                     gcagTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGACTTCGAACTGCCCCACCAGCTTc                                                              <  5:291511/149‑1 (MQ=255)
                                                                                      cagTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGACTTCGAACTGCCCCACCAGCTTCt                                                             <  8:152972/149‑1 (MQ=255)
                                                                                       agTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGACTTCGAACTGCCCCACCAGCTTCTg                                                            <  7:52957/149‑1 (MQ=255)
                                                                                           aTCTGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGATACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGACTTCGCACTGCCCCACCAGCTTCTGCAAc                                                        >  5:222023/1‑149 (MQ=255)
                                                                                           aTCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGACTTCGAACTGCCCCACCAGCTTCTGCaa                                                         <  3:147952/148‑1 (MQ=255)
                                                                                            tCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGACTTCGAACTGCCCCACCAGCTTCTGCAACt                                                       >  8:205033/1‑149 (MQ=255)
                                                                                                ggctggctTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGACTTCGAACTGCCCCACCAGCTTCTGCAACTTGcc                                                   <  3:262280/149‑1 (MQ=255)
                                                                                                 gctggctTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGACTTCGAACTGCCCCACCAGCTTCTGCAACTTGCCa                                                  <  6:132839/149‑1 (MQ=255)
                                                                                                  ctggctTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGACTTCGAACTGCCCCACCAGCTTCTGCAACTTGCCAg                                                 >  5:7650/1‑149 (MQ=255)
                                                                                                   tggctTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGACTTCGAACTGCCCCACCAGCTTCTGCAACTTGCCAg                                                 >  6:265407/1‑148 (MQ=255)
                                                                                                     gctTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACAGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGACTTCGAACTGCCCCACCAGCTTCTGCAACTTGCCAGCa                                               <  3:301387/148‑1 (MQ=255)
                                                                                                       ttGCCCCTGATCATTGGCAGGTACAGCCAGCCGGATCATGCCGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGACTTCGAACTGCCCCACCAGCTTCTGCAACTTGCCAGCAt                                              <  7:131443/147‑1 (MQ=255)
                                                                                                       ttGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGACTTCGAACTGCCCCACCAGCTTCTGCAACTTGCCAGCATcc                                            <  8:129462/149‑1 (MQ=255)
                                                                                                       ttGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGGTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGACTTCGCACTGCCCCACCCGCTTCTCCACCTGGCCAGcct                                              >  8:131443/1‑145 (MQ=255)
                                                                                                         gCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGACTTCGAACTGCCCCACCAGCTTCTGCAACTTGCCAGCATCca                                           <  5:265407/148‑1 (MQ=255)
                                                                                                          cACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGACTTCGAACTGCCCCACCAGCTTCTGCAACTTGCCAGCATCcac                                          <  3:139299/148‑1 (MQ=255)
                                                                                                             cTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGa                                                                                                             >  8:120890/1‑78 (MQ=255)
                                                                                                             cTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGa                                                                                                             <  7:120890/78‑1 (MQ=255)
                                                                                                              tGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGACTTCGAACTGCCCCACCAGCTTCTGCAACTTGCCAGCATCcaccttca                                     >  4:292645/1‑149 (MQ=255)
                                                                                                                  cATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATAAATGACTTCGAACTGCCCCACCAGCTTCTGCAACTTGCCAGCATCcaccttcacat                                  >  3:89155/1‑146 (MQ=255)
                                                                                                                                  cagccGGATCATGACAGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGACTTCGAACTGCCCCACCAGCTTCTGCAACTTGCCAGCATCCACCTTCACCTTCTGCGCACTGCTTTgc                 <  1:246282/149‑1 (MQ=255)
                                                                                                                                   agccGGATCATGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGACTTCGAACTGCCCCACCAGCTTCTGCAACTTGCCAGCATCCACCTTCACCTTCTGCGCACTGCTTTgca                <  7:171584/149‑1 (MQ=255)
                                                                                                                                   agccGGATCATGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGACTTCGAACTGCCCCACCAGCTTCTGCAACTTGCCAGCATCCACCTTCACCTTCTGCGCACTGCTTTgca                <  8:85179/149‑1 (MQ=255)
                                                                                                                                   agccGGATCATGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGACTTCGAACTGCCCCACCAGCTTCTGCAACTTGCCAGCATCCACCTTCACCTTCTGCGCACTGCTTTg                  >  3:264818/1‑147 (MQ=255)
                                                                                                                                            aTGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGACTTCGAACTGCCCCACCAGCTTCTGCAACTTGCCAGCATCCACCTTCACCTTCTGCGCACTGCTTTGCAGCACCAgg        <  6:128800/148‑1 (MQ=255)
                                                                                                                                            aTGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGACTTCGAACTGCCCCACCAGCTTCTGCAACTTGCCAGCATCCACCTTCACCTTCTGCGCACTGCTTTGCAGCACCAgg        >  5:128800/1‑148 (MQ=255)
                                                                                                                                             tGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGACTTCGAACTGCCCCACCAGCTTCTGCAACTTGCCAGCATCCACCTTCACCTTCTGCGCACTGCTTTGCAGCACCAGGGc      >  4:65610/1‑149 (MQ=255)
                                                                                                                                               aCGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGACTTCGAACTGCCCCACCAGCTTCTGCAACTTGCCAGCATCCACCTTCACCTTCTGCGCACTGCTTTGCAGCACCAGGGCg     <  6:7650/148‑1 (MQ=255)
                                                                                                                                                 gggCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGACTTCGAACTGCCCCACCAGCTTCTGCAACTTGCCAGCATCCACCTTCACCTTCTGCGCACTGCTTTGCAGCACCAGGGCGGTc  >  5:204801/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CCAACACTGCCACAGACATTGAGCCATTGCGCGGCTTTTTCGGCTGGGCAAACCCGCTCTCCGTCCGCTACCGGCCTGGAAGTCAGCAGTCATCGGGGCTGGCTTGCACCTGATCATTGGCAGGTACAGCCAGCCGGATCATGACGGGCTCCCGGGTTACCGGACAGGTGCGTTACTCAAGCAGCGACGCCTCGTTCGCCATCATCACTGACTTCGAACTGCCCCACCAGCTTCTGCAACTTGCCAGCATCCACCTTCACCTTCTGCGCACTGCTTTGCAGCACCAGGGCGGTC  >  NC_002947/4454631‑4454924

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: ‑