Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 1,299,499 Δ1,199 bp intergenic (+254/‑485) ychE → / → oppA UPF0056 family inner membrane protein/oligopeptide ABC transporter periplasmic binding protein

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1299499–1300696 1300697 2–1199 67 [0] [0] 65 ychE/oppA UPF0056 family inner membrane protein/oligopeptide ABC transporter periplasmic binding protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 12994980 (0.000)62 (0.940) 57/290 0.3 100% intergenic (+253/‑1684) ychE/oppA UPF0056 family inner membrane protein/oligopeptide ABC transporter periplasmic binding protein
?NC_000913 1300698 = 0 (0.000)intergenic (+1453/‑484) ychE/oppA UPF0056 family inner membrane protein/oligopeptide ABC transporter periplasmic binding protein

AACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATA  >  NC_000913/1300549‑1300845
                                                                                                                                                     |                                                                                                                                                   
cttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaat                                                                                                                                                    >  1:356657‑M2/150‑151 (MQ=255)
 ttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaat                                                                                                                                                    >  1:441333‑M2/149‑150 (MQ=255)
    tcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaataat                                                                                                                                                 >  1:628192‑M2/146‑150 (MQ=255)
     caaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaataat                                                                                                                                                 <  1:83440‑M2/5‑1 (MQ=255)
     caaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaataa                                                                                                                                                  <  2:645685‑M2/4‑1 (MQ=255)
        attacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATc                                                                                                                                             >  2:474154‑M2/142‑150 (MQ=255)
          tacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCaa                                                                                                                                           <  2:761204‑M2/11‑1 (MQ=255)
          tacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAAt                                                                                                                                          >  1:256065‑M2/140‑151 (MQ=255)
           acagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAAtt                                                                                                                                         <  1:324232‑M2/13‑1 (MQ=255)
             agttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGt                                                                                                                                       >  1:392285‑M2/137‑151 (MQ=255)
             agttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGt                                                                                                                                       <  2:405834‑M2/15‑1 (MQ=255)
                 agttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTaaa                                                                                                                                   <  1:848628‑M2/19‑1 (MQ=255)
                 agttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTaaa                                                                                                                                   <  2:904822‑M2/19‑1 (MQ=255)
                 agttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTaa                                                                                                                                    <  2:406241‑M2/18‑1 (MQ=255)
                 agttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTaa                                                                                                                                    <  1:139670‑M2/18‑1 (MQ=255)
                           atttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGc                                                                                                                          >  2:568214‑M2/123‑150 (MQ=255)
                            tttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCAt                                                                                                                        >  2:670640‑M2/122‑151 (MQ=255)
                             ttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCAt                                                                                                                        <  2:524250‑M2/30‑1 (MQ=255)
                               ccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactggctttggtataggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCAttt                                                                                                                      >  2:819130‑M2/119‑150 (MQ=255)
                               ccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTc                                                                                                                     <  1:219028‑M2/33‑1 (MQ=255)
                                          tctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAAc                                                                                                          >  1:223849‑M2/108‑151 (MQ=255)
                                             cctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACt                                                                                                         <  1:643438‑M2/45‑1 (MQ=255)
                                                 accatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAAt                                                                                                    >  2:806924‑M2/101‑150 (MQ=255)
                                                         gttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGaaa                                                                                              >  1:934556‑M2/93‑148 (MQ=255)
                                                           tattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGAt                                                                                          >  1:379312‑M2/91‑150 (MQ=255)
                                                               ttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGAc                                                                                      <  2:36407‑M2/64‑1 (MQ=255)
                                                                 cgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATg                                                                                   <  2:356646‑M2/67‑1 (MQ=255)
                                                                    cacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCtt                                                                                <  1:955769‑M2/70‑1 (MQ=255)
                                                                     acatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTa                                                                               <  2:713106‑M2/71‑1 (MQ=255)
                                                                      catcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTa                                                                               <  2:21268‑M2/71‑1 (MQ=255)
                                                                            aatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAAtc                                                                        <  2:733623‑M2/78‑1 (MQ=255)
                                                                            aatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAAtc                                                                        >  2:48051‑M2/74‑151 (MQ=255)
                                                                            aatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAAt                                                                         <  1:1059891‑M2/77‑1 (MQ=255)
                                                                               cctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAAtctc                                                                      <  1:277967‑M2/80‑1 (MQ=255)
                                                                                 tgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTAt                                                                   >  1:331812‑M2/69‑151 (MQ=255)
                                                                                 tgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTAt                                                                   >  2:829490‑M2/69‑151 (MQ=255)
                                                                                    gcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATttg                                                                >  1:457021‑M2/66‑151 (MQ=255)
                                                                                    gcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATttg                                                                <  2:1051367‑M2/86‑1 (MQ=255)
                                                                                     cttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATttg                                                                >  1:1017276‑M2/65‑150 (MQ=255)
                                                                                        gctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATttgattg                                                            <  2:370276‑M2/90‑1 (MQ=255)
                                                                                            aagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGattgaatt                                                        <  2:797928‑M2/94‑1 (MQ=255)
                                                                                                ataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGaaa                                                    <  2:331801‑M2/98‑1 (MQ=255)
                                                                                                    ttgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGAtg                                                >  2:109704‑M2/50‑151 (MQ=255)
                                                                                                         atgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGAtgtttgt                                           <  2:252897‑M2/107‑1 (MQ=255)
                                                                                                            atattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTaa                                        <  1:692796‑M2/110‑1 (MQ=255)
                                                                                                               ttattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGc                                     >  2:291019‑M2/39‑151 (MQ=255)
                                                                                                                tattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGc                                     <  1:275788‑M2/113‑1 (MQ=255)
                                                                                                                tattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGc                                     >  2:580818‑M2/38‑150 (MQ=255)
                                                                                                                      ttccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATg                              >  2:262038‑M2/32‑151 (MQ=255)
                                                                                                                         cactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCa                            >  1:581848‑M2/29‑150 (MQ=255)
                                                                                                                          actgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAGGATGTTTGTTAAGGCATGGATGCAAg                          <  2:412882‑M2/124‑1 (MQ=255)
                                                                                                                             gcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCt                        <  1:532981‑M2/126‑1 (MQ=255)
                                                                                                                                tttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCtat                      >  2:67403‑M2/22‑149 (MQ=255)
                                                                                                                                 ttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCtata                     <  2:452805‑M2/129‑1 (MQ=255)
                                                                                                                                   ggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGAt                  <  1:896006‑M2/132‑1 (MQ=255)
                                                                                                                                    gtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTc                >  1:939832‑M2/18‑151 (MQ=255)
                                                                                                                                    gtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTc                >  2:376905‑M2/18‑151 (MQ=255)
                                                                                                                                    gtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAAATGTTTGTTAAGGCATGGATGCAAGCTATAGATTc                <  2:176829‑M2/134‑1 (MQ=255)
                                                                                                                                         ggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATa           <  1:727789‑M2/139‑1 (MQ=255)
                                                                                                                                           aaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACg         >  2:494172‑M2/11‑151 (MQ=255)
                                                                                                                                            aagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACgg        >  1:737002‑M2/10‑151 (MQ=255)
                                                                                                                                                gctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCaa    <  2:508232‑M2/146‑1 (MQ=255)
                                                                                                                                                   aaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATa  <  2:671543‑M2/148‑1 (MQ=255)
                                                                                                                                                   aaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATa  <  2:763264‑M2/148‑1 (MQ=255)
                                                                                                                                                   aaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATa  <  2:101628‑M2/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATA  >  NC_000913/1300549‑1300845

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