Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,050,680 T→C 6.6% V238A (GTG→GCG)  oprB‑III → carbohydrate‑selective porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,050,6800TC 6.6% 173.5 / 2.5 60V238A (GTG→GCG) oprB‑IIIcarbohydrate‑selective porin
Reads supporting (aligned to +/- strand):  ref base T (27/29);  new base C (2/2);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.84e-01

CCCACCTCGGCGTGCGCGTGAAGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGCCATTACCTGCTGG  >  NC_002947/4050542‑4050820
                                                                                                                                          |                                                                                                                                            
cccACCTCGGCGTGCGCGTGAAGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGTACCACGGTGCCGGtg                                                                                                                                    <  6:367858/149‑2 (MQ=255)
 ccACCTCGGCGTGCGCGTGAAGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAg                                                                                                                                   <  3:128009/149‑1 (MQ=255)
   aCCTCGGCGTGCGCGTGAAGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGaa                                                                                                                                 >  5:142182/1‑149 (MQ=255)
   aCCTCGGCGTGCGCGTGAAGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGaa                                                                                                                                 >  2:51973/1‑149 (MQ=255)
     cTCGGCGTGCGCGTGAAGTACGACCTGTCAGCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACt                                                                                                                               <  5:40389/149‑1 (MQ=255)
              cgcgTGAAGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATAc                                                                                                                        >  5:59091/1‑147 (MQ=255)
                cgTGAAGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTAc                                                                                                                                         >  4:280999/1‑128 (MQ=255)
                cgTGAAGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTAc                                                                                                                                         <  3:280999/128‑1 (MQ=255)
                 gTGAAGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAg                                                                                                                     >  1:99749/1‑147 (MQ=255)
                 gTGAAGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAg                                                                                                                     <  2:99749/147‑1 (MQ=255)
                     aGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAg                                                                                                               <  5:183304/149‑1 (MQ=255)
                      gTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAgg                                                                                                              <  6:142182/149‑1 (MQ=255)
                       tACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGt                                                                                                                            <  6:228743/134‑1 (MQ=255)
                       tACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGt                                                                                                                            >  5:228743/1‑134 (MQ=255)
                         cGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTg                                                                                                           <  4:147787/149‑1 (MQ=255)
                           aCCTGTCAGCTGCCTGGCAGTTGCGCTTGGCGTCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGaa                                                                                                         <  7:159840/149‑1 (MQ=255)
                                  aCCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAg                                                                                                  >  2:68282/1‑149 (MQ=255)
                                  aCCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAg                                                                                                  >  8:168119/1‑149 (MQ=255)
                                   cctccccgGCAGTGGCGCCGGGCGGCGTCCAACCGCGACCCCGACCGCAACGGCACCTCCCGCCGCCCCTGGCACCTCGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTCCCGGGCGAACTGCCAgg                                                                                                 <  8:252647/149‑1 (MQ=255)
                                       cctgGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAg                                                                                             <  3:53849/149‑1 (MQ=255)
                                           gCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTAc                                                                                          >  7:22162/1‑148 (MQ=255)
                                            cAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACaa                                                                                        <  6:224589/149‑1 (MQ=255)
                                                  cgcgTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTg                                                                                   <  1:51973/148‑1 (MQ=255)
                                                            cgTTCAACGTCGACCCCGAAAGCAACGGCAACTCGAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTtagtactac                                                                        <  3:138961/149‑1 (MQ=255)
                                                             gTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTgg                                                                                                                             >  5:178728/1‑95 (MQ=255)
                                                             gTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTgg                                                                                                                             <  6:178728/95‑1 (MQ=255)
                                                                cAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGa                                                                      >  7:317349/1‑147 (MQ=255)
                                                                  aCGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTtac                                                                              <  5:36217/137‑1 (MQ=255)
                                                                  aCGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTtac                                                                              >  6:36217/1‑137 (MQ=255)
                                                                    gTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGctc                                                                <  8:22162/149‑1 (MQ=255)
                                                                        aCCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGAt                                                            <  1:68282/149‑1 (MQ=255)
                                                                           ccGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGAtgt                                                          >  6:77020/1‑148 (MQ=255)
                                                                                         aaCTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAg                                           <  7:168119/149‑1 (MQ=255)
                                                                                              cAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCa                                            <  4:341013/143‑1 (MQ=255)
                                                                                              cAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCa                                            >  3:341013/1‑143 (MQ=255)
                                                                                              cAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCgat                                       <  5:77020/148‑1 (MQ=255)
                                                                                               aGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCgatga                                     >  2:5209/1‑149 (MQ=255)
                                                                                                       cTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAgg                                                                                                 <  5:129897/81‑1 (MQ=255)
                                                                                                       cTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAgg                                                                                                 >  6:129897/1‑81 (MQ=255)
                                                                                                       cTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAgtgt                             >  2:228997/1‑149 (MQ=255)
                                                                                                        tGGCAGCTGGGGCCCAAGCACACCACCGGGACCGTGGTACCGGTAGAAGTGGTATACAAGCTGCAGTCTGAAGTGCCAGGCGAGTAGAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATTAAGTGTc                            <  6:121598/149‑1 (MQ=255)
                                                                                                          gCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCg                          >  3:167885/1‑149 (MQ=255)
                                                                                                                 gggCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAgg                   >  4:80054/1‑149 (MQ=255)
                                                                                                                    cccAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTgg                >  5:10919/1‑149 (MQ=255)
                                                                                                                     ccAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGc               >  6:225916/1‑149 (MQ=255)
                                                                                                                           acacCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGCCATTAc         <  8:317349/149‑1 (MQ=255)
                                                                                                                             accaccGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTg                                                                                   >  7:289290/1‑73 (MQ=255)
                                                                                                                             accaccGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTg                                                                                   <  8:289290/73‑1 (MQ=255)
                                                                                                                             accaccGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGCCAtt           >  6:213330/1‑145 (MQ=255)
                                                                                                                             accaccGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGCCAtt           <  5:213330/145‑1 (MQ=255)
                                                                                                                             accaccGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGCCATTACct       >  1:1718/1‑149 (MQ=255)
                                                                                                                               taccGGTACCACGGTGCCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGCCATTACctgc     >  8:10004/2‑149 (MQ=255)
                                                                                                                               caccGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGaaa                                                      <  8:65655/100‑1 (MQ=255)
                                                                                                                               caccGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGaaa                                                      >  7:65655/1‑100 (MQ=255)
                                                                                                                               caccGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGACGATGAAGTGTCCGGCCGAGGTGGCCATTACctgc     >  8:111023/1‑149 (MQ=255)
                                                                                                                                accGGTACCACGGTGCCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGa                                                                      >  8:94064/1‑83 (MQ=255)
                                                                                                                                accGGTACCACGGTGCCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGa                                                                      <  7:94064/83‑1 (MQ=255)
                                                                                                                                  cGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGCCATTACCTGCTgg  <  3:80054/149‑1 (MQ=255)
                                                                                                                                        cGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGCCAtt           <  4:186502/134‑1 (MQ=255)
                                                                                                                                        cGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGCCAtt           >  3:186502/1‑134 (MQ=255)
                                                                                                                                          tGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGCCa             <  8:317328/130‑1 (MQ=255)
                                                                                                                                          tGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGCCa             >  7:317328/1‑130 (MQ=255)
                                                                                                                                          |                                                                                                                                            
CCCACCTCGGCGTGCGCGTGAAGTACGACCTGTCACCTGCCTGGCAGTTGCGCGTGGCGGCGTTCAACGTCGACCCCGAAAGCAACGGCAACTCCAGCCGCGCCTGGCACCTGGGGCCCAAGCACACCACCGGCACCGTGGTACCGGTAGAACTGGTATACAAGCTGCAGGGTGAACTGCCAGGCGAGTACAAGCTTGGTTACTACTACGACAGCTCTGATGTGAAACGCATCGGCAGCGATGATGAAGTGTCCGGCCGAGGTGGCCATTACCTGCTGG  >  NC_002947/4050542‑4050820

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: