Predicted mutation
evidence seq id position mutation annotation gene description
RA 87 19,512 G→A intergenic (–/–)  /  –/–

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*8719,5120GA83.9% 220.0 / 32.6 93intergenic (–/–)–/––/–
Reads supporting (aligned to +/- strand):  ref base G (6/9);  new base A (47/31);  total (53/40)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

CAAGGACTGACGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGAAGACGGGTTTCTTCTCTGC‑‑‑‑‑‑‑‑‑‑  >  87/19412‑19524
                                                                                                    |                      
cAAGGACTGACGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGAAGACggg                        <  2:1326787/101‑1 (MQ=255)
  accgaaaggagactgccacgcataatttcttcagtcctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGAtt                      <  1:483648‑M2/49‑1 (MQ=255)
  accgaaaggagactgccacgcataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGAtt                      >  1:1128387‑M2/53‑101 (MQ=255)
  accgaaaggagactgccacgcataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGAtt                      <  1:1304569‑M2/49‑1 (MQ=255)
   ggACTGACGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGAAGACGGGttt                     <  2:922821/101‑1 (MQ=255)
   ccgaaaggagactgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGAttt                     <  2:1486429‑M2/52‑1 (MQ=255)
   ccgaaaggagactgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGAttt                     <  2:39195‑M2/52‑1 (MQ=255)
   ccgaaaggagactgccacgcataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGAttt                     >  2:871277‑M2/52‑101 (MQ=255)
    cgaaaggagactgccacgcataatttctttagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATttc                    >  1:697420‑M2/51‑101 (MQ=255)
    cgaaaggagactgccacgcataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATttc                    >  2:792193‑M2/51‑101 (MQ=255)
    cgaaaggagactgccacgcataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATttc                    <  1:870625‑M2/51‑1 (MQ=255)
    cgaaaggagactgccacgcataatttctccagccctttccattcaaacgtGAAATCAATAAGTACATCAGGGACATCGCAAACAATATGAAGACGGATttc                    <  2:533406‑M2/51‑1 (MQ=255)
    cagtaaaattaacacacattgtgaattttcagtcctttccattcaaacgtgaaatcaataagttcatcagggacatcgcAAACAATATGAATACGGATttc                    <  1:1481573‑M2/22‑1 (MQ=255)
    cagtaaaattaacacacattgtgaattttcagtcctttccattcaaacgtgaaatcaataagttcatcagggacatcgcAAACAATATGAATACGGATttc                    <  2:660616‑M2/22‑1 (MQ=255)
     gaaaggagactgccacgcataatttctttagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATttct                   <  2:1473867‑M2/52‑1 (MQ=255)
     gaaaggagactgccacgcataatttcttcagtcctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATttct                   >  1:866784‑M2/50‑101 (MQ=255)
     gaaaggagactgccacgcataatttcttcagtcctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATttct                   <  1:225780‑M2/52‑1 (MQ=255)
     gaaaggagactgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGATttct                   >  2:629765‑M2/48‑101 (MQ=255)
     agtaaaattaacacacattgtgaattttcagtcctttccattcaaacgtgaaatcaataagttcatcagggacatcgcAAACAATATGAATACGGATttct                   >  2:1442157‑M2/79‑101 (MQ=255)
      gtaaaattaacacacattgtgaattttcagtcctttccattcaaacgtgaaatcaataagttcatcagggacatcgcAAACAATATGAATACGGATttctt                  <  1:267975‑M2/24‑1 (MQ=255)
      cTGACGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGAAGACGGGTttctt                  >  1:752943/1‑101 (MQ=255)
      cTGACGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGAAGACGGGTttctt                  <  1:521406/101‑1 (MQ=255)
      aaaggagactgccacgcataatttctttagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATttctt                  >  1:1077777‑M2/49‑101 (MQ=255)
      aaaggagactgccacgcataatttcttcagtcctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATttctt                  <  2:209049‑M2/53‑1 (MQ=255)
       tGACGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGAAGACGGGTttcttc                 >  1:173128/1‑101 (MQ=255)
       tGACGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGAAGACGGGTttcttc                 <  2:1512896/101‑1 (MQ=255)
       aaggagactgccacgcataatttctttagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATttcttc                 >  1:280914‑M2/48‑101 (MQ=255)
       aaggagactgccacgcataatttcttcagtcctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATttcttc                 >  1:1374260‑M2/48‑101 (MQ=255)
       aaggagactgccacgcataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATttcttc                 >  1:774242‑M2/48‑101 (MQ=255)
       aaggagactgccacgcataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATttcttc                 <  1:1035360‑M2/54‑1 (MQ=255)
       aaggagactgccacgcataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATttcttc                 >  2:607614‑M2/48‑101 (MQ=255)
        aggagactgccacgcataatttctttagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTtcttct                <  1:1156362‑M2/55‑1 (MQ=255)
        aggagactgccacgcataatttctttagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTtcttct                <  1:206541‑M2/55‑1 (MQ=255)
        aggagactgccacgcataatttctttagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTtcttct                >  2:1304349‑M2/47‑101 (MQ=255)
        aggagactgccacgcataatttcttcagtcctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTtcttct                <  1:1665166‑M2/55‑1 (MQ=255)
        aggagactgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGATTtcttct                >  1:457247‑M2/45‑101 (MQ=255)
        aggagactgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGATTtcttct                >  2:1150829‑M2/45‑101 (MQ=255)
        aggagactgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGAGAAGAAGGatttcttat                >  2:1641390‑M2/45‑99 (MQ=255)
        aggagactgccacgcataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTtcttct                <  1:1536446‑M2/55‑1 (MQ=255)
        aaaattaacacacattgtgaattttcagtcctttccattcaaacgtgaaatcaataagttcatcagggacatcgcAAACAATATGAATACGGATTtcttct                <  1:47529‑M2/26‑1 (MQ=255)
        aaaattaacacacattgtgaattttcagtcctttccattcaaacgtgaaatcaataagttcatcagggacatcgcAAACAATATGAATACGGATTtcttct                >  2:1349707‑M2/76‑101 (MQ=255)
          gaggctgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTGATCAGGGACATCGCGAACAATGTCAAGACGGATTTCttctcc              <  2:1339081‑M2/59‑2 (MQ=255)
          gagactgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGATTTCTTctct              <  2:950773‑M2/59‑1 (MQ=255)
          cGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGAAGACGGGTTTCTTctct              <  2:1612132/101‑1 (MQ=255)
          cGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGAAGACGGGTTTCTTctct              >  1:1013473/1‑101 (MQ=255)
          cGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGAAGACGGGTTTCTTctct              <  1:322428/101‑1 (MQ=255)
          aattaacacacattgtgaattttcagtcctttccattcaaacgtgaaatcaataagttcatcagggacatcgcAAACAATATGAATACGGATTTCTTctct              >  1:183308‑M2/74‑101 (MQ=255)
           ngactgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGATTTCTTCTCTg             >  1:861138‑M2/42‑101 (MQ=255)
           gTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGAAGACGGGTTTCTTCTCTg             <  1:1003762/101‑1 (MQ=255)
           gTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGAAGACGGGTTTCTTCTCTg             <  1:636762/101‑1 (MQ=255)
           attaacacacattgtgaattttcagtcctttccattcaaacgtgaaatcaataagttcatcagggacatcgcAAACAATATGAATACGGATTTCTTCTCTg             >  2:252154‑M2/73‑101 (MQ=255)
           agactgccacgcataatttctttagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTTCTTCTCTg             <  1:716778‑M2/58‑1 (MQ=255)
           agactgccacgcataatttcttcagtcctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTTCTTCTCTg             >  1:1501981‑M2/44‑101 (MQ=255)
           agactgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGATTTCTTCTCTg             <  1:241329‑M2/60‑1 (MQ=255)
           agactgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGATTTCTTCTCTg             >  1:1275765‑M2/42‑101 (MQ=255)
           agactgccacgcataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTTCTTCTCTg             <  1:1593146‑M2/58‑1 (MQ=255)
            ttaacacacattgtgaattttcagtcctttccattcaaacgttaaatcaataagttcatcaggtacatcgcAAACAATATGAATACGGATTTCTTCTCTGc            <  1:985547‑M2/30‑1 (MQ=255)
            gactgccacgcataatttcttcagtcctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTTCTTCTCTGc            >  2:163516‑M2/43‑101 (MQ=255)
            gactgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGATTTCTTCTCTGc            >  2:1650876‑M2/41‑101 (MQ=255)
             taacacacattgtgaattttcagtcctttccattcaaacgtgaaatcaataagttcatcagggacatcgcAAACAATATGAATACGGATTTCTTCTCTgcc           >  1:1233426‑M2/71‑100 (MQ=255)
             actgccacgcataatttcttcagtcctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTTCTTCTCTgcc           >  2:1651310‑M2/42‑100 (MQ=255)
             actgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGATTTCTTCTCTgcc           <  2:79660‑M2/62‑2 (MQ=255)
             actgccacgcataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTTCTTCTCTgcc           >  2:97019‑M2/42‑100 (MQ=255)
              atgccacgcataatttcttcagtcctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTTCTTCTCtgccg          <  1:523366‑M2/61‑3 (MQ=255)
              aGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGAAGACGGGTTTCTTCTCTGccg          >  2:1436001/1‑99 (MQ=255)
              aGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGAAGACGGGTTTCTTCTCTGccg          >  1:239988/1‑99 (MQ=255)
               tgccacgcataatttctttagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTTCTTCTctgccgt         >  2:205568‑M2/40‑98 (MQ=255)
               tgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGATTTCTTCTctgccgt         >  2:293309‑M2/38‑98 (MQ=255)
               tgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGATTTCTTCTctgccgt         >  2:364635‑M2/38‑98 (MQ=255)
               tgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGATTTCTTCTctgccgt         >  1:1032233‑M2/38‑98 (MQ=255)
               tgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGATTTCTTCTctgccgt         >  1:422316‑M2/38‑98 (MQ=255)
                gccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGATTTCTTCtctgccgtg        <  2:305946‑M2/65‑5 (MQ=255)
                 tAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGAAGACGGGTTTCTTCTCTGccgtga       <  2:1458945/101‑6 (MQ=255)
                 ccacgcataatttcttcagtcctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTTCTTctctgccgtga       >  1:1177281‑M2/38‑96 (MQ=255)
                 acacattgtgaattttcagtcctttccattcaaacgtgaaatcaataagttcatcagggacatcgcAAACAATATGAATACGGATTTCTTctctgccgtga       >  1:1402575‑M2/67‑96 (MQ=255)
                  cacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGATTTCTtctctgccgtgac      >  2:1238466‑M2/35‑95 (MQ=255)
                   gTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGAAGACGGGTTTCTTCTCTGccgtgact     >  1:992893/1‑94 (MQ=255)
                   gTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGAAGACGGGTTTCTTCTCTGccgtgact     <  1:184086/101‑8 (MQ=255)
                   acgcataatttctttagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTTCttctctgccgtgact     <  2:1035343‑M2/66‑8 (MQ=255)
                   acgcataatttcttcagtcctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTTCttctctgccgtgact     <  1:123463‑M2/66‑8 (MQ=255)
                   acgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGATTTCttctctgccgtgact     >  1:1201447‑M2/34‑94 (MQ=255)
                   acgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGATTTCttctctgccgtgact     >  1:1435330‑M2/34‑94 (MQ=255)
                   acattgtgaattttcagtcctttccattcaaacgtgaaatcaataagttcatcagggacatcgcAAACAATATGAATACGGATTTCttctctgccgtgact     >  2:25257‑M2/65‑94 (MQ=255)
                   acattgtgaattttcagtcctttccattcaaacgtgaaatcaataagttcatcagggacatcgcAAACAATATGAATACGGATTTCttctcggccgtgact     >  2:954269‑M2/65‑94 (MQ=255)
                    cgcataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTTcttctctgccgtgactc    <  1:1468830‑M2/67‑9 (MQ=255)
                     gcataatttctttagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTtcttctctgccgtgactct   >  2:565134‑M2/34‑92 (MQ=255)
                     gcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGATTtcttctctgccgtgactct   >  1:1241548‑M2/32‑92 (MQ=255)
                     gcataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTtcttctctgccgtgactct   >  1:1518419‑M2/34‑92 (MQ=255)
                      cataatttctttagccctttcctttcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATatcttctctgccgtgactctt  >  2:1653265‑M2/33‑91 (MQ=255)
                      cataatttctttagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATttcttctctgccgtgactctt  >  2:530912‑M2/33‑91 (MQ=255)
                      cataatttcttcagtcctttccattcaaacgtAAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATttcttctctgccgtgactctt  >  1:1314948‑M2/33‑91 (MQ=255)
                      cataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGATttcttctctgccgtgactctt  <  1:1564453‑M2/71‑11 (MQ=255)
                      cataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGATttcttctctgccgtgactctt  >  1:518901‑M2/31‑91 (MQ=255)
                      cataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATttcttctctgccgtgactctt  <  2:911084‑M2/69‑11 (MQ=255)
                      cataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATttcttctctgccgtgactctt  >  2:1156989‑M2/33‑91 (MQ=255)
                      cataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATttcttctctgccgtgactctt  >  1:1117617‑M2/33‑91 (MQ=255)
                                                                                                    |                      
CAAGGACTGACGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGAAGACGGGTTTCTTCTCTGC‑‑‑‑‑‑‑‑‑‑  >  87/19412‑19524

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