Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 556,492 T→C 7.1% K21R (AAA→AGA)  ybcI ← DUF457 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913556,4920TC 7.1% 165.3 / 3.0 56K21R (AAA→AGA) ybcIDUF457 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (29/23);  new base C (2/2);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

CAAACGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGT  >  NC_000913/556349‑556639
                                                                                                                                               |                                                                                                                                                   
caaACGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAAAAACAGACGCGGAGGGATGACGTTctca                                                                                                                                                >  3:10308/1‑146 (MQ=255)
   aCGCAAACAACAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGaaccc                                                                                                                                             <  1:300731/149‑1 (MQ=255)
             taGCGAATGGGTAAACCCGCGATTACAAAATTCATTGCCGTAAGCAACGCCAAATTTTAACGACAATACGTGGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAAc                                                                                                                                   <  4:10308/148‑1 (MQ=255)
              aGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAAc                                                                                                                                   >  3:210472/1‑148 (MQ=255)
                      ggTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTa                                                                                                                                      >  7:300769/1‑137 (MQ=255)
                      ggTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTa                                                                                                                                      <  8:300769/137‑1 (MQ=255)
                       gTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGGTGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAg                                                                                                                         >  8:141335/1‑149 (MQ=255)
                         aaaCCCGCGATGACCCAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGAGGACTTTCGAACCCAACCCTAAACCAATGCATAggg                                                                                                                       >  7:160996/1‑149 (MQ=255)
                          aaCCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAgggg                                                                                                                      <  1:73998/149‑1 (MQ=255)
                           aCCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAgggg                                                                                                                      >  4:175855/1‑148 (MQ=255)
                             ccGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAACCCAATGCATAGGGGAAc                                                                                                                   >  4:195843/1‑149 (MQ=255)
                             ccGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAAc                                                                                                                   >  4:100955/1‑149 (MQ=255)
                                      ccAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGgcgtg                                                                                                          >  1:101440/1‑149 (MQ=255)
                                       cAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACACAACCCTAAACCAATGCATAGGGGAACGGCGgcgtgc                                                                                                         >  2:228058/1‑149 (MQ=255)
                                        aaaaaCATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGgcgtgc                                                                                                         >  4:233080/1‑148 (MQ=255)
                                         aaaaCATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCa                                                                                                                            <  5:278125/128‑1 (MQ=255)
                                         aaaaCATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCa                                                                                                                            >  6:278125/1‑128 (MQ=255)
                                          aaaCATTGCCGTAAGCAACGCCAAATTTAAACGACAATCCGTCGGCGTCGGGCAGCATCGCCAGGATGATGCCGGCAACTAACAGACGCGGAGGGATGACGTTCGAACCCAACCCTAAACCAAGGCGGGGGGGACCGGGGgcgggcgt                                                                                                       >  4:228576/1‑148 (MQ=255)
                                             cATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGtaata                                                                                                   <  8:28229/149‑1 (MQ=255)
                                              aTTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGtaataa                                                                                                  >  5:282680/1‑149 (MQ=255)
                                                 gCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACg                                                                                                <  6:243608/148‑1 (MQ=255)
                                                    gTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTg                                                                                            <  8:160996/149‑1 (MQ=255)
                                                        gCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACttt                                                                                                                                                    >  5:224330/1‑89 (MQ=255)
                                                        gCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACttt                                                                                                                                                    <  6:224330/89‑1 (MQ=255)
                                                        gCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCat                                                                                        <  3:175855/149‑1 (MQ=255)
                                                          aaCGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCatga                                                                                      >  3:320397/1‑149 (MQ=255)
                                                           aCGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCatgat                                                                                     <  3:270761/149‑1 (MQ=255)
                                                             gCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGCTGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATgg                                                                                   <  7:141335/149‑1 (MQ=255)
                                                                      aaaCGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCgg                                                                          >  4:235144/1‑149 (MQ=255)
                                                                                  gtcggcgtcggGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGAt                                                              <  3:100955/149‑1 (MQ=255)
                                                                                  gtcggcgtcggGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTCTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGAt                                                              >  4:245970/1‑149 (MQ=255)
                                                                                  gtcggcgtcggGCAGCATCGCCAGGATGAGTCCGGCAAATAACAGACGCGGAGGGATGCCTCTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGTTCGCTTCTGCAAAATGTAGAt                                                              >  4:241373/1‑149 (MQ=255)
                                                                                   tcggcgtcggGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATg                                                             <  3:235144/149‑1 (MQ=255)
                                                                                   tcggcggcggGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATg                                                             <  4:210472/149‑1 (MQ=255)
                                                                                    cggcgtcggGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGc                                                            >  2:269758/1‑149 (MQ=255)
                                                                                         tcggGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCa                                                       <  3:195843/149‑1 (MQ=255)
                                                                                                  aTCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTa                                                          >  4:312698/1‑137 (MQ=255)
                                                                                                  aTCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTa                                                          <  3:312698/137‑1 (MQ=255)
                                                                                                   tCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGc                                                                               <  6:269337/115‑1 (MQ=255)
                                                                                                   tCGCAAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGc                                                                               >  5:269337/1‑115 (MQ=255)
                                                                                                         ggATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACgggg                                       >  5:329981/1‑149 (MQ=255)
                                                                                                              atTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGa                                                <  7:358153/135‑1 (MQ=255)
                                                                                                              atTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGa                                                >  8:358153/1‑135 (MQ=255)
                                                                                                                 ccGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTagag                                <  2:101440/148‑1 (MQ=255)
                                                                                                                  cGGCAAATAACTGACGAGGAGGGATGAATTTCGAACCCAACGCTAAACCATTGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACggg                                        <  4:187771/139‑1 (MQ=255)
                                                                                                                  cGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAggg                              >  6:249883/1‑149 (MQ=255)
                                                                                                                  cGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACCGCGGCGTGCGTAATAACGGATGGCAAGATGGTCGCTTCCGCAAAATGTCGAAGCTATCAGCATGGATGAACggg                                        >  3:187771/1‑139 (MQ=255)
                                                                                                                         tAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGt                       >  1:122601/1‑149 (MQ=255)
                                                                                                                          aaCAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTc                      >  4:53085/1‑149 (MQ=255)
                                                                                                                                 gcgGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGaa                  >  4:96015/1‑146 (MQ=255)
                                                                                                                                 gcgGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGaa                  >  4:91325/1‑146 (MQ=255)
                                                                                                                                  cgGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAaga              <  6:329981/149‑1 (MQ=255)
                                                                                                                                         aTGACTTTCGAACCCAACCCTAACCCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAATTATGAAAAGAGAACCgg       >  4:210193/1‑149 (MQ=255)
                                                                                                                                          tGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGaa           >  5:218375/1‑144 (MQ=255)
                                                                                                                                          tGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGaa           <  6:218375/144‑1 (MQ=255)
                                                                                                                                          tGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGAGGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGc      <  3:233080/149‑1 (MQ=255)
                                                                                                                                          tGACTCTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGc      <  3:245970/149‑1 (MQ=255)
                                                                                                                                           gACTCTCTAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGAGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGc      <  3:241373/148‑1 (MQ=255)
                                                                                                                                              tttCGAACCCACCCCTAAACCCATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGAGGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTTGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGt  >  2:225372/1‑149 (MQ=255)
                                                                                                                                              tttCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGc                                                                                          <  6:209884/61‑1 (MQ=255)
                                                                                                                                              tttCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGc                                                                                          >  5:209884/1‑61 (MQ=255)
                                                                                                                                              tttCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGAGTAATAACGCTTCGCATGATGGTCGCTTAGGCAAAATGTCGATGCTATCAGAATGGATGAAAGGGGCGTAGAGGGAAAAAGTCTGAAAAGATAACAGGCCAGt  >  3:262214/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CAAACGCAAACACCAGCGAATGGGTAAACCCGCGATGACCAAAAACATTGCCGTAAGCAACGCCAAATTTAAACGACAATACGTCGGCGTCGGGCAGCATCGCCAGGATGATTCCGGCAAATAACAGACGCGGAGGGATGACTTTCGAACCCAACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGT  >  NC_000913/556349‑556639

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