Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,130,853 C→T S9N (AGC→AAC) 
*406* (TAG→TAA
wcaD ←
wcaC ←
putative colanic acid polymerase
putative glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,130,8530CT100.0% 225.8 / NA 69S9N (AGC→AAC) 
*406* (TAG→TAA
wcaD
wcaC
putative colanic acid polymerase
putative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (33/36);  total (33/36)

TAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGATTTCCGGTTTGCTTAACTGCA  >  NC_000913/2130722‑2130986
                                                                                                                                   |                                                                                                                                     
tAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACa                                                                                                                                   <  1:709319/136‑1 (MQ=255)
tAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACa                                                                                                                                   >  2:709319/1‑136 (MQ=255)
tAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGAt                                                                                                                                <  2:806043/139‑1 (MQ=255)
     cTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGAtt                                                                                                                               <  2:1705231/135‑1 (MQ=255)
     cTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGAtt                                                                                                                               >  1:1705231/1‑135 (MQ=255)
         tGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGa                                                                                                                       <  1:90188/139‑1 (MQ=255)
               cTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGa                                                                                                                 <  2:788826/139‑1 (MQ=255)
               cTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGa                                                                                                                 >  1:844896/1‑139 (MQ=255)
                 aaaaaCAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACa                                                                                                               <  1:1165058/139‑1 (MQ=255)
                   aaaCAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATa                                                                                                             <  1:1442029/139‑1 (MQ=255)
                          agTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACtcctcc                                                                                                      <  2:646734/139‑1 (MQ=255)
                               cAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCAt                                                                                                 >  2:912760/1‑139 (MQ=255)
                                  gACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGa                                                                                                                       >  1:1391739/1‑114 (MQ=255)
                                  gACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGa                                                                                                                       <  2:1391739/114‑1 (MQ=255)
                                   aCAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCtgt                                                                                             >  1:500735/1‑139 (MQ=255)
                                      ggTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGAt                                                                                                                          <  2:1451717/107‑1 (MQ=255)
                                      ggTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGAt                                                                                                                          >  1:1451717/1‑107 (MQ=255)
                                      ggTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGa                                                                                                                       <  2:1183870/110‑1 (MQ=255)
                                      ggTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGa                                                                                                                       >  1:1183870/1‑110 (MQ=255)
                                      ggTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCtgttgt                                                                                          >  1:152722/1‑139 (MQ=255)
                                      ggTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCtgttgt                                                                                          <  2:421841/139‑1 (MQ=255)
                                       gTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGa                                                                                                                 >  1:994188/1‑115 (MQ=255)
                                       gTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGa                                                                                                                 <  2:994188/115‑1 (MQ=255)
                                       gTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTc                                                                                         >  2:1810968/1‑139 (MQ=255)
                                           aaaaGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACt                                                                                     <  1:1070628/139‑1 (MQ=255)
                                                 gACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATa                                                                                                             <  1:699510/109‑1 (MQ=255)
                                                 gACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATa                                                                                                             >  2:699510/1‑109 (MQ=255)
                                                  aCGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGGTTTTTACGGTAACCAGCAGATAAATCAGCGGCAGCAACAGGGAGTTACAGATTCTGATAGAAGTTGACCTACTCCTCCAGCATCTGTTGTCCACTGTAggcg                                                                              >  2:68915/1‑139 (MQ=255)
                                                  aCGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAggcg                                                                              >  1:379586/1‑139 (MQ=255)
                                                          gATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCt                                                                      >  1:1463244/1‑139 (MQ=255)
                                                            tGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGc                                                                                                   <  2:1620191/108‑1 (MQ=255)
                                                            tGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGc                                                                                                   >  1:1620191/1‑108 (MQ=255)
                                                              gggAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCtg                                                                                              <  2:1882375/111‑1 (MQ=255)
                                                              gggAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCtg                                                                                              >  1:1882375/1‑111 (MQ=255)
                                                                  aaCTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGc                                                              <  1:1810968/139‑1 (MQ=255)
                                                                   aCTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACCGGTAGTTACAGATTCTGATAg                                                                                                                        >  1:440777/1‑80 (MQ=255)
                                                                   aCTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAg                                                                                                                        <  2:440777/80‑1 (MQ=255)
                                                                    cTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTg                                                            >  2:853718/1‑139 (MQ=255)
                                                                     tttCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGa                                                           >  1:140677/1‑139 (MQ=255)
                                                                        cGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACt                                                        >  1:973876/1‑139 (MQ=255)
                                                                        cGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACt                                                        >  2:629108/1‑139 (MQ=255)
                                                                                 gggCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCg                                               <  2:152722/139‑1 (MQ=255)
                                                                                 gggCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCg                                               <  2:1580449/139‑1 (MQ=255)
                                                                                 gggCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCg                                               <  1:853718/139‑1 (MQ=255)
                                                                                  ggCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGc                                                              >  2:1843384/1‑123 (MQ=255)
                                                                                  ggCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGc                                                              <  1:1843384/123‑1 (MQ=255)
                                                                                   gCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTg                                             <  1:629108/139‑1 (MQ=255)
                                                                                   gCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTg                                             <  2:973876/139‑1 (MQ=255)
                                                                                       ttttGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTAc                                         <  2:1463244/139‑1 (MQ=255)
                                                                                                                 aGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCg                                               >  1:142228/1‑107 (MQ=255)
                                                                                                                 aGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCg                                               <  2:142228/107‑1 (MQ=255)
                                                                                                                 aGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGATTTCCGGt               >  1:151290/1‑139 (MQ=255)
                                                                                                                 aGCGGCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGATTTCCGGt               >  2:1079569/1‑139 (MQ=255)
                                                                                                                    ggCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCAc                                                                                      <  2:431784/65‑1 (MQ=255)
                                                                                                                    ggCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCAc                                                                                      >  1:431784/1‑65 (MQ=255)
                                                                                                                    ggCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAggcgg                                                                             >  2:487270/1‑74 (MQ=255)
                                                                                                                    ggCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAggcgg                                                                             >  2:1217550/1‑74 (MQ=255)
                                                                                                                    ggCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAggcgg                                                                             <  1:1217550/74‑1 (MQ=255)
                                                                                                                    ggCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAggcgg                                                                             <  1:487270/74‑1 (MQ=255)
                                                                                                                    ggCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCg                                               >  2:670041/1‑104 (MQ=255)
                                                                                                                    ggCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCg                                               <  1:670041/104‑1 (MQ=255)
                                                                                                                    ggCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGATTTCCGGTTTg            <  2:844896/139‑1 (MQ=255)
                                                                                                                    ggCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCCGAACTCAGCCAGCGTGGTAc                                         <  2:83898/110‑1 (MQ=255)
                                                                                                                    ggCAGCAACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCCGAACTCAGCCAGCGTGGTAc                                         >  1:83898/1‑110 (MQ=255)
                                                                                                                      cagcaACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCg                                               >  2:1340254/1‑102 (MQ=255)
                                                                                                                      cagcaACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCg                                               <  1:1340254/102‑1 (MQ=255)
                                                                                                                        gcaACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGgcggcg                                                                           <  1:987623/72‑1 (MQ=255)
                                                                                                                        gcaACAGGTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGgcggcg                                                                           >  2:987623/1‑72 (MQ=255)
                                                                                                                              ggTAGTTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGATTTCCGGTTTGCTTAACTGCa  <  2:500735/139‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                     
TAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAGATAAATCAGCGGCAGCAACAGGTAGCTACAGATTCTGATAGAAGTTGACATACTCCTCCAGCATCTGTTGTCCACTGTAGGCGGCGCGGCTGCGTTGGCTGAACTCAGCCAGCGTGGTACCAAATATCGCCTGCGCGATTTCCGGTTTGCTTAACTGCA  >  NC_000913/2130722‑2130986

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