Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 3,547,518 G→A E359K (GAA→AAA)  malT → mal regulon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,547,5180GA100.0% 192.1 / NA 59E359K (GAA→AAA) malTmal regulon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/30);  total (29/30)

ACCGGCGAGTGGTTCTGCTATCACCCGCTGTTTGGTAACTTCCTGCGCCAGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCGAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCC  >  CP009273/3547383‑3547645
                                                                                                                                       |                                                                                                                               
aCCGGCGAGTGGTTCTGCTATCACCCGCTGTTTGGTAACTTCCTGCGCCAGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCaaa                                                                                                                               >  1:4290702/1‑138 (MQ=255)
aCCGGCGAGTGGTTCTGCTATCACCCGCTGTTTGGTAACTTCCTGCGCCAGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCaaa                                                                                                                               >  1:172743/1‑138 (MQ=255)
    gCGAGTGGTTCTGCTATCACCCGCTGTTTGGTAACTTCCTGCGCCAGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAAt                                                                                                                          >  1:560252/1‑139 (MQ=255)
     cGAGTGGTTCTGCTATCACCCGCTGTTTGGTAACTTCCTGCGCCAGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAAt                                                                                                                          >  1:4692105/1‑138 (MQ=255)
         tGGTTCTGCTATCACCCGCTGTTTGGTAACTTCCTGCGCCAGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATtcatc                                                                                                                     >  1:3689443/1‑139 (MQ=255)
              ccgCTATCACCCGCTGTTTGGTACCTTCCTGCGCCAGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATgcg                                                                                                                <  1:4245696/137‑1 (MQ=255)
              cTGCTATCACCCGCTGTTTGGTAACTTCCTGCGCCAGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATgc                                                                                                                 >  1:2184731/1‑138 (MQ=255)
                    tCACCCGCTGTTTGGTAACTTCCTGCGCCAGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTggcg                                                                                                          >  1:2446553/1‑139 (MQ=255)
                      aCCCGCTGTTTGGTAACTTCCTGCGCCAGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCATCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCAGCATGCTCTtgcggc                                                                                                        >  1:1783360/1‑139 (MQ=255)
                      aCCCGCTGTTTGGTAACTTCCTGCGCCAGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGCTTTCCCAGCAAAGCAATTCATCATGCGCTggcggc                                                                                                        >  1:1298155/1‑139 (MQ=255)
                      aCCCGCTGTTTGGTAACTTCCTGCGCCAGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTggcggc                                                                                                        >  1:3645160/1‑139 (MQ=255)
                      aCCCGCTGTTTGGTAACTTCCTGCGCCAGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTggcggc                                                                                                        >  1:3008692/1‑139 (MQ=255)
                      aCCCGCTGTTTGGTAACTTCCTGCGCCAGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTggcggc                                                                                                        >  1:341961/1‑139 (MQ=255)
                      aCCCGCTGTTTGGTAACTTCCTGCGCCAGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTggcggc                                                                                                        >  1:2526699/1‑139 (MQ=255)
                         cGCTGTTTGGTAACTTCCTGCGCCAGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAgg                                                                                                     <  1:4418294/139‑1 (MQ=255)
                           cTGTTTGGTAACTTCCTGCGCCAGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCg                                                                                                   >  1:2656646/1‑139 (MQ=255)
                           cTGTTTGGTAACTTCCTGCGCCAGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCg                                                                                                   >  1:2339307/1‑139 (MQ=255)
                               ttGGTAACTTCCTGCGCCAGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATg                                                                                                <  1:3142592/138‑1 (MQ=255)
                                     aCTTCCTGCGCCAGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTggcggc                                                                                                        <  1:2944062/124‑1 (MQ=255)
                                      cTTCCTGCGCCAGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCTATGCGCTGAtg                                                                                        <  1:1918581/139‑1 (MQ=255)
                                        tCCTGCGCCAGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGTGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGAtgct                                                                                      <  1:1654387/139‑1 (MQ=255)
                                        tCCTGCGCCAGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGAtgct                                                                                      <  1:3361552/139‑1 (MQ=255)
                                        tCCTGCGCCAGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGAt                                                                                         >  1:2314730/1‑136 (MQ=255)
                                        tCCTGCGCCAGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGAt                                                                                         >  1:3061719/1‑136 (MQ=255)
                                        tCCTGCGCCAGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGAt                                                                                         >  1:1377339/1‑136 (MQ=255)
                                           tGCGCCAGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTgcg                                                                                   <  1:1740661/139‑1 (MQ=255)
                                           tGCGCCAGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTgcg                                                                                   <  1:3879433/139‑1 (MQ=255)
                                            gcgcCAGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTgcgc                                                                                  <  1:3448479/139‑1 (MQ=255)
                                               ccAGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATgcg                                                                                              >  1:361447/1‑124 (MQ=255)
                                                cAGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGata                                                                              <  1:496893/139‑1 (MQ=255)
                                                 aGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGatat                                                                             >  1:1050240/1‑139 (MQ=255)
                                                   cgcTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGatat                                                                             >  1:2474761/1‑137 (MQ=255)
                                                    gcTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATAtt                                                                            <  1:3802986/137‑1 (MQ=255)
                                                             gggAACTGTCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAAt                                                                   <  1:3848985/137‑1 (MQ=255)
                                                             gggAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAAt                                                                   <  1:2755642/137‑1 (MQ=255)
                                                                  cTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGccc                                                            <  1:2112485/139‑2 (MQ=255)
                                                                    ggcggcGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTgg                                                          >  1:3548106/1‑139 (MQ=255)
                                                                      cggcggAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGa                                                         <  1:73575/138‑1 (MQ=255)
                                                                      cggcggAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGa                                                         <  1:1710294/138‑1 (MQ=255)
                                                                      cggcggAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGa                                                         <  1:1287423/138‑1 (MQ=255)
                                                                       ggcggAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAg                                                        <  1:980619/138‑1 (MQ=255)
                                                                            aGCTGCCGGAAATCCCCCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGGAGGCGATTCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGt                                                   <  1:426189/138‑1 (MQ=255)
                                                                               tGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCa                                                >  1:958396/1‑138 (MQ=255)
                                                                               tGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCa                                                >  1:4531021/1‑138 (MQ=255)
                                                                                 ccGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCa                                                                 <  1:685762/119‑1 (MQ=255)
                                                                                     aaaTCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCagcaaagcaa                                                                                                                           <  1:3116674/57‑1 (MQ=255)
                                                                                       aTCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGa                                                                                                  <  1:1268194/80‑1 (MQ=255)
                                                                                           aCCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGATGCGCGATATTCTGCTTAATCACGCCTGGAGt                                                       >  1:3267930/1‑119 (MQ=255)
                                                                                             cGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGt                                 <  1:383132/139‑1 (MQ=255)
                                                                                               tGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATgc                                                                                                                 >  1:534299/1‑57 (MQ=255)
                                                                                               tGCCGCCGCAGAAAGCTGGATGGCCCAGGGAATTCACAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGCTGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCCACCATAGCGAACTGTc                                >  1:2927438/1‑138 (MQ=255)
                                                                                                 ccgccgCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGAtg                                                                                        >  1:3849518/1‑80 (MQ=255)
                                                                                                      gcAGAAAGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTgg                        <  1:2532603/139‑1 (MQ=255)
                                                                                                           aaGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAagag                   <  1:4764193/139‑1 (MQ=255)
                                                                                                           aaGCTGGATGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAag                     <  1:5105506/137‑1 (MQ=255)
                                                                                                                   tGGCCCAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTa            <  1:2413387/138‑1 (MQ=255)
                                                                                                                      cccAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAgg         >  1:4992369/1‑138 (MQ=255)
                                                                                                                      cccAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGAGGCGCTGATGCTGCGCGATATTCTGCTTAATCGCGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAgg         <  1:1951042/138‑1 (MQ=255)
                                                                                                                        cAGTGATTTCCGAGGAAAGCAATTCACTATGCGCTGGCGGTAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCTCGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGccc      <  1:4363653/139‑1 (MQ=255)
                                                                                                                        cAGGGATTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGc        >  1:485444/1‑137 (MQ=255)
                                                                                                                             aTTTCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCGCGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGcc  <  1:4918401/138‑1 (MQ=255)
                                                                                                                                tCCCAGCAAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGACGAGTCGCTTAAgg         >  1:1971667/1‑128 (MQ=255)
                                                                                                                                       |                                                                                                                               
ACCGGCGAGTGGTTCTGCTATCACCCGCTGTTTGGTAACTTCCTGCGCCAGCGCTGCCAGTGGGAACTGGCGGCGGAGCTGCCGGAAATCCACCGTGCCGCCGCAGAAAGCTGGATGGCCCAGGGATTTCCCAGCGAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCC  >  CP009273/3547383‑3547645

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: