Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,353,731 T→C Q107R (CAG→CGG)  yhcC ← putative Fe‑S oxidoreductase, Radical SAM superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,353,7310TC100.0% 278.3 / NA 73Q107R (CAG→CGG) yhcCputative Fe‑S oxidoreductase, Radical SAM superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (35/38);  total (35/38)

GGCGGTTTGTAGCCCCAGCTCCAGCCACACTTCGTAACCCTGGTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCTGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTGCTGGTATACGCCTGAAAGTAGGCCAGATAGCGTTTAGCGCGGTTAACTAAATTCGCCTGGTGCGCCAGTTGCTC  >  NC_000913/3353604‑3353858
                                                                                                                               |                                                                                                                               
ggCGGTTTGTAGCCCCAGCTCCAGCCACACTTCGTAACCCTGGTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGcc                                                                                                                      >  2:1481774/1‑139 (MQ=255)
ggCGGTTTGTAGCCCCAGCTCCAGCCACACTTCGTAACCCTGGTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGcc                                                                                                                      >  2:628063/1‑139 (MQ=255)
ggCGGTTTGTAGCCCCAGCTCCAGCCACACTTCGTAACCCTGGTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGcc                                                                                                                      >  1:111353/1‑139 (MQ=255)
     tttGTAGCCCCAGCTCCAGCCACACTTCGTAACCCTGGTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCctgctg                                                                                                                 >  1:674402/1‑139 (MQ=255)
     tttGTAGCCCCAGCTCCAGCCACACTTCGTAACCCTGGTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCctgctg                                                                                                                 <  2:144801/139‑1 (MQ=255)
     tttGTAGCCCCAGCTCCAGCCACACTTCGTAACCCTGGTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCctgctg                                                                                                                 <  2:29217/139‑1 (MQ=255)
     tttGTACCCCCAGGCCCAGCCACACTTCTTAACCCTGCTCCTTATATTCGCAAAGCAGAACCGGCACCGCATCCAGCACGCAGTCCGGGCGGGTACCCACAAACAAACCGACAATATTGGCCCGGCTCACCGCcagctc                                                                                                                 <  2:363262/139‑2 (MQ=255)
       tGTAGCCCCAGCTCCAGCCACACTTCGTAACCCTGGTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGAt                                                                                                               <  2:10680/139‑1 (MQ=255)
               cAGCTCCAGCCACACTTCGTAACCCTGGTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGaa                                                                                                       >  1:1227316/1‑139 (MQ=255)
               cAGCTCCAGCCACACTTCGTAACCCTGGTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGaa                                                                                                       >  2:1444460/1‑139 (MQ=255)
                  cTCCAGCCACACTTCGTAACCCTGGTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACgc                                                                                                    <  2:1004652/139‑1 (MQ=255)
                       gCCACACTTCGTAACCCTGGTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGcc                                                                                                                      >  2:1054767/1‑116 (MQ=255)
                       gCCACACTTCGTAACCCTGGTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGcc                                                                                                                      <  1:1054767/116‑1 (MQ=255)
                            acTTCGTAACCCTGGTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGaa                                                                                          >  1:953495/1‑139 (MQ=255)
                                cGTAACCCTGGTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTc                                                                                      >  2:259081/1‑139 (MQ=255)
                                     cccTGGTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCaa                                                                                 <  2:1061425/139‑1 (MQ=255)
                                      ccTGGTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCaaa                                                                                >  2:1409726/1‑139 (MQ=255)
                                        tGGTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAgg                                                                              <  1:336111/139‑1 (MQ=255)
                                          gTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGcc                                                                                                                      >  2:1463303/1‑97 (MQ=255)
                                          gTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGcc                                                                                                                      <  1:1463303/97‑1 (MQ=255)
                                          gTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCctg                                                                                                                    <  1:404275/99‑1 (MQ=255)
                                          gTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCctg                                                                                                                    >  2:404275/1‑99 (MQ=255)
                                          gTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCaa                                                                                 >  2:590357/1‑134 (MQ=255)
                                          gTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCaa                                                                                 <  1:743179/134‑1 (MQ=255)
                                          gTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCaa                                                                                 <  1:590357/134‑1 (MQ=255)
                                          gTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCaa                                                                                 >  2:743179/1‑134 (MQ=255)
                                          gTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTg                                                                            >  1:1052611/1‑139 (MQ=255)
                                          gTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTg                                                                            >  2:222350/1‑139 (MQ=255)
                                           tCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTGc                                                                           <  2:1057011/139‑1 (MQ=255)
                                                  atTCGCAAATCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTGCTGGtata                                                                    <  2:538794/139‑1 (MQ=255)
                                                  atTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTGCTGGtata                                                                    <  1:1409726/139‑1 (MQ=255)
                                                      gCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTg                                                                            <  1:140747/127‑1 (MQ=255)
                                                      gCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTg                                                                            >  2:140747/1‑127 (MQ=255)
                                                         aaGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTGCTGGTATACGCCTGa                                                             <  2:111353/139‑1 (MQ=255)
                                                           gCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTGCTGGTATACGCCTGaaa                                                           <  2:674402/139‑1 (MQ=255)
                                                            cAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTGCTGGTATACGCCTGAAAg                                                          <  2:205907/139‑1 (MQ=255)
                                                             aGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTGCTGGTATACGCCTGAAAGt                                                         <  1:549002/139‑1 (MQ=255)
                                                                  cAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTGCTGGTATACGCCTGAAAGTAGGcc                                                    <  1:466776/139‑1 (MQ=255)
                                                                    gCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTGCTGGTATACGCCTGAAAGTAGGCCAg                                                  >  2:712577/1‑139 (MQ=255)
                                                                      aCCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTGCTGGTATACGCCTGAAAGTAGGCCAGAt                                                <  1:1444460/139‑1 (MQ=255)
                                                                        cGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGTCAATATTGGCCCGGCTCACCGcc                                                                                                                      >  1:1008498/1‑67 (MQ=255)
                                                                        cGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGcc                                                                                                                      <  2:1008498/67‑1 (MQ=255)
                                                                         gCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCttg                                                                                                                    <  2:1413512/68‑1 (MQ=255)
                                                                         gCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCttg                                                                                                                    >  1:1413512/1‑68 (MQ=255)
                                                                         gCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACgca                                                                                                   <  2:477786/85‑1 (MQ=255)
                                                                         gCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACgca                                                                                                   >  1:477786/1‑85 (MQ=255)
                                                                          cATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTGCTGGTATACGCCTGAAAGTAGGCCAGATAGCg                                            >  1:719942/1‑139 (MQ=255)
                                                                          cATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTGCTGGTATACGCCTGAAAGTAGGCCAGATAGCg                                            >  1:806136/1‑139 (MQ=255)
                                                                                 cACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCctg                                                                                                                    >  2:1147936/1‑60 (MQ=255)
                                                                                 cACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCctg                                                                                                                    <  1:1147936/60‑1 (MQ=255)
                                                                                  aCGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTGCTGGTATACGCCTGAAAGTAGGCCAGATAGCGTTTAgcgc                                    <  1:259081/139‑1 (MQ=255)
                                                                                   cGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTGCTGGTATACGCCTGAAAGTAGGCCAGATAGCGTTTAGcgcg                                   <  1:222350/139‑1 (MQ=255)
                                                                                    gCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTGCTGGTATACGCCTGAAAGTAGGCCAGATAGCGTTTAGCGCgg                                  <  2:771283/139‑1 (MQ=255)
                                                                                      aGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTGCTGGTATACGCCTGAAAGTAGGCCAGATAGCGTTTAGCGCGGtt                                >  2:800709/1‑139 (MQ=255)
                                                                                                  tACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGAt                                                                                                               >  2:1106104/1‑48 (MQ=255)
                                                                                                  tACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGAt                                                                                                               <  1:1106104/48‑1 (MQ=255)
                                                                                                  tACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTGCTGGtata                                                                    >  2:206714/1‑91 (MQ=255)
                                                                                                  tACCCACACACAAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTGCTGGtata                                                                    <  1:206714/91‑1 (MQ=255)
                                                                                                       acacacaAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCaa                                                                                 <  2:127489/73‑1 (MQ=255)
                                                                                                       acacacaAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCaa                                                                                 >  1:127489/1‑73 (MQ=255)
                                                                                                           acaAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTGCTGGTATACGCCTGAAAGTAGGCCAGATAGCGTTTAGCGCGGTTAACTAAATTCg                     >  2:1469480/1‑129 (MQ=255)
                                                                                                           acaAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTGCTGGTATACGCCTGAAAGTAGGCCAGATAGCGTTTAGCGCGGTTAACTAAATTCg                     <  1:1469480/129‑1 (MQ=255)
                                                                                                           acaAACCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTGCTGGTATACGCCTGAAAGTAGGCCAGATAGCGTTTAGCGCGGTTAACTAAATTCGCCTGGTgcgc           >  2:215682/1‑139 (MQ=255)
                                                                                                               aCCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACgc                                                                                                    >  1:1501583/1‑46 (MQ=255)
                                                                                                               aCCGACAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACgc                                                                                                    <  2:1501583/46‑1 (MQ=255)
                                                                                                                   aCAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCa                                                                                  >  2:377262/1‑60 (MQ=255)
                                                                                                                   aCAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCa                                                                                  <  1:377262/60‑1 (MQ=255)
                                                                                                                   aCAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTGCTGGTATACGCCTGAAAGTAGGCCAGATAGCGTTTAGCGCGGTTAACTAAATTCGCCTGGTGCGCCAGTTGCt   <  2:277329/139‑1 (MQ=255)
                                                                                                                   aCAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTGCTGGTATACGCCTGAAAGTAGGCCAGATAGCGTTTAGCGCGGTTAACTAAATTCGCCTGGTGCGCCAGTTGCt   >  2:180443/1‑139 (MQ=255)
                                                                                                                   aCAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTGCTGGTATACGCCTGAAAGTAGGCCAGATAGCGTTTAGCGCGGTTAACTAAATTCGCCTGGTGCGCCAGTTGCt   >  1:277329/1‑139 (MQ=255)
                                                                                                                    cAATATTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTGCTGGTATACGCCTGAAAGTAGGCCAGATAGCGTTTAGCGCGGTTAACTAAATTCGCCTGGTGCGCCAGTTGCTc  <  2:719942/139‑1 (MQ=255)
                                                                                                                       tatTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTGCTGGTATACGCCt                                                               <  2:321907/75‑1 (MQ=255)
                                                                                                                       tatTGGCCCGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTGCTGGTATACGCCt                                                               >  1:321907/1‑75 (MQ=255)
                                                                                                                               |                                                                                                                               
GGCGGTTTGTAGCCCCAGCTCCAGCCACACTTCGTAACCCTGGTCCTTATATTCGCAAAGCAGATCCAGCACCGCATCCGGCACGCAGTCCGGGCGGGTACCCACACACAAACCGACAATATTGGCCTGGCTCACCGCCTGCTGATACATAGAACGCAGCACCTGAACTTCCGCAAAGGTGCTGGTATACGCCTGAAAGTAGGCCAGATAGCGTTTAGCGCGGTTAACTAAATTCGCCTGGTGCGCCAGTTGCTC  >  NC_000913/3353604‑3353858

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