Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,496,758 Δ13,333 bp insC6insO insC6, insD6, yjgZ, insG, nanX, nanY, ythA, insO, insI3, insO, insO

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4496758 4510090 13333 39 [0] [0] 39 insC6–insO insC6,insD6,yjgZ,insG,nanX,nanY,ythA,insO,insI3,insO,insO

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 44967570 (0.000)39 (0.570) 39/494 0.5 100% intergenic (+268/‑1515) leuX/insC6 tRNA‑Leu/IS2 insertion element repressor InsA
?NC_000913 4510091 = 0 (0.000)intergenic (+298/+599) insO/fecE IS911B regulator fragment/ferric citrate ABC transporter ATP binding subunit

CTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGGCGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGATTGAGTGGCATGGCACCAATAAA  >  NC_000913/4496512‑4497004
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       
cTGAGGACTTTTTTTATGCCTGCAACCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTAACAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCACCGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagat                                                                                                                                                                                                                                                     >  2:420193‑M1/1‑246 (MQ=255)
  gAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagatat                                                                                                                                                                                                                                                   <  2:537655‑M1/250‑7 (MQ=255)
   aGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagatatca                                                                                                                                                                                                                                                 <  1:123547‑M1/251‑9 (MQ=255)
    ggTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagatatcat                                                                                                                                                                                                                                                <  1:196197‑M1/251‑10 (MQ=255)
     gTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagatatcatt                                                                                                                                                                                                                                               <  1:532662‑M1/251‑11 (MQ=255)
         ttttttATGCCGGAAAGCCAGTGTTTATCTATCTTGCCCGCTATATTAACTCTCTCAAGGTCAACCGATAGCAACGTCCATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGATTTTACCATCAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagatatcattgat                                                                                                                                                                                                                                            <  2:330767‑M1/250‑14 (MQ=255)
          tttttATGACTGAAATCCAGGGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACAGATATCAACGTACATCTACCAACATATGTTGGTACAGTTGATGGTATTTCCGGTTCGATAATGCTTGAACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGCCAGAGCAGCCAAGCCAACGGATACGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagatatcatttata                                                                                                                                                                                                                                           >  2:20172‑M1/1‑236 (MQ=255)
               aTGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagatatcattgataacaatc                                                                                                                                                                                                                                     >  1:266996‑M1/1‑231 (MQ=255)
               aTGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagatatcattgataacaatc                                                                                                                                                                                                                                     <  2:424025‑M1/251‑21 (MQ=255)
                   cTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagatatcattgataacaatcgatc                                                                                                                                                                                                                                 >  1:510428‑M1/1‑227 (MQ=255)
                      aaaTCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagatatcattgataacaatcgatccat                                                                                                                                                                                                                              >  2:214650‑M1/1‑224 (MQ=255)
                                 ttATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagatatcattgataacaatcgatccatcaactggcctg                                                                                                                                                                                                                   >  2:437898‑M1/1‑213 (MQ=255)
                                           ttCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagatatcattgataacaatcgatccatcaactggcctgtataac                                                                                                                                                                                                             <  2:457427‑M1/247‑45 (MQ=255)
                                                     atTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTTGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagatatcattgataacaatcgatccatcaactggcctgtataacgtggctgccggtg                                                                                                                                                                                                <  2:451823‑M1/250‑58 (MQ=255)
                                                     atTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagatatcattgataacaatcgatccatcaactggcctgtataacgtggctgccggtg                                                                                                                                                                                                >  2:446195‑M1/1‑193 (MQ=255)
                                                        aaCTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagatatcattgataacaatcgatccatcaactggcctgtataacgtggctgccggtgatt                                                                                                                                                                                             <  1:435257‑M1/250‑61 (MQ=255)
                                                                           gATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagatatcattgataacaatcgatccatcaactggcctgtataacgtggctgccggtgattaagaaagctgcacctacct                                                                                                                                                                          >  2:149052‑M1/1‑171 (MQ=255)
                                                                            atatCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagatatcattgataacaatcgatccatcaactggcctgtataacgtggctgccggtgattaagaaagctgcacctaccta                                                                                                                                                                         >  1:239920‑M1/1‑170 (MQ=255)
                                                                                   cGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagatatcattgataacaatcgatccatcaactggcctgtataacgtggctgccggtgattaagaaagctgcacctacctaagtagtag                                                                                                                                                                 <  1:152666‑M1/251‑89 (MQ=255)
                                                                                                    atGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagatatcattgataacaatcgatccatcaactggcctgtataacgtggctgccggtgattaagaaagctgcacctacctaagtagtagcaaacgcacacttttt                                                                                                                                                 >  1:148943‑M1/1‑146 (MQ=255)
                                                                                                      gTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagatatcattgataacaatcgatccatcaactggcctgtataacgtggctgccggtgattaagaaagctgcacctacctaagtagtagcaaacgcacactttttag                                                                                                                                               <  1:492009‑M1/250‑107 (MQ=255)
                                                                                                      gTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagatatcattgataacaatcgatccatcaactggcctgtataacgtggctgccggtgattaagaaagctgcacctacctaagtagtagcaaacgcacacttttt                                                                                                                                                 >  1:107243‑M1/1‑144 (MQ=255)
                                                                                                                   atgGAATTACCGGCTCGATAATGCTTGTACCAACAGGGAGGGAATACGATGGAATTAACAGATATCAACGTCAGAGCAGCCAAGCCAACGGATAAGCACTATTAGCTGACTGATGGTGGCGGTATGCATCTagatatcattgataacaatcgatccatcaactggcctgtataacgtggctgccggtgattaagaaagctgcacctacctaagtagaagcaaacgcaaactttttagaaaaatcgatggtc                                                                                                                                 >  2:370668‑M1/1‑131 (MQ=255)
                                                                                                                      gTATTTCCGGTTCGATAATGCTTGTACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagatatcattgataacaatcgatccatcaactggcctgtataacgtggctgccggtgattaagaaagctgcacctacctaagtagtagcaaacgcacactttttagaaaaatcgatggtc                                                                                                                                 >  1:314248‑M1/1‑128 (MQ=255)
                                                                                                                                  cGATAATGCTTGTACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagatatcattgataacaatcgatccatcaactggcctgtataacgtggctgccggtgattaagaaagctgcacctacctaagtagtagcaaacgcactctttttagaaaaatcgatggtcagaaactggattagc                                                                                                                  <  2:314248‑M1/251‑136 (MQ=255)
                                                                                                                                   gATAATGCTTGTACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagatatcattgataacaatcgatccatcaactggcctgtataacgtggctgccggtgattaagaaagctgcacctacctaagtagtagcaaacgcacactttttagaaaaatcgatggtcagaaactggattagca                                                                                                                 <  2:239920‑M1/251‑137 (MQ=255)
                                                                                                                                         gCTTGTACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagatatcattgataacaatcgatccatcaactggcctgtataacgtggctgccggtgattaagaaagctgcacctacctaagtagtagcaaacgcacactttttagaaaaatcgatggtcagaaactggattagcaattcc                                                                                                            <  1:496875‑M1/250‑142 (MQ=255)
                                                                                                                                          cTTGTACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagatatcattgataacaatcgatccatcaactggcctgtataacgtggctgccggtgattaagaaagctgcacctacctaagtagtagcaaacgcacactttttagaaaaatcgatggtcagaaactggattagcaattccgt                                                                                                          >  2:526566‑M1/1‑108 (MQ=255)
                                                                                                                                                     agggagggAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagatatcattgataacaatcgatccatcaactggcctgtataacgtggctgccggtgattaagaaagctgcacctacctaagtagtagcaaacgcacactttttagaaaaatcgatggtcagaaactggattagcaattccgttccatggttgc                                                                                               >  1:493843‑M1/1‑97 (MQ=255)
                                                                                                                                                              aTACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagatatcattgataacaatcgatccatcaactggcctgtataacgtggctgccggtgattaagaaagctgcacctacctaagtagtagcaaacgcacactttgtagaaaaatcgatggtcagaaactggattagcaattccgttccatggttgcttttgattt                                                                                      <  2:491931‑M1/251‑164 (MQ=255)
                                                                                                                                                                     ggCAGTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagatatcattgataacaatcgatccatcaactggcctgtataacgtggctgccggtgattaagaaagctgcacctacctaagtagtagcaaacgcacactttttagaaaaatcgatgggcagaaactggattagcaattccgttccatggttgcttttgatttacgttg                                                                                <  1:420193‑M1/250‑170 (MQ=255)
                                                                                                                                                                           aaCAGATATCAAAGTCATAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTTATGGTGGCGGTATGCATATagatatcattgataacaatcgttccatcaactggcctgtataacgtggctgccggtgattaagaaagctgcacctacctaagtagtagcaaacgcacacgttttagaaaaatcgatggtcagaaactggattagcgattccgtgccatggttgcttgtgattttcgatcgcgtct                                                                          <  2:489340‑M1/250‑176 (MQ=255)
                                                                                                                                                                                  atCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagatatcattgataacaatcgatccatcaactggcctgtataacgtggctgccggtgattaagaaagctgcacctacctaagtagtagcaaacgcacactttttagaaaaatcgatggtcagaaactggattagcaattccgttccatggttgcttttgatttacgttggcgtctgatcatt                                                                   <  2:462147‑M1/250‑183 (MQ=255)
                                                                                                                                                                                     aaaGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagatatcattgataacaatcgatccatcaactggcctgtataacgtggctgccggtgattaagaaagctgcacctacctaagtagtagcaaacgcacactttttagaaaaatcgatggtcagaaactggattagcaattccgttccatggttgcttttgatttacgttggcgtctgatcattgatt                                                               <  1:526566‑M1/251‑187 (MQ=255)
                                                                                                                                                                                      aaGTCAGAGCAGCCAAGCCAACGGCTAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagaaatcattgataacaatcgatccatcaactggcctgtataacgtggctgccggtgattaagaaagctgcacctacctaagtagtcgcaaacgcacactttttagaaaaatcgatggtcagaaactggattagcaattccgttccatggttgcttttgatttacgttggcgcctgatcattgattt                                                              >  2:306962‑M1/1‑64 (MQ=255)
                                                                                                                                                                                                 gccaagccaaCGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTagatatcattgataacaatcgatccatcaactggcctgtataacgtggctgccggtgattaagaaagctgcacctacctaagtagtagcaaacgcacactttttagaaaaatcgatggtcagaaactggattagcaattccgttccatggttgcttttgatttacgttggcgtctgatcattgatttatcctcaaa                                                     >  1:513320‑M1/1‑53 (MQ=255)
                                                                                                                                                                                                                        ataGCTTACTGATGGTGGCGGTATGCATCTagatatcattgataacaatcgatccatcaactggcctgtataacgtggctgccggtgattaagaaagctgcacctacctaagtagtagcaaacgcacactttttagaaaaatcgatggtcagaaactggattagcaattccgttccatggttgcttttgatttacgttggcgtctgatcattgatttatcctcaaaagcccaacctcattggtaatgaacc                            <  2:320237‑M1/250‑222 (MQ=255)
                                                                                                                                                                                                                                            gTATGCATCTagatatcattgataacaatcgatccatcaactggcctgtataacgtggctgccggtgattaagaaagctgcacctacctaagtagtagcaaacgcacactttttagaaaaatcgatggtcagaaactggattagcaattccgttccatggttgcttttgatttacgttggcgtctgatcattgatttatcctcaaaagcccaacctcattggtaatgaaccagctccgtgaatgtccgctc        >  1:231669‑M1/1‑10 (MQ=255)
                                                                                                                                                                                                                                                  aTCTagatatcgttgataacaatcgatccatcaactggcctgtataacgtggctgccggtgattaagaaagctgcacccacctaagcagtagcaaacgcacactttttagaaaaatcgatggtcagaaactggattagcaattccgttccatggttgctttttatttacgttggcgtctgatcattgatttatcctcaacagcccaacctcattggtaatgaacaagctccgtgaatgtccgctctggcac  >  2:333642‑M1/1‑4 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       
CTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTACCAACAGGGAGGGAATACGATGGCATTAACAGATATCAAAGTCAGAGCAGCCAAGCCAACGGATAAGCAATATTAGCTGACTGATGGTGGCGGTATGCATCTGCTTGTCCATCCAAATGGTTCTAAGTACTGGCGTTTGCAGTACCGTTATGAGGGAAAGCAAAAAATGCTGGCACTTGGGGTTTATCCTGAAATCACACTAGCGGATGCCAGAGTACGTCGTGACGAGGCGCGTAAGCTGCTTGCGAATGGCGTCGATCCGGGAGACAAAAAGAAAAATGATAAGGTTGAACAGAGTAAAGCACGAACCTTTAAAGAAGTCGCGATTGAGTGGCATGGCACCAATAAA  >  NC_000913/4496512‑4497004

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