Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,774,247 C→G 11.5% A317G (GCC→GGC)  PP_5743 → TonB‑dependent receptor protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,774,2470CG11.5% 142.3 / 10.2 52A317G (GCC→GGC) PP_5743TonB‑dependent receptor protein
Reads supporting (aligned to +/- strand):  ref base C (20/26);  new base G (2/4);  total (22/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.57e-01

CAGATAACGCTTGCCGGCTTCGGTCAGCGCGATGCGTCGGGTGGTGCGGTTGAGCAGGCGAATGTGCAGGTGAGCCTCCAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGC  >  NC_002947/3774099‑3774395
                                                                                                                                                    |                                                                                                                                                    
cAGATCACGCTTGCCGGCTTCGGTCCTCGCGATGCGTCGGGTGGTGCGGTTGAGCAGGCTAATGTGCAGGTGAGCCTCCAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGc                                                                                                                                                      <  7:228363/149‑1 (MQ=255)
cAGATAACGCTTGCCGGCTTCGGTCAGCGCGATGCGTCGGGTGGTGCGGTTGAGCAGGCGAATGTGCAGGTGAGCCTCCAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGc                                                                                                                                                      <  8:204001/149‑1 (MQ=255)
cAGATAACGCTTGCCGGCTTCGGTCAGCGCGATGCGTCGGGTGGTGCGGTTGAGCAGGCGAATGTGCAGGTGAGCCTCCAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGc                                                                                                                                                      <  4:216684/149‑1 (MQ=255)
 aGATAACGCTTGCCGGCTTCGGTCAGCGCGATGCGTCGGGTGGTGCGGTTGAGCAGGCGAATGTGCAGGTGAGCCTCCAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGcc                                                                                                                                                     <  7:181780/149‑1 (MQ=255)
  gATAACGCTTGCCGGCTTCGGTCAGCGCGATGCGTCGGGTGGTGCGGTTGAGCAGGCGAATGTGCAGGTGAGCCTCCAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCa                                                                                                                                                    <  6:22795/149‑1 (MQ=255)
  gATAACGCTTGCCGGCTTCGGTCAGCGCGATGCGTCGGGTGGTGCGGTTGAGCAGGCGAATGTGCAGGTGAGCCTCCAGGATCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCa                                                                                                                                                    <  2:247649/149‑1 (MQ=255)
      aCGCTTGCCGGCTTCGGTCAGCGCGATGCGTCGGGTGGTGCGGTTGAGCAGGCGAATGTGCAGGTGAGCCTCCAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGcc                                                                                                                                                 <  1:222475/148‑1 (MQ=255)
         cTTGCCGGCTTCGGTCAGCGCGATGCGTCGGGTGGTGCGGTTGAGCAGGCGAATGTGCAGGTGAGCCTCCAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGGCAGCCTGac                                                                                                                                             <  8:316216/149‑2 (MQ=255)
         cTTGCCGGCTTCGGTCAGCGCGATGCGTCGGGTGGTGCGGTTGAGCAGGCGAATGTGCAGGTGAGCCTCCAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGGCAGCCTGac                                                                                                                                             <  2:38766/149‑2 (MQ=255)
         cTTGCCGGCTTCGGTCAGCGCGATGCGTCGGGTGGTGCGGTTGAGCAGGCGAATGTGCAGGTGAGCCTCCAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGaa                                                                                                                                             <  1:286147/149‑1 (MQ=255)
          ttGCCGGCTTCGGTCAGCGCGATGCGTCGGGTGGTGCGGTTGAGCAGGCGAATGTGCAGGTGAGCCTCCAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGGCAGCCTGacg                                                                                                                                            <  8:68192/149‑3 (MQ=255)
              cGGCTTCGGTCAGCGCGATGCGTCGGGTGGTGCGGTTGAGCAGGCGAATGTGCAGGTGAGCCTCCAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACg                                                                                                                                        >  7:164062/1‑149 (MQ=255)
               ggCTTCGGTCAGCGCGATGCGTCGGGTGGTGCGGTTGAGCAGGCGAATGTGCAGGTGAGCCTCCAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGc                                                                                                                                       <  8:164062/149‑1 (MQ=255)
                gCTTCGGTCAGCGCGATGCGTCGGGTGGTGCGGTTGAGCAGGCGAATGTGCAGGTGAGCCTCCAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACg                                                                                                                                        <  7:241058/147‑1 (MQ=255)
                gCTTCGGTCAGCGCGATGCGTCGGGTGGTGCGGTTGAGCAGGCGAATGTGCAGGTGAGCCTCCAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACg                                                                                                                                        >  8:241058/1‑147 (MQ=255)
                gCTTCGGTCAGCGCGATGCGTCGGGTGGTGCGGTTGAGCAGGCGAATGTGCAGGTGAGCCTCCAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCa                                                                                                                                      <  1:268808/149‑1 (MQ=255)
                       tCAGCGCGATGCGTCGGGTGGTGCGGTTGAGCAGGCGAATGTGCAGGTGAGCCTCCAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAAc                                                                                                                               <  7:229805/149‑1 (MQ=255)
                                      gggTGGTGCGGTTGAGCAGGCGAATGTGCAGGTGAGCCTCCAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTg                                                                                                                 >  3:217055/1‑148 (MQ=255)
                                       ggtggtGCGGTTGAGCAGGCGAATGTGCAGGTGAGCCTCCAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGaa                                                                                                               <  7:252694/149‑1 (MQ=255)
                                           gtgCGGTTGAGCAGGCGAATGTGCAGGTGAGCCTCCAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTa                                                                                                            <  4:164031/148‑1 (MQ=255)
                                             gcagTTGAGCAGGCGAATGTGCAGGTGAGCCTCCAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTAcat                                                                                                         >  1:56767/4‑149 (MQ=255)
                                                        ggCGAATGTGCAGGTGAGCCTCCAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGt                                                                                              <  2:56767/149‑1 (MQ=255)
                                                           gAATGTGCAGGTGAGCCTCCAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAAc                                                                                           >  3:325920/1‑149 (MQ=255)
                                                             aTGTGCAGGTGAGCCTCCAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCt                                                                                         <  5:175301/149‑1 (MQ=255)
                                                                 gCAGGTGAGCCTCCAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCAt                                                                                      <  3:51399/148‑1 (MQ=255)
                                                                  cAGGTGAGCCTCCAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGt                                                                                    <  1:51054/149‑1 (MQ=255)
                                                                     gTGAGCCTCCAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGaa                                                                                                               <  7:252680/119‑1 (MQ=255)
                                                                     gTGAGCCTCCAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTg                                                                                  >  6:210341/1‑148 (MQ=255)
                                                                     gTGAGCCTCCAGGTTCGACAGCATCCCGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGaa                                                                                                               >  8:252680/1‑119 (MQ=255)
                                                                      tGAGCCTCCAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGGCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGa                                                                                <  6:248623/149‑1 (MQ=255)
                                                                           cTCCAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTACCCTCATGTTGGAGCTTg                                                                           >  2:311057/1‑149 (MQ=255)
                                                                            tCCAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTgg                                                                          >  3:249852/1‑149 (MQ=255)
                                                                             ccAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTggg                                                                         >  1:236983/1‑149 (MQ=255)
                                                                             ccAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTggg                                                                         >  3:190503/1‑149 (MQ=255)
                                                                             ccAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTgg                                                                          >  7:286038/1‑148 (MQ=255)
                                                                                  ttCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAg                                                                    <  2:236983/149‑1 (MQ=255)
                                                                                     gACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGt                                                                 <  4:217055/149‑1 (MQ=255)
                                                                                        aGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGc                                                               >  8:171555/1‑148 (MQ=255)
                                                                                           aTCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCACGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAg                                                            >  1:43958/1‑148 (MQ=255)
                                                                                            tCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGa                                                          <  4:190503/149‑1 (MQ=255)
                                                                                               gggAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGAcgg                                                       <  5:210341/149‑1 (MQ=255)
                                                                                               gggAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGAcgg                                                       >  4:152010/1‑149 (MQ=255)
                                                                                                ggAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTAcatccatc                                                                                                    >  7:109537/1‑103 (MQ=255)
                                                                                                ggAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTAcatccatc                                                                                                    <  8:109537/103‑1 (MQ=255)
                                                                                                          gCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCAc                                            >  6:251691/1‑149 (MQ=255)
                                                                                                             gtggtggtgTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGt                                         <  4:325920/149‑1 (MQ=255)
                                                                                                                   gtgtCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGa                                                          >  6:116389/1‑126 (MQ=255)
                                                                                                                   gtgtCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGa                                                          <  5:116389/126‑1 (MQ=255)
                                                                                                                        cAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCAt                              >  6:254808/1‑149 (MQ=255)
                                                                                                                          ggTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATgg                            <  8:286038/149‑1 (MQ=255)
                                                                                                                            tATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGtt                          <  7:171555/149‑1 (MQ=255)
                                                                                                                                      cgtggttCAGGCTGGCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCAt                              <  4:201281/129‑1 (MQ=255)
                                                                                                                                      cgtggttCAGGCTGGCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCAt                              >  3:201281/7‑135 (MQ=255)
                                                                                                                                         ggttCAGGCTGGCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCg             >  7:169359/4‑149 (MQ=255)
                                                                                                                                          cGTCAGGCTGCCAGCCTGAACCACGCAGTTGCACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCAACGTCAACTGATCATgg                            >  6:239861/1‑133 (MQ=255)
                                                                                                                                          cGTCAGGCTGCCAGCCTGAACCACGCAGTTGAAGGCGCGGATGGCTTGAAGTACTTCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGTCCACCGTCAAATGATCAttg                            <  5:239861/133‑3 (MQ=255)
                                                                                                                                                    ccAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAgcgc  >  1:228519/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAGATAACGCTTGCCGGCTTCGGTCAGCGCGATGCGTCGGGTGGTGCGGTTGAGCAGGCGAATGTGCAGGTGAGCCTCCAGGTTCGACAGCATCCGGGAAACGCTGGCCGTGGTGGTGTCCAGGTATTCGGCCGCACCCGTCAGGCTGCCAGCCTGAACCACGCAGTTGAACGCGCGCATGGCTTGAAGTACATCCATCGTCGGTAACCTCATGTTGGAGCTTGGGCAAAGGGTGCAAGGACGGCGGCCGGCCACCGTCAACTGATCATGGTTTAGCCGGGCACCGGCAGCCAGCGC  >  NC_002947/3774099‑3774395

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