Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 807,696 C→T D72N (GAT→AAT)  ybhB ← kinase inhibitor homolog, UPF0098 family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913807,6960CT100.0% 221.7 / NA 57D72N (GAT→AAT) ybhBkinase inhibitor homolog, UPF0098 family
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (31/26);  total (31/26)

GAACGGTAAAAATGTAGCGATGAGTTTCGCCTTTCGGCGGTGCTGCGCCATCGTACCCGGTTTTACCAAAGTCGGTACGCGTCTGCAAAACGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATCAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTAGCAGGTGACAACAAAACTTTTCGTTCCCGCAGGAACATCATCCCACGCCAGATGCGGTGAAATATTATCGCCATCGTAACCCAT  >  NC_000913/807549‑807843
                                                                                                                                                   |                                                                                                                                                   
gAACGGTAAAAATGTAGCGATGAGTTTCGCCTTTCGGCGGTGCTGCGCCATCGTACCCGGTTTTACCAAAGTCGGTACGCGTCTGCAAAACGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGc                                                                                                                                                  <  1:1070133/151‑1 (MQ=255)
 aaCGGTAAAAATGTAGCGATGAGTTTCGCCTTTCGGCGGTGCTGCGCCATCGTACCCGGTTTTACCAAAGTCGGTACGCGTCTGCAAAACGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCg                                                                                                                                                 >  1:632783/1‑151 (MQ=255)
    ggTAAAAATGTAGCGATGAGTTTCGCCTTTCGGCGGCGCTTCGCCATCGTACCCGGTTTTACCAAAGCCGGAACGCCACTGCAAAACGCCGTCCGGCATTGCTACCAAACCAGAGCCAAAACCCTGCGGTAATACGCGGGTATTAGCgg                                                                                                                                                >  2:56191/1‑149 (MQ=255)
     gTAAAAATGTAGCGATGAGTTTCGCCTTTCGGCGGTGCTGCGCCATCGTACCCGGTTTTACCAAAGTCGGTACGCGTCTGCAAAACGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCggg                                                                                                                                               >  2:930763/1‑149 (MQ=255)
     gTAAAAATGTAGCGATGAGTTTCGCCTTTCGGCGGTGCTGCGCCATCGTACCCGGTTTTACCAAAGTCGGTACGCGTCTGCAAAACGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTa                                                                                                                                             <  1:1050158/151‑1 (MQ=255)
       aaaaaTGTAGCGATGAGTTTCGCCTTTCGGCGGTGCTGCGCCATCGTACCCGGTTTTACCAAAGTCGGTACGCGTCTGCAAAACGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTaa                                                                                                                                            <  2:689736/150‑1 (MQ=255)
          aaTGTAGCGATGAGTTTCGCCTTTCGGCGGTGCTGCGCCATCGTACCCGGTTTTACCAAAGTCGGTACGCGTCTGCAAAACGCCGTCCGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATAAGCGGGGATTAGCGGGTAAGtt                                                                                                                                         >  1:1108917/1‑150 (MQ=255)
          aaTGTAGCGATGAGTTTCGCCTTTCGGCGGTGCTGCGCCATAGTACCCGGTTTTACCAAAGTCGGTACGCGTCTGCAAAACGCCGTCTGGCATTGCTACCCGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAACGGGTAAgat                                                                                                                                         >  1:803150/1‑148 (MQ=255)
                gCGATGAGTTTCGCCTTTCGGCGGTGCTGCGCCATCGTACCCGGTTTTACCAAAGTCGGTACGCGTCTGCAAAACGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACt                                                                                                                                  >  2:365755/1‑151 (MQ=255)
                   aTGAGTTTCGCCTTTCGGCGGTGCTGCGCCATCGTACCCGGTTTTACCAAAGTCGGTACGCGTCTGCAAAACGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTAcc                                                                                                                               <  2:360487/151‑1 (MQ=255)
                     gAGTTTCGCCTTTCGGCGGTGCTGCGCCATCGTACCCGGTTTTACCAAAGTCGGTACGCGTCTGCAAAACGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTAccc                                                                                                                              <  2:990076/150‑1 (MQ=255)
                     gAGTTTCGCCTTTCGGCGGTGCTGCGCCATCGTACCCGGTTTTACCAAAGTCGGTACGCGTCTGCAAAACGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCa                                                                                                                             <  1:64404/151‑1 (MQ=255)
                          tCGCCTTTCGGCGGTGCTGCGCCATCGTACCCGGTTTTACCAAAGTCGGTACGCGTCTGCAAAACGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGGATTAGCGGGTAAGTTAACAACTACCCAGTGcc                                                                                                                        >  2:757758/1‑151 (MQ=255)
                           cGCCTTTCGGCGGTGCTGCGCCATCGTACCCGGTTTTACCAAAGTCGGTACGCGTCTTCAAAACGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGGATTAGCGGGTAAGTTAACAAATACCCAGTGcc                                                                                                                        >  1:208594/1‑150 (MQ=255)
                           cGCCTTTCGGCGGTGCTGCGCCATCGTACCCGGTTTTACCAAAGTCGGTACGCGTCTGCAAAACGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGcca                                                                                                                       <  1:722930/151‑1 (MQ=255)
                            gCCTTTCGGCGGTGCCGCGCCATCGTACCCGGTTTTACCCAAGACGGTACGCGTCTTCAAAACGCCGTCTGGCATTGCTACCAGACCAAAGCCAAACCCCTGCGGTACTACGCGGGTTTTCACGGGCAAGTTAACAACTAACCAGGGcc                                                                                                                        >  1:785155/1‑149 (MQ=255)
                               tttCGGCGGTGCTGCGCCATCGTACCCGGTTTTACCAAAGTCGGTACGCGTCTGCAAAACGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGccacc                                                                                                                     >  1:838354/1‑149 (MQ=255)
                                ttCGGCGGTGCTGCGCCATCGTACCCGGTTTTACCAAAGTCGGTACGCGTCTGCAAAACGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGccacca                                                                                                                    <  1:873049/149‑1 (MQ=255)
                                    gcggTGCTGCGCCATCGTACCCGGTTTTACCAAAGTCGGTACGCGTCTGCAAAACGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCagccgga                                                                                                              <  2:151086/151‑1 (MQ=255)
                                            gcgcCATCGTACCCGGTTTTACCAAAGTCGGTACGCGTCTGCAAAACGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTg                                                                                                      >  1:476627/1‑151 (MQ=255)
                                              gcCATCGTACCCGGTTTTACCAAAGTCGGTACGCGTCTGCAAAACGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGgc                                                                                                    >  1:157993/1‑151 (MQ=255)
                                                   cGTACCCGGTTTTACCAAAGTCGGTACGCGTCTGCAAAACGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATc                                                                                               >  2:756686/1‑151 (MQ=255)
                                                        ccGGTTTTACCAAAGTCGGTACGCGTCTGCAAAACGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGt                                                                                          <  1:1167948/151‑1 (MQ=255)
                                                          ggTTTTACCAAAGTCGGTACGCGTCTGCAAAACGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCg                                                                                        <  1:477351/151‑1 (MQ=255)
                                                           gTTTTACCAAAGTCGGTACGCGTCTGCAAAACGCCGTCTGGCATTGCTACCAGACCATAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGt                                                                                       >  2:764533/1‑151 (MQ=255)
                                                             tttACCAAAGTCGGTACGCGTCTGCAAAACGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTa                                                                                      <  1:883796/150‑1 (MQ=255)
                                                             tttACCAAAGTCGGTACGCGTCTGCAAAACGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTa                                                                                      <  2:94022/150‑1 (MQ=255)
                                                                  cAAAGTCGGTACGCGTCTGCAAAACGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTAGCAgg                                                                                 >  2:896652/1‑150 (MQ=255)
                                                                   aaaGTCGGTACGCGTCTGCAAAACGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTAGCAGGTg                                                                               <  1:753121/151‑1 (MQ=255)
                                                                   aaaGTCGGTACGCGTCTGCAAAACGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTAGCAGGTg                                                                               <  2:300618/151‑1 (MQ=255)
                                                                    aaGTCGGTACGCGTCTGCAAAACGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTAGCAGGTGa                                                                              <  1:701341/151‑1 (MQ=255)
                                                                      gTCGGTACGCGTCTGCAAAACGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTAGCAGGTGaca                                                                            <  1:979004/151‑1 (MQ=255)
                                                                            aCGCGTCTGCAAAACGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTAGCAGGTGACAACaaa                                                                       >  1:909345/1‑150 (MQ=255)
                                                                                   tGCAAAACGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTAGCAGGTGACAACAAAACTTTTCg                                                               >  1:595192/1‑151 (MQ=255)
                                                                                    gCAAAACGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTAGCAGGTGACAACAAAACTTTTCg                                                               <  2:461624/150‑1 (MQ=255)
                                                                                    gCAAAACGCCGTCTGGCATGGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTAGCAGGTGACAACAAAACtttt                                                                 <  1:602388/148‑1 (MQ=255)
                                                                                         aCGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTAGCAGGTGGCAACAAAACTTTTCGTTccc                                                          >  1:270840/1‑150 (MQ=255)
                                                                                          cGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTAGCAGGTGACAACAAAACTTTTCGTTCCCGc                                                        >  2:62185/1‑151 (MQ=255)
                                                                                               tCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTAGCAGGTGACAACAAAACTTTTCGTTCCCGCAGGa                                                    >  2:1127053/1‑150 (MQ=255)
                                                                                               tCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGACACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTAGCAGGTGACAACAAAACTTTTCGTTCCCGCAGGa                                                    >  1:821403/1‑150 (MQ=255)
                                                                                                         cTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTAGCAGGTGACAACAAAACTTTTCGTTCCCGCAGGAACATCATccc                                          >  1:155086/1‑150 (MQ=255)
                                                                                                          tACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTAGCAGGTGACAACAAAACTTTTCGTTCCCGCAGGAACATCATCCCAc                                        >  2:110479/1‑151 (MQ=255)
                                                                                                                 ccagaGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTAGCAGGTGACAACAAAACTTTTCGTTCCCGCAGGAACATCATCCCACGcc                                     >  1:1099338/1‑147 (MQ=255)
                                                                                                                 ccagaGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTAGCAGGTGACAACAAAACTTTTCGTTCCCGCAGGAACATCATCCCACGCCAGAt                                 <  2:907773/151‑1 (MQ=255)
                                                                                                                 ccagaGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTAGCAGGTGACAACAAAACTTTTCGTTCCCGCAGGAACATCATCCCACGCCAGAt                                 <  1:778111/151‑1 (MQ=255)
                                                                                                                  cagaGCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTAGCAGGTGACAACAAAACTTTTCGTTCCCGCAGGAACATCATCCCACGCCAGATg                                >  1:786580/1‑151 (MQ=255)
                                                                                                                   agagCCAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTAGCAGGTGACAACAAAACTTTTCGTTCCCGCAGGAACATCATCCCACGCCAGATGc                               >  2:619499/1‑151 (MQ=255)
                                                                                                                       ccAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTAGCAGGTGACAACAAAACTTTTCGTTCCCGCAGGAACATCATCCCACGCCAGATGCGGTg                           >  2:368474/1‑151 (MQ=255)
                                                                                                                        cAAACCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTAGCAGGTGACAACAAAACTTTTCGTTCCCGCAGGAACATCATCCCACGCCAGATGCGGTg                           >  1:1139604/1‑150 (MQ=255)
                                                                                                                         aaaCCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTAGCAGGTGACAACAAAACTTTTCGTTCCCGCAGGAACATCATCCCACGCCAGATGCGGTGaa                         <  2:776978/151‑1 (MQ=255)
                                                                                                                          aaCCCTTGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTAGCAGGTGACAACAAAACTTTTCGTTCCCGCAGGAACATCATCCCACGCCAGATGCGGTGaaa                        <  2:674852/151‑1 (MQ=255)
                                                                                                                                tGCGGTAATACGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTAGCAGGTGACAACAAAACTTTTCGTTCCCGCAGGAACATCATCCCACGCCAGATGCGGTGAAAtatt                    <  1:505371/149‑1 (MQ=255)
                                                                                                                                         aCGCGGGTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTAGCAGGTGACAACAAAACTTTTCGTTCCCGCAGGAACATCATCCCACGCCAGATGCGGTGAAATATTATCGCCATCg          <  1:971735/150‑1 (MQ=255)
                                                                                                                                             gggTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTAGCAGGTGACAACAAAACTTTTCGTTCCCGCAGGAACATCATCCCACGCCAGATGCGGTGAAATATTATCGCCATCGTAAcc     >  1:977039/1‑151 (MQ=255)
                                                                                                                                              ggTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTAGCAGGTGACAACAAAACTTTTCGTTCCCGCAGGAACATCATCCCACGCCAGATGCGGTGAAATATTATCGCCATCGTAAccc    <  2:267844/151‑1 (MQ=255)
                                                                                                                                               gTATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTAGCAGGTGACAACAAAACTTTTCGTTCCCGCAGGAACATCATCCCACGCCAGATGCGGTGAAATATTATCGCCATCGTAAccc    >  2:1181750/1‑150 (MQ=255)
                                                                                                                                                tATTAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTAGCAGGTGACAACAAAACTTTTCGTTCCCGCAGGAACATCATCCCACGCCAGATGCGGTGAAATATTATCGCCATCGTAACCCAt  >  2:378139/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GAACGGTAAAAATGTAGCGATGAGTTTCGCCTTTCGGCGGTGCTGCGCCATCGTACCCGGTTTTACCAAAGTCGGTACGCGTCTGCAAAACGCCGTCTGGCATTGCTACCAGACCAGAGCCAAACCCTTGCGGTAATACGCGGGTATCAGCGGGTAAGTTAACAACTACCCAGTGCCACCAGCCGGAGCCGGTTGGCGCATCCGGGTCGTAGCAGGTGACAACAAAACTTTTCGTTCCCGCAGGAACATCATCCCACGCCAGATGCGGTGAAATATTATCGCCATCGTAACCCAT  >  NC_000913/807549‑807843

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