Predicted mutation
evidence seq id position mutation annotation gene description
MC JC CP009273 4,145,501 Δ320 bp coding (573‑892/1005 nt) argC → N‑acetyl‑gamma‑glutamylphosphate reductase, NAD(P)‑binding

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ CP009273 4145501 4145820 320 61 [0] [0] 61 argC N‑acetyl‑gamma‑glutamylphosphate reductase, NAD(P)‑binding

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? CP009273 = 41455000 (0.000)57 (1.090) 45/278 0.2 100% coding (572/1005 nt) argC N‑acetyl‑gamma‑glutamylphosphate reductase, NAD(P)‑binding
?CP009273 4145821 = 0 (0.000)coding (893/1005 nt) argC N‑acetyl‑gamma‑glutamylphosphate reductase, NAD(P)‑binding

AAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATCGCCACACACCTCGGTGCTGACGTTATCTTCACCCCACATCTGGGCAATTTCCCGCGCGGCATTCT  >  CP009273/4145352‑4145647
                                                                                                                                                    |                                                                                                                                                   
aaaTTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACTAGCGGGGTGAGCGGTGCAGGGCG                                                                                                                                                     >  1:904142‑M1/1‑149 (MQ=255)
 aaTTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGa                                                                                                                                                    <  1:870824‑M1/149‑2 (MQ=255)
   ttAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGGTGTCTACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagc                                                                                                                                                  >  1:630208‑M1/1‑146 (MQ=255)
    tAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagca                                                                                                                                                 >  1:129003‑M1/1‑145 (MQ=255)
        agaagCGGATTGGAGGGCGGGGCCGGGGTGGTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAAGCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagc                                                                                                                                                  <  2:797298‑M1/144‑4 (MQ=255)
             cGAATTTGATTGCGGTGCCGGGCTGTTACCAGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattat                                                                                                                                        <  1:863528‑M1/149‑14 (MQ=255)
                 tttGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtg                                                                                                                                    <  1:539012‑M1/149‑18 (MQ=255)
                 tttGATTCCGGTGCCGGGCTGTTATCCGACGGTGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtg                                                                                                                                    <  1:474031‑M1/149‑18 (MQ=255)
                    gATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgcttattg                                                                                                                                      <  2:542065‑M1/144‑16 (MQ=255)
                      ttGCGGTGCCAGGCTGTAAACCGAAAGCGGAACAGCTGGAGCTAAAACCATTGATTGATGCCGATCTTCTTGACCTCACTCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCTGTGCAGGGCGagcatctgattattgtggcgac                                                                                                                               <  1:886973‑M1/149‑23 (MQ=255)
                             gCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccg                                                                                                                             <  2:314064‑M1/144‑25 (MQ=255)
                              tcgGGCTGTTATCGGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagaca                                                                                                                       <  1:53500‑M1/148‑31 (MQ=255)
                                 ggCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagacaact                                                                                                                    <  1:128456‑M1/149‑34 (MQ=255)
                                 ggCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaaga                                                                                                                         <  2:315023‑M1/144‑29 (MQ=255)
                                  gCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagac                                                                                                                        <  2:263625‑M1/144‑30 (MQ=255)
                                   cTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGCCCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagacaaatta                                                                                                                  >  1:243885‑M1/1‑114 (MQ=255)
                                    tGTTATCCGACGGCGGCACAGCTGGCGCGGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCCACGCAACGAGCGGCATGCATGATCCAGGGCGagcatctgattattgtgtcgaccgaagacaacttac                                                                                                                 <  1:681010‑M1/149‑37 (MQ=255)
                                         tCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcgtctgagtattgtggcgaccgaaggcaacttac                                                                                                                 >  2:403826‑M1/1‑108 (MQ=255)
                                                 cgggcccgcgggcgctgaaaccGTTGATTGATGCCGATCTCCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCACCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagacaacttactgaaaggcgcggc                                                                                                    <  1:67874‑M1/139‑50 (MQ=255)
                                                   gcACAGCTGGCGCTGAAACCGTTGAATGATACCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCCACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagacaacttactgaaaggcgcggcgg                                                                                                  >  1:481619‑M1/1‑98 (MQ=255)
                                                    cacaGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagacaacttactgaaaggcgcggcggc                                                                                                 >  1:927018‑M1/1‑97 (MQ=255)
                                                      caGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCCGTGGCCGGTGATCAACGCTACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagacaacttactgaaaggcgcggc                                                                                                    >  2:820076‑M1/1‑95 (MQ=255)
                                                       aGCTGGCGCTGAAACCGTGGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagacaacttactgaaaggcgcggcg                                                                                                   <  2:793425‑M1/144‑51 (MQ=255)
                                                         cTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagacaacttactgaaaggcgcggcggc                                                                                                 >  2:212597‑M1/1‑92 (MQ=255)
                                                          tGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagacaacttactgaaaggcgcggcggcacaagc                                                                                           >  1:232664‑M1/1‑91 (MQ=255)
                                                               cTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagacaacttactgcaaggggcggcggcccaagc                                                                                           >  2:207386‑M1/1‑86 (MQ=255)
                                                                tGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCCGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagacaacttactgaaaggcgcggcggcacaagcg                                                                                          <  2:840553‑M1/144‑60 (MQ=255)
                                                                tGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagaaagtttttttgttggcgcggcggcacaagcggtaca                                                                                     <  1:225307‑M1/149‑65 (MQ=255)
                                                                tGAAACCGTTGATCGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagacaacttactgaaaggcgcggcggcacaagcggtaca                                                                                     <  1:268003‑M1/149‑65 (MQ=255)
                                                                cgaAACCGCTGATTGATGCCTATCTTCTCGCCCTCAATCAGTGGCCGGTGACCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgcacggtctcaaccgactgaaaggcgcggcggcacaagcggtaca                                                                                     <  1:224156‑M1/148‑65 (MQ=255)
                                                                    aCCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagacaacttactgaaaggcgcggcggcacaagcggtacagtgc                                                                                 >  1:462053‑M1/1‑81 (MQ=255)
                                                                             tgatGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagacaacttactgaaaggcgcggcggcacaagcggtacagtgcgccaatatt                                                                        >  1:925896‑M1/1‑72 (MQ=255)
                                                                             tgatGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagacaacttactgaaaggcgcggcggcacaagcggtacagtgcgccaatatt                                                                        >  1:608869‑M1/1‑72 (MQ=255)
                                                                                         tcttGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagacaacttactgaaaggcgcggcggcacaagcggtacagtgcgccaatattcgtttcg                                                                 <  2:418872‑M1/144‑85 (MQ=255)
                                                                                            tgactcctatcagtggcCGGTGATCAACCCCACCACGGGCGTGGGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagacaacttactgaaaggcgcggcggcacaagcggtacagtgcgccaatattcgtttcggct                                                              <  2:124483‑M1/136‑88 (MQ=255)
                                                                                               ccTCCATCAGGGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgatcattgtggcgaccgaagacaacttactgaaaggcgcggcggcacaagcggtacagtgcgccaatattcgtttcggctatg                                                           <  2:113610‑M1/144‑91 (MQ=255)
                                                                                                cTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaacgaagacaacttactgaaaggcgcggcggcaccagcggtgcagtgcgccaatatggcttgcggctatgg                                                          >  2:32895‑M1/1‑53 (MQ=255)
                                                                                                 tCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagacaacttactgaaaggcgcggcggcacaagcggtacagtgcgccaatattcgtttcggctatgcggaaac                                                    >  1:603677‑M1/1‑52 (MQ=255)
                                                                                                 tCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagacaacttactgaaaggagcggcggcacaaccggtacagccggccagtccccgtgtcggctatgcggggag                                                    >  1:398019‑M1/1‑52 (MQ=255)
                                                                                                 ccagtcggggtccggtgatcaacgCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagacaacttaccgaaaggcgcggcggcacaagcggtacagtgcgccaatattcgtttcggctatgcg                                                         <  2:537209‑M1/133‑93 (MQ=255)
                                                                                                  cAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagacaacttactgaaaggcgcggcggcacaagcggtacagtgcgccaatattcgtttcggctatgcgg                                                        <  2:472433‑M1/144‑94 (MQ=255)
                                                                                                         tGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagacaacgtagtgaagggcaggacggcacaagcggtacagtgcgccaatattcgtttcggctatgcggaaacgccgtctc                                            >  1:445666‑M1/1‑44 (MQ=255)
                                                                                                            ccGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagacaacttactgaaaggcgcggcggcacaagcggtacagtgcgccaatattcgtttcggctatgcggaaacgcagtctctta                                         >  1:537443‑M1/1‑41 (MQ=255)
                                                                                                             cGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagacaacttactgaaaggcgcggcggcacaagcggtacagtgcgccaacatgcggttcggctatgcggaaactcagtcg                                             >  2:472095‑M1/1‑40 (MQ=255)
                                                                                                                 gATCAAGGCGACCAGCGGGGTGAGCGGTGCGGGGCGagcgtctgattattgtggcgaccgaagacaacttactgaaaggcgcggcggcacaagcggtacagtgcgccaatattcgtttcggctatgcggaaacgcagtctctta                                         <  2:56506‑M1/144‑109 (MQ=255)
                                                                                                                  aTCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagacaacttactgaaaggcgcggcggcacaagcggtacagtgcgccaatattcgtttcggctatgcggaaacgcaccactccgc                                        <  2:46886‑M1/144‑110 (MQ=255)
                                                                                                                   tCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagacaacttactgaaaggcgcggcggcacaagcggtacagtgcgccaatattcgtttcggctatgcggaaacgcagtctcttatttaagg                                  >  1:855183‑M1/1‑34 (MQ=255)
                                                                                                                   tCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgaggcgaccgaagacaacttactgaaaggcgcggcggcacaagcggtacagtgcgccaatattcgtttcggctatgcggaaacgcagtctcttatttaagg                                  >  1:663489‑M1/1‑34 (MQ=255)
                                                                                                                     aaCGCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagacaacttactgaaaggcgcggcggcacaagcggtacagtgcgacaatactcggttcggctatgcggaaacgcagcctcctactga                                     >  2:506430‑M1/1‑32 (MQ=255)
                                                                                                                     aaCGCCACCAGCGGCGTGAGCGGTGAAGGGCGagcatctgattattgtggcgaccgaagacaacttactgaaaggcgcggcggcacaagcggtacagtgcgccaatattcgttgcggctctgcggaaccgcagccccttctttaagggt                                >  1:554164‑M1/1‑32 (MQ=255)
                                                                                                                        gCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagacaacttactgaaaggcgcggcggcacacgcggtgcagtgcgccaatattcgtttcggctatgaggaaaagcagtctcttatttaagggggca                             >  1:13679‑M1/1‑29 (MQ=255)
                                                                                                                        gCCACCAGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagacaacttactgaaaggcgcggcggcacaagcggtacagtgcgccaatattcgtttcggctatgcggaaacgcagtctcttatttaagggtgga                             >  1:273686‑M1/1‑29 (MQ=255)
                                                                                                                         ccaccaGCGGCGTGAGCGGTGCAGGGCGagcatctgattattgtggcgaccgaagacaacttactgaaaggcgcggcggcacaagcggtacagtgcgccaatattcgtttcggctatgcggaaacgcagtctcttatttaagggtgcaa                            >  1:802408‑M1/1‑28 (MQ=255)
                                                                                                                                     tGAGCGGTTCAGGGCGaggatctgattattgtggcgaccgaagacaacttactgaaaggcgcggcggcacaagcggtacagtgcgccaatattcgtttcggctatgcggaaacgcagtctcttatttaagggtgcaatgatgaatccat                <  1:59077‑M1/149‑134 (MQ=255)
                                                                                                                                        gCGGTGCAGGGCGagcatctgattattgtggcgaccgaagacaacttactgaaaggcgcggcggcacaagcggtacagtgcgccaatcttcgtttcggctatgcggaaacgcagtctcttatttaagggtgcaatgatgaatccattaa             <  1:872792‑M1/149‑137 (MQ=255)
                                                                                                                                         cGGTGCAGGGCGagcatctgattattgtggcgaccgaagacaacttactgaaaggcgcggcggcacaagcggtacagtgcgccaatattcgtttcggctatgcggaaacgcagtctcttatttaagggtgcaatgatgaatacc                 >  2:368220‑M1/1‑12 (MQ=255)
                                                                                                                                          ggTGCAGGGCGagcatctgattattgtggcgaccgaagacaacttactgaaaggcgcggcggcacaagcggtacagtgcgccaatattcgtttcggctatgcggaaacgcagtctcttatttaagggtgcaatgatgaatccat                <  2:521307‑M1/144‑134 (MQ=255)
                                                                                                                                            tGCAGGGCGagcatctgattattgtggcgaccgaagacaacttactgaaaggcgcggcggcacaagcggtacagtgcgccaatattcgtttcggctatgcggcaacgcagtctcttatttaagggtgcaatgatgaatccattaagtat         >  1:767613‑M1/1‑9 (MQ=255)
                                                                                                                                              cAGGGCGagcatctgattattgtggcgaccgaagacaacttactgaaaggcgcggcggcacaagcggtacagtgcgccaatattcgtttcggctatgcggaaacgcagtctcttatttaagggtgcaatgatgaatccattaattatca       >  1:529659‑M1/1‑7 (MQ=255)
                                                                                                                                                 ggCGggcatctgattattgtggcgaccgaagacaatttacagaaaggcgcggcggcacaagcggtacagtgcgccaatattcgtttcggctatgcggaaacgcagtctcttatttaagggtgcaatgatgaatccattaattat         <  2:88714‑M1/144‑141 (MQ=255)
                                                                                                                                                   cGagcatctgattattgtggcgaccgaagacaacttactgaaaggcgcggcggcacaagcggtacagtgcgccaatattcgtttcgtctatgcggaaacgcagtctcttatttaagggtgcaatgatgaatcccttaattatcaaactg  >  1:627512‑M1/1‑2 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATCGCCACACACCTCGGTGCTGACGTTATCTTCACCCCACATCTGGGCAATTTCCCGCGCGGCATTCT  >  CP009273/4145352‑4145647

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