Predicted mutation
evidence seq id position mutation freq annotation gene description
MC JC NC_000913 257,908 Δ776 bp 100% [crl] [crl]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 257908 258683 776 61 [0] [0] 62 [crl] [crl]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 2579070 (0.000)57 (0.900) 49/246 NT 100% intergenic (+8/‑769) crl/crl pseudogene, sigma factor‑binding protein, RNA polymerase holoenzyme formation stimulator;regulator; Surface structures; transcriptional regulator of cryptic csgA gene for curli surface fibers/pseudogene, sigma factor‑binding protein, RNA polymerase holoenzyme formation stimulator;regulator; Surface structures; transcriptional regulator of cryptic csgA gene for curli surface fibers
?NC_000913 258684 = 0 (0.000)pseudogene (9/331 nt) crl pseudogene, sigma factor‑binding protein, RNA polymerase holoenzyme formation stimulator;regulator; Surface structures; transcriptional regulator of cryptic csgA gene for curli surface fibers

CCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTC  >  NC_000913/258544‑258832
                                                                                                                                            |                                                                                                                                                    
ttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCaaa                                                                                                                                               <  5:103177‑M2/8‑1 (MQ=255)
    aaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactgggcccgtatattcgtgaaggtaAGTGCAAAGATaa                                                                                                                                          >  7:219585‑M2/137‑149 (MQ=255)
     aatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAAt                                                                                                                                         <  2:239253‑M2/14‑1 (MQ=255)
         tgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGAt                                                                                                                                     >  4:208593‑M2/132‑149 (MQ=255)
                ttttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCttttt                                                                                                                              <  3:293831‑M2/25‑1 (MQ=255)
                ttttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCttttt                                                                                                                              <  3:208593‑M2/25‑1 (MQ=255)
                     ccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGAtt                                                                                                                         <  8:377727‑M2/30‑1 (MQ=255)
                         tttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCt                                                                                                                     >  2:287843‑M2/116‑149 (MQ=255)
                            tgtaaaacagttgcatcacaacaggagatatcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGTAAAGATAATCGATTCTTTTTCGATTGTCTgg                                                                                                                   <  7:368033‑M2/36‑1 (MQ=255)
                               aaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTg                                                                                                                >  6:14062‑M2/110‑148 (MQ=255)
                                aaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTa                                                                                                              >  4:103728‑M2/109‑149 (MQ=255)
                                  acagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATg                                                                                                            >  6:297325‑M2/107‑149 (MQ=255)
                                               aacatgagatagcaatgacttttccgagtggacacccgaatagcagatttatcaaaaaatttaccgcactaggccagtatattcgtgaatgtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGGCAAc                                                                                                     <  6:375693‑M2/50‑1 (MQ=255)
                                               aacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCTATTGTCTGGCTGTATGCTTCAAc                                                                                                     >  5:375693‑M2/94‑143 (MQ=255)
                                               aacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGaaa                                                                                               <  8:257198‑M2/56‑1 (MQ=255)
                                                    gagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCgg                                                                                           >  8:239128‑M2/89‑148 (MQ=255)
                                                       atagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCttttt                                                                                                                              >  7:157875‑M2/86‑110 (MQ=255)
                                                       atagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCttttt                                                                                                                              <  8:157875‑M2/25‑1 (MQ=255)
                                                        tagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatatttgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCg                                                                                      >  6:342506‑M2/85‑149 (MQ=255)
                                                        tagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaagggaAGTGAAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCg                                                                                      >  5:288171‑M2/85‑149 (MQ=255)
                                                          gcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCgg                                                                                     >  4:337457‑M2/83‑148 (MQ=255)
                                                           caatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGaa                                                                                   <  7:234491‑M2/68‑1 (MQ=255)
                                                            aatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattagtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAg                                                                                  >  5:52356‑M2/81‑149 (MQ=255)
                                                                 cgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGt                                                                             <  3:103728‑M2/74‑1 (MQ=255)
                                                                 cgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGt                                                                             <  7:239128‑M2/74‑1 (MQ=255)
                                                                   ttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTa                                                                                                              >  2:63000‑M2/74‑114 (MQ=255)
                                                                   ttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTa                                                                                                              <  1:63000‑M2/41‑1 (MQ=255)
                                                                     accgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAAt                                                                         >  6:206503‑M2/72‑149 (MQ=255)
                                                                      ccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGAtt                                                                                                                                    <  8:280962‑M2/19‑1 (MQ=255)
                                                                       cgagtggacacccgaagagcagattgatcaaaaaatttaccgcactgggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTc                                                                       <  8:219585‑M2/80‑1 (MQ=255)
                                                                          gtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCt                                                                      >  6:332810‑M2/67‑147 (MQ=255)
                                                                          gtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCt                                                                      <  5:332810‑M2/81‑1 (MQ=255)
                                                                          gtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTgg                                                                    <  3:91682‑M2/83‑1 (MQ=255)
                                                                             gacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTc                                                                                                                      <  3:135170‑M2/33‑1 (MQ=255)
                                                                             gacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTc                                                                                                                      >  4:135170‑M2/64‑96 (MQ=255)
                                                                                  ccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTGTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCtggtg                                                            >  4:227011‑M2/59‑149 (MQ=255)
                                                                                  ccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCGCCGGAAGTGCGTGAATTCTGGGGCtggtg                                                            >  1:316755‑M2/59‑149 (MQ=255)
                                                                                  ccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCtggtg                                                            >  3:183537‑M2/59‑149 (MQ=255)
                                                                                  ccgaagagcagattgatcaaaaaatttaccgcactaggcccgtagattcgtgaaggtaAGTGCAAAGATAATCGCTTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCtggtg                                                            >  3:272739‑M2/59‑149 (MQ=255)
                                                                                          cagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGc                                                    >  6:207563‑M2/51‑149 (MQ=255)
                                                                                             attgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTgg                                                 <  3:37855‑M2/102‑1 (MQ=255)
                                                                                                gatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGaa                                               <  5:14062‑M2/104‑1 (MQ=255)
                                                                                                 atcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAgc                                             >  8:113901‑M2/44‑149 (MQ=255)
                                                                                                  tcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAgc                                             <  5:207563‑M2/106‑1 (MQ=255)
                                                                                                    aaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCa                                          <  1:257899‑M2/109‑1 (MQ=255)
                                                                                                     aaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTAt                                                                                                             >  2:262824‑M2/40‑81 (MQ=255)
                                                                                                     aaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTAt                                                                                                             <  1:262824‑M2/42‑1 (MQ=255)
                                                                                                     aaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGGTAGTAGCGCAg                                         <  5:297325‑M2/110‑1 (MQ=255)
                                                                                                          tttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATc                                    <  4:183537‑M2/115‑1 (MQ=255)
                                                                                                              ccgcacgaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGt                                >  2:256570‑M2/31‑149 (MQ=255)
                                                                                                                   ctaggcccgtatatttgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTAc                           <  5:342506‑M2/124‑1 (MQ=255)
                                                                                                                      ggcccgtatattcgtgaaggtaACTGCAAAGAAAATCGATTCTTTTTCGATTGTCTGGCTGTATGAGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGAAGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTa                        >  7:219426‑M2/23‑149 (MQ=255)
                                                                                                                       gcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGATTTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTAc                           <  2:16115‑M2/124‑1 (MQ=255)
                                                                                                                       gcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTAc                           >  1:16115‑M2/22‑145 (MQ=255)
                                                                                                                        cccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACa                      <  7:113901‑M2/129‑1 (MQ=255)
                                                                                                                          cgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGt                    >  1:279851‑M2/19‑149 (MQ=255)
                                                                                                                            tatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTa                  >  7:319934‑M2/17‑149 (MQ=255)
                                                                                                                               attcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTAc                 >  6:318582‑M2/14‑147 (MQ=255)
                                                                                                                                       aggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGtttt                             <  6:104085‑M2/122‑1 (MQ=255)
                                                                                                                                       aggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGtttt                             >  5:104085‑M2/6‑127 (MQ=255)
                                                                                                                                       aggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTc       <  6:52356‑M2/144‑1 (MQ=255)
                                                                                                                                            aGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTc  <  1:287843/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
CCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTC  >  NC_000913/258544‑258832

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