Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,903 G→A E523K (GAA→AAA)  thrA → fused aspartate kinase/homoserine dehydrogenase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,9030GA94.5% 233.1 / ‑3.0 73E523K (GAA→AAA) thrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base G (1/3);  new base A (35/34);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCG  >  NC_000913/1660‑2138
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            
gTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACGTGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACt                                                                                                                                                                                                                                      <  1:152025/251‑1 (MQ=255)
     aaaTAACGATGATTCGACCACTTGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTTC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTTGCAGAAAG‑AACTTgcg                                                                                                                                                                                                                                  <  1:517500/250‑1 (MQ=255)
     aaaTAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATCCCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGgcg                                                                                                                                                                                                                                  <  1:517518/250‑1 (MQ=255)
     aaaTAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGGGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGgcgc                                                                                                                                                                                                                                 >  1:232546/1‑251 (MQ=255)
     aaaTAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGGTATCGAAGTGTTTGTGATTGGCATCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGCAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTTGAAAACTGGCAGAAAG‑AACCGgcgc                                                                                                                                                                                                                                 >  2:149342/1‑251 (MQ=255)
         aaCGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAg                                                                                                                                                                                                                              <  2:131127/250‑1 (MQ=255)
               gatgCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCaa                                                                                                                                                                                                                          >  1:303440/1‑248 (MQ=255)
                    gACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGt                                                                                                                                                                                                                  >  2:160730/1‑251 (MQ=255)
                    gACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGAGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGGGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGt                                                                                                                                                                                                                  >  1:478130/1‑251 (MQ=255)
                        aCTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTAAGTGTCTGCGGTGTTGC‑CACATCGAAGGCTCTGCTCCCCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTa                                                                                                                                                                                                               >  1:520833/1‑250 (MQ=255)
                         cTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAAt                                                                                                                                                                                                             >  1:515102/1‑251 (MQ=255)
                           ggCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTCTCGAAGTGTTTGTGATTGGCGTCGGCGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACCTATCGACTTACGCGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTAAACAATGTACATGGCCTTAATATGGCAAACTGGAAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAAtct                                                                                                                                                                                                           >  1:395515/1‑251 (MQ=255)
                           ggCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGATCCGTTTAAt                                                                                                                                                                                                             <  2:36130/249‑1 (MQ=255)
                               ggcgcGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCggg                                                                                                                                                                                                       <  2:205281/250‑1 (MQ=255)
                                gcgcGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGATCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGgc                                                                                                                                                                                                      <  1:46700/251‑1 (MQ=255)
                                    gTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCtt                                                                                                                                                                                                  >  1:313570/1‑251 (MQ=255)
                                                 tGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTg                                                                                                                                                                                     >  1:162148/1‑251 (MQ=255)
                                                            cccGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGGTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCAAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATatca                                                                                                                                                                          <  2:406985/250‑1 (MQ=255)
                                                             ccGAGCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAAGCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATatcat                                                                                                                                                                         <  1:69576/251‑1 (MQ=255)
                                                              cGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCACAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATatcatc                                                                                                                                                                        >  1:188113/1‑251 (MQ=255)
                                                               gatccgggtATCAAAGTGTTTGTGATTGGCAACGTGGGCGTTGGCGGGACGCTGCTGGAGCAATTGAAGCGTCAGTAAAGGTGGCTGCTGAATTAACATATAGCATTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCCGTTCCCCAATGTACATGGCCTTAATCCGGACAACTGGCAGGACG‑CACTGGCGCACACCAACGACCCGTTTAATCCCCGGCGCTTAATTCGCCTAGTGGAAGAATatcatcg                                                                                                                                                                       >  2:214980/9‑250 (MQ=255)
                                                                aTCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCTCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATctg                                                                                                                                                                      >  2:104450/1‑251 (MQ=255)
                                                                  cAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATctgct                                                                                                                                                                    <  1:523986/251‑1 (MQ=255)
                                                                    ggTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGa                                                                                                                                                                  <  1:59802/251‑1 (MQ=255)
                                                                          cGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCgg                                                                                                                                                            <  2:313570/251‑1 (MQ=255)
                                                                             aGTGTTTTTGATTGGCGGCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCa                                                                                                                                                         <  2:310309/251‑1 (MQ=255)
                                                                               tgtTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGAAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCAtt                                                                                                                                                       >  2:292577/1‑251 (MQ=255)
                                                                                  ttGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCAttgtt                                                                                                                                                    >  1:146155/1‑251 (MQ=255)
                                                                                  ttGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACAGGGCCTGAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCGCGGGCGTTTAATGCGCCTCGTGAAAGAATATCATctgcggtactcgtgcattgt                                                                                                                                                     <  2:252172/250‑17 (MQ=255)
                                                                                            cGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCGTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACtt                                                                                                                                          <  2:104779/251‑1 (MQ=255)
                                                                                               cGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCCGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGCTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGACAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTcc                                                                                                                                        >  1:539450/1‑250 (MQ=255)
                                                                                                 gTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTccagc                                                                                                                                     <  1:11806/251‑1 (MQ=255)
                                                                                                 gTGGCGTTGGCCGTGAGCTGCTGGCGCCACTGAAGCGTCAGCAACGCTGGCTGAAGACTAAACATCTCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCGGCTAACCAATGTACAAGGCCTTAATCTGGAACACTGGAAGATAG‑AACTGGCGCAAGCCAAAGAGCAGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATATTATGCTGAACCCGGTCAATGTTGATTGCACTTccagc                                                                                                                                     >  2:36979/1‑251 (MQ=255)
                                                                                                  tGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTccagcc                                                                                                                                    >  2:166620/1‑251 (MQ=255)
                                                                                                  tGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTccagcc                                                                                                                                    >  1:110167/1‑251 (MQ=255)
                                                                                                    gCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTccagccag                                                                                                                                  <  1:289812/251‑1 (MQ=255)
                                                                                                     cGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGACTTAATCTGGAAAACTGGCAGAAAG‑ACCTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCACCCAgg                                                                                                                                 >  2:158629/1‑251 (MQ=255)
                                                                                                      ggtGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTTCGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAATCTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGTGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGc                                                                                                                                <  1:556671/249‑1 (MQ=255)
                                                                                                      gTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGACTCGTGAAAGCATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGc                                                                                                                                >  1:533408/1‑251 (MQ=255)
                                                                                                         ggCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGt                                                                                                                             <  1:160730/251‑1 (MQ=255)
                                                                                                                cgcTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCggtt                                                                                                                      >  1:340811/1‑249 (MQ=255)
                                                                                                                       ggAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTTATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAatat                                                                                                                <  2:546192/250‑1 (MQ=255)
                                                                                                                       ggAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATg                                                                                                               <  1:173143/251‑1 (MQ=255)
                                                                                                                           cAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGa                                                                                                           <  1:521955/251‑1 (MQ=255)
                                                                                                                                gAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATATGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTcc                                                                                                      >  1:293012/1‑251 (MQ=255)
                                                                                                                                    cgacAGCAAGGCTGGCGGAAGAATAAACATGTCGAGTATCGGGTCGGCGGTGT‑GCTCACCTAGCAGACTCTGCTCACCAATGTTCATGGCCGTAATCTGGAAAACTGGC‑TGAAGTAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTGATTCGCCTCGTGAAAGAAGATCGTCTGCTGAACCCGGTCATGGTTGGTTGCTCTTCCTGATAGGCAGTGGGGGATCAAGATGCCGACGTCCTgcg                                                                                                  <  2:566229/248‑1 (MQ=255)
                                                                                                                                         gCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAg                                                                                             <  2:75802/251‑1 (MQ=255)
                                                                                                                                           aaaGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGt                                                                                           >  1:246813/1‑251 (MQ=255)
                                                                                                                                               gttgctgAAGGATAGACATTTCGACTTACGTGTCGGCGGTGTTGC‑CGATTTGAAGGCGCGGGTCACCAATGTACATGTCCGTAATCTGGAAAACTGGGAGGAAG‑AACTGGCGCAAGACAAAGTGCCGTTTAAGCTCGGGCGCTTAATTCGCCTCGTGAAAGAAGATCAGCTGCGGAGCCCGGGCATTTTTGTGTGGGCTTCCAGGCAGGCAGTGGCGGATCAATATGCGGACGTCCTGCGGGAAGGTTTcc                                                                                       <  2:73754/248‑1 (MQ=255)
                                                                                                                                                 ggctgAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTc                                                                                                                                                                                             <  1:368347/147‑1 (MQ=255)
                                                                                                                                                 ggctgAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTc                                                                                                                                                                                             >  2:368347/1‑147 (MQ=255)
                                                                                                                                                  gctgAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGA‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCAGTCATTGTTGACTGCACTTCCAGCCAGGACGTGGCGGATCAACGTGCCGACTTCCCGCGCGACGGTTTCAACg                                                                                    >  2:487461/1‑251 (MQ=255)
                                                                                                                                                        gaaTAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCa                                                                              <  1:491539/251‑1 (MQ=255)
                                                                                                                                                                atatCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGTAGCAATATGCCGACTTCCTGCTCGAAGGTTTCCTCGTTGTCACGCCGAAc                                                                      <  1:553626/251‑1 (MQ=255)
                                                                                                                                                                        ttACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGc                                                              >  2:444018/1‑251 (MQ=255)
                                                                                                                                                                              gtCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAaca                                                         <  1:66191/250‑1 (MQ=255)
                                                                                                                                                                              gtCTGCGGGGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCAGGAAAACTGGCAGAACG‑AAATCGCGCAAGACAAAGAGCCGTTTCATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGCTCAATATGGCGCCTTCCTGCGCGAAGGCTTCCACGGTGTCACGCCGACAAAAAAGGGCAacac                                                        >  2:337963/1‑251 (MQ=255)
                                                                                                                                                                                    ggTGTTGC‑CAACTCGAAGGCTCTGCTCACCAAGGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCtcgtc                                                  <  2:232546/251‑1 (MQ=255)
                                                                                                                                                                                          gC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGCTCAATTTGCCGACTTCCTGCGCGAAGGGTTCCACGGTGTCACGCCGAACAAAAAGGCCAACACATCGTCGATgg                                             >  1:50502/1‑250 (MQ=255)
                                                                                                                                                                                              aaCTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGACGCCTTCCTGCGCGAAGGGTTCCACGGTGTCAAGCCGCACAAAAAGGCCAACACATCGTCGATggctt                                          >  2:480865/1‑247 (MQ=255)
                                                                                                                                                                                               aCTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATtac                                        <  1:374403/251‑1 (MQ=255)
                                                                                                                                                                                                         tctGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGt                              >  1:453373/1‑251 (MQ=255)
                                                                                                                                                                                                          ctGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATTCCGGCTTCCTGCGCGAAGGTTTCCACGTTGTTCCGCCGACCAAAAAGGCCAACACCTCGTCGATGGATTAATCCCATCaagt                             >  2:413162/1‑248 (MQ=255)
                                                                                                                                                                                                            gCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGtt                             <  1:443712/249‑1 (MQ=255)
                                                                                                                                                                                                               cACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGtt                        >  1:525782/1‑251 (MQ=255)
                                                                                                                                                                                                                 ccAATGTACATGGCCTTAATCTGGAAAACCGGCAGGTAT‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCACTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGAGCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTTCCAGCAGTTGCGTTag                      <  2:514197/251‑2 (MQ=255)
                                                                                                                                                                                                                  cAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATg                     >  1:342296/1‑251 (MQ=255)
                                                                                                                                                                                                                          aTGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGATTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATAAATATGCCGACTTCCTGCGCGAAGGGTTCCACGTTGTCACGACGAACAAAAAGGCCAACACCTCGTCGCTGGATTACTTCCATCCGTTGCGTTATGCGGCGGa              >  2:428681/1‑250 (MQ=255)
                                                                                                                                                                                                                           tGGCCTTAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGGATTCCTGCGCGAAGGTTTCCCCGTTGTCACGCCGAACAAAAAGGCCCACACCTCGTCGATGGCTTTCTACCATCCGTTGCGTTATGCAGCGGaa             >  2:3577/1‑250 (MQ=255)
                                                                                                                                                                                                                            ggCCTGAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATAGCATCTGCTGAACCCGGTCATTGTTGACTGCACGTCCAGCCAGGCAGTGGCGGATCAAGATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGGTGCGTTATGCGGCGGaaa            <  2:515102/250‑1 (MQ=255)
                                                                                                                                                                                                                                 tAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATcg       <  2:146155/250‑1 (MQ=255)
                                                                                                                                                                                                                                 tAATCTGGAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATc        <  1:464213/249‑1 (MQ=255)
                                                                                                                                                                                                                                       ggAAAACTGGCAGAAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCgcggcg  <  1:544117/249‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            
GTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGC‑CAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAG‑AACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCG  >  NC_000913/1660‑2138

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