Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 467,101 A→G 10.0% intergenic (‑31/+22) PP_5445 ← / ← PP_t71 hypothetical protein/tRNA‑Met

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947467,1010AG10.0% 160.3 / 7.5 60intergenic (‑31/+22)PP_5445/PP_t71hypothetical protein/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base A (24/30);  new base G (3/3);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

GGTGACCGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCTGTGCCATC  >  NC_002947/466963‑467245
                                                                                                                                          |                                                                                                                                                
ggTGACCGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCt                                                                                                                                       >  2:854486/1‑150 (MQ=255)
  tGACCGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCGATTGCGCGGCttt                                                                                                                                     <  1:946505/150‑1 (MQ=255)
  tGACCGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTAAGCGCCCCTGACCAGTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCttt                                                                                                                                     >  1:977447/1‑150 (MQ=255)
   gACCGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCGATTGCGCGGCtttt                                                                                                                                    <  1:152426/150‑2 (MQ=255)
     ccGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTc                                                                                                                                  <  2:498212/150‑1 (MQ=255)
     ccGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTc                                                                                                                                  <  2:341643/150‑1 (MQ=255)
        gCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTa                                                                                                                               >  2:945212/1‑150 (MQ=255)
          ttcttGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATc                                                                                                                             <  1:854486/150‑1 (MQ=255)
                      cTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCcac                                                                                                                 >  1:257266/1‑150 (MQ=255)
                           ttCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGa                                                                                                            <  1:486936/150‑1 (MQ=255)
                           ttCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGa                                                                                                            <  2:257266/150‑1 (MQ=255)
                             cGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACt                                                                                                          <  2:544931/150‑1 (MQ=255)
                               cTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTAGGCCCACACGGACTCg                                                                                                        >  1:1036796/1‑150 (MQ=255)
                                  agaTACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAAc                                                                                                     <  1:107324/150‑1 (MQ=255)
                                      aCCGTTTTAGCGCTACTTCCATGTTCATTCGCTCGACGCCGCTGCGCTTGATATCCTCGCCCCTGACCATTTTTTTCGTTCGTGCTTGATAAATCCGCTCAATCTCGCAGCTTTCTCGTATCTCGTGGGCCCACGCGGACTCAAACCGTg                                                                                                 <  2:508707/150‑1 (MQ=255)
                                      aCCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGAGCGTCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTg                                                                                                 >  1:508707/1‑150 (MQ=255)
                                       ccGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTgg                                                                                                >  2:802768/1‑150 (MQ=255)
                                          ttttAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATcgcgcg                                                                                                                                          >  2:750265/1‑105 (MQ=255)
                                          ttttAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATcgcgcg                                                                                                                                          <  1:750265/105‑1 (MQ=255)
                                             tAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCaaa                                                                                          <  1:789071/150‑1 (MQ=255)
                                              aGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAg                                                                                         <  2:21250/150‑1 (MQ=255)
                                               gcgcTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTc                                                                                                                                  >  2:1109976/1‑108 (MQ=255)
                                               gcgcTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTc                                                                                                                                  <  1:1109976/108‑1 (MQ=255)
                                                   tACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTa                                                                                    <  1:712510/150‑1 (MQ=255)
                                                      ttCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGCTGCCTGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGa                                                                                 <  1:336165/150‑1 (MQ=255)
                                                      ttCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGa                                                                                 >  2:529056/1‑150 (MQ=255)
                                                        ccAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGt                                                                               >  2:848204/1‑150 (MQ=255)
                                                         cAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCGATTGCGCGGCtt                                                                                                                                      >  1:177639/1‑94 (MQ=255)
                                                         cAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCGATTGCGCGGCtt                                                                                                                                      <  2:177639/94‑1 (MQ=255)
                                                            ggTCATGCGCTCGGCGCCGCTACGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCctc                                                                           >  1:1037103/1‑150 (MQ=255)
                                                                     cTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAAc                                                                  <  1:513354/150‑1 (MQ=255)
                                                                         gcgcCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCaa                                                               >  1:620583/1‑149 (MQ=255)
                                                                         gcgcCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCaa                                                               <  2:620583/149‑1 (MQ=255)
                                                                            ccgcTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGa                                                           <  1:848204/150‑1 (MQ=255)
                                                                             cgctgcgctgGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAg                                                          <  1:529056/150‑1 (MQ=255)
                                                                              gctgcgctgGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGc                                                         >  2:1053258/1‑150 (MQ=255)
                                                                                   gctgGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCAGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTAt                                                                                   <  2:1052891/119‑1 (MQ=255)
                                                                                   gctgGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCAGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTAt                                                                                   >  1:1052891/1‑119 (MQ=255)
                                                                                     tgGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCAAACCGTCGCCCAAATGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGc                                                  >  1:63174/1‑150 (MQ=255)
                                                                                       gATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCaaa                                                                                          <  1:104827/108‑1 (MQ=255)
                                                                                       gATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCaaa                                                                                          >  2:104827/1‑108 (MQ=255)
                                                                                       gATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAg                                                          <  1:687789/140‑1 (MQ=255)
                                                                                       gATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAg                                                          >  2:687789/1‑140 (MQ=255)
                                                                                             cgcgCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCtggtgg                                                                                                                       <  1:859162/73‑1 (MQ=255)
                                                                                             cgcgCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCtggtgg                                                                                                                       >  2:859162/1‑73 (MQ=255)
                                                                                                gcCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGt                                       <  1:680577/150‑1 (MQ=255)
                                                                                                   ccTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTggg                                    >  2:362959/1‑150 (MQ=255)
                                                                                                        ccctttatttCTTGCGGGCATGAAAAAGCCGCGCAATAACCTGTCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTTTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGAt                               <  1:851345/147‑1 (MQ=255)
                                                                                                         cATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGAtt                              >  2:1006843/1‑150 (MQ=255)
                                                                                                                 tCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAg                      <  1:802768/150‑1 (MQ=255)
                                                                                                                       gggCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTAGCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTc                <  1:560333/150‑1 (MQ=255)
                                                                                                                           aTGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCAc            <  1:945212/150‑1 (MQ=255)
                                                                                                                                aaaGCCGCGCGATTGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAAAGAAGGTTTTGCACCtata       >  2:237022/1‑147 (MQ=255)
                                                                                                                                   gCCGCGCGATTGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCTGTGCCa    >  2:839469/1‑150 (MQ=255)
                                                                                                                                   gCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGtt                  >  2:728036/1‑136 (MQ=255)
                                                                                                                                   gCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGtt                  >  2:728150/1‑136 (MQ=255)
                                                                                                                                   gCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGtt                  <  1:728036/136‑1 (MQ=255)
                                                                                                                                   gCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGtt                  <  1:728150/136‑1 (MQ=255)
                                                                                                                                    ccGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCTGTGCCAt   <  2:234882/150‑1 (MQ=255)
                                                                                                                                     cgcgCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCTGTGCCATc  >  2:772313/1‑150 (MQ=255)
                                                                                                                                          |                                                                                                                                                
GGTGACCGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCTGTGCCATC  >  NC_002947/466963‑467245

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