Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 913,349 G→A L156L (CTG→TTG)  hcp ← hybrid‑cluster [4Fe‑2S‑2O] subunit of anaerobic terminal reductases

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913913,3490GA100.0% 231.2 / NA 71L156L (CTG→TTG) hcphybrid‑cluster [4Fe‑2S‑2O] subunit of anaerobic terminal reductases
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (37/34);  total (37/34)

AGAAGCGCGTTCATATCGGCAGGCCAGGTCCCCAGCCACGCCATGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGTCAGCCATCGGGTTA  >  NC_000913/913213‑913479
                                                                                                                                        |                                                                                                                                  
agaagCGCGTTCATATCGGCAGGCCAGGTCCCCAGCCACGCCATGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGc                                                                                                                                  <  1:971423/139‑1 (MQ=255)
       cgTTCATATCGGCAGGCCAGGTCCCCAGCCACGCCATGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCAcagcag                                                                                                                           >  2:485955/1‑139 (MQ=255)
         ttCATATCGGCAGGCCAGGTCCCCAGCCACGCCATGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGAc                                                                                                                         >  1:1054410/1‑139 (MQ=255)
               tCGGCAGGCCAGGTCCCCAGCCACGCCATGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGc                                                                                                                   <  2:470999/139‑1 (MQ=255)
                    aGGCCAGGTCCCCAGCCACGCCATGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCgaga                                                                                                              <  1:1660329/139‑1 (MQ=255)
                      gCCAGGTCCCCAGCCACGCCATGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAAt                                                                                                            >  2:892227/1‑139 (MQ=255)
                                cAGCCACGCCATGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGc                                                                                                                                  <  2:888125/107‑1 (MQ=255)
                                     aCGCCATGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATcgc                                                                                             >  1:1848054/1‑139 (MQ=255)
                                              ttttATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATc                                                                                    <  1:892227/139‑1 (MQ=255)
                                                 tATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGca                                                                                                                                 >  1:1293386/1‑91 (MQ=255)
                                                 tATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGca                                                                                                                                 <  2:1293386/91‑1 (MQ=255)
                                                 tATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTa                                                                                      >  1:362712/1‑134 (MQ=255)
                                                 tATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTa                                                                                      <  2:362712/134‑1 (MQ=255)
                                                     gTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTa                                                                             <  1:1862464/139‑1 (MQ=255)
                                                     gTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTa                                                                             <  1:485955/139‑1 (MQ=255)
                                                        cTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATc                                                                                               >  1:329265/1‑118 (MQ=255)
                                                        cTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATc                                                                                               <  2:329265/118‑1 (MQ=255)
                                                        cTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTa                                                                             >  1:1811493/1‑136 (MQ=255)
                                                        cTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTa                                                                             <  2:1811493/136‑1 (MQ=255)
                                                         tGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGa                                                                                                                          <  2:676783/90‑1 (MQ=255)
                                                         tGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGa                                                                                                                          >  1:676783/1‑90 (MQ=255)
                                                                aaaTATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTa                                                                             >  1:1521570/1‑128 (MQ=255)
                                                                aaaTATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTa                                                                             <  2:1521570/128‑1 (MQ=255)
                                                                     tCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTctgctg                                                             >  2:1038195/1‑139 (MQ=255)
                                                                       gTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGtt                                                                              <  1:351793/120‑1 (MQ=255)
                                                                       gTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGtt                                                                              >  2:351793/1‑120 (MQ=255)
                                                                       gTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGcc                                                           >  2:724870/1‑139 (MQ=255)
                                                                         tGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTg                                                         <  2:1848054/139‑1 (MQ=255)
                                                                          gTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGa                                                        >  2:666404/1‑139 (MQ=255)
                                                                           tCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCt                                                                                         >  2:911382/1‑105 (MQ=255)
                                                                           tCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCt                                                                                         <  1:911382/105‑1 (MQ=255)
                                                                                  gACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATc                                                                                               <  2:1574064/92‑1 (MQ=255)
                                                                                  gACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATc                                                                                               >  1:1574064/1‑92 (MQ=255)
                                                                                  gACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAg                                                <  2:1054410/139‑1 (MQ=255)
                                                                                   aCCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGc                                               >  2:1328710/1‑139 (MQ=255)
                                                                                        gaACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGt                                                                               <  2:1450282/102‑1 (MQ=255)
                                                                                        gaACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGt                                                                               >  1:1450282/1‑102 (MQ=255)
                                                                                             tGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGAt                                     >  2:1684158/1‑139 (MQ=255)
                                                                                             tGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAACACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGAt                                     >  1:1334014/1‑139 (MQ=255)
                                                                                                 cgTGTTCCATATAGGCTGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGcc                                                                                                                  <  1:1016358/58‑1 (MQ=255)
                                                                                                 cgTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGcc                                                                                                                  >  2:1016358/1‑58 (MQ=255)
                                                                                                 cgTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGAtcgtc                                 >  1:864322/1‑139 (MQ=255)
                                                                                                   tgtTCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGc                               >  2:965123/1‑139 (MQ=255)
                                                                                                     ttCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGc                               >  2:365421/1‑137 (MQ=255)
                                                                                                     ttCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGc                               <  1:365421/137‑1 (MQ=255)
                                                                                                     ttCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGc                               >  2:1137626/1‑137 (MQ=255)
                                                                                                     ttCCATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGc                               <  1:1137626/137‑1 (MQ=255)
                                                                                                       ccATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCttt                                                                                       >  1:639576/1‑79 (MQ=255)
                                                                                                       ccATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCttt                                                                                       <  2:639576/79‑1 (MQ=255)
                                                                                                        cATATAGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCAcc                          <  1:666404/139‑1 (MQ=255)
                                                                                                          tataGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGtt                                                                              >  2:1778701/1‑85 (MQ=255)
                                                                                                          tataGGCCGCCGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGtt                                                                              <  1:1778701/85‑1 (MQ=255)
                                                                                                              ggCCGCCGCACCTTTCAGACCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAg                                                >  2:1621788/1‑111 (MQ=255)
                                                                                                              ggCCGCCGCACCTTTCAGACCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAg                                                <  1:1621788/111‑1 (MQ=255)
                                                                                                               gccgccGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCa                   >  2:1062712/1‑139 (MQ=255)
                                                                                                               gccgccGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCa                   >  1:146963/1‑139 (MQ=255)
                                                                                                               gccgccGCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCa                   >  1:305123/1‑139 (MQ=255)
                                                                                                                ccgccgCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGt                                                                        >  2:366645/1‑85 (MQ=255)
                                                                                                                ccgccgCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGt                                                                        <  1:366645/85‑1 (MQ=255)
                                                                                                                ccgccgCACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAg                  >  2:654669/1‑139 (MQ=255)
                                                                                                                    cgcACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATcg                                                                                              <  2:542413/59‑1 (MQ=255)
                                                                                                                    cgcACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATcg                                                                                              >  1:542413/1‑59 (MQ=255)
                                                                                                                    cgcACCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGTCa              >  1:416601/1‑139 (MQ=255)
                                                                                                                       aCCTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGTCAGcc           <  2:416601/139‑1 (MQ=255)
                                                                                                                         cTTTCAGGCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGTCAGCCAt         <  1:1684158/139‑1 (MQ=255)
                                                                                                                                gCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCt                                                                                         <  1:287274/52‑1 (MQ=255)
                                                                                                                                gCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCt                                                                                         >  2:287274/1‑52 (MQ=255)
                                                                                                                                gCCATACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGTCAGCCATCGGGTTa  <  1:654669/139‑1 (MQ=255)
                                                                                                                                   aTACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCt                                                                                         >  1:888123/1‑49 (MQ=255)
                                                                                                                                   aTACAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCt                                                                                         <  2:888123/49‑1 (MQ=255)
                                                                                                                                     aCAAGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCtgctgc                                                            <  1:661685/76‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                  
AGAAGCGCGTTCATATCGGCAGGCCAGGTCCCCAGCCACGCCATGATTTTATGGTACTGGGCATAAATATCGTTGTCGTATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTGCTGCCTGACGTTGCAGCTCGCCGAGATCGTCGCTCACCAGTTGCAGGTCAGCCATCGGGTTA  >  NC_000913/913213‑913479

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