Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,184,585 G→C 18.0% intergenic (‑201/‑211) mltF ← / → purL membrane‑bound lytic murein transglycosylase F/phosphoribosylformylglycinamidine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,184,5850GC18.0% 120.1 / 20.1 50intergenic (‑201/‑211)mltF/purLmembrane‑bound lytic murein transglycosylase F/phosphoribosylformylglycinamidine synthase
Reads supporting (aligned to +/- strand):  ref base G (17/24);  new base C (5/4);  total (22/28)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.24e-01

TTGTCTGGTAGCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGA  >  NC_002947/1184438‑1184730
                                                                                                                                                   |                                                                                                                                                  
ttGTCTGGTAGCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCcgcgc                                                                                                                                                   >  3:62605/1‑149 (MQ=255)
     tGGTAGCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCgcgggcg                                                                                                                                              >  5:127751/1‑149 (MQ=255)
          gCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGgcgc                                                                                                                                          <  1:88365/148‑1 (MQ=255)
           cgggtcgggtGGTTGGGGGAGGGGGGATGGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                        <  6:115559/149‑1 (MQ=255)
              gTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGAt                                                                                                                                     <  1:233871/149‑1 (MQ=255)
                 gtgtgGTTGGGGGAGGGGGGATGGCATCAGGGCTGCCAGCTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCGTAAGATATTAAGCGCCTAAGAGAAAGAGCGCCGCGCGGGCGGCgctcgctcg                                                                                                                                   >  6:100926/1‑147 (MQ=37)
                    tgGTTGGGGGAGGGGGGATTGCATCAGGGCCGCCAGGTGGCCGGCTTGAGATATTAAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                                                 >  3:136538/1‑147 (MQ=255)
                      gTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAgg                                                                                                                             >  8:192760/1‑149 (MQ=255)
                               ggggggATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGACATTCTGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCTGAAta                                                                                                                    <  8:88275/149‑2 (MQ=21)
                                gggggATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAAtgt                                                                                                                   <  6:176786/149‑1 (MQ=255)
                                         aaTCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGATATCGAGCGCCGCCAGGTGTACGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCAAAGGCGCAGAATGTGTCGAGcag                                                                                                          >  5:151964/2‑147 (MQ=255)
                                                   gCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTgcg                                                                                               <  4:55835/149‑1 (MQ=255)
                                                    ccAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTgcgg                                                                                              >  7:188095/1‑149 (MQ=255)
                                                    ccAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTgcgg                                                                                              >  4:202323/1‑149 (MQ=255)
                                                        gcggCCGGCTTGAGATATTCCGCGCCGGTGAGCTCGAGCGCCGCCAGGTGGCCGCCTTAAGATATGCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCt                                                                                          <  5:366711/147‑1 (MQ=38)
                                                          ggCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCg                                                                                        <  3:74320/149‑1 (MQ=255)
                                                           gCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGACGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGa                                                                                       >  7:294052/1‑149 (MQ=255)
                                                               gCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCgcgggcgg                                                                                                                                             <  7:130635/92‑1 (MQ=255)
                                                               gCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCgcgggcgg                                                                                                                                             >  8:130635/1‑92 (MQ=255)
                                                                cTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGacac                                                                                                     <  6:170931/131‑1 (MQ=255)
                                                                cTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGacac                                                                                                     >  5:170931/1‑131 (MQ=255)
                                                                   gagaTATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGtct                                                                               >  3:353192/1‑149 (MQ=255)
                                                                   gagaTATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCAATCTCGt                                                                                 <  6:340372/147‑1 (MQ=255)
                                                                         ttCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCccgcgcgg                                                                                                                                             <  8:120469/82‑1 (MQ=255)
                                                                         ttCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCccgcgcgg                                                                                                                                             >  7:120469/1‑82 (MQ=255)
                                                                           cAGCGCCTGTTAGATCGAGCGCCGCCAGGTGTCCTCCTTAAGATTTTCAGCGCCTTTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCg                                                                        <  4:62605/148‑1 (MQ=255)
                                                                           cAGCGACTGTGAGTGCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATGGCGCGCGGCGCGGGCGTCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCg                                                                        <  5:100926/148‑1 (MQ=38)
                                                                            aGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGa                                                                      >  5:277138/1‑149 (MQ=255)
                                                                               gcCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGacac                                                                   <  4:136538/149‑1 (MQ=255)
                                                                                 cTGTGAGATCGAGCGACTCCAGGTGTCCGCATTAAGATATTCAGCTCCTATGAGATCTAGCGCCGCGCGGGCTGCGCTCGATCTCACAGGCGCATAATGTGTCGAGCCTGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACct                                                                 <  6:151964/149‑1 (MQ=38)
                                                                                           gAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCATGCCGGACACCTCTGAAAACTg                                                       >  7:122805/1‑149 (MQ=255)
                                                                                                   ccAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCg                                                 >  8:32779/1‑147 (MQ=255)
                                                                                                          tCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGt                                        <  6:277138/149‑1 (MQ=255)
                                                                                                                ttAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCCCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTccgcc                                  >  3:355176/1‑149 (MQ=255)
                                                                                                                    gATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGccc                              <  6:127751/149‑1 (MQ=255)
                                                                                                                         tCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAg                         >  4:273664/1‑149 (MQ=255)
                                                                                                                           aGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGc                       >  8:89346/1‑149 (MQ=255)
                                                                                                                           aGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGc                       >  4:270990/1‑149 (MQ=255)
                                                                                                                             cgcCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGt                     <  7:32779/149‑1 (MQ=255)
                                                                                                                               ccTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTgg                   <  3:273664/149‑1 (MQ=255)
                                                                                                                                cTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGc                  <  3:202323/149‑1 (MQ=255)
                                                                                                                                  gtGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCTGAATATGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCgg                            <  4:125146/136‑1 (MQ=25)
                                                                                                                                  gtGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCTGAATATGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCgg                            >  3:125146/2‑137 (MQ=25)
                                                                                                                                  gtGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCTGAATATCTTAAGGCGGACACCTGGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCt                 >  2:229915/2‑149 (MQ=17)
                                                                                                                                    gagaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTa               <  8:122805/148‑1 (MQ=255)
                                                                                                                                     agaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCCCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTaga             <  4:355176/149‑1 (MQ=255)
                                                                                                                                     agaTCGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCTGAATATCTTAAGGCGGACACCTGGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTa               <  1:229915/148‑1 (MQ=18)
                                                                                                                                      gaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTccg                                    >  1:225553/1‑125 (MQ=255)
                                                                                                                                      gaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTccg                                    <  2:225553/125‑1 (MQ=255)
                                                                                                                                        tCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCa          <  7:89346/149‑1 (MQ=255)
                                                                                                                                        tCGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCTGAATATGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGt                     <  7:296532/138‑1 (MQ=37)
                                                                                                                                        tCGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCTGAATATGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGt                     >  8:296532/1‑138 (MQ=37)
                                                                                                                                         cGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCTGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAg         >  4:273466/1‑149 (MQ=255)
                                                                                                                                          gAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCTGAATATCTTAAGGCGGACACCTGGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAg         >  3:133054/1‑149 (MQ=18)
                                                                                                                                           aGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGt       <  7:192760/149‑1 (MQ=255)
                                                                                                                                             cgccgcGCTGTCGGAGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGc     <  8:323308/149‑1 (MQ=255)
                                                                                                                                                 gcgcgGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCAGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGa  >  1:275566/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTGTCTGGTAGCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGA  >  NC_002947/1184438‑1184730

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: