Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 3,361,226 Δ9,811 bp [gltF][nanE] [gltF], yhcA, yhcD, yhcE, insH10, yhcE, yhcF, yhcG, nanQ, nanK, [nanE]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3361226 3371036 9811 36 [0] [0] 38 [gltF]–[nanE] [gltF],yhcA,yhcD,yhcE,insH10,yhcE,yhcF,yhcG,nanQ,nanK,[nanE]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 33612250 (0.000)36 (0.550) 35/490 0.5 100% coding (50/765 nt) gltF periplasmic protein GltF
?NC_000913 3371037 = 0 (0.000)intergenic (‑1/+47) nanE/nanT putative N‑acetylmannosamine‑6‑phosphate 2‑epimerase/N‑acetylneuraminate:H(+) symporter

CACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGAACGTCAGACTGAACGAGAGCGATGCCAGCGCAGTACCCA  >  NC_000913/3370842‑3371283
                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
cgtgactataaaattcattaaataaatacatcccagcaagttggctatttctgaacgagttatttatttacattcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACttt                                                                                                                                                                                                  <  1:156194‑M2/55‑1 (MQ=255)
   aactataaaattcattaaataaatacatcccagcaagttggctatttctgaacgagttatttatttacattcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTgg                                                                                                                                                                                               <  2:286652‑M2/58‑1 (MQ=255)
         aaaattcattaaagaaatacatcccagcaagttggctatttctgaacgagttatttatttacattcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcAACTCTTCCCTTAGCGACAGGCACGGTACATAGACCGGGCAACAGGATTTACTTCTGGTTTTGAc                                                                                                                                                                                        >  2:533022‑M2/187‑251 (MQ=255)
                                   gcaagttggctatttctgaacgagttatttatttacattcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACATTTGGTTTTGACTACATCGTTTTTGGCGTTGCCACacg                                                                                                                                                              >  1:530302‑M2/161‑251 (MQ=255)
                                              atttctgaacgagttatttatttacattcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCg                                                                                                                                                   >  1:465202‑M2/150‑251 (MQ=255)
                                                                   tttcattcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCGTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCTCTGCCAAACGGCGCGGCACCGCTGAATGGTTTACCGTCGAt                                                                                                                               <  2:60742‑M2/122‑1 (MQ=255)
                                                                      cattcacaaactgatatatatactgcttcagctttttattttaaaaaataaatttgatattaaggatatgttcatatgtttttcaaaaagaaccacacaacagctgctatttgtgcagcactttcACACCTTCCCTTAGCGAAAGCCCCGGTACATAGACCGGGCAACAGGATTATCTTTTGGTTTTGACTAAATCGTTGTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTaccgtcgttatcg                                                                                                                           <  2:139745‑M2/126‑7 (MQ=255)
                                                                            caaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATg                                                                                                                      >  1:28982‑M2/120‑250 (MQ=255)
                                                                                     atatatactgcttcaactttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagccgctatttgtgcagctctttcTCAGCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAAc                                                                                                            <  2:94902‑M2/141‑1 (MQ=255)
                                                                                       atatactgcttcagctttttattttaaaaaataaatttgctattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGc                                                                                                          >  1:364821‑M2/109‑251 (MQ=255)
                                                                                                  cagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGCCTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGgcgc                                                                                                >  2:436670‑M2/98‑250 (MQ=255)
                                                                                                           attttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGATGCCAAACAGCACGGCACCGCTGAATGGTTTACAGTCGATAGCGTCATGAGTACGCACCGCTTACGGGCGCAACAAACGCt                                                                                      >  2:98988‑M2/89‑251 (MQ=255)
                                                                                                             tttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGcaacca                                                                                           <  1:352455‑M2/158‑1 (MQ=255)
                                                                                                             tttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGcaacca                                                                                           >  2:352455‑M2/87‑244 (MQ=255)
                                                                                                                       aaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCgaga                                                                            >  2:146020‑M2/77‑249 (MQ=255)
                                                                                                                        aatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAgg                                                                         >  2:4677‑M2/76‑251 (MQ=255)
                                                                                                                          tttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACg                                                                                        <  2:34165‑M2/161‑1 (MQ=255)
                                                                                                                          tttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACg                                                                                        >  1:34165‑M2/74‑234 (MQ=255)
                                                                                                                               tattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGTTAGCGTCATGGGCACGCAACCCTTCCGGGGGCAACCAACGCTGAACGCGAGAAGGCATATcc                                                                  >  2:492277‑M2/69‑251 (MQ=255)
                                                                                                                                   aaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGcc                                                              >  1:63306‑M2/65‑251 (MQ=255)
                                                                                                                                       atatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAAt                                                          <  1:342863‑M2/191‑1 (MQ=255)
                                                                                                                                          tgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATc                                                         >  2:30539‑M2/58‑249 (MQ=255)
                                                                                                                                                  tgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATc                                                >  2:532214‑M2/50‑250 (MQ=255)
                                                                                                                                                     ttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATcacc                                             <  1:265351‑M2/204‑1 (MQ=255)
                                                                                                                                                     ttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATcacc                                             >  1:269298‑M2/47‑250 (MQ=255)
                                                                                                                                                                 acctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGAACGTCAGAc                                >  1:122764‑M2/35‑251 (MQ=255)
                                                                                                                                                                        aacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACTCGAGAAGGCATATCCAGCCCAAGCAGCAGGATCACCACGAACGTCAGACTGAACg                          >  2:368321‑M2/28‑250 (MQ=255)
                                                                                                                                                                              tgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAg                                                                                                                             <  2:18245‑M2/124‑1 (MQ=255)
                                                                                                                                                                              tgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAg                                                                                                                             >  1:18245‑M2/22‑145 (MQ=255)
                                                                                                                                                                               gctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGAACGTCAGACTGAACGAGAGCGAt                  <  1:146020‑M2/231‑1 (MQ=255)
                                                                                                                                                                                  atttgtgcagcccttcaacatCTTCCCCTATAGAAAGGCCCGGTCCATCGACCGGGCAACAGGTTTATCTTTTTGTTTTGGCCAAACCGTTTATGGCGCTGCCAAACGGCACGGAACCGCTGAATGGTTTACCGTCGATAGCGTCAAGAGTACGCAACGCTGCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGTATCACCACGATCGTCAGCCTGAACGAGAGCGATGc                <  2:1034/233‑1 (MQ=255)
                                                                                                                                                                                      gtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGAACGTCAGACTGAACGAGAGCGATGCCAgcg           >  1:344585‑M2/14‑251 (MQ=255)
                                                                                                                                                                                            cactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAAACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGAACGTCAGACTGAACGAGAGCGATGCCAGCGCAGTAc     >  1:282324‑M2/8‑251 (MQ=255)
                                                                                                                                                                                             actttcACATCTTCCCTTAGCGAAAGGCCCGGTACATCGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAACCGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCCATCAGCAGGATCACCACGAACGTCAGACTGAACGAGAGCGATGCCAGCGCAGTAcc    >  1:244597‑M2/7‑251 (MQ=255)
                                                                                                                                                                                              ctttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCGGGATCACCACGAACGTCAGACTGAACGAGAGCGATGCCAGCGCAGTAccc   >  1:385602‑M2/6‑251 (MQ=255)
                                                                                                                                                                                              ctttcACATCTACCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAGCTTTTGGGTTTGACTAAATCGTTTTTTGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGATTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATACAGCCCAATCAGCAGGATCACCACTAACGTCAGACTGAACGAGAGCGATGCCAGCGCAGTAc     <  2:164034‑M2/244‑1 (MQ=255)
                                                                                                                                                                                                ttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGAACGTCAGACTGAACGAGAGCGATGCCAGCGCAGTACCCa  <  1:30539‑M2/247‑1 (MQ=255)
                                                                                                                                                                                                   acatCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGAACGTCAGACTGAACGAGAGCGATGCCAGCGCAGTAccc   >  1:335184/1‑246 (MQ=255)
                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
CACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGAACGTCAGACTGAACGAGAGCGATGCCAGCGCAGTACCCA  >  NC_000913/3370842‑3371283

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