Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 705,916 T→C 7.0% intergenic (+26/‑177) nagE → / → glnS N‑acetyl glucosamine specific PTS enzyme IIC, IIB, and IIA components/glutamyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913705,9160TC 7.0% 176.7 / 3.9 57intergenic (+26/‑177)nagE/glnSN‑acetyl glucosamine specific PTS enzyme IIC, IIB, and IIA components/glutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (31/22);  new base C (2/2);  total (33/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.08e-01

CGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATACGCCGTTA  >  NC_000913/705770‑706063
                                                                                                                                                  |                                                                                                                                                   
cGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTGATCTGCTTTATGCCTGATGCGACGCTTGa                                                                                                                                                   <  4:20280/149‑1 (MQ=255)
     gCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTc                                                                                                                                              <  1:581930/149‑1 (MQ=255)
          atgatTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTCAAGCGTCGCATc                                                                                                                                         <  4:414679/149‑1 (MQ=255)
           tgatTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTCAAGCGTCGCATca                                                                                                                                        <  4:243477/149‑2 (MQ=255)
           tgatTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGGCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTTAGCGTCGCCTcc                                                                                                                                        >  2:569520/1‑149 (MQ=255)
            gagtAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTTTGAAATCAAAAAGTAATTTGCTTTTTGCCTGATGCGACGCTTGAGCGTCGCAttcc                                                                                                                                       >  2:604755/4‑146 (MQ=255)
             atTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCcaa                                                                                                                                      <  3:210120/149‑1 (MQ=255)
              ttAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCcaac                                                                                                                                     >  1:728000/1‑149 (MQ=255)
              ttAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCcaac                                                                                                                                     >  2:483426/1‑149 (MQ=255)
                aGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCcaacaa                                                                                                                                   <  3:529634/149‑1 (MQ=255)
                 gCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAAt                                                                                                                                  <  1:742788/149‑1 (MQ=255)
                   ccGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGa                                                                                                                                <  4:422525/149‑1 (MQ=255)
                      gtggtTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACaa                                                                                                                             >  4:470042/1‑149 (MQ=255)
                              cagcaATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGga                                                                                                                     <  3:702079/149‑1 (MQ=255)
                                      tCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGAtcttct                                                                                                             <  2:584824/149‑1 (MQ=255)
                                        gacgaTTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTttct                                                                                                           >  2:5461/1‑149 (MQ=255)
                                         acgaTTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTg                                                                                                          <  4:546428/149‑1 (MQ=255)
                                         acgaTTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCAACCAACAATGACAAGCGGTGGAGGTCTTCTCTg                                                                                                          >  4:6118/1‑149 (MQ=255)
                                           gaTTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGcc                                                                                                        >  2:580806/1‑149 (MQ=255)
                                                   tGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCtttttt                                                                                                >  3:533828/1‑149 (MQ=255)
                                                    ggCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCttttttt                                                                                               <  1:254326/149‑1 (MQ=255)
                                                     gCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCtttttttt                                                                                              <  4:305050/149‑1 (MQ=255)
                                                     gCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCtttttttt                                                                                              >  1:366234/1‑149 (MQ=255)
                                                      cTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCttttttttt                                                                                             <  2:245030/149‑1 (MQ=255)
                                                      cTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGGAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCttttttttt                                                                                             >  2:37223/1‑149 (MQ=255)
                                                         gATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTtcat                                                                                          >  4:600019/1‑149 (MQ=255)
                                                          aTCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTtcatc                                                                                         >  1:717677/1‑149 (MQ=255)
                                                            cATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCCACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCatcaa                                                                                       >  4:110801/1‑149 (MQ=255)
                                                             aTTAAAGCTCAGGGCCATATTGTGGCTGGTCAAACACCGCTGTATGAAAACAAAAAGTAAATTGCTTTTTGCCCGATGCCGCCCTTGAGCGTCCCATCCAACAATGACAAGCGGTTGAGAGCTTCCCCGCCGCCTTTTTTTTCaacaat                                                                                      >  4:233343/1‑149 (MQ=255)
                                                                       aGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATa                                                                            >  3:716334/1‑149 (MQ=255)
                                                                         ggCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAAt                                                                          <  3:600019/149‑1 (MQ=255)
                                                                             atatTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCtt                                                                      <  4:533828/149‑1 (MQ=255)
                                                                              tatTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCTTCCCCATAATCCTTg                                                                     >  4:709047/1‑149 (MQ=255)
                                                                                  gtgGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCgg                                                                                                                         >  1:704753/1‑93 (MQ=255)
                                                                                  gtgGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCgg                                                                                                                         <  2:704753/93‑1 (MQ=255)
                                                                                       gggTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTAt                                                            >  3:496708/1‑149 (MQ=255)
                                                                                           cAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAAt                                                        <  3:470042/149‑1 (MQ=255)
                                                                                              acacCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAAtttt                                                     >  2:706794/1‑149 (MQ=255)
                                                                                              acacCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAAtttt                                                     >  1:672370/1‑149 (MQ=255)
                                                                                                  cGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTaaaa                                                 >  4:650612/1‑149 (MQ=255)
                                                                                                  cGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTaaaa                                                 >  4:573632/1‑149 (MQ=255)
                                                                                                          gAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACa                                         <  4:144765/149‑1 (MQ=255)
                                                                                                          gAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACa                                         >  3:304832/1‑149 (MQ=255)
                                                                                                                aaaaaGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTc                                   >  2:154595/1‑149 (MQ=255)
                                                                                                                 aaaaGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCa                                  <  3:709047/149‑1 (MQ=255)
                                                                                                                         tCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTcc                          >  3:721226/1‑149 (MQ=255)
                                                                                                                             cTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTCTCCCGCt                      <  2:419146/149‑1 (MQ=255)
                                                                                                                                 ttgttggATGCGACGCTCAAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTtata                  >  2:209581/7‑149 (MQ=255)
                                                                                                                                  tGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACTAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTCTCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATaa                 <  1:604755/149‑1 (MQ=255)
                                                                                                                                  tGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATaa                 <  3:573632/149‑1 (MQ=255)
                                                                                                                                   gCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACaa                                                                                                                             >  3:710421/1‑40 (MQ=255)
                                                                                                                                   gCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACaa                                                                                                                             <  4:710421/40‑1 (MQ=255)
                                                                                                                                     cTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGAt              >  4:58651/1‑149 (MQ=255)
                                                                                                                                         tGCGACGCTCAAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCCTCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCCTa          >  4:686026/1‑149 (MQ=255)
                                                                                                                                            gACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATACGc       >  1:762591/1‑149 (MQ=255)
                                                                                                                                             aCGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATACGcc      >  4:253240/1‑149 (MQ=255)
                                                                                                                                             aCGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATACGcc      >  3:268949/1‑149 (MQ=255)
                                                                                                                                                 ttGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATACGCCGTta  >  4:127412/1‑149 (MQ=255)
                                                                                                                                                 ttGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATACGCCGTta  <  1:37223/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CGCCCGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATACGCCGTTA  >  NC_000913/705770‑706063

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