Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 2,558,699 Δ6,790 bp intZ[eutA] intZ, yffL, yffM, yffN, yffO, yffP, yffQ, yffR, yffS, [eutA]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 2558699 2565488 6790 80 [0] [6] 86 intZ–[eutA] intZ,yffL,yffM,yffN,yffO,yffP,yffQ,yffR,yffS,[eutA]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 25586980 (0.000)79 (1.290) 71/468 0.1 100% intergenic (‑19/‑160) eutB/intZ ethanolamine ammonia‑lyase subunit alpha/putative phage integrase IntZ
?NC_000913 2565489 = 0 (0.000)coding (1396/1404 nt) eutA ethanolamine ammonia‑lyase reactivase EutA

TGAATCAGCCGCGTCACGCAGTCATCTTCATAGGCAATCACCGGATTATTGCGGATGTCCGCTACGGTCATTTCCGACAACACCTGCTTTGCCGCCACGCGCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTGGCATTTTTCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATT  >  NC_000913/2558450‑2558933
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                           
tGAATCAGCCGCGTCACGCAGTCATCTTCATAGGCAATCACCGGATTATTGCGGATGTCCGCTACGGTCATTTCCGACGACACCTGCTTTGCCGCCACGCGCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTACCTCAGGAAGg                                                                                                                                                                                                                                            >  1:117396‑M1/1‑249 (MQ=255)
  aaTCAGCCGCGTCACGCAGTCATCTTCATAGGCAATCACCGGATTATTGCGGATGTCCGCTACGGTCATTTCCGACAACACCTGCTTTGCCGCCACGCGCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaa                                                                                                                                                                                                                                          >  2:624397‑M1/1‑247 (MQ=255)
   aTCAGCCGCGTAACGCAGTCATCTTCATAGGCAATCACCGGATTATTGCGGCTGTCCGCTACGGTCATTTCCGACAACACCTGCTTTGCCGCCACGCGCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaat                                                                                                                                                                                                                                        >  2:309914‑M1/1‑246 (MQ=255)
     cAGCCGCGTCACGCAGTCATCTTCATAGGCAATCACCGGATTATTGCGGATGTCCGCTACGGTCATTTCCGACAACACCTGCTTTGCCGCCACGCGCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgc                                                                                                                                                                                                                                      >  1:347586‑M1/1‑244 (MQ=255)
     cAGCCGCGTCACGCAGTCATCTTCATAGGCAATCACCGGATTATTGCGGATGTCCGCTACGGTCATTTCCGACAACACCTGCTTTGCCGCCACGCGCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgc                                                                                                                                                                                                                                      >  2:231376‑M1/1‑244 (MQ=255)
       ccgtcgccACGCAGTCATCTTATTAGGCAAACACCGTATTTTTGCGGACGTGCGCTACGTTCATTCCCGATAACCCCTGCTTTGCCCCCCCGCCCTCTGGTGACCTGGCCGCTGCGCCGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGgaaggaaatgcga                                                                                                                                                                                                                                    <  2:132841/244‑10 (MQ=255)
           cgTCACTCAGTCATCTTCATAGGCAATCACCGGATTATTGCGGATGTCCGCTACGGTCATTTCCGACAACACCTGCTTTGCCGCCACGCGCTCCTGTGAGCTTGCCGCTGCGACACCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACATCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtg                                                                                                                                                                                                                                 >  1:523858‑M1/1‑238 (MQ=255)
            gTCACGCAGTCATCTTCATAGGCAATCACCGGATTATTGCGGATGTCCGCTACGGTCATTTCCGACAACACCTGCTTTGCCGCCACGCGCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGAACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGGCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgat                                                                                                                                                                                                                               >  1:141349‑M1/1‑237 (MQ=255)
                  cAGTCATCTTCATAGGCATCCCCCGTATTATGGCGGATGCCCGCTACGTTCATTTCCAACAACACCTGCTTTGCCCCCACGCCCTCCTGTGAGCTTGCCGCTGCGCCGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcacc                                                                                                                                                                                                                         <  1:309914‑M1/251‑21 (MQ=255)
                    gTCATCTTCATAGGCAATCACCGGATTATTGCGGATGTCCGCTACGGTCATTTCCGACAACACCTGCTTTGCCGCCACGCGCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgt                                                                                                                                                                                                                       >  1:635535‑M1/1‑229 (MQ=255)
                     tCATCTTCATAGGCAATCACCGGATTATTGCGGATGTCCGCTACGGTCATTTCCGACAACACCTGCTTTGCCGCCACGCGCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtc                                                                                                                                                                                                                      >  1:332534‑M1/1‑228 (MQ=255)
                     tCATCTTCATAGGCAATCACCGGATTATTGCGGATGTCCGCTACGGTCATTTCCGACAACACCTGCTTTGCCGCCACGAGCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatgtcaccgt                                                                                                                                                                                                                       >  1:88298‑M1/1‑228 (MQ=255)
                      cATCTTCATAGCCACTCACCGGATTTTTGCGGATGTCCGCTCCGGTCTTTTCCGACAACCCCTGCTTTGCCCCCCCGCGCTCCTGTGAGCTTGCCGCTGCGCCGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtca                                                                                                                                                                                                                     <  2:564988‑M1/251‑25 (MQ=255)
                         cTTCATAGGCAATCACCGGATTATTGCGGATGTCCGCTACGGTCATTTCCGACAACACCTGCTTTGCCGCCACGCGCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccg                                                                                                                                                                                                                  >  1:191308‑M1/1‑224 (MQ=255)
                             aTAGGCAATCACCGGATTATTGCGGATGTCCGCTACGGTCATTTCCGACAACACCTGATTTGCCGCCACGCGCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcac                                                                                                                                                                                                              >  2:38030‑M1/1‑220 (MQ=255)
                             aTAGCCAATCACCGGTTTATTGCTGATGTCCGCTACGGTCATTTCCGACAACGCCGGCTTTGCCGCCCCGCGCCCCTGTGATCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggca                                                                                                                                                                                                               <  1:186598‑M1/250‑31 (MQ=255)
                               aGGCAATCACCGGATTATTGCGGATGTCCGCTACGGTCATTTCCGACAACACCTGCTTTGCCGCCACGCGCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcaca                                                                                                                                                                                                             >  1:5352‑M1/1‑218 (MQ=255)
                                 gCAATCACCGGATTATTGCGGTTGTCCGCTCCGGTCCTTTCCGACACCCCCTTCTTTGCCGCCCCGCGCTCCTTTGAGCTGGCCCCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaacc                                                                                                                                                                                                          <  1:403822‑M1/251‑36 (MQ=255)
                                      cACCGGATTATTGCGGATGTCCGCTACGGTCATTTCCGACAACACCTGCTTTGCCGCCACGCGCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatcc                                                                                                                                                                                                     >  1:356053‑M1/1‑211 (MQ=255)
                                               attGCGGATGTCCGCTACGGTCATTTCCGACAACACCTGCTTTGCCGCCACGCGCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaa                                                                                                                                                                                              >  2:467343‑M1/1‑202 (MQ=255)
                                                  gcagatgTCCGCTACGGTCATTTCCGACAACCCCTGCTTAGCCGCCCCGCGCTCCTGGGAGCTGGCCGGTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagag                                                                                                                                                                                          <  1:172744‑M1/247‑52 (MQ=255)
                                                   gggttgtccGCTAGGGTCTTTTCCAACAACACCTGTTTCGCCCCCACGCGCTCCGGGGAGCTTGCCCCTGCCACACCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggc                                                                                                                                                                                        <  1:6230‑M1/247‑54 (MQ=255)
                                                      aTGTCCGCTACGGTCATTTCCGACAACACCTGCTTTGCCGCCACGCGCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgt                                                                                                                                                                                      >  1:132704‑M1/1‑195 (MQ=255)
                                                           cGCTACGGTCATTTCCGACAACACCTGCTTTGCCGCCACGCGCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaat                                                                                                                                                                                >  2:526749‑M1/1‑190 (MQ=255)
                                                               ccgGTCATTTCCGACAACCCCGGCTTTGCCCCCACGCCCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcg                                                                                                                                                                            <  1:544422‑M1/250‑66 (MQ=255)
                                                                gggtCATTTCCGACACCCCCTGCTTGGCCGCCACGCGTTCCTGGGGGCTGCCCGCTGCGCCGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCGCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcga                                                                                                                                                                           <  1:602562‑M1/250‑67 (MQ=255)
                                                                     aTTTCCGACAACACCTGCTTTGCCGCCACGCGCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccg                                                                                                                                                                                                    >  1:203464‑M1/1‑180 (MQ=255)
                                                                          cGACAACACCTGCTTTGCCGCCACGCGCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccc                                                                                                                                                                 >  2:332894‑M1/1‑175 (MQ=255)
                                                                             caacaCCTGCTTTGCCGCCACGCGCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgc                                                                                                                                                              >  2:176806‑M1/1‑172 (MQ=255)
                                                                               acacCTGCTTTGCCGCCACGCGCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgc                                                                                                                                                            >  1:262054‑M1/1‑170 (MQ=255)
                                                                                        ttGCCGCCCCGCGCTCCTGTGAGCTTGCCGCTGCGCCGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaat                                                                                                                                                     <  1:566438‑M1/249‑89 (MQ=255)
                                                                                                    gcTCCTGTGAGCTTGCCGCTGCGACACCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatacgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatg                                                                                                                                        >  2:419702‑M1/1‑149 (MQ=255)
                                                                                                          gtggtgttgccgctgggtccaccgccagcacatcccccgaacgcagtTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgcca                                                                                                                                 <  2:396804‑M1/230‑109 (MQ=255)
                                                                                                            gAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcaca                                                                                                                                                                                                             >  2:42626‑M1/1‑141 (MQ=255)
                                                                                                            gAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcaca                                                                                                                                                                                                             <  1:42626‑M1/173‑33 (MQ=255)
                                                                                                             aGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacg                                                                                                                              >  1:455559‑M1/1‑140 (MQ=255)
                                                                                                               cTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacggg                                                                                                                            >  2:240880‑M1/1‑138 (MQ=255)
                                                                                                                       tGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggagttgttg                                                                                                                    >  1:105538‑M1/1‑130 (MQ=255)
                                                                                                                          gccgcccgcctgCACATTCCCCGAACGCAGTTCGTTTGCTTTAGCCAGCACTTCTTTTACATCCTTAAACTTATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgtgtgatttcaccgtcaccggcacatccgatcctccaaaaataagcgtaccaatgtcaatatagtcccccgcgcggacaatcacttcctcaatgactgccaacgggttttgttg                                                                                                                    <  2:576326‑M1/246‑122 (MQ=255)
                                                                                                                          gACGCACGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaacgaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgtctgccaacgggagttgttgt                                                                                                                   >  2:263059‑M1/1‑127 (MQ=255)
                                                                                                                           ccgcccgccaGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggagttgttgtagc                                                                                                                <  2:560220‑M1/250‑126 (MQ=255)
                                                                                                                                 gccAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggagttgttgtagctgtggg                                                                                                          >  1:4348‑M1/1‑120 (MQ=255)
                                                                                                                                     gCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggagttgttgtagctgtgggcgca                                                                                                      <  2:312822‑M1/251‑136 (MQ=255)
                                                                                                                                       acaTCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggagttgttgtagctgtgggcgcaac                                                                                                    <  1:546571‑M1/251‑138 (MQ=255)
                                                                                                                                               cGAACGCAGATCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACGCTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGAAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggagtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggagttgttgtagctgtgggcgcaacaacatgcc                                                                                            >  2:296173‑M1/1‑106 (MQ=255)
                                                                                                                                                gAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggagttgttgtagctgtgggcgcaacaacatgccc                                                                                           >  1:51384‑M1/1‑105 (MQ=255)
                                                                                                                                                    gCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggagttgttgtagctgtgggcgcaacaacatgcccag                                                                                         >  1:146396‑M1/1‑101 (MQ=255)
                                                                                                                                                     cAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggagttgttgtagctgtgggcgcaacaacatgcccagagc                                                                                      <  2:4348‑M1/251‑152 (MQ=255)
                                                                                                                                                     cAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatata                                                                                                                                                                       >  1:555247‑M1/1‑100 (MQ=255)
                                                                                                                                                     cAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatata                                                                                                                                                                       <  2:555247‑M1/170‑71 (MQ=255)
                                                                                                                                                       gTTCGTTGGCTTTAGCCAGCACCTCTTTTGCATCCTTAAACTGATATACATGGCCGAACAATTTGGTCTTTAGTTTCAGAAGTCGTTCCCTCAGGAAGgaaatgcgattgatttcaccgtcacctgcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgtctgccaacgggagttgttgtagctgtgggcgcaacaacatgcccagagctt                                                                                    <  2:130107‑M1/251‑154 (MQ=255)
                                                                                                                                                       gTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgcaaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggagttgttgtagctgtgggcgaaacaacatgcccagagctt                                                                                    >  1:173453‑M1/1‑98 (MQ=255)
                                                                                                                                                       gTTCGTGGGATTTAGCCAGCAACTCTTTTCCATCCTTAAACGGAGATACATGGGTGAATAAGGTGGTCTTGAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagggagttcaccgtcaccggcaaaaccgatctgccaaaaagaggcgtaccaatttcgatatagtcccccgcgcggacaatcacttcgtcaatgactcccaacgggagttgttgtagctgtgggcgcaacaacatgcccagagctt                                                                                    <  1:609943‑M1/251‑154 (MQ=255)
                                                                                                                                                           gTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggagttgttgtagctgtgggcgcaacaacatgcccagagctttacc                                                                                >  1:523963‑M1/1‑94 (MQ=255)
                                                                                                                                                               gCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgctatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggagttgttgtagctgtgggcgcaacaccatgcccagagctttaccaa                                                                              >  2:5281‑M1/1‑90 (MQ=255)
                                                                                                                                                                      ccAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggagttgttgtagctgtgggcgcaacaacatgcccagagctttaccaaagtcctgcc                                                                     >  1:45400‑M1/1‑83 (MQ=255)
                                                                                                                                                                          cACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggagttgttgtagctgtgggcgcaacaacatgcccagagctttaccaaagtcctgcccggc                                                                 >  2:554496‑M1/1‑79 (MQ=255)
                                                                                                                                                                           aCCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggagttgttgtagctgtgggcgcaacaacatgcccagagctttaccaaagtcctgcccggcc                                                                >  2:348717‑M1/1‑78 (MQ=255)
                                                                                                                                                                                    aCATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggagttgttgtagctgtgggcgcaacaacatgcccagagctttaccaaagtcctgcccggccaccaccagc                                                       >  1:459989‑M1/1‑69 (MQ=255)
                                                                                                                                                                                      aTCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatttcgatatagtcccccgcgctgacaatcccttcgtcaatgactgccaccgggagttgttgtagctgtgggctcaacaccatgcccagagctttaccaaagtcctgcccgtccgccgccagcag                                                     <  2:88298‑M1/251‑185 (MQ=255)
                                                                                                                                                                                              cctGATATACATTGCCGCACAATGTGGTCTTTAGTTTCTTAAGTCGTTCCCTCAGGAAGgaaatgcgcgtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccattgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactcccaacgggagttgttgtagctgggggcgcaactccatgcccagagctttaccaaagtcctgcccgtccaccaccagcaccggtggcg                                             >  2:9627‑M1/2‑59 (MQ=255)
                                                                                                                                                                                               cTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggagttgttgtagctgtgggcgcaacaacatgcccagagctttaccaaagtcctgcccggccaccaccagcaggggatgcgg                                            >  1:114864‑M1/1‑58 (MQ=255)
                                                                                                                                                                                                 gATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatggctgccaacgggagttgttgtagctgtgggcgcaacaacatgcccagagctttaccaaagtcctgcccggccaccaccagcagtggatgcgg                                            <  2:73885‑M1/249‑194 (MQ=255)
                                                                                                                                                                                                    ataCATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggagttgttgtagctgtgggcgcaacaacatgcccagagctttaccaaagtcctgcccggccaccaccagcaggggatgcggattcg                                       <  1:624397‑M1/251‑199 (MQ=255)
                                                                                                                                                                                                     taCATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggagttgttgtagctgtgggcgcaacaacatgcccagagctttaccaaagtcctgcccggccaccaccagcaggggatgcggattcgg                                      >  1:121539‑M1/1‑52 (MQ=255)
                                                                                                                                                                                                      aCATTGCCGAACAATGTGGTCTTTAGTTTCAGAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggagttgttgtagctgtgggcgcaacaacatgcccagagctttaccaaagtcctgcccggccaccaccagcaggggatgcggattcaga                                     <  1:616717‑M1/251‑201 (MQ=255)
                                                                                                                                                                                                       cATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggagttgttgtagctgtgggcgcaacaacatgcccagagctttaccaaagtcctgcccggccaccaccagcaggggatgcggattcggaa                                    <  1:296173‑M1/251‑202 (MQ=255)
                                                                                                                                                                                                        aTTNCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggagttgttgtagctgtgggcgcaacaacatgcccagagcgttaccaaagtcctgcccggccaccaccagcaggggatgcggattcggaaa                                   <  2:601705‑M1/251‑203 (MQ=255)
                                                                                                                                                                                                        aTTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggagttgttgtagctgtgggcgcaacaacatgcccagagctttaccaaagtcctgcccggccaccaccagcaggggatgcggattcggaaa                                   <  1:654257‑M1/251‑203 (MQ=255)
                                                                                                                                                                                                            ccGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggagttgttgtagctgtgggcgcaacaacatgcccagagctttaccaaagtcctgcccggccaccaccagcaggggatgcggattcggaaaacgc                               >  1:529604‑M1/1‑45 (MQ=255)
                                                                                                                                                                                                              gAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggagttgttgtagctgtgggcgcaacaacatgcccagagctttaccaaagtcctgcccggccaccaccagcaggggatgcggcttcggaaaacgcgc                             >  1:170651‑M1/1‑43 (MQ=255)
                                                                                                                                                                                                                aCAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggagttgttgtagctgtgggcgcaacaacatgcccagagctttaccaaagtcctgcccggccaccaccagcaggggatgcggattcggaaaacgcgcg                            >  2:559292‑M1/1‑41 (MQ=255)
                                                                                                                                                                                                                 cAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgtctgccaacgggagttgttgtagctgtgggcgcaacaacatgcccagagctttaccaaagtcctgcccggccaccaccagcaggggatgcggattcggaaaacgcgcgac                          <  1:263059‑M1/251‑212 (MQ=255)
                                                                                                                                                                                                                       ggtctttagatacctcagtcgttccctcaggaaGgaaatgcgagagattacaccgtcaccggcccaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggcgttgttgtagctgtgggcgcaacaacatgcccacagctttagcaaagtcctgcccggccaacacaagcaggggaggcggattcggaaaacgcgcgacg                         >  2:649342‑M1/17‑34 (MQ=255)
                                                                                                                                                                                                                          cTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggagttgttgtagctgtgggcgcaacaacatgcccagagctttaccaaagtcctgcccggccaccaccagcaggggatgcggattcggaaaacgcgcgacgaaatcgac                 >  1:131308‑M1/1‑31 (MQ=255)
                                                                                                                                                                                                                               gTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggagttgttgtagctgtgggcgcaacaacatgcccagagctttaccaaagtcctgcccggccaccaccagcaggggatg                                               <  1:308426‑M1/216‑191 (MQ=255)
                                                                                                                                                                                                                               gTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggagttgttgtagctgtgggcgcaacaacatgcccagagctttaccaaagtcctgcccggccaccaccagcaggggatg                                               >  2:308426‑M1/1‑26 (MQ=255)
                                                                                                                                                                                                                                      aaGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggagttgttgtagctgtgggcgcaacaacatgcccagagctttaccaaagtcctgcccggccaccaccagcaggggatgcggattcggaaaacgcgcgacgaaatcgaccaacgcgttgat     >  1:159429‑M1/1‑19 (MQ=255)
                                                                                                                                                                                                                                        gTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggagttgttgtagctgtgggcgcaacaacatgcccagagctttaccaaagtcctgcccggccaccaccagccggggatgcggattcggaaaacgcgcgacgaaatcgaccaacgcgttgatg    >  2:495378‑M1/1‑17 (MQ=255)
                                                                                                                                                                                                                                          cGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaagaggcgtaccaatgtcgatatagtcccccgcgcggacaatcacttcgtcaatgactgccaacgggagttgttgtagctgtgggcgcaacaacatgcccagagctttaccaaagtcctgcccggccaccaccagcaggggatgcggattcggaaaacgcgcgacgaaatcgaccaacgcgttgatgac  <  1:231376‑M1/250‑236 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                           
TGAATCAGCCGCGTCACGCAGTCATCTTCATAGGCAATCACCGGATTATTGCGGATGTCCGCTACGGTCATTTCCGACAACACCTGCTTTGCCGCCACGCGCTCCTGTGAGCTTGCCGCTGCGACGCCCGCCAGCACATCCCCCGAACGCAGTTCGTTGGCTTTAGCCAGCACCTCTTTTACATCCTTAAACTGATATACATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTGGCATTTTTCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATT  >  NC_000913/2558450‑2558933

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