Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,385,274:1 +G 8.8% intergenic (+19/‑44) lldP → / → lldD L‑lactate permease/L‑lactate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,385,2741.G 8.8% 212.4 / 15.8 57intergenic (+19/‑44)lldP/lldDL‑lactate permease/L‑lactate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (25/27);  new base G (2/3);  total (27/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 = 9.19e-01

CGGCCTGGTCGGCAAGGAATCCGACCTGTTCCGCTTCACCGTCAAGCACAGCCTGTTCTTCGCCACCATCGTCGGCCTGATCACCCTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGCCTCTACCGACTATCGCGCCGCGGCCCAACGCAAGCTGCCCCCCTTCCTGTTCCACTACGCTGACGGCGGCGCCTACGCCGAGCACAC  >  NC_002947/5385126‑5385418
                                                                                                                                                     |                                                                                                                                                
cGGCCTGGTCGGCAAGGAATCCGACCTGTTCCGCTTCACCGTCAAGCACAGCCTGTTCTTCGCCACCATCGTCGGCCTGATCACCCTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGcc                                                                                                                                                 >  2:330911/1‑150 (MQ=255)
   ccTGGTCGGCAAGGAATCCGACCTGTTCCGCTTCACCGTCAAGCACAGTCTGTTCTTCGCCACCATCGTCGGCCTGATCACCCTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGgcgcgcga                                                                                                                                               >  1:342683/1‑146 (MQ=255)
       gTCGGCAAGGAATCCGACCTGTTCCGCTTCACCGTCAAGCACAGCCTGTTCTTCGCCACCATCGTCGGCCTGATCACCCTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCTGGCGCCGTcgcgc                                                                                                                                          <  1:1257499/150‑1 (MQ=255)
          ggCAAGGAATCCGACCTGTTCCGCTTCACCGTCAAGCACAGCCTGTTCTTCGCCACCATCGTCGGCCTGATCACCCTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGGCGGGGGCCTTCGCGCGcc                                                                                                                                       >  2:277956/1‑150 (MQ=255)
             aaGGAATCCGACCTGTTCCGCTTCACCGTCAAGCACAGCCTGTTCTTCGCCACCATCGTCGGCCTGATCACCCTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCgg                                                                                                                                    >  1:334985/1‑150 (MQ=255)
             aaGGAATCCGACCTGTTCCGCTTCACCGTCAAGCACAGCCTGTTCTTCGCCACCATCGTCGGCCTGATCACCCTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCgg                                                                                                                                    >  2:580604/1‑150 (MQ=255)
                gAATCCGACCTGTTCCGCTTCACCGTCAAGCACAGCCTGTTCTTCGCCACCATCGTCGGCCTGATCACCCTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCAc                                                                                                                                 >  1:430723/1‑150 (MQ=255)
                  aTCCGACCTGTTCCGCTTCACCGTCAAGCACAGTCTGTTCTTCGCCACCATCGTCGGCCTGATCACCCTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCGCGACGGCGCCCGGcccct                                                                                                                               <  2:342683/150‑5 (MQ=255)
                       aCCTGTTCCGCTTCACCGTCAAGCACAGCCTGTTCTTCGCCACCATCGTCGGCCTGATCACCCTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCcaca                                                                                                                          <  2:1027029/150‑1 (MQ=255)
                           gTTCCGCTTCACCGTCAAGCACAGCCTGTTCTTCGCCACCATCGTCGGCCTGATCACCCTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCGCGACGGCGCCCGGcccctttc                                                                                                                            <  1:1272037/144‑8 (MQ=255)
                           gTTCCGCTTCACCGTCAAGCACAGCCTGTTCTTCGCCACCATCGTCGGCCTGATCACCCTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCGCGACGGCGCCCGGcccctttc                                                                                                                            >  2:1272037/1‑137 (MQ=255)
                           gTTCCGCTTCACCGTCAAGCACAGCCTGTTCTTCGCCACCATCGTCGGCCTGATCACCCTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCa                                                                                                                      >  2:858483/1‑150 (MQ=255)
                                  ttCACCGTCAAGCACAGCCTGTTCTTCGCCACCATCGTCGGCCTGATCACCCTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGc                                                                                                               <  1:226243/150‑1 (MQ=255)
                                    cACCGTCAAGCACAGCCTGTTCTTCGCCACCATCGTCGGCCTGATCACCCTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCga                                                                                                             <  2:253066/150‑1 (MQ=255)
                                        gTCAAGCACAGCCTGTTCTTCGCCACCATCGTCGGCCTGATCACCCTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAAt                                                                                                         >  2:770182/1‑150 (MQ=255)
                                             gCACAGCCTGTTCTTCGCCACCATCGTCGGCCTGATCACCCTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCGCGACGGCGCCCGGc                                                                                                                                   <  1:17653/119‑1 (MQ=255)
                                             gCACAGCCTGTTCTTCGCCACCATCGTCGGCCTGATCACCCTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCGCGACGGCGCCCGGc                                                                                                                                   >  2:17653/1‑119 (MQ=255)
                                                  gCCTGTTCTTCGCCACCATCGTCGGCCTGATCACCCTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCACGCTGCGAGAATCCATGATCAtt                                                                                              <  2:665218/150‑1 (MQ=255)
                                                   ccTGTTCTTCGCCACCATCGTCGGCCTGATCACCCTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCAtt                                                                                              >  2:420742/1‑150 (MQ=255)
                                                   ccTGTTCTTCGCCACCATCGTCGGCCTGATCACCCTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCAtt                                                                                              >  1:775231/1‑150 (MQ=255)
                                                        tcttgGCCACCATCGTCGGCCTGATCACCCTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGc                                                                                         <  1:218285/150‑1 (MQ=255)
                                                        tcttcGCCACCATCGTCGGCCTGATCACCCTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGc                                                                                         <  1:556729/150‑1 (MQ=255)
                                                        tcttcGCCACCATCGTCGGCCTGATCACCCTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGc                                                                                         <  1:330911/150‑1 (MQ=255)
                                                              ccaccaTCGTCGGCCTGATCACCCTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGCCTCTAc                                                                                   >  1:1050467/1‑150 (MQ=255)
                                                               caccaTCGTCGGCCTGATCACCCTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGCCTCTAcc                                                                                  >  2:345975/1‑150 (MQ=255)
                                                               caccaTCGTCGGCCTGATCACCCTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGCCTCTAcc                                                                                  <  2:1050467/150‑1 (MQ=255)
                                                                  caTCGTCGGCCTGATCACCCTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATTGATGATCATTTCTGCCTCTACCGAc                                                                               <  1:495107/150‑1 (MQ=255)
                                                                          gCCTGATCACCCTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGCCTCTACCGACTATCgcgc                                                                       >  1:468575/1‑150 (MQ=255)
                                                                                tCACCCTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGCCTCTACCGACTATCGCGCCGCGGc                                                                 <  2:1064107/150‑1 (MQ=255)
                                                                                    ccTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTcc                                                                                                                              <  1:134630/85‑1 (MQ=255)
                                                                                    ccTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTcc                                                                                                                              >  2:134630/1‑85 (MQ=255)
                                                                                     cTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGCCTCTACCGACTATCGCGCCGCGGCCCAAc                                                            <  2:430723/150‑1 (MQ=255)
                                                                                     cTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGCCTCTACCGACTATCGCGCCGCGGCCCAAc                                                            <  1:420742/150‑1 (MQ=255)
                                                                                         tCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGCCTCTACCGACTATCGCGCCGCGGCCCAACGCaa                                                        <  1:1088303/150‑1 (MQ=255)
                                                                                         tCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCAAGATCATTTCTGCCTCTACCGACTATCGCGCCGCGGCCCAACGCaa                                                        >  1:771945/1‑150 (MQ=255)
                                                                                            aGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGCCTCTACCGACTATCGCGCCGCGGCCCAACGCAAg                                                       <  2:385575/148‑1 (MQ=255)
                                                                                            aGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGCCTCTACCGACTATCGCGCCGCGGCCCAACGCAAg                                                       >  1:385575/1‑148 (MQ=255)
                                                                                              gCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGCCTCTACCGACTATCGCGCCGCGGCCCAACGCAAGCTGc                                                   >  1:971608/1‑150 (MQ=255)
                                                                                                ttACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGCCTCTACCGACTATCGCGCCGCGGCCCAACGCAAGCTGccc                                                 <  2:558827/150‑1 (MQ=255)
                                                                                                ttACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGCCTCTACCGACTATCGCGCCGCGGCCCAACGCAAGCTGccc                                                 <  2:334985/150‑1 (MQ=255)
                                                                                                          aCCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCGCTCCCACAAGCAACGCTGCGAGAATCCa                                                                                                      <  2:405344/87‑1 (MQ=255)
                                                                                                          aCCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCa                                                                                                      >  1:405344/1‑87 (MQ=255)
                                                                                                            cGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGCCTCTACCGACTATCGCGCCGCGGCCCAACGCAAGCTGCCCCCCTTCCTGTTc                                     >  1:966421/1‑150 (MQ=255)
                                                                                                              gCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGAGCGGCTTGCTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGCCTCTACCGACTATCGCGCCGCGGCCCAACGCAAGCTGCCCCCCTTCCTGTTCCa                                   <  1:719840/150‑1 (MQ=255)
                                                                                                               cATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGCCTCTACCGACTATCgcgc                                                                       <  2:154537/113‑1 (MQ=255)
                                                                                                               cATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGCCTCTACCGACTATCgcgc                                                                       >  1:154537/1‑113 (MQ=255)
                                                                                                                aTGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGCCTCTACCGACTATCGCGCCGCGGCCCAACGCAAGCTGCCCCCCTTCCTGTTCCACt                                 >  2:628791/1‑150 (MQ=255)
                                                                                                                  gCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGCCTCTACCGACTATCGCGCCGCGGCCCAACGCAAGCTGCCCCCCTTCCTGTTCCACTAc                               <  1:628791/150‑1 (MQ=255)
                                                                                                                    tGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGAGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGCCTCTACCGACTATCGCGCCGCGGCCCAACGCAAGCTGCCCCCCTTCCTGTTCCACTACGc                             >  1:1222585/1‑150 (MQ=255)
                                                                                                                       ttCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGCCTCTACCGACTATCGCGCCGCGGCCCAACGCAAGCTGCCCCCTTTCCTGTTCCACTACGCTGa                          >  2:1144435/1‑150 (MQ=255)
                                                                                                                               tAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCTGCACTCCCACAAACACCACTGCGAGAATCCATGATCATTTCTGCCTCTACCGACTATCGCGCCGCGGCCCAACGCAAGCTGCCCCCCTTCCTGTTCCACTACGCTGACGgcggcg                  >  1:157732/1‑150 (MQ=255)
                                                                                                                                 aaGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGCCTCTACCGACTATCGCGCCGCGGCCCAACGCAAGCTGCCCCCCTTcc                                          <  1:1001685/124‑1 (MQ=255)
                                                                                                                                 aaGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGCCTCTACCGACTATCGCGCCGCGGCCCAACGCAAGCTGCCCCCCTTcc                                          >  2:1001685/1‑124 (MQ=255)
                                                                                                                                   gTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGCCTCTACCGACTATCGCGCCGCGGCCCAACGCAAGCTGCCCCCCTTCCTGTTCCACTACGCTGACGGCGGCGCCTa              <  1:345975/150‑1 (MQ=255)
                                                                                                                                         ggggCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGCCTCTACCGACTATCGCGCCGCGGCCCAACGCAAGCTGCCCCCCTTCCTGTTCCACTACGCTGACGGCGGCGCCTACGCCGa        <  2:775231/150‑1 (MQ=255)
                                                                                                                                          gggCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGCCTCTACCGACTATCGCGCCGCGGCCCAACGCAAGCTGCCCCCCTTCCTGTTCCACTACGCTGACGGCGGCGCCTACGCCGAg       >  1:184412/1‑150 (MQ=255)
                                                                                                                                            gCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGCCTCTACCGa                                                                                >  2:561224/1‑75 (MQ=255)
                                                                                                                                            gCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGCCTCTACCGa                                                                                <  1:561224/75‑1 (MQ=255)
                                                                                                                                               gggCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGCCTCTACCGACTATCGCGCCGCGGCCCAACGCAAGCTGCCCCCCTTCCTGTTCCACTACGCTGACGGCGGCGCCTACGCCGAGCacac  <  1:770182/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
CGGCCTGGTCGGCAAGGAATCCGACCTGTTCCGCTTCACCGTCAAGCACAGCCTGTTCTTCGCCACCATCGTCGGCCTGATCACCCTGGTCCAGGCTTACTGGCTGACCGGCATGCTGGTTCATCACTAAAGTGAAAGGGGCCGGGCGCCGTCGCGCGCCCGGCACTCCCACAAGCAACGCTGCGAGAATCCATGATCATTTCTGCCTCTACCGACTATCGCGCCGCGGCCCAACGCAAGCTGCCCCCCTTCCTGTTCCACTACGCTGACGGCGGCGCCTACGCCGAGCACAC  >  NC_002947/5385126‑5385418

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: