Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,385,797:1 +A 7.1% coding (81/576 nt) PP_2093 ← two‑component system response regulator NasT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,385,7971.A 7.1% 221.8 / 10.2 56coding (81/576 nt)PP_2093two‑component system response regulator NasT
Reads supporting (aligned to +/- strand):  ref base . (24/28);  new base A (2/2);  total (26/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 = 4.11e-01

GATGGGCCGCGGCTGGTCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCATGCAAGCTGAAA  >  NC_002947/2385650‑2385940
                                                                                                                                                     |                                                                                                                                               
gATGGGCCGCGGCTGGTCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGAtct                                                                                                                                                >  7:157243/1‑149 (MQ=255)
     gCCGCGGCTGGTCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGAGATCATCGAAg                                                                                                                                           <  7:147192/149‑1 (MQ=255)
     gCCGCGGCTGGTCACGGCTGACCAGCACTACCTGTTACATCACATCGCGGTCTGGCGAGTCGGTATCGATAAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAg                                                                                                                                           >  7:258712/1‑149 (MQ=255)
        gcgGCTGGTCACGGCTGCCCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCg                                                                                                                                        >  5:337555/1‑149 (MQ=255)
         cgGCTGGTCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCgg                                                                                                                                       <  8:293651/149‑1 (MQ=255)
                    ggCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAAc                                                                                                                            <  7:339772/149‑1 (MQ=255)
                      cTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGAGCTCGACGCCGGCCTCACGCACCGc                                                                                                                          >  6:115496/1‑149 (MQ=255)
                       tGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACg                                                                                                                           >  2:43235/1‑147 (MQ=255)
                         aCCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAg                                                                                                                        <  4:329361/148‑1 (MQ=255)
                         aCCAGCACTACCTGCTCCATCAGATCGCGGTCTGGCGAGTCGGAATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCa                                                                                                                         <  7:48994/147‑1 (MQ=255)
                           cAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCg                                                                                                                                        >  7:318897/1‑130 (MQ=255)
                           cAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCg                                                                                                                                        <  8:318897/130‑1 (MQ=255)
                             gCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTc                                                                                                                   <  6:370351/149‑1 (MQ=255)
                                 tACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGc                                                                                                               >  7:226719/1‑149 (MQ=255)
                                   ccTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGAtc                                                                                                                                                  >  5:53146/1‑113 (MQ=255)
                                   ccTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGAtc                                                                                                                                                  <  6:53146/113‑1 (MQ=255)
                                   ccTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGa                                                                                                             >  2:101940/1‑149 (MQ=255)
                                      gTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAgg                                                                                                                <  2:353482/143‑1 (MQ=255)
                                      gTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAgg                                                                                                                >  1:353482/1‑143 (MQ=255)
                                       ttCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGAcc                                                                                                           >  1:184879/1‑147 (MQ=255)
                                           aTCACGTCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACAAACCCCTCAGACGACAGATCtt                                                                                                     >  6:290193/1‑149 (MQ=255)
                                           aTCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCtt                                                                                                     <  2:253279/149‑1 (MQ=255)
                                            tCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTc                                                                                                                   <  4:263893/134‑1 (MQ=255)
                                            tCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTc                                                                                                                   >  3:263893/1‑134 (MQ=255)
                                                aTCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGctct                                                                                                >  6:17345/1‑149 (MQ=255)
                                                   gcgGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTg                                                                                               <  1:101940/148‑1 (MQ=255)
                                                   gcgGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTgg                                                                                              <  8:226719/148‑1 (MQ=255)
                                                       tCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTgtc                                                                                         <  2:184879/149‑1 (MQ=255)
                                                         gggCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTgtcgt                                                                                       <  4:283320/148‑1 (MQ=255)
                                                                   gTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATcagcag                                                                             >  6:67168/1‑149 (MQ=255)
                                                                   gTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGAGATCATCGAAGcccca                                                                                                                                      >  2:410172/1‑89 (MQ=38)
                                                                   gTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGAGATCATCGAAGcccca                                                                                                                                      <  1:410172/92‑4 (MQ=38)
                                                                        gATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATgc                                                                        <  8:157243/149‑1 (MQ=255)
                                                                            gtatccacGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCaa                                                                    >  4:130028/5‑149 (MQ=255)
                                                                             gCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCaa                                                                    >  1:41065/1‑148 (MQ=255)
                                                                             gCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAAc                                                                   <  2:41065/149‑1 (MQ=255)
                                                                               accacGTCCGGGCGCACCGTTGCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCTGCCTCAGTCAACGGAGTGCTCAGGCGACTGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACAt                                                                 <  5:290193/149‑1 (MQ=255)
                                                                                ccacGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATg                                                                >  7:42103/1‑149 (MQ=255)
                                                                                 cacGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGc                                                                                                       <  4:307760/109‑1 (MQ=255)
                                                                                 cacGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGc                                                                                                       >  3:307760/1‑109 (MQ=255)
                                                                                 cacGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTgg                                                                                              <  4:290735/118‑1 (MQ=255)
                                                                                 cacGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTgg                                                                                              >  3:290735/1‑118 (MQ=255)
                                                                                  acGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGtt                                                              <  5:67168/149‑1 (MQ=255)
                                                                                    gTCCGGGCGCACCGTTTCGACGCGGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCg                                                            >  5:227499/1‑149 (MQ=255)
                                                                                         ggCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTcccc                                                       <  8:42103/149‑1 (MQ=255)
                                                                                                  tttCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACgcgc                                              <  8:258712/149‑1 (MQ=255)
                                                                                                      gACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGac                                          >  4:323582/1‑149 (MQ=255)
                                                                                                                   aGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCaa                                                                    <  2:347159/110‑1 (MQ=255)
                                                                                                                   aGGTCGATGGTCAGGGCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCGCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGc                             <  6:337555/149‑1 (MQ=255)
                                                                                                                             tCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTccccc                                                      <  8:79978/114‑1 (MQ=255)
                                                                                                                             tCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTccccc                                                      >  7:79978/1‑114 (MQ=255)
                                                                                                                             tCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGccc                    >  8:174854/1‑148 (MQ=255)
                                                                                                                             tCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCACCGGGCAACTCGccc                    >  3:176959/1‑148 (MQ=255)
                                                                                                                               aGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAg                 <  1:43235/149‑1 (MQ=255)
                                                                                                                                   ggccggCTTCGAGATCATCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCAt             >  7:224520/5‑149 (MQ=255)
                                                                                                                                    tGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCATg            >  6:205841/1‑149 (MQ=255)
                                                                                                                                      gggCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCATGc           <  3:130028/148‑1 (MQ=255)
                                                                                                                                        gCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCATGCAAg        <  5:17345/149‑1 (MQ=255)
                                                                                                                                              atgatCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTGCGTCCCCCTCACGCGCCGACACGGGCAGCGGGCAACTCGCCCAAGGCATGCAAGCTGaaa  <  5:115496/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
GATGGGCCGCGGCTGGTCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATCTCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCATGCAAGCTGAAA  >  NC_002947/2385650‑2385940

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: