Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 466,148 G→A P389S (CCA→TCA)  ybaE ← uncharacterized protein YbaE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913466,1480GA97.1% 104.1 / NA 35P389S (CCA→TCA) ybaEuncharacterized protein YbaE
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (14/20);  total (15/20)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTCCGGTGCTTCGCCAATTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGACCTGGCTGACCCGTTGCAGCTCCTCCGGTTTGCCGATGGTGATTTGCACGGGATGCCGACAACTGGTTCCCAAATCTTT  >  NC_000913/465905‑466391
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        
atTCCGGTGCTTCGCCAATTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAAt                                                                                                                                                                                                                                                   >  1:399309/1‑251 (MQ=255)
                 attaatctcgtgccaactaggAGGTCTGCGTGTGTCTGTTGGGTCGTGTCTCCCTGGTTTTTTTCGGTATGAAAAATAATTGTGGGTTCCCAGCCTCCTGCTGCCATTGGCGCCTGTAGGCGTTCTGCCATGGTATGAATTTCTATCGGTAGGTGTTAACCCAGCGTCAAGTTTTTCGGTAGTTTTACCTTATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGAt                                                                                                                                                                                                                                 <  2:184612/231‑1 (MQ=255)
                    aaTCTGTCGCCCATCAGGAGGTTTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCa                                                                                                                                                                                                                               <  2:60190/250‑1 (MQ=255)
                           cGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCTTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTGCCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTtc                                                                                                                                                                                                                       <  1:412345/250‑1 (MQ=255)
                            gCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTtctc                                                                                                                                                                                                                     >  1:75391/1‑251 (MQ=255)
                             cccATCATGAGGTCTGCGTGTGCCTGTGGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCAATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTCCCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTcc                                                                                                                                                                                                                    <  2:279355/251‑1 (MQ=255)
                              ccATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGGTCAGGGTTTCTccc                                                                                                                                                                                                                   >  1:514234/1‑250 (MQ=255)
                                        gTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAACGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGGTTTCTCCGACTTCTAACg                                                                                                                                                                                                         >  1:369276/1‑251 (MQ=255)
                                             cGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCCTGGTATGAAGTTCTATCGGTAGGGGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCACTGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGAGATCAGGTTTTTTCAGACTTCTAACGGTTGt                                                                                                                                                                                                    >  2:385999/1‑251 (MQ=255)
                                                         gggTCGTGTAGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAGTAGTCCAGCCTGACAGTAATGCATGACTGGAGGGTATCAGACTTTCTCCGCCTTTTAACCATTGTAATAAACCgggt                                                                                                                                                                                       >  2:562919/1‑249 (MQ=255)
                                                          ggTCGTGTCGCCCTATTTTTTTGCGTTATGAAAAATATTTGGGAGTTCTCATCCTTACGTGGTCATTGTCGCCTGTAGGCGTTTGGCCATGTTTTGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTATTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGAACATGTTTTCGCCGACTTCTAACGTTTGTAATAAACCGGAtt                                                                                                                                                                                      <  2:365090/251‑1 (MQ=255)
                                                                       cctggtttttGCGTTATTAAAAATAATTGTTAGTTCACAGCCTACCACTGCATGTGTCGCGTGTAGGCGATCTGCCATGGTAGTAATGTCAATCGGTAGGAGATAAACTAACTTCAAGGTGTTCGTTAGTTTGACTTCATCCGGTTCCTGCCAATGCTGAATAGTCCAGCCTTACAGTAATGCATGACTGTCGGTGATCAGGTTTTCTCCGAATTCTAACGTTTGTAATAAACCGGATTGATGAATAATgg                                                                                                                                                                          <  1:540186/246‑1 (MQ=255)
                                                                        catttttttttggtATGAAAAATAATTGGTAGTTCCAAGCCTTCCGCTGCCAGTTTCGCCTGTAGGCGTTCTGCCAGTGTATGAAGTTCTTTCGGTAGGGGATAAACCAGCGTCAAGGTTTTCGGTATTTTGATTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGga                                                                                                                                                                         <  2:399309/238‑1 (MQ=255)
                                                                           ttttttGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATgggaat                                                                                                                                                                      >  1:385413/1‑247 (MQ=255)
                                                                              tgtgtgttttgaaaaagaatttggcgTTCACCGCCGTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCAATGGTATGAAGTGCTATCGGTAGGTGATAAACCAACGTCACGGTTTTCGGTAGTTTGACTTCTTCCTTTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATTAACTGGATTGATGAATAATGGAGATCa                                                                                                                                                                    <  2:369276/225‑1 (MQ=255)
                                                                                   gTTATGAGAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGCGGGGGTCAGGGTTTCTACGCCTTCTAACGGTTGTAATAAACCGGATTGTTGGATAATGGGGATCCCttttt                                                                                                                                                              >  2:29589/1‑250 (MQ=255)
                                                                                            aataatTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAATTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGcc                                                                                                                                                     <  1:248378/251‑1 (MQ=255)
                                                                                             ataatTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCa                                                                                                                                                    <  1:5561/251‑1 (MQ=255)
                                                                                             ataatTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGAGATCAGGTTTTCTCCCTCTTTTTACGTTTGTGATAAACCAGGTTAATGAATAATGGAGCACACttttatcgacctccc                                                                                                                                                    >  2:124116/1‑240 (MQ=255)
                                                                                                      aGTTCACAGCCTTCCGCTGCCGGGTTCCCTTGTTGGGGTTCTTCCGTTGGATGAAGTTCTATGGTGTGGTGTTAAACCACGGTCAAGGTTTTAGGTAGTTTGACTTCATCCGGTTCTTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGagag                                                                                                                                           <  2:431221/251‑1 (MQ=255)
                                                                                                         tCACAGCCTTCCGCTGCCCGTGTAGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGGACCGGGTTTTTCCTCCTTCTCACGTGTTTTATAACACCGATTTGTGAATAATTGAGATCACTTTTTGAGCCTGcccggggtagaggcg                                                                                                                                        >  2:362161/1‑238 (MQ=255)
                                                                                                             aGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGa                                                                                                                                    <  1:447472/251‑1 (MQ=255)
                                                                                                              cccTTCCGCGGCCAGTGCCGCCTGTAGGCGTTCTGCCATGGTATGAATTTCTAGCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTCCCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGAc                                                                                                                                   <  2:403633/250‑1 (MQ=255)
                                                                                                                              gtCGCCTGTAGGCTTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGAGAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAAGAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCaa                                                                                                                  <  1:123022/251‑1 (MQ=255)
                                                                                                                                                  cATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGGGATCAGGGTTTTTACCGCTTCTAACGTTTGGTATAAACCGGGTTGACGGATAATTGAGATCACTTTTCGCGCCTTGCAGACGGAGAAGCGGGGGCTTTTGGGCAACATCCCAATAACAAAAACCAAACTgg                                                                                             >  2:70398/1‑250 (MQ=255)
                                                                                                                                                    gctcgtccgtactgtatcgatgggggATAAACCAGCGTCAAGGTTTTCGGTAGTTTGATTGCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCGCCGACTTCTAACGTTTGGAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGGCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGAt                                                                                            <  2:281740/226‑1 (MQ=255)
                                                                                                                                                                cTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGGTCAGGGTTTCTCCGACTTCTCACGTTTTTAATAAACCGGATTGATGAATAATGGTGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGGCTTTTTAGCCACGGCCAATAGGCAAAAACTAAACTGaatagtcgggaaggg                                                                               >  1:281740/1‑237 (MQ=255)
                                                                                                                                                                  aTCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGGCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGAc                                                                                <  1:169280/248‑1 (MQ=255)
                                                                                                                                                                           gTGATAACCCATCGTAAAGGTTGTCGTTAGTTGTACTGCATCCGGAACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGACCTGGCTGAccc                                                                     <  2:514234/251‑1 (MQ=255)
                                                                                                                                                                             taaaaaCCAGTGCCAAGGTTTTGTGTAGTTGAACTTCACCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGaata                                                                                                                                                                              <  2:507730/143‑1 (MQ=255)
                                                                                                                                                                             gATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAAACGGATTGATGaata                                                                                                                                                                              >  1:507730/1‑146 (MQ=255)
                                                                                                                                                                                        gTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTACCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGTTCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGGTTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCCCAACGGCAAATAACAAAAACCTAAACTGATGCCGCTACTGACCCTGCTGACCCGTTGCAGcccccc                                                        >  1:474681/1‑246 (MQ=255)
                                                                                                                                                                                                 ttCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGACCTGGCTGACCCGTTGCAGCTCCTCCGGTTTGc                                                <  1:29589/250‑1 (MQ=255)
                                                                                                                                                                                                                                         aaGTCCAGCTTGACAGTAATGCATGACTGGCGGTGAACAGGTTTGCCCCGACTTTTAACGTTTGTAATAAACCGGATTGATGAAGAATGGAGAGCACGTTTCGCGCCTGCCAGAGGGAGAGGCGGGGACTTTTGCGCATCGTCAAATAGTAAAAACCTAAACTGATGCCGCTACTGACCTGGCTGTCCCGTTGCAGCTCCTCCGGTTTGCCGATGGTGATTTGCACGGGATGCCGACAACTGGTTCCCaaa       <  2:463637/250‑1 (MQ=255)
                                                                                                                                                                                                                                              cAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGGTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGGTGACGCTACTGACCTTGCTGCACCGTTTCACCTCCTCCGGTTTGCCGGCGGTGATTTGCACGGGAAGCCCGCAACAGGTTCCCCAAACttt  >  1:412980/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        
ATTCCGGTGCTTCGCCAATTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGACCTGGCTGACCCGTTGCAGCTCCTCCGGTTTGCCGATGGTGATTTGCACGGGATGCCGACAACTGGTTCCCAAATCTTT  >  NC_000913/465905‑466391

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