Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 1,357,803 Δ9,014 bp [puuP]puuE [puuP], ymjE, puuA, puuD, puuR, puuC, puuB, puuE

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1357803 1366816 9014 39 [0] [0] 39 [puuP]–puuE [puuP],ymjE,puuA,puuD,puuR,puuC,puuB,puuE

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 13578020 (0.000)39 (0.600) 35/474 0.4 100% coding (1386/1386 nt) puuP putrescine:H(+) symporter PuuP
?NC_000913 1366817 = 0 (0.000)intergenic (+2/+118) puuE/pspF 4‑aminobutyrate aminotransferase PuuE/DNA‑binding transcriptional dual regulator PspF

GTTCCGGCCATTGCTGCGGTACGCGGCCTGGGGTCGATGATTGCGGTAGAGTTTAACGATCCGCAAACGGGCGAGCCGTCAGCGGCGATTGCACAGAAAATCCAGCAACGCGCGCTGGCGCAGGGGCTGCTCCTGCTGACCTGTGGCGCATACGGCAACGTGATTCGCTTCCTGTATCCGCTGACCATCCCGGATGCGCAATTCGATGCGGCAATGAAAATTTTGCAGGATGCGCTGAGCGATTAATAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGCTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTCCCTTGTTGCAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATC  >  NC_000913/1366570‑1367064
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
aactggaagtagtgggccagcggtgtgcgttcagcctctgcgacaggctctgcgcattcggtcgcgtactcatcggcaatatcaaaatattttagcgggatatcgtggttcattatttgcccctggggttaatgcgctaagcataagggcggatgatagcactcgtgatccgggtcatacagctttttatagttaatggcgcggcgcgttaccctcaggcaggataatgcgccgcgcatccgactatAGc                                                                                                                                                                                                                                                       <  1:303581‑M2/3‑1 (MQ=255)
     gaagtagtgggccagcggtgtgcgttcagcctctgcgacaggctctgcgcattcggtcgcgtactcatcggcaatatcaaaatattttagcgggatatcgtggttcattatttgcccctggggttaatgcgctaagcataagggcggatgatagcactcgtgatccgggtcatacatctttttatagttaatggcgcggcgcgttaccctcaggcaggataatgcgccgcgcatccgactatAGCCCAACg                                                                                                                                                                                                                                                 <  1:471953‑M2/9‑1 (MQ=255)
        gtagtgggccagcggtgtgcgttcagcctctgcgacaggctctgcgcattcggtcgcgtactcatcggcaatatcaaaatattttagcgggatatcgtggttcattatttgcccctggggttaatgcgctaagcataagggcggatgatagcactcgtgatccgggtcatacatctttttatagttaatggcgcggcgcgttaccctcaggcaggataatgcgccgcgcatccgactatAGCCCAACGCaa                                                                                                                                                                                                                                              <  1:70018‑M2/12‑1 (MQ=255)
                         ggcgtttagcctctgcgacaggctcggcgcattcggtcgcgtactcatcggcaatatcaaaatattttagcgggatatcggggttcattatttgcccctggtgttaaggcgctaagcataagggcggatgatagcactcgtgatccgggtcatacatctttttatagttaatggcgcggcgcgttaccctcaggcaggataatgcgccgcgcatccgactatAGCCCAACGCAATAATGTCTGAtgcgctg                                                                                                                                                                                                                             <  2:455459‑M2/29‑1 (MQ=255)
                                                           ggtcgcgtactcatcggcaatatcaaaatattttagcgggatatcgtggttcattatttgcccctggggttaatgcgctaagcataagggcggatgatagcactcgtgatccggttcatacatctttttatagttaatggcgcggcgcgttaccctcaggcaggataatgcgccgcgcatccgactatAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATtat                                                                                                                                                                                           <  1:369995‑M2/63‑1 (MQ=255)
                                                           ggtcgcgtactcatcggcaatatcaaaatatnttagcgggatatcgtggttcattatttgcccctggggttaatgcgctaagcataagggcggatgatagcactcgtgatccgggtcatacatctttttatagttaatggcgcggcgcgttaccctcaggcaggataatgcgccgcgcatccgactatAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATtat                                                                                                                                                                                           <  2:241086‑M2/63‑1 (MQ=255)
                                                                   actcatcggcaatatcaaaatattttagcgggatatcgtggttcattatttgcccctggggttaatgcgctaagcataagggcggatgatagcactcgtgatccgggtcatacatctttttatagttaatggcgcggcgcgttaccctcaggcaggataatgcgccgcgcatccgactatAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGcc                                                                                                                                                                                    <  1:457962‑M2/70‑1 (MQ=255)
                                                                      catcggcaatatcaaaatattttagcgggatatcgtggttcattatttgcccctggggttaatgcgctaagcataagggcggatgatagcactcgtgatccgggtcatacatctttttatagttaatggcgcggcgcgttaccctcaggcaggataatgcgccgcgcatccgactatAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCaaa                                                                                                                                                                                                          >  2:205102‑M2/178‑225 (MQ=255)
                                                                      catcggcaatatcaaaatattttagcgggatatcgtggttcattatttgcccctggggttaatgcgctaagcataagggcggatgatagcactcgtgatccgggtcatacatctttttatagttaatggcgcggcgcgttaccctcaggcaggataatgcgccgcgcatccgactatAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCaaa                                                                                                                                                                                                          <  1:205102‑M2/48‑1 (MQ=255)
                                                                                atcaaaatattttagcgggatatcgtggttcattatttgcccctggggttaatgcgctaagcataagggcggatgatagcactcgtgatccgggtcatacatctttttatagttaatggcgcggcgcgttaccctcaggcaggataatgcgccgcgcatccgactatAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTc                                                                                                                                                                       <  1:120997‑M2/83‑1 (MQ=255)
                                                                                              gcgggatatcgtggttcattatttgccactggggttaatgcgctaagcataagggcggatgatagcactcgtgctccgggtcatacatctttttatagttaatggcgcggcgcgttaccctcaggcaggataatgcgccgcgcatccgactatAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGc                                                                                                                                                         >  2:493822‑M2/154‑250 (MQ=255)
                                                                                                 ggatatcgtggttcattatttgcccctggggttaatgcgctaagcataagggcggatgatagcactcgtgatccgggtcatacatctttttatagttaatggcgcggcgcgttaccctcaggcaggataatgcgccgcgcatccgactatAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGt                                                                                                                                                     >  1:361058‑M2/151‑251 (MQ=255)
                                                                                                 ggatatcgtggttcattatttgcccctggggttaatgcgctaagcataagggcggatgatagcactcgtgatccgggtcatacatctttttatagttaatggcgcggcgcgttaccctcaggcaggataatgcgccgcgcatccgactatAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGt                                                                                                                                                     <  1:126966‑M2/101‑1 (MQ=255)
                                                                                                     atcgtggttcattatttgcccctggggttaatgcgctaagcataagggcggacgatagcactcgtgatccgggtcatacatctttttatagttaatggcgcggcgcgttaccctcaggcaggataatgcgccgcgcatccgactatAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTa                                                                                                                                                 <  1:86285‑M2/105‑1 (MQ=255)
                                                                                                      tcgtggttcattatttgcccctggggttaatgcgctaagcataagggcggatgatagcactcgtgatccgggtcatacatctttttatagttaatggcgcggcgcgttaccctcaggcaggataatgcgccgcgcatccgactatAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTAt                                                                                                                                                >  1:488843‑M2/146‑251 (MQ=255)
                                                                                                        gtggttcattatttgcccctggggttaatgcgctaagcataagggcggctgatagcactcgtgatccgggtcatacatctttttatagttaatggcgcggcgcgttaccctcaggcaggataatgcgccgcgcatccgactatAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTattt                                                                                                                                              >  2:206412‑M2/144‑250 (MQ=255)
                                                                                                              cattatttgcccctggggttaatgcgctaagcataagggcggatgatagcactcgtgatccgggtcatacatctttttatagttaatggcgcggcgcgttaccctcaggcaggataatgcgccgcgcatccgactatAGCCCAACGCCATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAAc                                                                                                                                        >  2:267055‑M2/138‑251 (MQ=255)
                                                                                                                            ggggttaatgcgctaagcataagggcggatgatatcactcgtgatccgggtcatacatctttttatagttaatggcgcggcgcgttaccctcaggcaggataatgcgccgcgcatccgacaatAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGCTAAATCt                                                                                                                            <  2:14658‑M2/126‑1 (MQ=255)
                                                                                                                                  aatgcgctaagcataagggcggatgatagcactcgtgatccgggtcatacatctttttatagttaatggcgcggcgcgttaccctcaggcaggataatgcgccgcgcatccgactatAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGCTAAATCTGGTGCttt                                                                                                                    >  2:132675‑M2/118‑251 (MQ=255)
                                                                                                                                            gcataagggcggatgatagcactcgtgatccgggtcatacatctttttatagttaatggcgcggcgcgttaccctcaggcaggataatgcgccgcgcatccgactatAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGCTAAATCTGGTGCTTTTTcaacaacc                                                                                                          >  2:151512‑M2/108‑250 (MQ=255)
                                                                                                                                              ataagggcggatgatagcactcgtgatccgggtcatacatctttttatagttaatggcgcggcgcgttaccctcaggcaggataatgcgccgcgcatccgactatAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGCTAAATCTGGTGCTTTTTCAACAAcgcg                                                                                                        >  1:333877‑M2/106‑251 (MQ=255)
                                                                                                                                                aaggtcgtatgatagcactcgggatccgggtcatacatttttttatagttaagggcgcggcgcgttcccctcaggcaggataatgcgccgcgcatccgacaatAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAATGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGCTAAATCTGGTGCTTTTTCAACAAcgcgc                                                                                                       <  1:453488‑M2/147‑1 (MQ=255)
                                                                                                                                                                    gtgatccgggtcatacatctttttatagttaatggcgcggcgcgttaccctcaggcaggataatgcgccgcgcatccgactatAGCCCAACGCAATAATGTCTGAt                                                                                                                                                                                                                                   <  1:290670‑M2/23‑1 (MQ=255)
                                                                                                                                                                    gtgatccgggtcatacatctttttatagttaatggcgcggcgcgttaccctcaggcaggataatgcgccgcgcatccgactatAGCCCAACGCAATAATGTCTGAt                                                                                                                                                                                                                                   >  2:290670‑M2/84‑106 (MQ=255)
                                                                                                                                                                                atccatctttttatagttaatggcgcggcgcgttaccctcaggcaggataatgcgccgcgcatccgactatAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGCTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAg                                                                      <  2:2238‑M2/180‑1 (MQ=255)
                                                                                                                                                                                    atctttttatagttaatggcgcggcgcgttaccctcaggcaggataatgcgccgcgcatccgactatAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCCAAAGACGTATTGATTAAGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGATGTATTGCTCAACTTCGCTAAATCTGGTGCTTTTTCAACAACCCGCGAACCTGATGACAAGCTAACCCCAGTAATTCAGAcgc                                                                  >  2:483336‑M2/68‑251 (MQ=255)
                                                                                                                                                                                              agttaatggcgcggcgcgttaccctcaggcaggataatgcgccgcgcatccgactatAGCCCAACGCAAGAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGCTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTg                                                         <  1:163049‑M2/193‑1 (MQ=255)
                                                                                                                                                                                                 taatggcggggcgcgttgccctcaggcaggataatgcgccgcgcatccgactatAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGCTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGtt                                                      <  1:475600‑M2/196‑1 (MQ=255)
                                                                                                                                                                                                 taatggcgcggggcgttaccctcgggaaggatcatgccccgcgcagccggctatAGCCCAACGCAATAAGGGCGGATGCGCTTCGTTTAACAGGCCTGCAAACGCCTGATTGATTATGTATGCCGAATAAAGCATGCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGCTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGtt                                                      <  2:476777‑M2/196‑1 (MQ=255)
                                                                                                                                                                                                     ggggcggcgcgttacccccaggcaggataatgcgccgcgcatccgactatAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGCTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATAtt                                                 <  2:56191‑M2/201‑1 (MQ=255)
                                                                                                                                                                                                        gcggcgcggtaccctcaggcagggtaatgcgccgcgcacccgacgatAGCCAAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGAGGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGCTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTcc                                              <  2:89020‑M2/204‑1 (MQ=255)
                                                                                                                                                                                                                   ccctcaggcaggataatgcgccgcgcatccgactatAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGCTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTCCCTTGTTGCaa                                    >  1:102610‑M2/37‑250 (MQ=255)
                                                                                                                                                                                                                     ctcaggcaggataatgcgccgcgcatccgactatAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGCTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTCCCTTGTTGCAAACt                                 >  1:253848‑M2/35‑251 (MQ=255)
                                                                                                                                                                                                                       caggctggataatgcgccgcgcatccgactatAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGTCGTATTTATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGCTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTCCCTTGTtgcaaaatg                                <  1:299442‑M2/218‑4 (MQ=255)
                                                                                                                                                                                                                                  agtccccgcgctaccgactatAGCACATCGCAATAATGTCTGATGCGCTGCGCTTGTCAGGCCTGCAAACGACGTATGGATTATGTATGCCAAAGAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGCTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTCCCTTGTTGCAAACTGAGttgcagcatc                    <  2:63625‑M2/230‑3 (MQ=255)
                                                                                                                                                                                                                                        cgcgcatccgactatAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGCTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTCCCTTGTTGCAAACTGAGTTGCAGCAACTCtttt              >  1:335317‑M2/16‑251 (MQ=255)
                                                                                                                                                                                                                                               ccgactatAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTTGCTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTCCCTTGTTGCAAACTGAGTTGCAGCAACTCTTTTTCctgc        <  1:96961‑M2/242‑1 (MQ=255)
                                                                                                                                                                                                                                                 gactatAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGCTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTCCCTTGTTGCAAACTGAGTTGCAGCAACTCTTTTTCCtgctgc     >  1:376023‑M2/7‑251 (MQ=255)
                                                                                                                                                                                                                                                    tatAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGCTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTCCCTTGTTGCAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATc  >  2:31663‑M2/4‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
GTTCCGGCCATTGCTGCGGTACGCGGCCTGGGGTCGATGATTGCGGTAGAGTTTAACGATCCGCAAACGGGCGAGCCGTCAGCGGCGATTGCACAGAAAATCCAGCAACGCGCGCTGGCGCAGGGGCTGCTCCTGCTGACCTGTGGCGCATACGGCAACGTGATTCGCTTCCTGTATCCGCTGACCATCCCGGATGCGCAATTCGATGCGGCAATGAAAATTTTGCAGGATGCGCTGAGCGATTAATAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGCTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTCCCTTGTTGCAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATC  >  NC_000913/1366570‑1367064

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