Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,584,495 C→G 6.8% P173P (CCG→CCC)  PP_1388 ← EmrB/QacA family drug resistance transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,584,4950CG 6.8% 184.0 / 4.1 59P173P (CCG→CCC) PP_1388EmrB/QacA family drug resistance transporter
Reads supporting (aligned to +/- strand):  ref base C (25/29);  new base G (2/2);  total (28/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 = 3.98e-01

CTGGCCCACCAAGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGCACGAAGTCGC  >  NC_002947/1584348‑1584640
                                                                                                                                                   |                                                                                                                                                 
ctggcccACCAAGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCg                                                                                                                                                  >  5:43396/1‑149 (MQ=255)
   gcccACCAAGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCa                                                                                                                                               <  8:223545/149‑1 (MQ=255)
   gcccACCAAGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCa                                                                                                                                               <  6:214482/149‑1 (MQ=255)
          aaGGTGATGCCCAGCAGCATGATGCCAATGCCGATAATCAACATCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCTGCCGCTGGGGCTGGTTgct                                                                                                                                        <  2:126714/149‑3 (MQ=255)
          aaGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGAt                                                                                                                                        <  1:145275/149‑1 (MQ=255)
           aGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATc                                                                                                                                       <  3:128872/149‑1 (MQ=255)
           aGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATc                                                                                                                                       <  3:133673/149‑1 (MQ=255)
           aGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATc                                                                                                                                       <  5:180672/149‑1 (MQ=255)
             gTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGTCCCAGCGGCAGGTTGATCCa                                                                                                                                     <  2:17208/149‑1 (MQ=255)
              tGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCg                                                                                                                                                  <  6:132661/135‑1 (MQ=255)
              tGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAg                                                                                                                                    <  5:305006/149‑1 (MQ=255)
              tGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGAATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCg                                                                                                                                                  >  5:132661/1‑135 (MQ=255)
               gATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAg                                                                                                                                    >  8:41200/1‑148 (MQ=255)
                aTGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGaa                                                                                                                                  <  3:299011/149‑1 (MQ=255)
                     cagcagcagGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACAcc                                                                                                                              >  5:285148/1‑148 (MQ=255)
                         agcagGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAgcg                                                                                                                         >  3:199738/1‑149 (MQ=255)
                            agGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCa                                                                                                                      <  4:199738/149‑1 (MQ=255)
                                   cAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCa                                                                                                                                               <  7:59243/117‑1 (MQ=255)
                                   cAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCa                                                                                                                                               >  8:59243/1‑117 (MQ=255)
                                      ggCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTAc                                                                                                            >  7:124038/1‑149 (MQ=255)
                                         cGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCg                                                                                                         >  2:288364/1‑149 (MQ=255)
                                                 aaCAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAAc                                                                                                   >  1:151708/1‑147 (MQ=255)
                                                  aCAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCTGCCGCTGGGGCTGGTTgc                                                                                                                                         <  7:15043/108‑2 (MQ=255)
                                                  aCAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCTGCCGCTGGGGCTGGTTgc                                                                                                                                         >  8:15043/1‑107 (MQ=255)
                                                  aCAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAg                                                                                                <  8:124038/149‑1 (MQ=255)
                                                       aTGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAgccgcc                                                                                           <  7:41200/149‑1 (MQ=255)
                                                            gcCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGCCAGGAACTCGGTCAACCAGCCGCCCAGCa                                                                                      >  8:235532/1‑149 (MQ=255)
                                                                     gTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCg                                                                              <  6:285148/148‑1 (MQ=255)
                                                                        cACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGACCCGGCCa                                                                          >  3:207736/1‑149 (MQ=255)
                                                                           cTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAAc                                                                                                   >  5:295977/1‑121 (MQ=255)
                                                                           cTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAAc                                                                                                   <  6:295977/121‑1 (MQ=255)
                                                                             gtgCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTg                                                                     >  1:312417/1‑149 (MQ=255)
                                                                              tgCCTCCCGCCGCCGCACTGGCATCGCACCGAGCGCACGGCGAATGGCCCACACAGCAACCAGCCCCAGCGGCAGGGTGATCCAGAACACCCCCCGCCACGACCGGTACTCGGGCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTgg                                                                    <  6:230976/149‑1 (MQ=255)
                                                                                     cgacgCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGAGCCAGAACACCCAGCGCCTCGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCgcggc                                                             <  1:303512/149‑1 (MQ=255)
                                                                                      gacgCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGc                                                                                 <  5:173617/128‑1 (MQ=255)
                                                                                      gacgCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGc                                                                                 >  6:173617/1‑128 (MQ=255)
                                                                                      gacgCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCgcggcg                                                            >  2:252451/1‑149 (MQ=255)
                                                                                      gacgCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCgcggcg                                                            <  2:151708/149‑1 (MQ=255)
                                                                                         gCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTa                                                          >  2:20602/1‑148 (MQ=255)
                                                                                             cACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACAtg                                                      <  6:282062/148‑1 (MQ=255)
                                                                                                     tGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAg                                             <  1:288364/149‑1 (MQ=255)
                                                                                                         agcgagcgCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGc                                         >  1:5147/1‑149 (MQ=255)
                                                                                                          gcgagcgCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGcc                                        >  7:60273/1‑149 (MQ=255)
                                                                                                                 cACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTAc                                                         >  7:277727/1‑125 (MQ=255)
                                                                                                                 cACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTAc                                                         <  8:277727/125‑1 (MQ=255)
                                                                                                                          tGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGt                        >  2:99601/1‑149 (MQ=255)
                                                                                                                            gCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTc                      >  2:302645/1‑149 (MQ=255)
                                                                                                                             cccACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTcg                     >  8:54078/1‑149 (MQ=255)
                                                                                                                              ccACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTcgc                    <  1:20602/149‑1 (MQ=255)
                                                                                                                               cacaAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCg                                                                              >  2:254884/1‑90 (MQ=255)
                                                                                                                               cacaAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCg                                                                              <  1:254884/90‑1 (MQ=255)
                                                                                                                                  ggatcAACCTGCCGCTGGGGCTGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCgg                >  7:179095/5‑149 (MQ=255)
                                                                                                                                   gatcAACCTGCCGCTGAGGCTGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCggg               >  8:101860/4‑149 (MQ=255)
                                                                                                                                   aaGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGaa                                              <  4:107593/118‑1 (MQ=255)
                                                                                                                                   aaGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGaa                                              >  3:107593/1‑118 (MQ=255)
                                                                                                                                    aGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCgggg              <  1:214118/149‑1 (MQ=255)
                                                                                                                                      cAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGc            <  1:99601/149‑1 (MQ=255)
                                                                                                                                            gCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGCACGAAg      <  2:5147/149‑1 (MQ=255)
                                                                                                                                                cAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGCACGAGGTCGc  >  8:108331/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CTGGCCCACCAAGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGCACGAAGTCGC  >  NC_002947/1584348‑1584640

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: ‑