Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 2,525,917 Δ9 bp intergenic (+26/+5) yfeN → / ← yfeR putative outer membrane protein/transcriptional regulator of yefH

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 2525917 2525925 9 80 [0] [0] 70 yfeN/yfeR putative outer membrane protein/transcriptional regulator of yefH

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 25259160 (0.000)70 (1.180) 54/210 0.1 100% intergenic (+25/+14) yfeN/yfeR putative outer membrane protein/transcriptional regulator of yefH
?NC_000913 2525926 = 0 (0.000)intergenic (+35/+4) yfeN/yfeR putative outer membrane protein/transcriptional regulator of yefH

TGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTT  >  NC_000913/2525778‑2526037
                                                                                                                                          |                                                                                                                         
tggAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCa                                                                                                                           <  2:679520/139‑1 (MQ=255)
tggAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCa                                                                                                                           <  2:261985/139‑1 (MQ=255)
  gAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAca                                                                                                                         <  1:469098‑M1/139‑3 (MQ=255)
       cAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCa                                                                                                                           <  1:64105/132‑1 (MQ=255)
       cAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCa                                                                                                                           >  2:64105/1‑132 (MQ=255)
                    tactacGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCa                                                                                                                           >  2:16023/1‑119 (MQ=255)
                    tactacGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCa                                                                                                                           <  1:16023/119‑1 (MQ=255)
                     actacGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacg                                                                                                      <  2:611140‑M1/139‑22 (MQ=255)
                         cGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatc                                                                                                  >  1:56461‑M1/1‑114 (MQ=255)
                            ggATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcact                                                                                                                       <  1:1099992‑M1/115‑5 (MQ=255)
                            ggATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcact                                                                                                                       >  2:1099992‑M1/1‑111 (MQ=255)
                              aTTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatcccctt                                                                                             >  1:915403‑M1/1‑109 (MQ=255)
                              aTTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatcccctt                                                                                             >  2:350902‑M1/1‑109 (MQ=255)
                                   gTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttata                                                                                                                  <  1:651103‑M1/113‑10 (MQ=255)
                                   gTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttata                                                                                                                  >  2:651103‑M1/1‑104 (MQ=255)
                                    tATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatc                                                                                                                >  2:443501‑M1/1‑103 (MQ=255)
                                    tATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatc                                                                                                                <  1:443501‑M1/114‑12 (MQ=255)
                                      tCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgccca                                                                                     <  2:915403‑M1/139‑39 (MQ=255)
                                        gATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAACGTAACTTCTGATTGATTACTGCCGGTGCTTCCGGCAcactttatatctgatacaacggatccccttcccgcccagc                                                                                   >  1:1332360‑M1/1‑99 (MQ=255)
                                         atatGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgcccagcc                                                                                  >  1:647829‑M1/1‑98 (MQ=255)
                                             gCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgcccagccatta                                                                              <  1:89569‑M1/139‑46 (MQ=255)
                                              cGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgcccagccattag                                                                             >  2:890095‑M1/1‑93 (MQ=255)
                                              cGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttat                                                                                                                   <  2:1261191‑M1/101‑9 (MQ=255)
                                              cGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttat                                                                                                                   >  1:1261191‑M1/1‑93 (MQ=255)
                                                  tGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTGCTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgcccagccattagcgca                                                                         <  1:890095‑M1/139‑51 (MQ=255)
                                                   gACAAACTGGTCGATGATTTTTACCAGGATGCTATAATCTACTCCATCAAATTTAACTTCTCATTGATTACTGCCGCGGCTCCCGCCAcactttatatctgatacaccggatccccttcccgcccagccattagcgcat                                                                        <  2:104810‑M1/139‑52 (MQ=255)
                                                    aCAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCa                                                                                                                           >  2:728743/1‑87 (MQ=255)
                                                    aCAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCa                                                                                                                           <  1:728743/87‑1 (MQ=255)
                                                        cctGGGCTGAGATTTCTACCAGGATGCGATAATCTACTCCAACAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatccgaaacaacggagccccttcccgcccagccattagcgcattcccc                                                                   <  2:1332360‑M1/138‑57 (MQ=255)
                                                        aCTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgcccagccattagcgcattgccc                                                                   <  1:1085704‑M1/139‑57 (MQ=255)
                                                          tGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgcccagccattagcgcattgccctg                                                                 >  1:104810‑M1/1‑81 (MQ=255)
                                                           gggCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgcccagccattagcgcattgccctgg                                                                >  1:199218‑M1/1‑80 (MQ=255)
                                                             gCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatc                                                                                                                >  1:575454‑M1/1‑78 (MQ=255)
                                                             gCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatc                                                                                                                <  2:575454‑M1/89‑12 (MQ=255)
                                                             gCGAAGATTTTTACCAGGATGCGATAATCTACTCCAGCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgcccagccattagcgcattgccctggtc                                                              >  2:431376‑M1/1‑78 (MQ=255)
                                                              cGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcacttt                                                                                                                     >  1:732724‑M1/1‑77 (MQ=255)
                                                              cGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcacttt                                                                                                                     <  2:732724‑M1/83‑7 (MQ=255)
                                                                     ttttACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgcccagccattagcgcattgccctggtctcgtacca                                                      <  1:350902‑M1/139‑70 (MQ=255)
                                                                             ggATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgccc                                                                                      <  2:227168‑M1/99‑38 (MQ=255)
                                                                             ggATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgccc                                                                                      >  1:227168‑M1/1‑62 (MQ=255)
                                                                              gATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgcccagccattagcgcattgccctggtctcgtaccacgtcccata                                             >  2:1140493‑M1/1‑61 (MQ=255)
                                                                               aTGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgcccagccattagcgcattgccctggtctcgtaccacgtcccataa                                            <  1:1027636‑M1/139‑80 (MQ=255)
                                                                               aTGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgcccagccattagcgcattgccctggtctcgtaccacgtcccataa                                            >  2:1027636‑M1/1‑60 (MQ=255)
                                                                               aTGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatc                                                                                                  <  1:849738‑M1/85‑26 (MQ=255)
                                                                               aTGCGATAATCTACTCCATCAAATTTAACTTCAGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatc                                                                                                  >  2:849738‑M1/1‑60 (MQ=255)
                                                                                 gCGATAATCTACTCCATCGAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgccca                                                                                     >  2:822847‑M1/1‑58 (MQ=255)
                                                                                 gCGATAATCTACTCCATCGAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgccca                                                                                     <  1:822847‑M1/96‑39 (MQ=255)
                                                                                            cTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCa                                                                                                                           >  2:64107/1‑47 (MQ=255)
                                                                                            cTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCa                                                                                                                           <  1:64107/47‑1 (MQ=255)
                                                                                               cATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttttatctgataccacgggtccccttcccgcccagccattaccgcattgccctggtcccgtccctcgtccccttcagcttcaggggcggga                            >  2:1020594‑M1/1‑44 (MQ=255)
                                                                                               cATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgcccagccattagcgcattgccctggtctcgtaccacgtcccataaagcttcagcggcggtg                            >  1:1242145‑M1/1‑44 (MQ=255)
                                                                                                 tCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgcccagccattagcgcattgccctggtctcgtaccacgtcccataaagcttcagcggcggtgg                           >  2:1306077‑M1/1‑42 (MQ=255)
                                                                                                 tCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgcccagccattagcgcattgccctggtctcgtaccacgtcccataaagcttcagcggcggtgg                           <  1:1306077‑M1/138‑97 (MQ=255)
                                                                                                 tCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgcccagccattagcgcattgccc                                                                   <  1:522039‑M1/98‑57 (MQ=255)
                                                                                                 tCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgcccagccattagcgcattgccc                                                                   >  2:522039‑M1/1‑42 (MQ=255)
                                                                                                 tCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatcaccttaccgcccagccattagcgcattgccctggtctcgtaccacgtcccatcaagcttcagcggcggtgga                          >  2:1000015‑M1/1‑42 (MQ=255)
                                                                                                      tttAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgcccagccattagcgcattgccctggtctcgtaccacgtcccataaagcttcagcggcggtggaaagtg                     <  2:1242145‑M1/139‑103 (MQ=255)
                                                                                                      tttAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgcccagccattagcgcattgccctggtc                                                              >  1:389689‑M1/1‑37 (MQ=255)
                                                                                                      tttAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgcccagccattagcgcattgccctggtc                                                              <  2:389689‑M1/98‑62 (MQ=255)
                                                                                                         aaCTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgc                                                                                        >  1:404157‑M1/1‑34 (MQ=255)
                                                                                                         aaCTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgc                                                                                        <  2:404157‑M1/69‑36 (MQ=255)
                                                                                                           cTTCTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgcccagccattagcgcattgccctggtctcgtaccacgtcccataaagcttcagcggcggtggaaagtgaacgg                >  2:482342‑M1/1‑32 (MQ=255)
                                                                                                              cTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgcccagccattagcgcattgccctggtctcgtaccacgtcccataaagcttcagcggcggtggaaagtg                     <  2:1049275‑M1/131‑103 (MQ=255)
                                                                                                              cTGATTGATTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgcccagccattagcgcattgccctggtctcgtaccacgtcccataaagcttcagcggcggtggaaagtg                     >  1:1049275‑M1/1‑29 (MQ=255)
                                                                                                               tgattgatTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgcccagccattagcgcattgccctggtctcgtaccacgtcccataaagcttcagcggcggtggaaagtgaacggttct            <  1:431376‑M1/139‑112 (MQ=255)
                                                                                                               tgattgatTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgcc                                                                                       <  1:1246945‑M1/64‑37 (MQ=255)
                                                                                                               tgattgatTACTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgcc                                                                                       >  2:1246945‑M1/1‑28 (MQ=255)
                                                                                                                        aCTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgcccagccattagcgcattgccctggtctcgtaccacgtcccataaagcttcagcggcggtggaaagtgaac                  <  1:244514‑M1/124‑106 (MQ=255)
                                                                                                                        aCTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgcccagccattagcgcattgccctggtctcgtaccacgtcccataaagcttcagcggcggtggaaagtgaac                  >  2:244514‑M1/1‑19 (MQ=255)
                                                                                                                         cTGCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgcccagccattagcgcattgccctggtctcgtaccacgtcccataaagcttcagcggcggtggaaagtgaacggttcttcctgcgcac  <  2:647829‑M1/139‑122 (MQ=255)
                                                                                                                           gCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgcccagccattagcgcattgcc                                                                    >  2:28017‑M1/1‑16 (MQ=255)
                                                                                                                           gCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttcccgcccagccattagcgcattgcc                                                                    <  1:28017‑M1/71‑56 (MQ=255)
                                                                                                                           gCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttc                                                                                            >  1:1318012‑M1/1‑16 (MQ=255)
                                                                                                                           gCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttc                                                                                            >  1:282433‑M1/1‑16 (MQ=255)
                                                                                                                           gCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttc                                                                                            <  2:282433‑M1/47‑32 (MQ=255)
                                                                                                                           gCCGGGGCTCCCGGCAcactttatatctgatacaacggatccccttc                                                                                            <  2:1318012‑M1/47‑32 (MQ=255)
                                                                                                                               gggCTCCCGGCAcactttatatctgatacaacggatccccttcccgcccagccattagcgcattgccctggtctcgtaccacgtcccataaagcttcagcggcggtgga                          >  2:727823‑M1/1‑12 (MQ=255)
                                                                                                                               gggCTCCCGGCAcactttatatctgatacaacggatccccttcccgcccagccattagcgcattgccctggtctcgtaccacgtcccataaagcttcagcggcggtgga                          <  1:727823‑M1/109‑98 (MQ=255)
                                                                                                                                    cccGGCAcactttatatctgatacaacggatccccttcccgcccagccattagcgcattgccctggtctcgtaccacgtcccataaagcttcagcggcggtggaaagtga                    >  1:1168142‑M1/1‑7 (MQ=255)
                                                                                                                                    cccGGCAcactttatatcggatacaacggatccccttcccgcccagccattagcgcattgccctggtctcgtaccacgtcccataaagcttcagcggcggtggaaagtga                    <  2:1168142‑M1/110‑104 (MQ=255)
                                                                                                                                          |                                                                                                                         
TGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTT  >  NC_000913/2525778‑2526037

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